# Integral conservation of mass species equation

### From Thermal-FluidsPedia

Line 5: | Line 5: | ||

for the <big><math>i</math><sup>th</sup></big> species with <big><math>\Phi = {m_i}</math></big> and <big><math>\phi = {\rho _i}/\rho </math></big> (where <math>{\rho _i}</math> is the concentration of the <big><math>i</math><sup>th</sup></big> component), i.e., | for the <big><math>i</math><sup>th</sup></big> species with <big><math>\Phi = {m_i}</math></big> and <big><math>\phi = {\rho _i}/\rho </math></big> (where <math>{\rho _i}</math> is the concentration of the <big><math>i</math><sup>th</sup></big> component), i.e., | ||

- | <center><math>{\left. {\frac{{d{m_{k,i}}}}{{dt}}} \right|_{system}} = \frac{\partial }{{\partial t}}\int_V {{\rho _{k,i}}dV + \int_A {{\rho _{k,i}}({{\mathbf{V}}_{k,i,rel}} \cdot {{\mathbf{n}}_k})dA} } \qquad \qquad(1)</math></center> | + | <center><math>{\left. {\frac{{d{m_{k,i}}}}{{dt}}} \right|_{system}} = \frac{\partial }{{\partial t}}\int_V {{\rho _{k,i}}dV + \int_A {{\rho _{k,i}}({{\mathbf{V}}_{k,i,rel}} \cdot {{\mathbf{n}}_k})dA} } \qquad \qquad(1)</math></center> |

If there is no chemical reaction, the total mass of the <big><math>i</math><sup>th</sup></big> species for a closed system remains constant. Chemical reactions, on the other hand, will result in the production or consumption of the <big><math>i</math><sup>th</sup></big> species, which can be modeled as a mass source or sink for the <big><math>i</math><sup>th</sup></big> species, i.e., | If there is no chemical reaction, the total mass of the <big><math>i</math><sup>th</sup></big> species for a closed system remains constant. Chemical reactions, on the other hand, will result in the production or consumption of the <big><math>i</math><sup>th</sup></big> species, which can be modeled as a mass source or sink for the <big><math>i</math><sup>th</sup></big> species, i.e., | ||

- | <center><math>{\left. {\frac{{d{m_i}}}{{dt}}} \right|_{system}} = \int_V {{{\dot m'''}_i}dV} \qquad \qquad(2)</math></center> | + | <center><math>{\left. {\frac{{d{m_i}}}{{dt}}} \right|_{system}} = \int_V {{{\dot m'''}_i}dV} \qquad \qquad(2)</math></center> |

where <math>{m_i^{'''}}</math> is the mass production rate for the <math>i</math><sup>th</sup> species, which can be either positive (mass production) or negative (mass consumption). For the cases without chemical reaction, <math>{m_i^{'''}}=0</math>. The conservation of species mass for a control volume is therefore | where <math>{m_i^{'''}}</math> is the mass production rate for the <math>i</math><sup>th</sup> species, which can be either positive (mass production) or negative (mass consumption). For the cases without chemical reaction, <math>{m_i^{'''}}=0</math>. The conservation of species mass for a control volume is therefore | ||

- | <center><math>\frac{\partial }{{\partial t}}\int_V {{\rho _i}dV + \int_A {{\rho _i}({{\mathbf{V}}_{i,rel}} \cdot {\mathbf{n}})dA} } = \int_V {{{\dot m'''}_i}dV} \qquad \qquad(3) </math></center> | + | <center><math>\frac{\partial }{{\partial t}}\int_V {{\rho _i}dV + \int_A {{\rho _i}({{\mathbf{V}}_{i,rel}} \cdot {\mathbf{n}})dA} } = \int_V {{{\dot m'''}_i}dV} \qquad \qquad(3) </math></center> |

which is valid for each component in the control volume. If the total number of species is <math>N</math>, summation of the conservation of mass of all species results in | which is valid for each component in the control volume. If the total number of species is <math>N</math>, summation of the conservation of mass of all species results in | ||

