Homogeneous model

(Difference between revisions)
 Revision as of 16:09, 7 July 2010 (view source)← Older edit Revision as of 18:28, 7 July 2010 (view source)Newer edit → Line 6: Line 6: -
$\frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}} \right)+\nabla \cdot \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}}=\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk}}}}\qquad \qquad( )$
+
$\frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}} \right)+\nabla \cdot \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}}=\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk}}}}\qquad \qquad(1)$
(2.385) (2.385) - The right-hand side of equation (2.385) must be zero because the total mass of all phases produced by phase change must equal the total mass of all phases consumed by phase change. Considering this fact and eq. (2.330), the continuity equation becomes + The right-hand side of equation (1) must be zero because the total mass of all phases produced by phase change must equal the total mass of all phases consumed by phase change. Considering this fact and eq. (2.330), the continuity equation becomes -
$\frac{\partial \left\langle \rho \right\rangle }{\partial t}+\nabla \cdot \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}}=0\qquad \qquad( )$
+
$\frac{\partial \left\langle \rho \right\rangle }{\partial t}+\nabla \cdot \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}}=0\qquad \qquad(2)$
(2.386) (2.386) Line 20: Line 20: -
$\mathbf{\tilde{V}}=\frac{1}{\left\langle \rho \right\rangle }\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}}\qquad \qquad( )$
+
$\mathbf{\tilde{V}}=\frac{1}{\left\langle \rho \right\rangle }\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}}\qquad \qquad(3)$
(2.387) (2.387) - Substituting eq. (2.387) into eq. (2.386), the final form of the continuity equation for a multiphase mixture is + Substituting eq. (3) into eq. (2), the final form of the continuity equation for a multiphase mixture is -
$\frac{\partial \left\langle \rho \right\rangle }{\partial t}+\nabla \cdot \left( \left\langle \rho \right\rangle \mathbf{\tilde{V}} \right)=0\qquad \qquad( )$
+
$\frac{\partial \left\langle \rho \right\rangle }{\partial t}+\nabla \cdot \left( \left\langle \rho \right\rangle \mathbf{\tilde{V}} \right)=0\qquad \qquad(4)$
(2.388) (2.388) - It can be seen that eq. (2.388) has the same form as the local continuity equation (2.51), except that the volume-averaged density and velocity are used in eq. (2.388), where $\left\langle \rho \right\rangle =\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}.$ + It can be seen that eq. (4) has the same form as the local continuity equation (2.51), except that the volume-averaged density and velocity are used in eq. (4), where $\left\langle \rho \right\rangle =\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}.$ Line 42: Line 42: & \frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right) \\ & \frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right) \\ & =\nabla \cdot \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{{\mathbf{{\tau }'}}}_{k}} \right\rangle }^{k}}}+\sum\limits_{k=1}^{\Pi }{\left( {{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}} \right)}\mathbf{X}+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{\left( \left\langle {{\mathbf{F}}_{jk}} \right\rangle +\left\langle {{{{\dot{m}}'''}}_{jk}} \right\rangle {{\left\langle {{\mathbf{V}}_{k,I}} \right\rangle }^{k}} \right)}} \\ & =\nabla \cdot \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{{\mathbf{{\tau }'}}}_{k}} \right\rangle }^{k}}}+\sum\limits_{k=1}^{\Pi }{\left( {{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}} \right)}\mathbf{X}+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{\left( \left\langle {{\mathbf{F}}_{jk}} \right\rangle +\left\langle {{{{\dot{m}}'''}}_{jk}} \right\rangle {{\left\langle {{\mathbf{V}}_{k,I}} \right\rangle }^{k}} \right)}} \\ - \end{align}\qquad \qquad( ) [/itex] + \end{align}\qquad \qquad(5) [/itex] (2.389) (2.389) Line 49: Line 49: -
$\left\langle {\mathbf{{\tau }'}} \right\rangle =\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\mathbf{\tau }}_{k}} \right\rangle }^{k}}}=-\left\langle p \right\rangle \mathbf{I}+\mu \left[ \nabla \mathbf{\tilde{V}}+\nabla {{{\mathbf{\tilde{V}}}}_{{}}}^{T} \right]-\frac{2}{3}\mu (\nabla \cdot \mathbf{\tilde{V}})\mathbf{I}\qquad \qquad( )$
+
$\left\langle {\mathbf{{\tau }'}} \right\rangle =\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\mathbf{\tau }}_{k}} \right\rangle }^{k}}}=-\left\langle p \right\rangle \mathbf{I}+\mu \left[ \nabla \mathbf{\tilde{V}}+\nabla {{{\mathbf{\tilde{V}}}}_{{}}}^{T} \right]-\frac{2}{3}\mu (\nabla \cdot \mathbf{\tilde{V}})\mathbf{I}\qquad \qquad(6)$
(2.390) (2.390) - The summation of all interphase forces must be zero since $\left\langle {{\mathbf{F}}_{jk}} \right\rangle =-\left\langle {{\mathbf{F}}_{kj}} \right\rangle + The summation of all interphase forces must be zero since [itex]\left\langle {{\mathbf{F}}_{jk}} \right\rangle =-\left\langle {{\mathbf{F}}_{kj}} \right\rangle$, i.e., - [/itex], i.e., + -
$\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{}}\left\langle {{\mathbf{F}}_{jk}} \right\rangle }=0\qquad \qquad( )$
+
$\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{}}\left\langle {{\mathbf{F}}_{jk}} \right\rangle }=0\qquad \qquad(7)$
(2.391) (2.391) - Considering eqs. (2.387), (2.390) and (2.391), the momentum equation becomes + Considering eqs. (3), (6) and (7), the momentum equation becomes -
$\frac{\partial }{\partial t}\left( \left\langle \rho \right\rangle \mathbf{\tilde{V}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right)=\nabla \cdot \left\langle {\mathbf{{\tau }'}} \right\rangle +\left\langle \rho \right\rangle \mathbf{X}+{{\mathbf{{\dot{M}}'''}}_{I}}\qquad \qquad( )$
+
$\frac{\partial }{\partial t}\left( \left\langle \rho \right\rangle \mathbf{\tilde{V}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right)=\nabla \cdot \left\langle {\mathbf{{\tau }'}} \right\rangle +\left\langle \rho \right\rangle \mathbf{X}+{{\mathbf{{\dot{M}}'''}}_{I}}\qquad \qquad(8)$
(2.392) (2.392) Line 71: Line 70: -
${{\mathbf{{\dot{M}}'''}}_{I}}=\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{}}\left\langle {{{{\dot{m}}'''}}_{jk}} \right\rangle }{{\left\langle {{\mathbf{V}}_{k,I}} \right\rangle }^{k}}\qquad \qquad( )$
+
${{\mathbf{{\dot{M}}'''}}_{I}}=\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{}}\left\langle {{{{\dot{m}}'''}}_{jk}} \right\rangle }{{\left\langle {{\mathbf{V}}_{k,I}} \right\rangle }^{k}}\qquad \qquad(9)$
(2.393) (2.393) - Equation (2.393) represents the momentum production rate due to interaction between different phases along their separating interfaces. It must be specified according to the combination of phases in the multiphase system that is under consideration. + Equation (9) represents the momentum production rate due to interaction between different phases along their separating interfaces. It must be specified according to the combination of phases in the multiphase system that is under consideration. ====Energy Equation==== ====Energy Equation==== Line 85: Line 84: & \frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{h}_{k}} \right\rangle }^{k}}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{h}_{k}} \right\rangle }^{k}}} \right)=-\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{\left\langle {{{\mathbf{{q}''}}}_{k}} \right\rangle } \right)+\sum\limits_{k=1}^{\Pi }{\left\langle {{{{q}'''}}_{k}} \right\rangle } \\ & \frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{h}_{k}} \right\rangle }^{k}}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{h}_{k}} \right\rangle }^{k}}} \right)=-\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{\left\langle {{{\mathbf{{q}''}}}_{k}} \right\rangle } \right)+\sum\limits_{k=1}^{\Pi }{\left\langle {{{{q}'''}}_{k}} \right\rangle } \\ & +\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}\frac{D{{\left\langle {{p}_{k}} \right\rangle }^{k}}}{Dt}}+\sum\limits_{k=1}^{\Pi }{\nabla {{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}:{{\varepsilon }_{k}}{{\left\langle {{\mathbf{\tau }}_{k}} \right\rangle }^{k}}}+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{\left\langle {{{{q}'''}}_{jk}} \right\rangle }}+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk}}{{\left\langle {{h}_{k,I}} \right\rangle }^{k}}}} \\ & +\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}\frac{D{{\left\langle {{p}_{k}} \right\rangle }^{k}}}{Dt}}+\sum\limits_{k=1}^{\Pi }{\nabla {{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}:{{\varepsilon }_{k}}{{\left\langle {{\mathbf{\tau }}_{k}} \right\rangle }^{k}}}+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{\left\langle {{{{q}'''}}_{jk}} \right\rangle }}+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk}}{{\left\langle {{h}_{k,I}} \right\rangle }^{k}}}} \\ - \end{align}\qquad \qquad( ) [/itex] + \end{align}\qquad \qquad(10) [/itex] (2.394) (2.394) Line 91: Line 90: -
$\tilde{h}=\frac{1}{\left\langle \rho \right\rangle }\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{h}_{k}} \right\rangle }^{k}}}\qquad \qquad( )$
+
$\tilde{h}=\frac{1}{\left\langle \rho \right\rangle }\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{h}_{k}} \right\rangle }^{k}}}\qquad \qquad(11)$
(2.395) (2.395) - The fifth term on the right-hand side of eq. (2.394) is for summation of all interphase heat transfer and it must be zero. The last term on the right-hand side of eq. (2.394) accounts for contribution of interphase phase change energy flux due to phase change; it can be defined as + The fifth term on the right-hand side of eq. (10) is for summation of all interphase heat transfer and it must be zero. The last term on the right-hand side of eq. (10) accounts for contribution of interphase phase change energy flux due to phase change; it can be defined as -
$\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk}}{{\left\langle {{h}_{k,I}} \right\rangle }^{k}}}}={{{q}'''}_{I}}\qquad \qquad( )$
+
$\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk}}{{\left\langle {{h}_{k,I}} \right\rangle }^{k}}}}={{{q}'''}_{I}}\qquad \qquad(12)$
(2.396) (2.396) It is usually not zero although $\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{\dot{m}}'''}}}=0$. It is usually not zero although $\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{\dot{m}}'''}}}=0$. - Considering eqs. (2.395) and (2.396), the energy equation (2.394) becomes + Considering eqs. (11) and (12), the energy equation (10) becomes -
$\frac{\partial }{\partial t}\left( \left\langle \rho \right\rangle \tilde{h} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{h}_{k}} \right\rangle }^{k}}} \right)=-\nabla \cdot \left\langle {\mathbf{{q}''}} \right\rangle +\frac{D\left\langle p \right\rangle }{Dt}+\left\langle {{q}'''} \right\rangle +\nabla \mathbf{\tilde{V}}:\left\langle \mathbf{\tau } \right\rangle +{{{q}'''}_{I}}\qquad \qquad( )$
+ - (2.397) +
$\frac{\partial }{\partial t}\left( \left\langle \rho \right\rangle \tilde{h} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{h}_{k}} \right\rangle }^{k}}} \right)=-\nabla \cdot \left\langle {\mathbf{{q}''}} \right\rangle +\frac{D\left\langle p \right\rangle }{Dt}+\left\langle {{q}'''} \right\rangle +\nabla \mathbf{\tilde{V}}:\left\langle \mathbf{\tau } \right\rangle +{{{q}'''}_{I}}\qquad \qquad(13)$
(2.397) ====Species==== ====Species==== Line 116: Line 115: & \frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k,i}} \right\rangle }^{k}}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k,i}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right) \\ & \frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k,i}} \right\rangle }^{k}}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k,i}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right) \\ & =-\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{\left\langle {{\mathbf{J}}_{k,i}} \right\rangle } \right)+\sum\limits_{k=1}^{\Pi }{\left\langle {{{{\dot{m}}'''}}_{k,i}} \right\rangle }+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk,i}}}} \\ & =-\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{\left\langle {{\mathbf{J}}_{k,i}} \right\rangle } \right)+\sum\limits_{k=1}^{\Pi }{\left\langle {{{{\dot{m}}'''}}_{k,i}} \right\rangle }+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk,i}}}} \\ - \end{align}\qquad \qquad( ) [/itex] + \end{align}\qquad \qquad(14) [/itex] (2.398) (2.398) Line 122: Line 121: -
$\left\langle {{\rho }_{i}} \right\rangle =\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k,i}} \right\rangle }^{k}}}\qquad \qquad( )$
+
$\left\langle {{\rho }_{i}} \right\rangle =\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k,i}} \right\rangle }^{k}}}\qquad \qquad(15)$
(2.399) (2.399) Line 129: Line 128: -
$\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk,i}}}}=0\qquad \qquad( )$
+
$\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk,i}}}}=0\qquad \qquad(16)$
(2.400) (2.400) - Substituting eqs. (2.399) and (2.400) into eq. (2.398), and using the mass-averaged velocity defined in eq. (2.387), the conservation of species mass becomes + Substituting eqs. (15) and (16) into eq. (14), and using the mass-averaged velocity defined in eq. (3), the conservation of species mass becomes -
$\frac{\partial \left\langle {{\rho }_{i}} \right\rangle }{\partial t}+\nabla \cdot \left\langle {{\rho }_{i}} \right\rangle \mathbf{\tilde{V}}=-\nabla \cdot {{\mathbf{J}}_{i}}+{{{\dot{m}}'''}_{i}}\qquad \qquad( )$
+
$\frac{\partial \left\langle {{\rho }_{i}} \right\rangle }{\partial t}+\nabla \cdot \left\langle {{\rho }_{i}} \right\rangle \mathbf{\tilde{V}}=-\nabla \cdot {{\mathbf{J}}_{i}}+{{{\dot{m}}'''}_{i}}\qquad \qquad(17)$
(2.401) (2.401)

