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



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


  • Linux讲师 OpenFOAM讲师

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



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


  • Linux讲师 OpenFOAM讲师

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

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

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



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


  • 离散相副教授

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



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



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



  • @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节也有讨论。这个方程没人去解。自由度太多了。

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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2