CFD中文网

@李东岳 好的，谢谢李老师，我试一试

这些国内软件包都是直接fork国外github上面的OpenFOAM。因为有时候github访问不上去。就fork到国内。

下列地址可直接使用git clone，例如，打开终端，在终端直接将下面的复制进去：

git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-6.git

即可在本地创建ThirdParty-6文件夹。如果提示你没有git，那么输入下面的即可：

sudo apt install git

正常的话是下图这样：

捕获.JPG

相关git地址： git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2312.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2306.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2206.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2106.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v2006.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1912.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1906.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-v1712.git

第三方包地址：

git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2306.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2212.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2206.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2106.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v2006.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-v1906.git

在下载的时候，需要下载相应的OpenFOAM与相应版本的第三方包

这些国内软件包都是直接fork国外github上面的OpenFOAM。因为有时候github访问不上去。就fork到国内。

下列地址可直接使用git clone，例如，打开终端，在终端直接将下面的复制进去：

git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-6.git

即可在本地创建ThirdParty-6文件夹。如果提示你没有git，那么输入下面的即可：

sudo apt install git

正常的话是下图这样：

捕获.JPG

相关git地址： git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-11.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-10.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-9.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-8.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-7.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-6.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-5.x.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-4.x.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-3.0.x.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-2.3.0.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-2.2.0.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-2.1.1.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-2.1.0.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-2.0.0.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-1.7.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-1.6.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-1.5.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-1.4.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-1.3.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-1.2.git git clone https://e.coding.net/dyfluid/OpenFOAM/OpenFOAM-1.1.git

第三方包地址：

git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-11.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-10.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-9.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-8.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-7.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-6.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-5.x.git git clone https://e.coding.net/dyfluid/OpenFOAM/ThirdParty-3.0.x.git

在下载的时候，需要下载相应的OpenFOAM与相应版本的第三方包

如果有人需要scotch_6.0.3: scotch_6.0.3.tar.gz

@lwjetmann 老师您好，再请教一下您。

1. 用wallshearstress后处理算壁面粘性剪切力，然后在paraview中计算出在壁面各点法向和切向的投影； 2. 基于wallshearstress functionObject进行修改，输出投影后的结果

按照您提到的上述方案1，对于圆球绕流（球体界面是“noSlip”边界条件），我顺利的输出了粘性剪切力，结果也验证的上。

想请您再帮忙看看另外一种情况：绕流气泡，即将边界条件改为Symmetry：
因为wallshearstress只对wall类型的边界计算粘性剪应力，因此不能再通过【方案1】输出界面上的粘性剪应力结果。
关于设置Symmetry边界条件的原因是源于这里：Stokes流中的圆球绕流的阻力系数，怎么设置才能算准。

对此，针对绕流气泡的结果（边界条件为Symmetry）：

我尝试手动在结果中将Symmetry边界条件改为Wall，然后执行一次wallshearstress； 通过对比文献中气泡的压差阻力系数、法向粘性应力和切向粘性应力结果，输出的结果并不对

除此以外，是否只有像您提到的【方案二】基于wallshearstress functionObject进行修改，才能输出“绕流气泡”的粘性剪切力结果。
如果您方便的话，能否请您指点一些修改思路，谢谢老师。

我在算liutex，本来准备仍udm里 但是因为要矩阵运算不想写轮子，就导出数据到csv里在fluent外算的。本来准备仍tecplot 但缺少网格拓扑 导入之后只有一堆线。想问有没有其他可视化的思路

我觉得你这个挺有意思的。初学可以练手用。但是个人觉得算法应该没有问题，出现不一致的情况可能还是哪里出现了小bug。

@小泽同学 我此前还在想是否需要重新考虑从右下角的CD/DVD中加载一个适用于Ubuntu14.04版本的linux.iso，里面有vm-tools的安装包，相当于也是用里面写好的脚本重装VMware-tools。

emm... 总之，能用上就太好了，祝收敛！

已经解决了！谢谢各位

@BlancheDH 好嘞，我去试试，谢谢大佬:xinxin:

@sxz0823 可以加个联系方式交流一下吗，我用turbinesFoam和LES但是结果很不对我的qq是：1141101064

@李东岳 嗯嗯，好的，谢谢老师

更新到了5.9.1节

@李东岳 我用的v2006，Phi是速度势，用-grad(Phi)计算的速度和输出的U有些差异，尤其是在远场边界。

@lizhisongsjtu 这个图太酷炫了。 有doi不

有些大佬对liutex法的态度偏消极(审稿意见)。

是，我也看到了。为什么呢

@YuPeng 我也遇到了相同的问题，请问一下您现在解决了吗

@RQLI_CFD 感谢，我研究研究

我补充一句，最近也在写湍流模型UDF，需要用到壁面距离，发现C_WALL_DIST不能用之后不信邪，测试了很多（一个多小时测了十几二十个case），终于测试出来发现：
C_WALL_DIST只有在选择了fluent已有湍流模型（我选了k-omega）的时候才能用，而当选择了laminar模式，用自己的湍流模型UDF的时候，如果涉及到C_WALL_DIST，就会报错。我的fluent版本为2022R2