- | <center><math>\frac{\partial }{{\partial t}}\int_V {\left[ {\sum\limits_{i = 1}^N {{\rho _i}} } \right]dV + \int_A {\left[ {\sum\limits_{i = 1}^N {({\rho _i}{{\mathbf{V}}_{i,rel}}) \cdot {\mathbf{n}}} } \right]dA} } = \int_V {\left[ {\sum\limits_{i = 1}^N {{{\dot m'''}_i}} } \right]dV} \qquad \qquad(4) </math></center> | + | <center><math>\frac{\partial }{{\partial t}}\int_V {\left[ {\sum\limits_{i = 1}^N {{\rho _i}} } \right]dV + \int_A {\left[ {\sum\limits_{i = 1}^N {({\rho _i}{{\mathbf{V}}_{i,rel}}) \cdot {\mathbf{n}}} } \right]dA} } = \int_V {\left[ {\sum\limits_{i = 1}^N {{{\dot m'''}_i}} } \right]dV} \qquad \qquad(4) </math></center> |

The summation of the densities of the individual species is equal to the bulk density of the multicomponent substance: | The summation of the densities of the individual species is equal to the bulk density of the multicomponent substance: | ||

- | <center><math> \rho = \sum\limits_{i = 1}^N {{\rho _i}} \qquad \qquad(5) </math></center> | + | <center><math> \rho = \sum\limits_{i = 1}^N {{\rho _i}} \qquad \qquad(5) </math></center> |

The bulk velocity of the multicomponent substance is the mass-averaged velocity of the velocity of each individual species: | The bulk velocity of the multicomponent substance is the mass-averaged velocity of the velocity of each individual species: | ||

- | <center><math>\rho {{\mathbf{V}}_{rel}} = \sum\limits_{i = 1}^N {({\rho _i}{{\mathbf{V}}_{i,rel}})} \qquad \qquad(6) </math></center> | + | <center><math>\rho {{\mathbf{V}}_{rel}} = \sum\limits_{i = 1}^N {({\rho _i}{{\mathbf{V}}_{i,rel}})} \qquad \qquad(6) </math></center> |

The right-hand side of eq. (4) must be zero, because the total mass of species produced is equal to the total mass of species consumed, i.e., | The right-hand side of eq. (4) must be zero, because the total mass of species produced is equal to the total mass of species consumed, i.e., | ||

- | <center><math> \sum\limits_{i = 1}^N {{{\dot m'''}_i}} = 0 \qquad \qquad(7) </math></center> | + | <center><math> \sum\limits_{i = 1}^N {{{\dot m'''}_i}} = 0 \qquad \qquad(7) </math></center> |

Substituting eqs. (5) – (7) into eq. (4), the continuity equation <math>\int_V {\left( {\frac{{\partial \rho }}{{\partial t}} + \nabla \cdot \rho {{\mathbf{V}}_{rel}}} \right)dV = 0}</math> from [[Conservation of mass species equation]] is obtained. This means that we can write conservation of species mass for <math>N</math> species, but only <math>N-1</math> of these equations are independent. In practice, one can use these <math>N-1</math> equations for conservation of species mass in combination with the continuity equation, eq. <math>\int_V {\left( {\frac{{\partial \rho }}{{\partial t}} + \nabla \cdot \rho {{\mathbf{V}}_{rel}}} \right)dV = 0}</math>, to describe multicomponent systems. | Substituting eqs. (5) – (7) into eq. (4), the continuity equation <math>\int_V {\left( {\frac{{\partial \rho }}{{\partial t}} + \nabla \cdot \rho {{\mathbf{V}}_{rel}}} \right)dV = 0}</math> from [[Conservation of mass species equation]] is obtained. This means that we can write conservation of species mass for <math>N</math> species, but only <math>N-1</math> of these equations are independent. In practice, one can use these <math>N-1</math> equations for conservation of species mass in combination with the continuity equation, eq. <math>\int_V {\left( {\frac{{\partial \rho }}{{\partial t}} + \nabla \cdot \rho {{\mathbf{V}}_{rel}}} \right)dV = 0}</math>, to describe multicomponent systems. | ||

Line 35: | Line 35: | ||

