体积分区

当模拟双向耦合拉格朗日多相流或离散多相 (DMP) 流时，离散相的体积将减去连续流体所占的实际体积和计算通量的实际面积。对于液膜模拟，与液膜相邻的流体的可用体积小于网格单元体积。同样，在使用物理速度公式的多孔介质建模中，一部分网格单元体积被固相占据，因此减少流体的可用体积。对于二维模型，假定网格具有单位深度（以 SI 单位表示）。

在下图中，网格单元体积 $V$ 的体积分区可以由流体体积 $V^{F}$ 、多孔固体体积 $V^{S}$ 和离散相体积（此处为拉格朗日颗粒）体积 $V^{L}$ 之和组成，每个体积都可以由下标 $i$ 所指示的各个相组成：

总网格单元体积由以下公式给出：

V = \sum_{i}^{} V_{i}^{F} + \sum_{i}^{} V_{i}^{S} + \sum_{i}^{} V_{i}^{L}

(837)

流体 $η$ 的可用体积分数则定义为所有液相所占体积与总网格单元体积之比：

η = \frac{\sum_{i}^{} V_{i}^{F}}{V}

(838)

流体 $α_{i}$ 的一个相的体积分数定义为该相所占体积与总流体域体积之比：

α_{i} = \frac{V_{i}^{F}}{\sum_{i}^{} V_{i}^{F}}

(839)

其中， $\sum_{i}^{} α_{i} = 1$

对于体积分区，通过使用无量纲流体体积分数 $α_{i}$ 和可用体积分数 $η$ 来对控制流体流方程 Eqn. (623) 至 Eqn. (625) 进行比例缩放：

连续性方程: $\frac{\partial}{\partial t} \int_{V} α_{i} η ρ d V + \oint_{A} α_{i} η ρ v \cdot d a = \int_{V} S_{u} d V + \int_{V} S_{int} d V$
(840)
动量方程: $\begin{matrix} \frac{\partial}{\partial t} \int_{V} α_{i} η ρ v d V + \oint_{A}^{} α_{i} η ρ v \otimes v \cdot d a = \end{matrix} - \oint_{A} α_{i} η p I \cdot d a + \oint_{A} α_{i} η T \cdot d a + \int_{V} α_{i} η f_{b} d V + \int_{V} s_{u} d V + \int_{V} s_{int} d V$
(841)
能量方程: $\frac{\partial}{\partial t} \int_{V} α_{i} η ρ E d V + \oint_{A}^{} α_{i} η ρH v \cdot d a = - \oint_{A}^{} α_{i} η q \cdot d a + \oint_{A}^{} α_{i} η T \cdot v d a + \int_{V} α_{i} η f_{b} \cdot v d V + \int_{V} S_{u} d V + \int_{V} S_{int} d V$
(842); 其中 $S_{u}, s_{u}$ 是用户源项， $S_{int}$ 是由建模交互引起的内部源项。

如果要对组分或被动标量求解，则相应的传输方程 (Eqn. (1874) 和 Eqn. (1902)) 也会受到体积分区的影响：

组分传输: $\frac{\partial}{\partial t} \int_{V} α_{i} η ρ Y_{i} dV + \oint_{A} α_{i} η ρ Y_{i} v \cdot d a = \oint_{A} α_{i} η [J_{i} + \frac{μ_{t}}{σ_{t}} \nabla Y_{i}] \cdot d a + \int_{V} S_{Y_{i}} d V$

(843)
被动标量传输: $\frac{\partial}{\partial t} \int_{V} α_{i} ηρ ϕ_{i} d V + \oint_{A} α_{i} ηρ ϕ_{i} v \cdot d a = \oint_{A} α_{i} η j_{i} \cdot d a + \int_{V} S_{ϕ_{i}} d V$

(844)

注	源项不按可用体积分数进行缩放，而是基于几何网格单元体积进行缩放。