@流体中二少年

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:

不完全一致，某些关键指标一致就行。

https://ntnuopen.ntnu.no/ntnu-xmlui/bitstream/handle/11250/2496219/18290_FULLTEXT.pdf?sequence=1

@cresendo 哇 我还不知道 但是我在用v2106 那我去查查吧看这个版本有没有类似的 感谢感谢~

把这两个东西混合一下:yes:

https://mdpi-res.com/d_attachment/water/water-13-03225/article_deploy/water-13-03225-v2.pdf?version=1637055877

https://sv.20file.org/up1/916_0.pdf