@xiaolin 我是因爲安裝了不同版本的OpenFOAM而導致的問題.比如說OpenFOAM6和OpenFOAM18.06,這種不同編號的版本是不同組織創建的.所以裏面有一些東西不兼容.你需要下載相應的版本.

@bestucan 谢谢你的回复，有所帮助的，不是猜的，您谦虚了啊。不过，我这个方法是比较特殊的MPPIC，跟普通的粒子追踪有些差别。我又翻了两天，还是发现，那个文献里提到的particle distribution function是不被OpenFOAM进行求解的。只是一个描述方程。欢迎您多多留言啊。很有帮助！

@bestucan 不好意思，我没表达清楚。我看到的就是文献里多次出现的MPPIC对粒子处理方程。我是想问大家，有没有人在OpenFOAM里见过这个方程，我真的找了很久，都没发现。这个方程是专门应用在MPPIC方法里对粒子状态和性能的描述的。但是我没有找到这个方程的源文件，所以不知道它是怎么被求解的。希望能得到指点。谢谢！

@bestucan 谢谢，已经说的很好了。我之前想过rho=(alpha1rho1+alpha2rho2)/(alpha1+alpha2).这个我觉得是比较合适的。虽然这样定义，直接就说明了alphac=alpha1+alpha2，而不是从后面的1-theta来的，但是我觉得这不会引起什么问题，因为宏观上来说，alphac就是等于alpha1+alpha2. 这样也满足了alpha1+alpha2+theta=alphac+theta=1. 都是归一后的了。我也试了自带的算例。计算出来的结果与之前定义rho=alpha1rho1+alpha2rho2有差别。因为没有真实实验数据可以对比，所以不敢说哪个就是对的，但是后改的这个，更能理解。单独定义rho后，我也就可以重新考虑alphac*rho这个项。我还是觉得对流项中是要考虑alphac的。希望，有朋友可以指出错误，谢谢！

@bestucan 不好意思,我確實沒看懂,爲什麼一開始就說了alpha1和alpha2是未歸一化的值.

大家好，请问有没有人知道MPPIC这个方法的particle distribution function具体怎么运作的?我查找了很多文献，也在OpenFOAM上看了好久的codes，但是没有说要解这个方程。那么这个方程到底在这个方法里的贡献提现在哪里啊？希望有人可以解答，谢谢！

@东岳 好的,多謝回覆.我會接着再試試重新定義rho,等您有時間了,希望能多指導.

@zhe 不好意思,打錯了.最後的對流項應該是div(alphaRhoPhic, U).

考慮固體粒子也佔有一部分體積,並且將氣體和液體聯合求解,我推了一下氣液混合的密度公式,並不是rho=rho1*alpha1+rho2*alpha2.而是alphac*rho=alpha1*rho1+alpha2*rho2. 所以在原來creatfile.H裏的定義應該是有問題的,請查證.那也就是說,凡是涉及到通過密度求質量的公式,在MPPICInterFoam裏都應該用alphac*rho. 那麼在對流項,我覺得應該改成div(alphacRho, U). 有錯的話,希望指出,謝謝.

@yfclark 謝謝回覆.那麼高濃度粒子的情況呢?像MPPIC可以允許粒子佔體積65%,這樣的情況是不是會有很大的影響了?