还有很神奇的事情，就是440核的话我用20个节点，每个节点22核测是可以正常跑，但是换成22个节点每个节点20核跑也出现上述问题，不知道具体问题是出在内存不够呢还是说编译器并行方式不同呢？一个节点满核28核，上面的两个我都是节点上满核跑的，也就是504核用了18个节点，644核用了23个节点。

我没仔细看，你的代码思路上没啥问题，addparticle传入一个particletype类型的指针，就是你前面用new创建的指针，但是好像参数不太对吧，new sprayparcel（一些参数），应该不是传入this指针，你可以参考一下injection model里面入射粒子部分的写法

@星星星星晴 傻了 知道了

怎么没法写公式了 = = 是不是方法用错了

67069363-6745-40a9-8a27-6aa33e64e681-image.png
质量守恒很简单，但是根据code，child的速度和数量就没办法设置。。。

我再看看哪有方法解决吧

哈啊，喝 :xiabanle:

@东岳 在 在OF2.3.0版本如何编译低版本下的湍流模型？ 中说：

模型计算公式拔出来，自

懂了，，，谢谢东岳老师!

下面的图按照初始setField→计算一个时间步排列。
不知道为什么，经过计算后内部场的速度很快发生变化。 而内部场应该是利用势流理论估算的，但是拿N-S算出来的速度不应该跟势流估算的差这么大啊...

初始波峰在左侧
1Abefore.png
1Before.0000.png
1AAfter.png
1After.0001.png

波峰传递过来时
0A.0000.png
0.0000.png
0A.0001.png
0.0001.png

波峰在右侧
2.5A.0000.png
2.5.0000.png
2.5A.0001.png
2.5.0001.png

波峰离开
-5A.0000.png
-5.0000.png
-5A.0001.png
-5.0001.png

一阶格式更加稳定，在发散的情况下追求高阶格式会更加发散，为何不把可调节时间不长开启？

好像是因为你发散了想看连续性误差？这个没什么用。不能表现发散的本质信息

如果你想看的话，openfoam里面这个叫continuity error，有个sum，有个global，有个cumulative，可以参考sum和global

完结帖子，后来找到了，我的openmpi不是装在系统文件夹的，换一下bash文件里的export WM_MPLIB=OPENMPI就好了（如果openmpi-x.xx.x装在系统文件夹就是SYSTEMOPENMPI)

我们可以交流下的哦。我的QQ：2106936

关注点不太一样吧，我是做雾化模型的，主要是算的相对准，计算量要小，现在大佬们做雾化直接DNS，Quasi-DNS，SPH的@东岳

@Samuel-Tu 在 东岳老师CES论文中两相kEpsilon模型实现的疑问 中说：

这说明a相的粘度是靠b相的粘度计算来的，但是我没有找到相关代码，希望大家能告诉我一下文件路径，让我学习一下如何实现。

你可以参考下Gosman 1992年那篇文章

对于公式（32）中C3项（绿框部分），是否少了一个epsilon。因为源码C3项中含有epsilon：

应该是少了个epsilon

多相流中的混相问题，目前的研究问题是油气水三相，其中油为重烃、气体为CO2。CO2在一定温度压力下会与原油发生混相（CO2扩散进原油，导致原油粘度降低，最终两者充分混合形成一种流体），这个过程中要考虑流体的可压、溶解、相变等问题，这个过程应该基于现有OpenFoam中的哪个求解器进行模拟？需要改动哪些地方？流体性质随温度压力变化能否考虑？

@aiweimo 解决了！万分感谢啊！:baobao:

多谢！已更新

好像foam_extend里面已经有了Rusche这个代码？你看看？
这个算法我看过但没自己植入过。因此这些细节问题也不好说，没什么建议了

您好，我现在也在用openfoam做固液两相流，能求教您吗？

forAll循环在openfoam里面特别多，这种小操作不会减慢速度，

@东岳 感谢东岳老师，意思就是 在 P 文件中写成
inlet
type codedfixedvalue
value -ρ*（V-V'/δt）
这种格式就好么
东岳老师写的这个Δt是时间步长么，V-V(t-Δt)要怎么表述呢
感谢东岳老师，谢谢你的指导