谢谢李老师，那我是不是直接用turbulence->nut()就可以获得混合湍流运动粘度了？我是看到createFields.H定义了这样的指针，对应找到他的 incompressible::turbulenceModel 类，就发现mut和nut是一样的，对于这种多相流的，它实际用的是哪个类，在哪里体现了密度的影响？对openFOAM不太了解，谢谢老师指点
image.png

@xzdbjx 大佬，可以提供一下编译进Fluent里面的udf吗？ CFD小白，自己写的总是报错:zoule:

没事，找到错误就好。这个方法其实比较笨，因为每次你修改网格，就要重新处理一遍数据。
有的时候会麻烦一些。

祝收敛

倒数第二个图，这种提示都看不出来具体咋回事
倒数第一个图，我挂libtorch的时候用过那个，你那个挂libtorch了么

@上级 大佬，想问下cfMesh上手难吗？我在做相似的研究，想转成cfMesh试一下

@李东岳 李老师，我发现应该是dev2(T(……)）这一小段出的问题，删掉这个转置操作就能编译成功。因为我是在传热求解器的基础上进行修改的，T在createFields.H里声明了是温度场，所以是不是程序将T识别成了温度而进行不了转置操作了呢？

请教各位大佬一个边界条件的选择问题
我的研究区域是一个三维区域，长宽高分别为（6000m*6000m*600m），如下图（俯视图）所示：
08ae4b29-0f40-4ce5-b8b2-db687c3fc3fb-image.png

图中的红色点是在研究区域内的观测站点，每一个观测点位可以得到从地面到600m高度处的风速、气压、温度等数据，
请问应该选择什么边界条件来把这些观测点位的数据输入到openfoam中？

@WYing 在 fvSolution中PIMPLE算法的参数设置问题 中说：

@李东岳 谢谢李老师的回复！抱歉我的上个残差图绘制有问题，正确的应该是这样

Residuals_Ur5.png

从上到下依次是速度残差，压力残差和湍流统计量残差。

上面那个code是我随意找的，我目前使用的code是：

PIMPLE { correctPhi yes; nOuterCorrectors 2; nCorrectors 2; nNonOrthogonalCorrectors 1; turbOnFinalIterOnly false; } relaxationFactors { fields { p 0.3; } equations { "(U|k|omega)" 0.7; "(U|k|omega)Final" 1.0; } }

我主要做了以下修改：
（1）调整 nOuterCorrectors=2，结果几乎不变
（2）turbOnFinalIterOnly=false，即每次PIMPLE loop都计算湍流，结果几乎不变
（3）去掉全部的松弛因子，结果几乎不变。

在模拟中，CFL保持在0.5-1.5之间，但是和实验结果对的不好。李老师觉得还有可能是哪里出了问题呢？

请问同学，你的问题解决了吗

@李东岳 嗯对，感谢李老师！

2、经验教训
不知道说的对不对，前辈们可以建议一下，在新编译软件的时候需要注意，尤其是新加PATH的时候，编译后可能并不会直接出现问题，但是重启电脑将原有.bashrc重新source会引起错误

!感谢大神的回答！

@Shihang-Chen 哇！感谢，确实解决了！

绝对到位！嘎嘎一嘎子！

好的老师，我算完给您回复。

谢谢岳哥，成了，牛牛牛

好的，谢谢李老师，我试一下

0时间步的U，是用湍流生成方法直接生成的，可以计算得到div(U)。然后就可以计算div(phi)： U -> phi = fvc::flux(U) -> div(phi)。测试结果是，这种方式计算的div(U)和div(phi)结果基本一样的，区别是小数点后五六位。此时div(U)=div(phi)=0

那这个算法很屌。不错。

@李东岳 在 CodedFixedValue边界设置问题 中说：

http://dyfluid.com/code.html 来这查一下

好的，谢谢东岳哥

velocity.,txt
(
(0 (4 0 0))
(0.5 (4 0.0698 0))
(1 (4 0.0698 0))
(1.5 (4 0.0698 0))
(2 (4 0.0698 0))
(2.5 (4 0.0698 0))
(3 (4 0.0698 0))
(3.5 (4 0.0698 0))
(4 (4 0.1397 0))
(4.5 (4 0.2096 0))
(5 (4 0.2797 0))
(5.5 (4 0.3499 0))
(6 (4 0.4204 0))
(6.5 (4 0.4911 0))
(7 (4 0.5622 0))
(7.5 (4 0.6335 0))
(8 (4 0.7053 0))
(8.5 (4 0.7775 0))
(9 (4 0.8502 0))
(9.5 (4 0.9235 0))
(10 (4 0.9973 0))
(10.5 (4 1.0718 0))
(11 (4 0.9973 0))
(11.5 (4 0.9235 0))
(12 (4 0.8502 0))
(12.5 (4 0.7775 0))
(13 (4 0.7053 0))
(13.5 (4 0.6335 0))
(14 (4 0.5622 0))
(14.5 (4 0.4911 0))
(15 (4 0.4204 0))
(15.5 (4 0.3499 0))
(16 (4 0.2797 0))
(16.5 (4 0.2096 0))
(17 (4 0.1397 0))
(17.5 (4 0.0698 0))

)
入口设置
{
type uniformFixedValue;
uniformValue tableFile;
uniformValueCoeffs
{
file "velocity.txt";
outOfBounds repeat;
}
value $internalField;
}

流场内的模型（stl文件）已经创建好，想在模型外面创建一个矩形的外流场，类似snappyHexMesh先利用blockMesh生成矩形网格。cfmesh有这个功能吗？