感谢李老师！经过修改，已经可以成功运行了。:146:

请问在使用turbulentDFSEMInlet入口类型计算过程中，随着计算的进行，后半部分的风场慢慢的变得不是对数廓线了，请问是什么原因，有没有大佬知道。
我的初始边界信息：
计算域：xyz为1011（m）
U：

inlet { type turbulentDFSEMInlet; delta 1; U { type mappedFile; mapMethod nearest; } R uniform (0.8 0 0 0.4 0 0.4); L uniform 0.004; d 1; nCellPerEddy 1; scale 1; value $internalField; } outlet { type advective; value $internalField; } "(topWall|front|back)" { type zeroGradient; } "(bottomWall|Model)" { type noSlip; }

输入文件：
point：

( ( 0 0 0.01 ) ( 0 0 0.05 ) ( 0 0 0.1 ) ( 0 0 0.15 ) ( 0 0 0.2 ) ( 0 0 0.25 ) ( 0 0 0.3 ) ( 0 0 0.35 ) ( 0 0 0.4 ) ( 0 0 0.45 ) ( 0 0 0.5 ) ( 0 0 0.55 ) ( 0 0 0.6 ) ( 0 0 0.65 ) ( 0 0 0.7 ) ( 0 0 0.75 ) ( 0 0 0.8 ) ( 0 0 0.85 ) ( 0 0 0.9 ) ( 0 0 0.95 ) ( 0 0 0.99 ) )

U：

( ( 0 0 0 ) ( 1.048097372 0 0 ) ( 1.499488341 0 0 ) ( 1.763535131 0 0 ) ( 1.95087931 0 0 ) ( 2.096194744 0 0 ) ( 2.2149261 0 0 ) ( 2.31531203 0 0 ) ( 2.402270279 0 0 ) ( 2.47897289 0 0 ) ( 2.547585713 0 0 ) ( 2.609653562 0 0 ) ( 2.666317069 0 0 ) ( 2.718442442 0 0 ) ( 2.766702999 0 0 ) ( 2.811632503 0 0 ) ( 2.853661247 0 0 ) ( 2.893141184 0 0 ) ( 2.930363858 0 0 ) ( 2.965573488 0 0 ) ( 2.992431707 0 0 ) )

计算3600s后，平均风速廓线，U1U2U3U4U5U6U7U8U9的位置为0.1，2，4，5，6，7，8，9，9.9
e595f49f-347a-4180-82e2-b6295cd45a15-图片.png

lref设置为直径D=1，Aref=D×径向长度，改成2.68328就ok了，你使用什么软件导入的网格，fluent.mesh导入的时候会告诉你径向拉伸了多少，你再尝试一下吧https://blog.csdn.net/weixin_39124457/article/details/120187874
你可以看看这个网站设置的东西

@李东岳 就是单核跑也不行，现在解决了，应该是openfoam10里面的pitzDaily哪里要调一下，换成其他算例就能跑了，谢谢东哥

我目前想要计算一个后向台阶流动问题，实验中台阶上方为充分发展的湍流边界层，因此入口速度需要额外的处理。目前认为最准确的方法是单独计算一个槽道流的充分发展湍流并提取速度平面作为入口条件。然而这两个计算需要满足对应关系。Edmond Shehadi博士的博士“Large Eddy Simulation of Turbulent Flow over a Backward-Facing Step”中给出了通过后向台阶流动动量耗散雷诺数和台阶高度雷诺数得到对应槽道流参数的流程，然而根据他们的流程，我计算他们的台阶流参数得到的槽道流摩擦雷诺数为730.7073，他们最终的结果为388。我检查过公式没有问题，这说明有额外的变化在这里并没有交代，想请问各位老师和同学能否指点迷津，非常感谢。下方是论文中给出的换算流程。槽道流的半高他们说应该取为台阶流中边界层厚度δ。
屏幕截图 2024-04-20 212245.png
屏幕截图 2024-04-20 212331.png

我在做射流模拟的时候，使用了SST k-omega湍流模型，控制了y+～1，但模拟出来的边界层和自由剪切层与实验相比厚了两倍，没有使用壁面函数，这种情况下该从哪个方面着手去解决呢？湍流边界条件omega是根据文献中建议的值去给定的。

fluent 2D，重叠网格，圆柱入水，为什么这个样子呢，液面随着overset下降，很奇怪
8c2e5309-bef2-4c3f-911f-160729d98921-55d86309dbfea66064014c30d1894e1.png

烦请各位大佬指点一下，workbench里面mesh画的带角度的翼型的网格，为什么边界层厚度不一样啊？怎么处理一下呢？上面的感觉太窄了

e9d97488-fa0f-453d-9f2b-1884d5136695-image.png

root文件夹是不是没空间了？有一些系统root单独空间划片了。我看你跑在了root下。操作不是很常规