In each equation of the conservation of species mass, there are two new unknowns: mass density of the <math>i</math><sup>th</sup> species <math>{{\rho}_i}</math>, and the species mass velocity <math>{{\mathbf{V}}_{i,rel}}</math>. This yields more unknowns than the number of equations. So in order to have a properly-posed equation set, it is necessary to reduce the number of unknowns in each equation from two to one. The second term on the left-hand side of eq. (3) represents the species <math>i</math> mass flow across the surface of the control volume, which results from convection by bulk flow and diffusion relative to the bulk convection. | In each equation of the conservation of species mass, there are two new unknowns: mass density of the <math>i</math><sup>th</sup> species <math>{{\rho}_i}</math>, and the species mass velocity <math>{{\mathbf{V}}_{i,rel}}</math>. This yields more unknowns than the number of equations. So in order to have a properly-posed equation set, it is necessary to reduce the number of unknowns in each equation from two to one. The second term on the left-hand side of eq. (3) represents the species <math>i</math> mass flow across the surface of the control volume, which results from convection by bulk flow and diffusion relative to the bulk convection. | ||

- | <center><math>\int_A {{\rho _i}({{\mathbf{V}}_{i,rel}} \cdot {\mathbf{n}})dA} = \int_A {{\rho _i}({{\mathbf{V}}_{rel}} \cdot {\mathbf{n}})dA} + \int_A {{{\mathbf{J}}_i} \cdot {\mathbf{n}}dA} \qquad \qquad(8) </math></center> | + | <center><math>\int_A {{\rho _i}({{\mathbf{V}}_{i,rel}} \cdot {\mathbf{n}})dA} = \int_A {{\rho _i}({{\mathbf{V}}_{rel}} \cdot {\mathbf{n}})dA} + \int_A {{{\mathbf{J}}_i} \cdot {\mathbf{n}}dA} \qquad \qquad(8) </math></center> |

where <math>{{\mathbf{J}}_i}</math> is the diffusive mass flux vector of species <big><math>i</math></big>, which includes mass fluxes due to ordinary diffusion driven by the concentration gradient, and thermal (Soret) diffusion. | where <math>{{\mathbf{J}}_i}</math> is the diffusive mass flux vector of species <big><math>i</math></big>, which includes mass fluxes due to ordinary diffusion driven by the concentration gradient, and thermal (Soret) diffusion. | ||

Line 41: | Line 41: | ||

Substituting eq. (8) into eq. (3), one obtains an expression for the conservation of species mass that contains only one new additional variable, <big><math>{{\rho}_i}</math></big>: | Substituting eq. (8) into eq. (3), one obtains an expression for the conservation of species mass that contains only one new additional variable, <big><math>{{\rho}_i}</math></big>: | ||

- | <center><math>\frac{\partial }{{\partial t}}\int_V {{\rho _i}dV + \int_A {{\rho _i}({{\mathbf{V}}_{rel}} \cdot {\mathbf{n}})dA} } = - \int_A {{{\mathbf{J}}_i} \cdot {\mathbf{n}}dA} + \int_V {{{\dot m'''}_i}dV} \qquad \qquad(9) </math></center> | + | <center><math>\frac{\partial }{{\partial t}}\int_V {{\rho _i}dV + \int_A {{\rho _i}({{\mathbf{V}}_{rel}} \cdot {\mathbf{n}})dA} } = - \int_A {{{\mathbf{J}}_i} \cdot {\mathbf{n}}dA} + \int_V {{{\dot m'''}_i}dV} \qquad \qquad(9) </math></center> |

