我插一句：openfoam这个compressibleInterFoam的温度方程，应该是多相流里面最后一个尚未完善的地方了。温度并不是守恒变量。理想情况应该求解e方程。

@Rachel0096
找出可能出错的中间量，将这个量改为AUTO_WRITE，或者将这个量的初始值给的更合理一点

@nbyjn multiPhaseInterFoam

555c7ea9-b5fe-484b-ae06-64ec9cbd376d-image.png 这个是我想求解的方程形式，这段时间看了点相关东西，请问各位前辈，foam-extend里的blockMatrix能实现我想要的求解方式吗，然后就是曳力项应该OF里应该怎么写呢？是fvm : : Sp ( K, Ug ) - fvm : : Sp ( K, Ul )吗？希望各位大佬不吝赐教啊:mianmo:

https://www.cfd-online.com/Forums/openfoam-solving/181258-vof-model-mixture-model-ras-les-model.html

https://www.openfoam.com/documentation/guides/v2112/doc/

右上方有个小窗口，从那里搜相关关键词，会快一些定位错误

也可以编译本地的这种文档，比网页更快一些

@bestucan 非常感谢您百忙之中的回复:baobao:

有可能。非结构网格有的时候算出来是这样。这算的查的还不是特别多。纯粹是数值问题。

由于相变的存在，导致interPhaseChangeFoam的压力方程中额外添加了由相变产生的速度散度项，这一项是直接由空化模型得到的。如果我并没有使用空化模型而是使用别的方法的得到了每个时刻的alpha场，如何通过新的alpha与旧的alpha构造出压力方程中的速度散度呢

@李老师和各位老师：
LES 的 inlet 的初始条件，好像就是先运行 （同条件的） RANS，然后 把下游某处 map 给 LES 的 inlet

刚接触 LES，还有个问题要和老师们确认：
波浪的 LES，还有必要人为添加扰动吗？

谢谢

@Samuel-Tu 我在用拉格朗日方法计算颗粒的时候也出现了相似的问题，单核能算，但是并行算几步就不行了，想问一下你当初怎么解决了，你楼上的帖子说“往scalarField里面加数据的时候超过了定义的大小”我没看懂。
920e3c52-5310-4117-85d0-015244998d03-image.png

@李东岳
明白了，谢谢李老师。。我另外那问题，还要再麻烦您或其他老师帮指点下

@李东岳 好吧，谢谢大佬！

@Foamer 可以设置不同的边界，然后设置不同的润湿角？

@lys 你可以看下你的solver中有没有添加storeGlobalPositions这个函数，我对你使用的这个solver不是很了解，拉格朗日是你自己添加进去的吗

请问大家，一般对风力发电机方向的仿真一般基于哪个标准求解器？然后怎么做出过那种涡流线的结果图？

壁面的速度确实不能为0，壁面的surface是0. 从这个图来看，【猜测】是粒子直径太大了导致的。这个速度场不太正确。

@李东岳 嗷嗷，好的，谢谢李老师，我再调研一下。

看不出具体原因,只能一步一步排查,最后减小时间步长,可以运行.

@李东岳 好的好的，谢谢大家