multiphaseEulerFoam，颗粒分选，动理学模型

[xuyan]

@李东岳 感谢老师的建议！目前做了三组模拟，大颗粒的分压系数Bm分别选用1、5和10。效果如下图，增大分压系数会影响堆积形态，但对于颗粒分选现象并没有太大的影响。

[李东岳]

这一项长什么样？

[xuyan]

@李东岳 这个是摩擦压力的计算公式

[李东岳]

你植入的思路是什么，放在动量方程了么？

这个压力为什么会把大颗粒弄到上面，小颗粒弄到下面。从方程的角度怎么理解这个机制？

[xuyan]

@李东岳 这个力放在了动量方程中，具体修改直接在颗粒动理学模块摩擦压力返回值中加上了这个分压系数。图为动量方程的表达式。目前的理解是小颗粒承担的摩擦压力更小，坍塌速度更快，更早达到稳定的堆积形态，大颗粒坍塌速度更慢，在小颗粒上方坍塌，并且达到更远的坍塌距离。
论文中的模拟结果好像也表现为这种形式，小颗粒先坍塌聚集，大颗粒在小颗粒上方坍塌

[李东岳]

你那个$p_m^f$比较长，目前我也理解不上去这个的梯度跟$\bfU$的关系。$p_m^f$有简化模型么？

[xuyan]

@李东岳 老师目前用的这个公式是找的几个里面最简单的公式，近期我在读文献看看其他人的方法，有更简化的模型或者其他进展再和您交流。之前找了分压系数Bm的几个不同的表达式

[李东岳]

最好能简化一下。比如湍流分散力的模型都带$\nabla\alpha$，这就很好理解，alpha梯度大，湍流分散力大，可以把alpha的梯度抹平。你那个摩擦压力目前看不出来怎么解释颗粒分选。我对这个摩擦压力也不熟悉。那你再看看。

[xuyan]

@李东岳 嗯嗯，谢谢老师！

[李东岳]

后来这个怎么样了