这个问题解决了，主要是因为OpenFOAM在solarCalculator.C中有出纰漏
修改如下并编译就好了，可以模拟在useExternalBeam下给定太阳辐射值（随时间变化）的效果

void Foam::solarCalculator::correctSunDirection() { if (sunDirectionModel_ == mSunDirTracking) { calculateBetaTheta(); calculateSunDirection(); //directSolarRad_ = A_/exp(B_/sin(max(beta_, ROOTVSMALL))); if (sunLoadModel_ == mSunLoadTimeDependent) { directSolarRad_ = directSolarRads_->value(mesh_.time().value()); } else if (sunLoadModel_ == mSunLoadFairWeatherConditions) { directSolarRad_ = A_/exp(B_/sin(max(beta_, ROOTVSMALL))); } } } void Foam::solarCalculator::correctDirectSolarRad() { if (sunLoadModel_ == mSunLoadTimeDependent) { directSolarRad_ = directSolarRads_->value(mesh_.time().value()); } else if (sunLoadModel_ == mSunLoadFairWeatherConditions) { directSolarRad_ = A_/exp(B_/sin(max(beta_, ROOTVSMALL))); } }

@李东岳 好的，李老师，那我自己再修改修改定义步长的代码，感谢李老师，有什么参考的资料吗qaq

半圆柱没有

@李东岳 谢谢李老师的回答，我再找找原因！

@沐风 你上传的作者联系我删除你上传的pdf文件，我将在明天将其进行删除

@知行 ok，我的网格都是ICEM画的，我加密试试，谢谢

两相界面的处理是用contour，然后切alpha=0.5的分界面

@张德聪 https://cfd-china.com/ 用这个域名试一下

那个就是发散了。

@李东岳 李老师，就是感觉有时候，不能满足通过理论计算得到的边界层信息，比如不能满足边界层厚度啊，层数也达不到，这种一般怎么处理啊？感觉workbench的mesh虽然可以，但是我希望通过这个方式画网格，这样对我后续的计算方便一点。💋

一般泰勒格林涡的文章都会指定计算域长度与涡的大小啊。比如这个：https://www.jianguoyun.com/p/DblYnzoQ9s3ZBhil274FIAA 方程6.2

把log贴出来看一下

@zhangK 客气，祝顺利！

@李东岳 好滴，谢谢李老师~

发现这个参数匹配好麻烦啊，层数、第一层高度等等之类的，感觉按照计算得到的画一下基本不出来

:136:

d8ddb624-ae67-45a5-b91b-a37014bcda10-image.png
最近用tecplot处理fluent云图，发现在应用edge时，出现了图示裂痕，怀疑是多核计算造成的。有大佬知道怎么解决吗？我目前了解到的是用单核打开，再输出tecplot文件。不过好像也有概率问题，并且算例较大且多，每次都这样很麻烦。有比较好的方法吗？

CFDEM3.8及OF5.X是用这个reconstructPar -noLagrangian，或者直接在solver中把postpro改成true

热空气从左到右，一开始假设热空气中没有水蒸气，开启多组分模型（空气+水蒸气），中间多孔介质区域设置了一个水蒸气的源项，为什么中间多孔介质区域中水蒸气质量分数比干燥后的空气中水蒸气少？为什么右边水蒸气质量分数一直没啥变化？

cbc4c76d6a01019f59da4b5821b9841.png b88edb1d22b9badd6c01ade2964c0f8.png

https://www.bilibili.com/video/BV1or421x7Ut/?spm_id_from=333.999.0.0&vd_source=5dfca494f85875fa7f031a9f632b1a89
第9分钟介绍这种High Aspect Ratio对于瞬态、稳态的具体影响。

我尝试了olaFlow 来进行造波，水池深0.5米，波高0.024,周期1.26s，水池长12米，以下是各个位置的波高数据：
x= 0.png
x= 1.png
x= 2.png
x= 3.png
x= 4.png
x= 5.png
x= 6.png
x= 7.png
x= 8.png
x= 9.png

CPU型号： Intel(R) Xeon(R) Platinum 8475B 系统：ubuntu 22.04 版本：openfoam10 内存：128G 16G 3200MHz*8 48 51 32 60 16 91 8 117 4 226 2 410 1 848

@liujunCFD 最后用matlab怎么解决的呢，我目前也遇到这个问题了，没人带，只能自己摸索，所以大佬可以告知一下方法吗？

请问各位优秀的学术者，我现在已经用fluent算完了一个三维冷态的（即没进行点火燃烧）超音速燃烧仿真，现在想要去分析计算域里某一截面上燃料与空气的混合效率，但是这个 a 的值我没办法直接获取，因为这个截面上的燃料分布是不均匀的，如果用fluent里面自带的面加权平均来算的话，这个 a 的值就比实验值小很多，所以我想请教一下在fluent或者tecplot、cfd-post这样的后处理软件怎么去正确获取a的值呢？谢谢各位

混合效率.png