Revision as of 18:28, 7 July 2010

The multi-fluid model presented above is obtained by performing phase averaging as defined in eqs. (2.316) and (2.317). If spatial averaging is performed for all phases within a multiphase control volume, the homogeneous (or mixture) model can be obtained. The relationship between volume averaging and phase averaging is given in eq. (2.319), which indicates that the homogeneous model can be obtained by summing the individual phase equations of the multi-fluid model.

Contents

Continuity Equation

The continuity equation for phase k in the multifluid model is expressed by eq. (2.332). Summing the continuity equations for all Π phases together, one obtains

$\frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}} \right)+\nabla \cdot \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}}=\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk}}}}\qquad \qquad(1)$

(2.385)

The right-hand side of equation (1) must be zero because the total mass of all phases produced by phase change must equal the total mass of all phases consumed by phase change. Considering this fact and eq. (2.330), the continuity equation becomes

$\frac{\partial \left\langle \rho \right\rangle }{\partial t}+\nabla \cdot \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}}=0\qquad \qquad(2)$

(2.386)

The bulk velocity of the multiphase mixture is the mass-averaged velocity of all the individual phases:

$\mathbf{\tilde{V}}=\frac{1}{\left\langle \rho \right\rangle }\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}}\qquad \qquad(3)$

(2.387)