The above analysis is based on the assumption that the control volume contains only one phase. For a control volume containing <math>\Pi</math> phases, the conservation of species mass is [[#References|(Faghri and Zhang, 2006)]]. | The above analysis is based on the assumption that the control volume contains only one phase. For a control volume containing <math>\Pi</math> phases, the conservation of species mass is [[#References|(Faghri and Zhang, 2006)]]. | ||

Line 48: | Line 48: | ||

\sum\limits_{k = 1}^\Pi {\left[ {\frac{\partial }{{\partial t}}\int_{{V_k}(t)} {{\rho _{k,i}}dV} + \int_{{A_k}(t)} {{\rho _{k,i}}({{\mathbf{V}}_{k,rel}} \cdot {{\mathbf{n}}_k})dA} } \right]} \\ | \sum\limits_{k = 1}^\Pi {\left[ {\frac{\partial }{{\partial t}}\int_{{V_k}(t)} {{\rho _{k,i}}dV} + \int_{{A_k}(t)} {{\rho _{k,i}}({{\mathbf{V}}_{k,rel}} \cdot {{\mathbf{n}}_k})dA} } \right]} \\ | ||

= \sum\limits_{k = 1}^\Pi {\left[ { - \int_{{A_k}(t)} {({{\mathbf{J}}_{k,i}} \cdot {{\mathbf{n}}_k})dA} + \int_{{V_k}(t)} {{{\dot m'''}_{k,i}}dV} } \right]} \\ | = \sum\limits_{k = 1}^\Pi {\left[ { - \int_{{A_k}(t)} {({{\mathbf{J}}_{k,i}} \cdot {{\mathbf{n}}_k})dA} + \int_{{V_k}(t)} {{{\dot m'''}_{k,i}}dV} } \right]} \\ | ||

- | \end{array} \qquad \qquad(10) </math></center> | + | \end{array} \qquad \qquad(10) </math></center> |

==References== | ==References== |

## Revision as of 23:32, 18 November 2009

The continuity equation states that the total mass for a closed system is constant. For a system containing one phase but more than one component, the total mass of the system is composed of different species. If the concentrations of each of these species are not uniform, mass transfer occurs in a way that makes the concentrations more uniform. Therefore, it is necessary to track the individual components by applying the principle of conservation of species mass. For a system of multiple components, each component can have its own mass density and velocity. The conservation of mass for the *i*^{th} species in a single phase system for a fixed-mass system is obtained by applying eq.

for the *i*^{th} species with Φ = *m*_{i} and φ = ρ_{i} / ρ (where ρ_{i} is the concentration of the *i*^{th} component), i.e.,

If there is no chemical reaction, the total mass of the *i*^{th} species for a closed system remains constant. Chemical reactions, on the other hand, will result in the production or consumption of the *i*^{th} species, which can be modeled as a mass source or sink for the *i*^{th} species, i.e.,

where is the mass production rate for the *i*^{th} species, which can be either positive (mass production) or negative (mass consumption). For the cases without chemical reaction, . The conservation of species mass for a control volume is therefore

which is valid for each component in the control volume. If the total number of species is *N*, summation of the conservation of mass of all species results in

The summation of the densities of the individual species is equal to the bulk density of the multicomponent substance:

The bulk velocity of the multicomponent substance is the mass-averaged velocity of the velocity of each individual species:

The right-hand side of eq. (4) must be zero, because the total mass of species produced is equal to the total mass of species consumed, i.e.,

Substituting eqs. (5) – (7) into eq. (4), the continuity equation from Conservation of mass species equation is obtained. This means that we can write conservation of species mass for *N* species, but only *N* − 1 of these equations are independent. In practice, one can use these *N* − 1 equations for conservation of species mass in combination with the continuity equation, eq. , to describe multicomponent systems.

In each equation of the conservation of species mass, there are two new unknowns: mass density of the *i*^{th} species ρ_{i}, and the species mass velocity . This yields more unknowns than the number of equations. So in order to have a properly-posed equation set, it is necessary to reduce the number of unknowns in each equation from two to one. The second term on the left-hand side of eq. (3) represents the species *i* mass flow across the surface of the control volume, which results from convection by bulk flow and diffusion relative to the bulk convection.

where is the diffusive mass flux vector of species *i*, which includes mass fluxes due to ordinary diffusion driven by the concentration gradient, and thermal (Soret) diffusion.

Substituting eq. (8) into eq. (3), one obtains an expression for the conservation of species mass that contains only one new additional variable, ρ_{i}:

The above analysis is based on the assumption that the control volume contains only one phase. For a control volume containing Π phases, the conservation of species mass is (Faghri and Zhang, 2006).

## References

Faghri, A., and Zhang, Y., 2006, *Transport Phenomena in Multiphase Systems*, Burlington, MA.