编程人生
文章平均质量分 86
Kino Chan
这个作者很懒,什么都没留下…
展开
-
Fluent的msh格式网格学习
本文介绍Fluent的msh格式文件的一些内在含义,方便自编程。原创 2022-06-23 12:54:18 · 5611 阅读 · 1 评论 -
OpenFOAM的泊肃叶流动算例
泊肃叶流动是流体力学里面一个基本算例,排在库埃特流动后面。其应用很广(管道流动),有理论解,经常被用来验证算例、测试算法等。本文介绍如何用OpenFOAM复现泊肃叶流动。原创 2022-03-25 10:14:27 · 2483 阅读 · 0 评论 -
OpenFOAM的圆柱绕流算例记录
圆柱绕流是流体力学里最经典的算例之一,受到众多学者的“喜爱”。初学CFD的人除了空腔、管道流之外,最先接触的莫过于圆柱绕流了。圆柱绕流能反映流体流过障碍物时的各种情况:层流向湍流的转捩、流体的剪切与分离、边界层、卡门涡街现象等等,是深入理解流动机理的重要算例。因为经常要算圆柱绕流,在这里做一些记录,记录一些细节,以防以后忘掉。这次做Re=200时的圆柱绕流算例,画网格用ICEM,CFD软件用开源软件OpenFOAM,后处理用tecplot。1. 前处理:画网格用ICEM画网格。圆柱直径D=1m,计算原创 2021-09-08 23:58:13 · 12257 阅读 · 9 评论 -
“cannot be called for a calculatedFvPatchField“——OpenFOAM中实现ρU(rhoU)边界条件的源码学习
因研究需要,特写一篇非单一变量(ρU, rhoU)边界条件的实现过程。在解可压NS方程时,rhoCentralFoam(解析)对动量方程的ρU进行直接的插值求解。换句话说,以ρU作为一个守恒变量,在方程中先进行求解,再分别更新ρ和U来。我们知道, 有限体积法的边界条件就是计算边界面上的通量。在OpenFOAM中单一变量ρ和U的边界条件都必须要分别设置好,这意味着它们边界面上的通量如何计算将由我们来进行决定,其中最常用的有fixedValue(第一类边界条件)和zeroGradient(第二类边界条件)。*原创 2021-09-03 16:58:15 · 1522 阅读 · 2 评论 -
修改OpenFOAM中的pisoFoam使其按照最大库朗数修改最小时间步
前言在学习OpenFOAM中的LES模拟时,发现一开始的库朗数会比较大,我不得不将时间步调到1e-7秒。但当时间推进到1000步左右时,库朗数就降到0.2左右了(以我自己的算例来举例)。如果仍是以这样的时间步计算到最后,以来增加了计算的时间,二来库朗数太小又会影响计算结果(曾算过负的阻力系数)。因此,在计算的过程中将库朗数维持在某个值就显得尤为重要。要在OpenFOAM上实现这一功能,需要将solver(求解器)进行小改并重新编译。本文以pisoFoam为例将这一实现过程进行展示。网上能搜到关于这方面的原创 2020-12-19 17:51:20 · 4728 阅读 · 12 评论 -
用tecplot导出圆柱绕流中的表面平均压力系数
需要软件:tecplot15,matlab需要语音讲解请转到我录的b站视频,ID:带一司机,还上传了其他前处理(画网格),后处理视频等等,后续会继续更新,欢迎关注:)https://www.bilibili.com/video/BV1Mr4y1w71c/1. 导入算好的云图,这里是双圆柱绕流的算例,Re=2e6,使用DES湍流模型,经过一段时间计算之后,平均压力场如下:这里平均压力场(pMean)是用OpenFOAM的后处理工具导出的,这里不讨论其他CFD软件。1.5 计算压力系数原创 2020-11-08 16:53:20 · 9584 阅读 · 3 评论 -
持续不定期更新:CFD&C++之不可压库埃特流(Coutte)的数值解(隐式克兰克-尼科尔森方法)
之前一篇https://blog.csdn.net/weixin_39124457/article/details/101038218用到的是显式方法,这次尝试用隐式方法解决不可压库埃特流的问题。资料参考《计算流体力学基础及其应用》第9章第一节,电子书可自行搜索下载。代码自己打的,可下载:https://download.csdn.net/download/weixin_39124457/...原创 2020-01-13 19:32:11 · 1656 阅读 · 1 评论 -
持续不定期更新:CFD&C++之拟一维喷管流动的数值解(2)
这篇博文将剩下的问题解决完。第一篇在:https://mp.csdn.net/postedit/101038218在初始化步骤之后,就到了计算下一时间步的步骤了。计算之前先讲一讲这里用到的计算方法:显式麦考马克方法。简单地来讲,显式意味着可以根据当前时间步的流场量去算下一个时间步的流场量。相对的,隐式则要用到下一个时间步的流场量,换言之要求解方程组,计算及其复杂。麦考马克方法是比较早期的十...原创 2020-01-11 17:03:00 · 1648 阅读 · 3 评论 -
安装Openfoam的acoustic开源程序包libAcoustics
研究需要,从github上下载仅有的基于OpenFOAM的Acoustic程序,网址:https://github.com/unicfdlab/libAcoustics一开始我电脑上OpenFOAM的版本是v6的,编译程序包的时候终端显示缺少一些文件,于是问了问作者。对方回答说你的版本不对,于是我下载了他们最新的一个OpenFOAM-v1812的基础程序包OpenDTEF:https:...原创 2019-11-02 20:45:19 · 952 阅读 · 0 评论 -
持续不定期更新:CFD&C++之拟一维喷管流动的数值解(1)
前言:学习openfoam,或者准确来说,CFD,断断续续也有5个多月。3月底在本科学校写下的第一篇openfoam学习的博文,到现在在研一的学校里,尝试自己写代码计算cfd问题,慢慢一点一点地构建自己的自信心。其实说实话,学CFD很难,尤其是前期要战胜纯理论带来的「空虚」感,找到学习的落脚点。我不敢说现在的自己又足够的自信,但坚信「即使没有希望也要坚持下来」,终有一天会在学习中找到自我。...原创 2020-01-11 16:45:33 · 3915 阅读 · 2 评论 -
C++ 实现数字黑洞
昨晚看了一篇公众号的文章:https://mp.weixin.qq.com/s/idzZVuncx3XPcQwEZ4C3cg上面有介绍数字黑洞“6174”,论坛上也有很多,大概的历史由来不在赘述,随便找一个数如“1314”,看下图足矣:其中,流程分两步,一步是自己输入一个数字不全相同的四位数(1111,2222等等的不要),由程序计算出数字黑洞;然后由程序随机出100个四位数,分别算数...原创 2019-08-20 11:04:45 · 4130 阅读 · 0 评论