Substituting eq. (3) into eq. (2), the final form of the continuity equation for a multiphase mixture is

$\frac{\partial \left\langle \rho \right\rangle }{\partial t}+\nabla \cdot \left( \left\langle \rho \right\rangle \mathbf{\tilde{V}} \right)=0\qquad \qquad(4)$

(2.388)

It can be seen that eq. (4) has the same form as the local continuity equation (2.51), except that the volume-averaged density and velocity are used in eq. (4), where $\left\langle \rho \right\rangle =\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}.$

Momentum Equation

The momentum equation for phase k in the multi-fluid model is expressed in eq. (2.334). By adding together the momentum equations for all Π phases, one obtains

\begin{align} & \frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right) \\ & =\nabla \cdot \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{{\mathbf{{\tau }'}}}_{k}} \right\rangle }^{k}}}+\sum\limits_{k=1}^{\Pi }{\left( {{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}} \right)}\mathbf{X}+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{\left( \left\langle {{\mathbf{F}}_{jk}} \right\rangle +\left\langle {{{{\dot{m}}'''}}_{jk}} \right\rangle {{\left\langle {{\mathbf{V}}_{k,I}} \right\rangle }^{k}} \right)}} \\ \end{align}\qquad \qquad(5)

(2.389)

The stress tensor of the multiphase mixture is

