@冠竹 我也是新手，这是我个人理解。打开paraview，可以看到UPrime2Mean的XX、YY、ZZ三个对角线雷诺应力，代表了脉动速度u'、v'、w'的速度方差，这三个雷诺应力开根号对应就能得到u'、v'、w'的速度标准差。
1.png

@vbcwl 是的，

advective 边界条件有人用过嘛？求指导

@流体中二少年

forAll(crossDiff,i) { if (crossDiff[i] > 0.0) { arg1[i] = 0.125; } else { arg1[i] = 0.0; } }

这两项也可以写为：
$$
\frac{\rho T}{c_{v}} \left( \frac{\partial c_{v}}{\partial t} + \boldsymbol{U} \cdot \nabla c_{v} \right)
$$

@洱聿 好的，我先试一下。

@sibo 大佬，请教下，如何在这个buoyantPimpleFoam中添加两种流体，进而形成双扩散自然对流的情形哈？

你好，我想问一下，你成功编译interFluid了吗？有些问题想求助一下。

是的，如果不关注界面，不需要VOF

@jpzhao

看这个说明

https://www.openfoam.com/documentation/guides/latest/api/splitMeshRegions_8C_source.html

- cellZonesFileOnly behaves like -cellZonesOnly but reads the cellZones
from the specified file. This allows one to explicitly specify the region
distribution and still have multiple cellZones per region.

看看有没有能用的

如果没有我也不知道了:135:

@五好青年

DILU 预条件器出错 对应代码在

DILUPreconditioner.C 90行

for (label cell=0; cell<nCells; cell++) { 90 rDPtr[cell] = 1.0/rDPtr[cell]; }

更换预条件器 换成diagonal或者none一般就可以，或者换光顺求解器也可以。
一般不完全ILU分解做预条件器，是不能够完全分解成一个L和U矩阵的，还会剩下一个对角阵D，如果某行完全分解了，那么对角阵会有元素为0，由于这第90行代码会将D中的元素做除数，OpenFoam刚好又会捕捉浮点数异常，程序就退出了。

@zhanghan

有哪位大神知道为什么用reactingTwoPhaseEulerFoam求解器计算时不输出气相和液相速度的残差吗？

一般多相流不求解速度。单相流可可以调节，是求还是不求。多相流比较激进，索性直接不求了。没有开关。

reactingTwoPhaseEulerFoam求解器的U.Eqn中无该项，是这个原因吗？如果没solve，为什么会有速度场呢？

速度方程可以不求解。http://dyfluid.com/icoFoam.html 你可以看下这个。

太流弊了。屌屌屌。回想我10年前开始装OpenFOAM，竟然也是装了一个星期。

@cccrrryyy 好的，感谢老师

@coolhhh 感谢，我去看看

@bestucan 感谢

@pcwang 楼主思路十分清晰，祝顺利毕业:146:

@calf-z-dns 你好，我也在做DNS的管道传热，壁面是恒定热流密度，可是我的温度云图结果几乎和层流一样，可以请教下你这种情况下该怎么解决么

@李东岳 谢谢东岳老师，我这里温度设置是参考的hotroom，然后把壁面设置为了恒定热流密度，应该是没有扰动的。我先去学习下如何添加温度扰动，之后再来看看结果

@bestucan 完美解决，谢谢大佬:146: