# Conservation of mass species equation

### From Thermal-FluidsPedia

The integral conservation of mass species equation is

The surface integrals in the above equation may be converted to volume integrals by applying

as follows:

Substituting eqs. (1) and (2) into the integral conservation of mass species equation and considering , the entire left-hand side of the integral conservation of mass species equation is included in a single volume integral, i.e.,

The only condition that makes eq. (3) true regardless of the shape and size of the control volume is that the integrand must equal zero, i.e.,

The first term on the left-hand side is the rate of increase of mass of the species *i* per unit volume, and the second term is net rate of additions of mass of the *i*^{th} species per unit volume by convection. The terms on the right-hand side are the net rate of mass of *i*^{th} species per unit volume by diffusion, and rate of production of species *i* by chemical reaction.

Equation (4) is the equation of conservation of mass for species. If the total number of species is *N*, a total of *N* − 1 independent equations for conservation of species mass can be obtained.

After defining the mass fraction of species *i* as

eq. (4) can be rewritten as

According to the continuity equation, , the first bracket on the left-hand side of eq. (6) is zero. The second bracket on the left-hand side is the substantial derivative of the mass fraction. Therefore, the conservation of species mass in terms of mass fraction becomes

Assuming binary system of A and B, one can use Fick’s law in eq. (7) to yield

which is useful in determining the diffusion in dilute liquid solution at constant temperature and pressure. Equation (8), with , is similar to energy equation with no internal heat source and viscous dissipation:

and therefore it is used for analogy between heat and mass transfer analysis.

In the preceding discussion to develop eq. (4), the mass fraction and mass flux were used. The species equation can also be developed in term of molar concentration (or molar fraction) and molar flux. By following a similar procedure, the species equation is

where the molar flux relative to the stationary coordinate axes can be obtained from (see Introduction to Mass transfer):

i.e.,

Substituting eq. (10) into eq. (9), we have

where the first term on the left-hand side is rate of increase of mole of the species *i* per unit volume, and the second term is net rate of additions of mole of the *i*^{th} species per unit volume by convection. The terms on the right-hand side are net rate of mole of *i*^{th} species per unit volume by diffusion, and the molar rate of production of species *i* by chemical reaction.

For a binary system of components A and B with constant pressure, eq. (11) reduces to

where c is the mixture molar concentration and is molar-averaged velocity defined in Chapter 1. Equation (12) can be applied to low density gases with constant temperature and pressure.

The production rate of the *i*^{th} species, (or ), is still unknown at this point, but it can be obtained by analyzing the chemical reaction. If the number of chemical reactions taking place in the system is *N*_{c}, the mass production rate is (Kleijn, 1991; Mahajan, 1996)

where *a*_{ij} is the stoichiometric coefficient that describes the proportions of the mole numbers of reactants disappearing and mole numbers of products appearing as a result of the reaction process (see Section 2.3.3). The net reaction rate of the *j*^{th} chemical reaction is the difference between the forward and backward reactions, i.e.,

If chemical reactions take place in gas mixture, the forward and backward reaction rates are

where *N*_{r} and *N*_{p} are number of reactants and products, respectively. The two reaction rate constants, and , depend on the specific chemical reaction under consideration. These constants are related by

where *p*^{0} is the standard pressure and *K*_{j}(*T*) is the thermodynamic equilibrium constant of the *j*^{th} chemical reaction:

where is the standard Gibbs energy for the *j*^{th} chemical reaction; its value depends on the specific chemical reaction considered.

For the case in which only one chemical reaction (*N*_{c} = 1) leads to production of the *i*^{th} component, the molar production rate can be obtained by a simplified from

where the index *n* represents the order of the reaction and *k*'''_{n} is a rate constant (1/s^{n}) depends on the temperature.

As before, the mass flux, *J*_{i}, in eqs. (4) and (7) includes mass fluxes due to ordinary diffusion driven by the concentration gradient, pressure diffusion, body force diffusion, and thermal (Soret) diffusion [see eq. (1)]. For a binary mixture, the mass flux can be calculated by eq. (16).

Equation (7) is valid for the case that chemical reaction occurs in the entire volume – referred to as *homogeneous* reacting system. For the case where the chemical reaction takes place on a surface – referred to as *heterogeneous* reaction – the source term in eq. (7) will not appear, and the rate of production will be accounted as a boundary condition. More discussion about modeling of the heterogeneous reacting system can be found in Faghri and Zhang (2006).

## References

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

Faghri, A., Zhang, Y., and Howell, J. R., 2010, *Advanced Heat and Mass Transfer*, Global Digital Press, Columbia, MO.

Kleijn, C.R., 1991, “A Mathematical Model of the Hydrodynamics and Gas Phase Reaction in Silicon LPCVD in a Single Wafer Reactor,” *Journal of Electrochemical Society*, Vol. 138, pp. 2190-2200.

Mahajan, R.L., 1996, “Transport Phenomena in Chemical Vapor-Deposition Systems,” *Advances in Heat Transfer*, Academic Press, San Diego, CA.