$\left\langle {\mathbf{{\tau }'}} \right\rangle =\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\mathbf{\tau }}_{k}} \right\rangle }^{k}}}=-\left\langle p \right\rangle \mathbf{I}+\mu \left[ \nabla \mathbf{\tilde{V}}+\nabla {{{\mathbf{\tilde{V}}}}_{{}}}^{T} \right]-\frac{2}{3}\mu (\nabla \cdot \mathbf{\tilde{V}})\mathbf{I}\qquad \qquad(6)$

(2.390)

The summation of all interphase forces must be zero since $\left\langle {{\mathbf{F}}_{jk}} \right\rangle =-\left\langle {{\mathbf{F}}_{kj}} \right\rangle$, i.e.,

$\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{}}\left\langle {{\mathbf{F}}_{jk}} \right\rangle }=0\qquad \qquad(7)$

(2.391)

Considering eqs. (3), (6) and (7), the momentum equation becomes

$\frac{\partial }{\partial t}\left( \left\langle \rho \right\rangle \mathbf{\tilde{V}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right)=\nabla \cdot \left\langle {\mathbf{{\tau }'}} \right\rangle +\left\langle \rho \right\rangle \mathbf{X}+{{\mathbf{{\dot{M}}'''}}_{I}}\qquad \qquad(8)$

(2.392)

where

${{\mathbf{{\dot{M}}'''}}_{I}}=\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{}}\left\langle {{{{\dot{m}}'''}}_{jk}} \right\rangle }{{\left\langle {{\mathbf{V}}_{k,I}} \right\rangle }^{k}}\qquad \qquad(9)$

(2.393)

Equation (9) represents the momentum production rate due to interaction between different phases along their separating interfaces. It must be specified according to the combination of phases in the multiphase system that is under consideration.

Energy Equation

By summing the energy equations for all Π phases in the multifluid model, eq. (2.341), one obtains

\begin{align} & \frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{h}_{k}} \right\rangle }^{k}}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{h}_{k}} \right\rangle }^{k}}} \right)=-\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{\left\langle {{{\mathbf{{q}''}}}_{k}} \right\rangle } \right)+\sum\limits_{k=1}^{\Pi }{\left\langle {{{{q}'''}}_{k}} \right\rangle } \\ & +\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}\frac{D{{\left\langle {{p}_{k}} \right\rangle }^{k}}}{Dt}}+\sum\limits_{k=1}^{\Pi }{\nabla {{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}:{{\varepsilon }_{k}}{{\left\langle {{\mathbf{\tau }}_{k}} \right\rangle }^{k}}}+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{\left\langle {{{{q}'''}}_{jk}} \right\rangle }}+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk}}{{\left\langle {{h}_{k,I}} \right\rangle }^{k}}}} \\ \end{align}\qquad \qquad(10)

(2.394)

The mass average enthalpy of the multiphase mixture is

$\tilde{h}=\frac{1}{\left\langle \rho \right\rangle }\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{h}_{k}} \right\rangle }^{k}}}\qquad \qquad(11)$

(2.395)

The fifth term on the right-hand side of eq. (10) is for summation of all interphase heat transfer and it must be zero. The last term on the right-hand side of eq. (10) accounts for contribution of interphase phase change energy flux due to phase change; it can be defined as

$\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk}}{{\left\langle {{h}_{k,I}} \right\rangle }^{k}}}}={{{q}'''}_{I}}\qquad \qquad(12)$

(2.396)

It is usually not zero although $\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{\dot{m}}'''}}}=0$.

Considering eqs. (11) and (12), the energy equation (10) becomes

$\frac{\partial }{\partial t}\left( \left\langle \rho \right\rangle \tilde{h} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}}{{h}_{k}} \right\rangle }^{k}}} \right)=-\nabla \cdot \left\langle {\mathbf{{q}''}} \right\rangle +\frac{D\left\langle p \right\rangle }{Dt}+\left\langle {{q}'''} \right\rangle +\nabla \mathbf{\tilde{V}}:\left\langle \mathbf{\tau } \right\rangle +{{{q}'''}_{I}}\qquad \qquad(13)$
(2.397)

Species

Summing the equations for conservation of species mass, eq. (2.349), for all phases yields

\begin{align} & \frac{\partial }{\partial t}\left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k,i}} \right\rangle }^{k}}} \right)+\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k,i}} \right\rangle }^{k}}{{\left\langle {{\mathbf{V}}_{k}} \right\rangle }^{k}}} \right) \\ & =-\nabla \cdot \left( \sum\limits_{k=1}^{\Pi }{\left\langle {{\mathbf{J}}_{k,i}} \right\rangle } \right)+\sum\limits_{k=1}^{\Pi }{\left\langle {{{{\dot{m}}'''}}_{k,i}} \right\rangle }+\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk,i}}}} \\ \end{align}\qquad \qquad(14)

(2.398)

By applying eq. (2.330) to the mass density of the ith component, one obtains

$\left\langle {{\rho }_{i}} \right\rangle =\sum\limits_{k=1}^{\Pi }{{{\varepsilon }_{k}}{{\left\langle {{\rho }_{k,i}} \right\rangle }^{k}}}\qquad \qquad(15)$

(2.399)

In accordance with the conservation of mass, the mass source (or sink) of the ith component due to phase change in all phases must add up to zero, i.e.,

$\sum\limits_{k=1}^{\Pi }{\sum\limits_{j=1(j\ne k)}^{\Pi }{{{{{\dot{m}}'''}}_{jk,i}}}}=0\qquad \qquad(16)$

(2.400)

Substituting eqs. (15) and (16) into eq. (14), and using the mass-averaged velocity defined in eq. (3), the conservation of species mass becomes

$\frac{\partial \left\langle {{\rho }_{i}} \right\rangle }{\partial t}+\nabla \cdot \left\langle {{\rho }_{i}} \right\rangle \mathbf{\tilde{V}}=-\nabla \cdot {{\mathbf{J}}_{i}}+{{{\dot{m}}'''}_{i}}\qquad \qquad(17)$

(2.401)