@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有这个功能吗？

好的
我试试

@dxl 您好，在你的 https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/47 这篇文章中我也遇到了相同的情况可以讨论一下吗

@学流体的小明 感谢感谢:146: :146:

如题，在一个三维计算域中，均匀提取空间中各个位置的数据，除了使用probesDict文件添加监测点外，有没有其他方法实现？

@李东岳 谢谢李老师的回复，期待:xiexie:

多谢老师！

谢谢啦