计算科学 - 了解 OpenFOAM 中的 MP-PIC 实现 - 吾爱随笔录

多相粒子单元(MP-PIC) 方法的特点是将粒子属性从拉格朗日坐标映射到欧拉网格。然而，这种方法在 OpenFOAM 中的实现仅在一个“子模型”中求解欧拉网格上的插值方程，即隐式包装模型，对此有替代方案（无或显式包装模型）。这怎么可能？为什么它仍然被认为是 MP-PIC？

此外，OpenFOAM 似乎解决了这个方程以生成粒子速度校正：其中是粒子相体积分数（在欧拉网格上平均），是粒子间应力。但是，我在下面列出的两篇论文中都找不到这个等式。如何推导出方程，以及 OpenFOAM 中使用的速度校正？

\frac{\partial α}{\partial t} - Δ (τ^{'} (α) α) = 0,

$\frac{\partial\alpha}{\partial t} - \Delta\left(\tau^\prime(\alpha) \alpha\right) = 0,$

α

$\alpha$

τ

$\tau$

参考：

Andrews, MJ 和 O'Rourke, PJ (1996)。密集粒子流的多相细胞内粒子 (MP-PIC) 方法。国际多相流杂志，22（2）：379–402。
斯奈德，DM (2001)。用于致密粒子流的不可压缩三维多相粒子单元模型。计算物理学杂志，170：523–549。