@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:

@马乔 有没有什么办法输出这个误差呢，一般大家是怎么判断pisoFoam或者pimpleFoam是否收敛呢?
我发现误差降到１e-5以下就很困难了。。。所以我是设置的监测点，发现关注的变量变化不大了就认为收敛。。

你这是什么版本啊，为嘛那一通运算没有量纲检查。forAll这个只会遍历内场。

看论文发现有作者说fluent稳态计算残差降到1e-7则认为收敛停止计算，但是我用of的simpleFoam的时候发现残差根本降不到那么低，连降到１e-６都已经是极限了。我都不敢在论文里写残差判断标准了，一般你们怎么说收敛标准呢

@bestucan 好滴，我试一下

@xpqiu 好滴

@Samuel-Tu 用了，确实快，但是有时候分网不成功，而且网格质量很难保证

本来想说试试加深水槽。 静水流深嘛，一看是浅水方程……
如果边界条件设定为带反射的，那就应该不停波动。如果是自由边界，不停波动不正常

在cfd online 里逛了一下，有些收货。
似乎non-orthogonal corrector 是在组装矩阵Ax=b的时候用的，它会把扩散项中的正交部分放入系数矩阵Ａ中，而把非正交部分放入方程右边的源项b中，使得系数矩阵更稀疏且主对角线控制，使得计算更稳定？？cfd online

@cccrrryyy 受教受教，授我以渔，万分感谢~ 帖子链接里开篇就阐明了我的疑惑，也难怪k-epsilon可以在工程领域运用的如此广泛~ 我就先对照着openfoam tutorials和CFDEM案例里的k-epsilon取值作参考，计算后再做调整，再次感谢指点~

@东岳 非常感谢老师:xinxin: ！ :xiexie:
就是这样迂回的话好像有些难，我还要再思考思考，再尝试尝试 :144:

@东岳 好滴

Already solved.

@bestucan Thanks for your attention.

谢谢老师

@东岳 这个当时手一抖就点上了，发现没法删掉，真是尴尬:mihu:

@莫晓柒 编译之后得到的库和求解器程序都在 platforms 下面，所以只能重新安装了。

@东岳 我不算运营人员，只是参与了号主的一些活动，有时候帮帮忙。。。

@upc_ngh
1.是的，$\mathbf{U}$代表的是KinematicParcel中的$\mathbf{U_c}$，也就是连续相的在颗粒所处位置的速度插值，具体的插值方式OF里面有cell，cellpoint那些，$\mathbf{U}_\rm{d}$表示的就是颗粒的速度。
2.是的，系数乘以速度差。

新知识点