关于openfoam中使用的双流体模型

[anubis]

在看《Computational Fluid Dynamics of Dispersed Two-Phase Flows at High Phase Fractions》的时候碰到一个不太理解的地方。
一般的相间动量传递可以写成方程（3.10）

为了模拟phase inversion， 方程（3.11）加入了b相变成离散相时的曳力、升力等。
但对于方程(3.11)的曳力相，如果仍然是a相为离散相的情况，$\alpha_a$→0 ，$ f_a$→0，
那曳力相前面的部分不就变成$3/4\alpha_a\alpha_b(f_b\frac{C_{db}\rho_a}{d_b})$，和（3.10）里的不一样，b相变成分散相时候应该才是这个形式，虚拟质量力也是一样的情况。 升力也对不上。