请问，MPPIC这个方法的particle distribution function的运作

zhe

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

bestucan

至少把你的现有进展说一下，进行到哪一步了。particle distribution function 是代码里的还是文献里的，在代码哪里。你看了什么文献、哪个代码却没有找到。好歹让想帮你的人把你踩过的坑给避开了。

zhe

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

bestucan

@zhe DPMFoam.C 这个的108行在算粒子，下面就开始用算过的信息更新alphac了。

108行的定义在这里 KinematicClud.C 673行。这里用到的solve函数就在这个文件的91行。

然后这solve里一大堆追着函数名找，大概就是particle distribution function的过程了。我猜的

zhe

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

马乔

你要的MPPIC大概都在这里了 MPPIC

zhe

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

zhe

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

zhe

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

刘雄国

@zhe qtcreator line-by-line debuging

东岳

@zhe 在请问，MPPIC这个方法的particle distribution function的运作中说：

固相是由particle distribution function表示的，并且固相的相分数也是通过它进行求解的。那么这个方程在OpenFOAM里有体现吗？

没有，你说的这个方程是一个非常基本的控制方程，其实就是玻尔兹曼方程，很多模型都是从这个模型上展开对他进行封闭的。比如

双流体模型，就是粒子分布函数传输方程取2次矩的模型。
矩方法就是取高阶矩的模型。
MPPIC在这个方程上对碰撞进行BGK封闭。可以更简单的处理碰撞。看这个文献的方程2

A model for collisional exchangeing as /liquid/solidfluidizedbeds

《无痛苦N-S方程笔记》4,8节也有讨论。这个方程没人去解。自由度太多了。

可以彼此融会贯通左右逢源一下