@zhe 不好意思啊,之前看文献的时候前面看的快了,没留意.我的问题的答案是:在实际求解时,particle distribution function不被调用,而是通过直接对parcel进行追踪得到的.因为每个parcel中所包含的粒子具有所有相同的性质.所以用这个来计算出粒子的体积分数么的问题.希望能帮到跟我一样有疑问的人.谢谢.

@马乔 您好，不好意思，还是得麻烦您。我一直在研究MPPIC，想要把codes都捋出来。可是一直有个关卡过不去。这个在您的博客里也没提到，所以想请您帮忙解答一下。MPPIC中，固相是由particle distribution function表示的，并且固相的相分数也是通过它进行求解的。那么这个方程在OpenFOAM里有体现吗？我查看了很多遍，都没有找到哪个主程序需要调用这个方程，不知道他们是怎么用这个方程的。还有，对于固相粒子速度的修正是通过damping 和isotropy，那么这两个相是怎么被调用的啊？它们也应该在这个particle distribution function里体现啊，可是哪里有这个求解公式啊？希望您能抽时间为我解答一下，谢谢啦🙏

@东岳 稠密的.sediment transportation.

@马乔 您好,我并不是想要追踪精确的粒子运动轨迹,只是当网格饱和时,这个方法的处理方式确实会引起结果的不确定性.我再看看是否可以运用其他的方法,谢谢您!

@马乔 这个是固相应力以压力方式排斥粒子，不让它进入该网格。但是由于用的是normal stress，那么粒子返回的路径就是个问题了。我做了几个小tests，粒子都是直接被弹走，并且不能预测它去哪里。还有就是，被弹走的粒子也对原来饱和的网格内的粒子有反作用力，使得它们分散看来。不知道这个现象有没有什么办法可以规避？希望您多指导一下。

@马乔 我刚刚去看了，这就是我想弄懂的部分。太感谢了。还是您的原创！了不起！

大家好,我想问一下,在研究MPPIC这个方法时,大家应该都遇到了"固体粒子相分数极限"这个项.也就是控制网格内固体粒子的占比最大值.但是我的问题是,如果超过这个占比的粒子想要移动到这个网格里,是被阻挡了吗?是通过压力吗?那么被阻挡下来的粒子又是怎么安置的呢?是阻挡在了附近网格,还是在一个时间步长内,把原本装满粒子的网格挤出相对数量的粒子,来保证质量守恒?希望大家能多发言讨论,哪怕是猜测的,只要有道理,欢迎!!!

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

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

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