找到原因了，是粒子入注的时候打到并行分块的交界面上了，很容易进入Cloud:move的死循环，正在思考怎么修改

不知道你的源项具体是怎样的，我在压力方程中植入过源项，基本计算都不稳定，我的源项和相体积分数，压力温度相关，耦合十分严重，后来也是重新推导方程强行引入隐式的压力项求解才稳定的。另外reactingTwoPhaseEulerFoam的算例下有injection算例，通过fvOption植入的质量入注源项，你可以看一下

converge cfd提供的理论手册，可以作为openfoam和fluent理论手册的补充。包含Euler，Largraian多相流模型，燃烧模型，湍流模型等。个人用户下载的话需要注册账号并且等待官方核实，我也等了三天，这里分享给有需求的小伙伴
链接: https://pan.baidu.com/s/1LOfuYd6OVZK7R4Pnz7wUbg 提取码: 5q54

謝謝老哥，後來用的歪門邪道搞定的，你説的這個方法不錯 @bestucan

在slice云图上作出一条黑线怎么做想做一条X=0这样的线
(

模板类函数的继承在子类里面调用需要声明用的是父类的函数，所以你直接改最底层的的那个调用的时候是找不到的

速度平均还是不会做，不过发现湍流强度对计算结果影响很大，大家算两相流的注意一下

对的，这个不好说，但是SMD对得上，就想着速度也赶紧对上就好了

各位大佬好，最近在用自己开发的求解器LES算横向射流，最后统计粒子在某条线上的粒子的平均速度，发现和实验差别还蛮大的，需要像欧拉场一样做时间平均吗？有如何处理的idea吗？

@东岳 ok，明白了，谢谢东岳老师

@东岳 Ueqn没有处理，是通量的符号出了问题吗？

@东岳 对的，按照理论来讲只是方程的系数变了一个符号，相当于
solve(pEqn==0) 和solve（-pEqn==0)

最近修改方程发现
-fvc::div(phiHbyA)+fvm::laplacian(rAUf, p_rgh)
fvc::div(phiHbyA)-fvm::laplacian(rAUf, p_rgh)
求解的稳定性完全不一样前者稳定，后者发散了，觉得很奇怪，大佬们有何高见@东岳