Interfacial mass, momentum, energy, and species balances
From ThermalFluidsPedia
(10 intermediate revisions not shown)  
Line 5:  Line 5:  
At a liquidvapor interface, the mass balance is  At a liquidvapor interface, the mass balance is  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{{\dot{m}}''}_{\delta }}={{\rho }_{\ell }}({{\mathbf{V}}_{\ell }}{{\mathbf{V}}_{I}})\cdot \mathbf{n}={{\rho }_{v}}({{\mathbf{V}}_{v}}{{\mathbf{V}}_{I}})\cdot \mathbf{n}</math>  
+  </center>  
+  {{EquationRef(1)}}  
}  }  
  where <math>{{{  +  where <math>{{{\dot{m}}''}_{\delta }}</math> is mass flux at the interface due to phase change and VI is the velocity of the interface. For a threedimensional interface, there are three components of velocity: the normal direction, and two tangential directions, denoted by n, t<sub>1</sub>, and t<sub>2</sub>, respectively. Therefore, the velocity components should be defined according to these directions, as follows: 
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\mathbf{V}\cdot \mathbf{n}={{V}_{\mathbf{n}}}</math>  
+  </center>  
+  {{EquationRef(2)}}  
}  }  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\mathbf{V}\cdot {{\mathbf{t}}_{1}}={{V}_{{{\mathbf{t}}_{1}}}}</math>  
+  </center>  
+  {{EquationRef(3)}}  
}  }  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\mathbf{V}\cdot {{\mathbf{t}}_{2}}={{V}_{{{\mathbf{t}}_{2}}}}</math>  
+  </center>  
+  {{EquationRef(4)}}  
}  }  
Line 29:  Line 45:  
The interfacial mass balance can be rewritten in these terms by:  The interfacial mass balance can be rewritten in these terms by:  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{{\dot{m}}''}_{\delta }}={{\rho }_{\ell }}\left( {{V}_{\ell ,\mathbf{n}}}{{V}_{I,\mathbf{n}}} \right)={{\rho }_{v}}\left( {{V}_{v,\mathbf{n}}}{{V}_{I,\mathbf{n}}} \right)</math>  
+  </center>  
+  {{EquationRef(5)}}  
}  }  
Line 38:  Line 58:  
For a situation where the effects of surface tension and disjoining pressure are negligible, the momentum balance at the liquidvapor interface can be described by eq. (3.173), i.e.,  For a situation where the effects of surface tension and disjoining pressure are negligible, the momentum balance at the liquidvapor interface can be described by eq. (3.173), i.e.,  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\left( {{{\mathbf{{\tau }'}}}_{\ell }}{{{\mathbf{{\tau }'}}}_{v}} \right)\cdot \mathbf{n}={{{\dot{m}}''}_{\delta }}({{\mathbf{V}}_{\ell }}{{\mathbf{V}}_{v}})</math>  
+  </center>  
+  {{EquationRef(6)}}  
}  }  
  where n is the normal direction of the interface and points toward the vapor. The lefthand side of eq. (6) represents the force per unit area acting on the interface, while the righthand side represents the change of momentum across the interface. For applications involving thin film evaporation and condensation, the effects of surface tension and disjoining pressure will create additional forces on the interface; in these cases, the lefthand side of eq. (6) should be modified. The force per unit area created by the surface tension as indicated by eq. (  +  where n is the normal direction of the interface and points toward the vapor. The lefthand side of eq. (6) represents the force per unit area acting on the interface, while the righthand side represents the change of momentum across the interface. For applications involving thin film evaporation and condensation, the effects of surface tension and disjoining pressure will create additional forces on the interface; in these cases, the lefthand side of eq. (6) should be modified. The force per unit area created by the surface tension as indicated by eq. (51) is <math>\sigma (T)(1/{{R}_{I}}+1/{{R}_{II}})\mathbf{n}</math>. The force per unit area created by the disjoining pressure is <math>{{p}_{d}}\mathbf{n}.</math> For the situations where the interfacial temperature is not a constant, the contribution of the Marangoni effect, <math>(d\sigma /dT)\nabla {{T}_{\delta }}</math>, should also be included. Therefore, the momentum balance at the interface becomes 
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\begin{align}  
+  & \left( {{{\mathbf{{\tau }'}}}_{\ell }}{{{\mathbf{{\tau }'}}}_{v}} \right)\cdot \mathbf{n}+\sigma (T)\left( \frac{1}{{{R}_{I}}}+\frac{1}{{{R}_{II}}} \right)\mathbf{n}{{p}_{d}}\mathbf{n} \\  
+  & \left( \frac{d\sigma }{dT} \right)(\nabla {{T}_{\delta }}\cdot \mathbf{t})\mathbf{t}={{{{\dot{m}}''}}_{\delta }}({{\mathbf{V}}_{\ell }}{{\mathbf{V}}_{v}}) \\  
+  \end{align}</math>  
+  </center>  
+  {{EquationRef(7)}}  
}  }  
Line 52:  Line 83:  
On the lefthand side, the first term is the stress tensor, the second term is the capillary pressure, the third term is the disjoining pressure, and the fourth term is the Marangoni stress. The righthand side is the momentum transfer due to inertia. In this equation, the tangential direction, t, can either be t<sub>1</sub> or t<sub>2</sub>. The stress tensor is:  On the lefthand side, the first term is the stress tensor, the second term is the capillary pressure, the third term is the disjoining pressure, and the fourth term is the Marangoni stress. The righthand side is the momentum transfer due to inertia. In this equation, the tangential direction, t, can either be t<sub>1</sub> or t<sub>2</sub>. The stress tensor is:  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{\tau }'=p\mathbf{I}+2\mu \mathbf{D}\frac{2}{3}\mu \left( \nabla \cdot \mathbf{V} \right)\mathbf{I}</math>  
+  </center>  
+  {{EquationRef(8)}}  
}  }  
Line 59:  Line 94:  
The deformation tensor can be written for a reference frame that is adjusted to the interface:  The deformation tensor can be written for a reference frame that is adjusted to the interface:  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
  ==  +   width='100%' <center> 
  +  <math>\mathbf{D}=\frac{1}{2}\left[ \nabla \mathbf{V}+{{\left( \nabla \mathbf{V} \right)}^{T}} \right]=\left[ \begin{matrix}  
  +  \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{\mathbf{n}}}} & \frac{1}{2}\left( \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{1}}}}}{\partial {{x}_{\mathbf{n}}}} \right) & \frac{1}{2}\left( \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{2}}}}}{\partial {{x}_{\mathbf{n}}}} \right) \\  
  +  \frac{1}{2}\left( \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{1}}}}}{\partial {{x}_{\mathbf{n}}}} \right) & \frac{\partial {{V}_{{{\mathbf{t}}_{1}}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}} & \frac{1}{2}\left( \frac{\partial {{V}_{{{\mathbf{t}}_{1}}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{2}}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}} \right) \\  
  {  +  \frac{1}{2}\left( \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{2}}}}}{\partial {{x}_{\mathbf{n}}}} \right) & \frac{1}{2}\left( \frac{\partial {{V}_{{{\mathbf{t}}_{1}}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{2}}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}} \right) & \frac{\partial {{V}_{{{\mathbf{t}}_{2}}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}} \\ 
  +  \end{matrix} \right]</math>  
  +  </center>  
  +  {{EquationRef(9)}}  
  +  
  +  
  +  
  +  
  +  
}  }  
Line 79:  Line 109:  
The normal direction of the interface is [1 0 0], the first tangential direction is [0 1 0] and the second tangential direction is [0 0 1]. Therefore,  The normal direction of the interface is [1 0 0], the first tangential direction is [0 1 0] and the second tangential direction is [0 0 1]. Therefore,  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
  +   width='100%' <center>  
  +  <math>\begin{align}  
+  & {\tau }'\cdot \mathbf{n}=\left[ \begin{matrix}  
+  p & 0 & 0 \\  
+  \end{matrix} \right]+ \\  
+  & 2\mu \left[ \begin{matrix}  
+  \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{\mathbf{n}}}} & \frac{1}{2}\left( \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{1}}}}}{\partial {{x}_{\mathbf{n}}}} \right) & \frac{1}{2}\left( \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{2}}}}}{\partial {{x}_{\mathbf{n}}}} \right) \\  
+  \end{matrix} \right]\frac{2}{3}\mu \left[ \begin{matrix}  
+  \nabla \cdot \mathbf{V} & 0 & 0 \\  
+  \end{matrix} \right] \\  
+  \end{align}</math>  
+  </center>  
+  {{EquationRef(10)}}  
}  }  
Line 88:  Line 129:  
This can be reduced to the three components to obtain:  This can be reduced to the three components to obtain:  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{\tau }'\cdot \mathbf{n}\cdot \mathbf{n}=p+2\mu \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{\mathbf{n}}}}\frac{2}{3}\mu \nabla \cdot \mathbf{V}=p+\frac{4}{3}\mu \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{\mathbf{n}}}}\frac{2}{3}\mu \left( \frac{\partial {{V}_{{{\mathbf{t}}_{1}}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{2}}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}} \right)</math>  
+  </center>  
+  {{EquationRef(11)}}  
}  }  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{\tau }'\cdot \mathbf{n}\cdot {{\mathbf{t}}_{1}}=\mu \left( \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{1}}}}}{\partial {{x}_{\mathbf{n}}}} \right)</math>  
+  </center>  
+  {{EquationRef(12)}}  
}  }  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{\tau }'\cdot \mathbf{n}\cdot {{\mathbf{t}}_{2}}=\mu \left( \frac{\partial {{V}_{\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{{{\mathbf{t}}_{2}}}}}{\partial {{x}_{\mathbf{n}}}} \right)</math>  
+  </center>  
+  {{EquationRef(13)}}  
}  }  
Line 107:  Line 160:  
==Normal Direction==  ==Normal Direction==  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
  +   width='100%' <center>  
  +  <math>\begin{align}  
  +  & {{p}_{\ell }}+{{p}_{v}}+\frac{4}{3}\left( {{\mu }_{\ell }}\frac{\partial {{V}_{\ell ,\mathbf{n}}}}{\partial {{x}_{\ell ,\mathbf{n}}}}{{\mu }_{v}}\frac{\partial {{V}_{v,\mathbf{n}}}}{\partial {{x}_{v,\mathbf{n}}}} \right) \\  
  +  & \frac{2}{3}\left[ {{\mu }_{\ell }}\left( \frac{\partial {{V}_{\ell ,{{\mathbf{t}}_{1}}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{\ell ,{{\mathbf{t}}_{2}}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}} \right){{\mu }_{v}}\left( \frac{\partial {{V}_{v,{{\mathbf{t}}_{1}}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{v,{{\mathbf{t}}_{2}}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}} \right) \right] \\  
  {  +  & +\sigma \left( \frac{1}{{{R}_{I}}}+\frac{1}{{{R}_{II}}} \right){{p}_{d}}={{{{\dot{m}}''}}_{\delta }}\left( {{V}_{\ell ,\mathbf{n}}}{{V}_{v,\mathbf{n}}} \right) \\ 
  +  \end{align}</math>  
+  </center>  
+  {{EquationRef(14)}}  
}  }  
Line 120:  Line 175:  
Tangential 1  Tangential 1  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\mu }_{\ell }}\left( \frac{\partial {{V}_{\ell ,\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{\ell ,{{\mathbf{t}}_{1}}}}}{\partial {{x}_{\mathbf{n}}}} \right){{\mu }_{v}}\left( \frac{\partial {{V}_{v,\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{v,{{\mathbf{t}}_{1}}}}}{\partial {{x}_{\mathbf{n}}}} \right)\left( \frac{d\sigma }{dT} \right)\left( \frac{\partial {{T}_{\delta }}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}} \right)={{{\dot{m}}''}_{\delta }}\left( {{V}_{\ell ,{{\mathbf{t}}_{1}}}}{{V}_{v,{{\mathbf{t}}_{1}}}} \right)</math>  
+  </center>  
+  {{EquationRef(15)}}  
}  }  
Line 127:  Line 186:  
Tangential 2  Tangential 2  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\mu }_{\ell }}\left( \frac{\partial {{V}_{\ell ,\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{\ell ,{{\mathbf{t}}_{2}}}}}{\partial {{x}_{\mathbf{n}}}} \right){{\mu }_{v}}\left( \frac{\partial {{V}_{v,\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{v,{{\mathbf{t}}_{2}}}}}{\partial {{x}_{\mathbf{n}}}} \right)\left( \frac{d\sigma }{dT} \right)\left( \frac{\partial {{T}_{\delta }}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}} \right)={{{\dot{m}}''}_{\delta }}\left( {{V}_{\ell ,{{\mathbf{t}}_{2}}}}{{V}_{v,{{\mathbf{t}}_{2}}}} \right)</math>  
+  </center>  
+  {{EquationRef(16)}}  
}  }  
The nonslip condition at the liquidvapor interface requires that  The nonslip condition at the liquidvapor interface requires that  
  <math>{{V}_{  +  <math>{{V}_{\ell ,{{\mathbf{t}}_{1}}}}={{V}_{v,{{\mathbf{t}}_{1}}}}</math> and <math>{{V}_{\ell ,{{\mathbf{t}}_{2}}}}={{V}_{v,{{\mathbf{t}}_{2}}}}.</math> The momentum balance at the tangential directions becomes 
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\mu }_{\ell }}\left( \frac{\partial {{V}_{\ell ,\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{\ell ,{{\mathbf{t}}_{1}}}}}{\partial {{x}_{\mathbf{n}}}} \right)={{\mu }_{v}}\left( \frac{\partial {{V}_{v,\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{v,{{\mathbf{t}}_{1}}}}}{\partial {{x}_{\mathbf{n}}}} \right)+\left( \frac{d\sigma }{dT} \right)\left( \frac{\partial {{T}_{\delta }}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}} \right)</math>  
+  </center>  
+  {{EquationRef(17)}}  
}  }  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\mu }_{\ell }}\left( \frac{\partial {{V}_{\ell ,\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{\ell ,{{\mathbf{t}}_{2}}}}}{\partial {{x}_{\mathbf{n}}}} \right)={{\mu }_{v}}\left( \frac{\partial {{V}_{v,\mathbf{n}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}}+\frac{\partial {{V}_{v,{{\mathbf{t}}_{2}}}}}{\partial {{x}_{\mathbf{n}}}} \right)+\left( \frac{d\sigma }{dT} \right)\left( \frac{\partial {{T}_{\delta }}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}} \right)</math>  
+  </center>  
+  {{EquationRef(18)}}  
}  }  
For most applications, the evaporation or condensation rate –  For most applications, the evaporation or condensation rate –  
  <math>{{{  +  <math>{{{\dot{m}}''}_{\delta }}</math> – is not very high; therefore, it can be assumed that <math>{{{\dot{m}}''}_{\delta }}({{\mathbf{V}}_{\ell }}{{\mathbf{V}}_{v}})\doteq 0</math>. If the liquid and vapor phases are further assumed to be inviscid (<math>{{\mathbf{\tau }}_{\ell }}={{\mathbf{\tau }}_{v}}\doteq 0</math>), the momentum equation at the interface can be reduced to 
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{p}_{v}}{{p}_{\ell }}=\sigma (T)\left( \frac{1}{{{R}_{I}}}+\frac{1}{{{R}_{II}}} \right){{p}_{d}}</math>  
+  </center>  
+  {{EquationRef(19)}}  
}  }  
Line 157:  Line 232:  
The energy balance at the interface can be obtained from eq. (3.180), i.e.,  The energy balance at the interface can be obtained from eq. (3.180), i.e.,  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\left( {{{\mathbf{{q}''}}}_{\ell }}{{{\mathbf{{q}''}}}_{v}} \right)\cdot \mathbf{n}(\mathbf{n}\cdot {{\mathbf{{\tau }'}}_{\ell }})\cdot {{\mathbf{V}}_{\ell ,rel}}+(\mathbf{n}\cdot {{\mathbf{{\tau }'}}_{v}})\cdot {{\mathbf{V}}_{v,rel}}={{{\dot{m}}''}_{\delta }}\left[ \left( {{e}_{v}}+\frac{\mathbf{V}_{v,rel}^{2}}{2} \right)\left( {{e}_{\ell }}+\frac{\mathbf{V}_{\ell ,rel}^{2}}{2} \right) \right]</math>  
+  </center>  
+  {{EquationRef(20)}}  
}  }  
Line 164:  Line 243:  
where the work done by surface tension and disjoining pressure is neglected. If the velocity of the reference frame is taken as the interfacial velocity VI, eq. (20) can be rewritten as  where the work done by surface tension and disjoining pressure is neglected. If the velocity of the reference frame is taken as the interfacial velocity VI, eq. (20) can be rewritten as  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\begin{align}  
+  & \left( {{k}_{v}}\nabla {{T}_{v}}{{k}_{\ell }}\nabla {{T}_{\ell }} \right)\cdot \mathbf{n}(\mathbf{n}\cdot {{{\mathbf{{\tau }'}}}_{\ell }})\cdot ({{\mathbf{V}}_{\ell }}{{\mathbf{V}}_{I}})+(\mathbf{n}\cdot {{{\mathbf{{\tau }'}}}_{v}})\cdot ({{\mathbf{V}}_{v}}{{\mathbf{V}}_{I}}) \\  
+  & ={{{{\dot{m}}''}}_{\delta }}\left\{ \left[ {{e}_{v}}+\frac{{{({{\mathbf{V}}_{v}}{{\mathbf{V}}_{I}})}^{2}}}{2} \right]\left[ {{e}_{\ell }}+\frac{{{({{\mathbf{V}}_{\ell }}{{\mathbf{V}}_{I}})}^{2}}}{2} \right] \right\} \\  
+  \end{align}</math>  
+  </center>  
+  {{EquationRef(21)}}  
}  }  
where the heat fluxes in the liquid and vapor phases have been determined by Fourier’s law of conduction.  where the heat fluxes in the liquid and vapor phases have been determined by Fourier’s law of conduction.  
  Equation (21) can be rewritten in terms of enthalpy (<math>h=e+p/  +  
+  Equation (21) can be rewritten in terms of enthalpy (<math>h=e+p/\rho </math>), i.e.,  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\begin{align}  
+  & \left( {{k}_{v}}\nabla {{T}_{v}}{{k}_{\ell }}\nabla {{T}_{\ell }} \right)\cdot \mathbf{n}(\mathbf{n}\cdot {{{\mathbf{{\tau }'}}}_{\ell }})\cdot ({{\mathbf{V}}_{\ell }}{{\mathbf{V}}_{I}})+(\mathbf{n}\cdot {{{\mathbf{{\tau }'}}}_{v}})\cdot ({{\mathbf{V}}_{v}}{{\mathbf{V}}_{I}}) \\  
+  & ={{{{\dot{m}}''}}_{\delta }}\left[ {{h}_{\ell v}}\left( \frac{{{p}_{v}}}{{{\rho }_{v}}}\frac{{{p}_{\ell }}}{{{\rho }_{\ell }}} \right)+\frac{1}{2}\mathbf{V}_{v}^{2}{{\mathbf{V}}_{v}}\cdot {{\mathbf{V}}_{I}}\frac{1}{2}\mathbf{V}_{\ell }^{2}+{{\mathbf{V}}_{\ell }}\cdot {{\mathbf{V}}_{I}} \right] \\  
+  \end{align}</math>  
+  </center>  
+  {{EquationRef(22)}}  
}  }  
  where <math>{{h}_{  +  where <math>{{h}_{\ell v}}</math> is the difference between the enthalpy of vapor and liquid at saturation, i.e., the latent heat of vaporization. The stress tensor in eq. (22) can be expressed as 
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{\tau }'=p\mathbf{I}+2\mu \mathbf{D}\frac{2}{3}\mu \left( \nabla \cdot \mathbf{V} \right)\mathbf{I}=p\mathbf{I}+\tau </math>  
+  </center>  
+  {{EquationRef(23)}}  
}  }  
Line 186:  Line 284:  
Since the relative velocities at the interface satisfy and , we have  Since the relative velocities at the interface satisfy and , we have  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{p}_{\ell }}\left( {{\mathbf{V}}_{\ell }}{{\mathbf{V}}_{I}} \right)\cdot \mathbf{n}{{p}_{v}}\left( {{\mathbf{V}}_{v}}{{\mathbf{V}}_{I}} \right)\cdot \mathbf{n}={{{\dot{m}}''}_{\delta }}\left[ \left( \frac{{{p}_{v}}}{{{\rho }_{v}}}\frac{{{p}_{\ell }}}{{{\rho }_{\ell }}} \right) \right]</math>  
+  </center>  
+  {{EquationRef(24)}}  
}  }  
Line 193:  Line 295:  
Substituting eqs. (23) and (24) into eq. (22), the energy balance can be written as  Substituting eqs. (23) and (24) into eq. (22), the energy balance can be written as  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\begin{align}  
+  & \left( {{k}_{v}}\nabla {{T}_{v}}{{k}_{\ell }}\nabla {{T}_{\ell }} \right)\cdot \mathbf{n}\left( \mathbf{n}\cdot {{\tau }_{\ell }} \right)\left( {{\mathbf{V}}_{\ell }}{{\mathbf{V}}_{I}} \right)+\left( \mathbf{n}\cdot {{\tau }_{v}} \right)\left( {{\mathbf{V}}_{v}}{{\mathbf{V}}_{I}} \right) \\  
+  & ={{{{\dot{m}}''}}_{\delta }}\left[ {{h}_{\ell v}}+\frac{1}{2}\mathbf{V}_{v}^{2}{{\mathbf{V}}_{v}}\cdot {{\mathbf{V}}_{I}}\frac{1}{2}\mathbf{V}_{\ell }^{2}+{{\mathbf{V}}_{\ell }}\cdot {{\mathbf{V}}_{I}} \right] \\  
+  \end{align}</math>  
+  </center>  
+  {{EquationRef(25)}}  
}  }  
Line 200:  Line 309:  
To simplify the energy equation, the kinetic energy terms are considered negligible and noslip conditions are assumed at the interface,  To simplify the energy equation, the kinetic energy terms are considered negligible and noslip conditions are assumed at the interface,  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{V}_{\ell ,\mathbf{t}}}={{V}_{v,\mathbf{t}}}={{V}_{I,\mathbf{t}}}</math>  
+  </center>  
+  {{EquationRef(26)}}  
}  }  
Line 207:  Line 320:  
Therefore, the energy equation can be rewritten as  Therefore, the energy equation can be rewritten as  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
  +   width='100%' <center>  
  +  <math>\begin{align}  
  +  & \left( {{k}_{v}}\frac{\partial {{T}_{v}}}{\partial {{x}_{\mathbf{n}}}}{{k}_{\ell }}\frac{\partial {{T}_{\ell }}}{\partial {{x}_{\mathbf{n}}}} \right) \\  
  +  & ={{{{\dot{m}}''}}_{\delta }}\left[ {{h}_{\ell v}}+\frac{4}{3}{{\nu }_{\ell }}\frac{\partial {{V}_{\ell ,\mathbf{n}}}}{\partial {{x}_{\mathbf{n}}}}\frac{2}{3}{{\nu }_{\ell }}\left( \frac{\partial {{V}_{\ell ,{{\mathbf{t}}_{1}}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{\ell ,{{\mathbf{t}}_{2}}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}} \right) \right. \\  
  {  +  & \left. \text{ }\frac{4}{3}{{\nu }_{v}}\frac{\partial {{V}_{v,\mathbf{n}}}}{\partial {{x}_{\mathbf{n}}}}+\frac{2}{3}{{\nu }_{v}}\left( \frac{\partial {{V}_{v,{{\mathbf{t}}_{1}}}}}{\partial {{x}_{{{\mathbf{t}}_{1}}}}}+\frac{\partial {{V}_{v,{{\mathbf{t}}_{2}}}}}{\partial {{x}_{{{\mathbf{t}}_{2}}}}} \right) \right] \\ 
  +  \end{align}</math>  
+  </center>  
+  {{EquationRef(27)}}  
}  }  
Line 221:  Line 336:  
The energy balance at the interface can be simplified by assuming that the change in the kinetic energy across the interface is negligible, i.e.,  The energy balance at the interface can be simplified by assuming that the change in the kinetic energy across the interface is negligible, i.e.,  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>\left( {{k}_{v}}\nabla {{T}_{v}}{{k}_{\ell }}\nabla {{T}_{\ell }} \right)\cdot \mathbf{n}={{{\dot{m}}''}_{\delta }}{{h}_{\ell v}}</math>  
+  </center>  
+  {{EquationRef(28)}}  
}  }  
Line 233:  Line 352:  
Some dependent variables, such as temperature, are continuous across the phases; other variables, such as concentration, are discontinuous.  Some dependent variables, such as temperature, are continuous across the phases; other variables, such as concentration, are discontinuous.  
  For a general interface between phases k and j in a multicomponent system, a local balance in mass flux of species i must be upheld. The total species mass flux, <math>{{{  +  For a general interface between phases k and j in a multicomponent system, a local balance in mass flux of species i must be upheld. The total species mass flux, <math>{{{\dot{m}}''}_{i}}</math>, at an interface is: 
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{{\dot{m}}''}_{i}}={{\rho }_{k,i}}\left( {{\mathbf{V}}_{k,i}}{{\mathbf{V}}_{I}} \right)\cdot \mathbf{n}={{\rho }_{j,i}}\left( {{\mathbf{V}}_{j,i}}{{\mathbf{V}}_{I}} \right)\cdot \mathbf{n}</math>  
+  </center>  
+  {{EquationRef(29)}}  
}  }  
The velocity of species i in phase k and phase j is  The velocity of species i in phase k and phase j is  
  <math>{{  +  <math>{{\mathbf{V}}_{k,i}}</math> and <math>{{\mathbf{V}}_{j,i}}</math>, respectively. These velocities are defined as: 
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\rho }_{k,i}}{{\mathbf{V}}_{k,i}}={{\mathbf{J}}_{k,i}}+{{\omega }_{k,i}}{{\rho }_{k}}{{\mathbf{V}}_{k}}</math>  
+  </center>  
+  {{EquationRef(30)}}  
}  }  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\rho }_{j,i}}{{\mathbf{V}}_{j,i}}={{\mathbf{J}}_{j,i}}+{{\omega }_{j,i}}{{\rho }_{j}}{{\mathbf{V}}_{j}}</math>  
+  </center>  
+  {{EquationRef(31)}}  
}  }  
Line 255:  Line 386:  
Using the interfacial species mass balance, and substituting the definition of the species velocity, the interfacial species mass flux is:  Using the interfacial species mass balance, and substituting the definition of the species velocity, the interfacial species mass flux is:  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{{\dot{m}}''}_{i}}={{\mathbf{J}}_{k,i}}\cdot \mathbf{n}+{{\omega }_{k,i}}{{\rho }_{k}}\left( {{\mathbf{V}}_{k}}{{\mathbf{V}}_{I}} \right)={{\mathbf{J}}_{j,i}}\cdot \mathbf{n}+{{\omega }_{j,i}}{{\rho }_{j}}\left( {{\mathbf{V}}_{jk}}{{\mathbf{V}}_{I}} \right)</math>  
+  </center>  
+  {{EquationRef(32)}}  
}  }  
Line 262:  Line 397:  
Remembering the overall mass conservation at the interface,  Remembering the overall mass conservation at the interface,  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{\dot{m}}''={{\rho }_{k}}\left( {{\mathbf{V}}_{k}}{{\mathbf{V}}_{I}} \right)\cdot n={{\rho }_{j}}\left( {{\mathbf{V}}_{j}}{{\mathbf{V}}_{I}} \right)\cdot n</math>  
+  </center>  
+  {{EquationRef(33)}}  
}  }  
Line 269:  Line 408:  
The interfacial species mass flux is:  The interfacial species mass flux is:  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{{\dot{m}}''}_{i}}={{\mathbf{J}}_{k,i}}\cdot \mathbf{n}+{{\omega }_{k,i}}{\dot{m}}''={{\mathbf{J}}_{j,i}}\cdot \mathbf{n}+{{\omega }_{j,i}}{\dot{m}}''</math>  
+  </center>  
+  {{EquationRef(34)}}  
}  }  
Line 276:  Line 419:  
In some problems, the species mass flux will be specified, and the total mass flux is simply a sum of all the species mass fluxes.  In some problems, the species mass flux will be specified, and the total mass flux is simply a sum of all the species mass fluxes.  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{\dot{m}}''=\sum\limits_{i}^{{}}{{{{{\dot{m}}''}}_{i}}}</math>  
+  </center>  
+  {{EquationRef(35)}}  
}  }  
Line 283:  Line 430:  
If the species mass flux is not specified, the total mass flux at an interface can be calculated from the interfacial species mass flux equation, eq. (34):  If the species mass flux is not specified, the total mass flux at an interface can be calculated from the interfacial species mass flux equation, eq. (34):  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{\dot{m}}''=\frac{\left( {{\mathbf{J}}_{j,i}}{{\mathbf{J}}_{k,i}} \right)\cdot n}{{{\omega }_{k,i}}{{\omega }_{j,i}}}</math>  
+  </center>  
+  {{EquationRef(36)}}  
}  }  
Line 290:  Line 441:  
Since the species equation is second order in space, and there are two phases, there needs to be a secondary condition relating the species mass fraction in phases k and j. This can be done by expressing the mass fraction of species i in phase k as a function of the species mass fraction i in phase j, or vice versa.  Since the species equation is second order in space, and there are two phases, there needs to be a secondary condition relating the species mass fraction in phases k and j. This can be done by expressing the mass fraction of species i in phase k as a function of the species mass fraction i in phase j, or vice versa.  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\omega }_{k,i}}={{\omega }_{k,i}}\left( {{\omega }_{j,i}} \right)</math>  
+  </center>  
+  {{EquationRef(37)}}  
}  }  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\omega }_{j,i}}={{\omega }_{j,i}}\left( {{\omega }_{k,i}} \right)</math>  
+  </center>  
+  {{EquationRef(38)}}  
}  }  
Line 302:  Line 461:  
Note that, in general, the species mass fraction is not a continuous function, and there is almost always a jump condition at the interface when two or more species are present. To make this point clear, the simplest case of a binary mixture is considered. For a binary mixture, the species diffusion flux, J, can be calculated by Fick’s law.  Note that, in general, the species mass fraction is not a continuous function, and there is almost always a jump condition at the interface when two or more species are present. To make this point clear, the simplest case of a binary mixture is considered. For a binary mixture, the species diffusion flux, J, can be calculated by Fick’s law.  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\mathbf{J}}_{k,1}}\cdot n={{\rho }_{k}}{{D}_{k,12}}\nabla {{\omega }_{k,1}}\cdot \mathbf{n}</math>  
+  </center>  
+  {{EquationRef(39)}}  
}  }  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\mathbf{J}}_{j,1}}\cdot n={{\rho }_{j}}{{D}_{j,12}}\nabla {{\omega }_{j,1}}\cdot \mathbf{n}</math>  
+  </center>  
+  {{EquationRef(40)}}  
}  }  
Line 314:  Line 481:  
There are several examples of different phenomena of a binary mixture in the presence of a liquid/vapor, solid/vapor, or liquid/solid interface. In other examples, the species flux in one of the phases is unimportant; therefore, only one phase must be considered.  There are several examples of different phenomena of a binary mixture in the presence of a liquid/vapor, solid/vapor, or liquid/solid interface. In other examples, the species flux in one of the phases is unimportant; therefore, only one phase must be considered.  
  A classical example of species flux is the condition of zero mass flux due to an impermeable surface where species A does not diffuse to the stationary media <math>  +  A classical example of species flux is the condition of zero mass flux due to an impermeable surface where species A does not diffuse to the stationary media <math>\left( \nabla {{x}_{A}}\cdot n=0 \right)</math> or <math>\left( \partial {{x}_{A}}/\partial y=0 \right)</math>. The case of constant surface concentration is more typical; however, all types of mass transfer boundary conditions are demonstrated in three categories below. 
a. The mass transfer from solid or liquid to a gas stream  a. The mass transfer from solid or liquid to a gas stream  
  Evaporation and sublimation are typical examples of mass transfer from liquid or solid to a gas mixture, as shown in Fig.  +  Evaporation and sublimation are typical examples of mass transfer from liquid or solid to a gas mixture, as shown in Fig. 1, where 
  <math>{{x}_{A,  +  <math>{{x}_{A,\ell }}</math>, <math>{{x}_{A,s}}</math>, <math>{{x}_{A,g}}</math> are molar fractions in liquid, solid, and gas, respectively, and 
  <math>{{{  +  <math>{{{\dot{m}}''}_{A}}</math> is mass flux for species A. There are two cases shown in Fig. 1. In case I, a pure liquid/solid is evaporating/sublimating into a gaseous mixture and the evaporation /sublimation rate is controlled purely by the species gradient in the gas phase. In case II, the liquid/solid is not pure, therefore concentration gradients also exist in the liquid/solid phases. The mass transfer rate is limited by both phases. Also note that the concentration gradient in the solid [Fig. 1(b)] is steeper and does not penetrate as deeply when compared to the liquid concentration gradient [Fig. 1(a)]. Also, the liquid concentration gradient is steeper and does not penetrate as far compared to the gas. These trends are due to the ratio of the mass diffusivities in each of the phases. 
A simplified relationship for the mass concentration condition at the interface can be obtained assuming the gas mixture is approximated by the ideal gas and the solid or liquid has a high concentration of A. The condition within the gas stream is of interest, because the main resistance to mass transfer is within that region. With the preceding assumptions, the partial vapor pressure of A in the gas mixture at the interface can be approximated from Raoult’s Law:  A simplified relationship for the mass concentration condition at the interface can be obtained assuming the gas mixture is approximated by the ideal gas and the solid or liquid has a high concentration of A. The condition within the gas stream is of interest, because the main resistance to mass transfer is within that region. With the preceding assumptions, the partial vapor pressure of A in the gas mixture at the interface can be approximated from Raoult’s Law:  
+  
+  [[Image:Species concentration and mass transfer from solid and or liquid to a gas mixture.jpgthumb400 pxalt=Figure 1 Species concentration and mass transfer from solid and or liquid to a gas mixture. Figure 1 Species concentration and mass transfer from solid and or liquid to a gas mixture.]]  
  <math>{{  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\left. {{p}_{A}} \right}_{y={{0}^{}}}}={{\left. {{x}_{A,\ell }} \right}_{y={{0}^{+}}}}\left( {{p}_{A,sat}} \right)</math>  
+  </center>  
+  {{EquationRef(41)}}  
+  }  
+  
  <math>{{  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\left. {{p}_{A}} \right}_{y={{0}^{}}}}={{\left. {{x}_{A,s}} \right}_{y={{0}^{+}}}}\left( {{p}_{A,sat}} \right)</math>  
+  </center>  
+  {{EquationRef(42)}}  
+  }  
+  
  where pA is the partial vapor pressure of A in the gas mixture,  +  where pA is the partial vapor pressure of A in the gas mixture, <math>{{\left. {{x}_{A,\ell }} \right}_{y={{0}^{+}}}}</math> is the mole fraction of species A in liquid, <math>{{\left. {{x}_{A,s}} \right}_{y={{0}^{+}}}}</math> is the mole fraction of species in solid, and pA,sat is the normal saturation pressure of species A at the surface interface. Clearly, if the solid or liquid is made of pure species A, then 
  <math>{{  +  <math>{{\left. {{x}_{A,\ell }} \right}_{y=0}}{{\left. ={{x}_{A,s}} \right}_{y=0}}=1</math> and eqs. (41) and (42) reduce to 
  +  
  <math>{{  +  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\left. {{p}_{A}} \right}_{y={{0}^{}}}}={{p}_{A,sat}}</math>  
+  </center>  
+  {{EquationRef(43)}}  
}  }  
  This means that the partial vapor pressure of species A in the gas mixture at the interface is equal to the normal saturation pressure for species A, which is a function of interfacial temperature and can be obtained from thermodynamic tables in Appendix B. At the interface, species A and B in the gas phase must be in equilibrium with species A and B in the liquid or solid phase, except in extreme circumstances. Knowing <math>{{  +  This means that the partial vapor pressure of species A in the gas mixture at the interface is equal to the normal saturation pressure for species A, which is a function of interfacial temperature and can be obtained from thermodynamic tables in Appendix B. At the interface, species A and B in the gas phase must be in equilibrium with species A and B in the liquid or solid phase, except in extreme circumstances. Knowing <math>{{\left. {{p}_{A}} \right}_{y={{0}^{}}}}</math> and total pressure p, one can easily calculate the mole fraction <math>{{x}_{A,g}}</math> and mass fraction <math>{{\omega }_{A,g}}</math> of species A at the interface on the gas side by the following relation 
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{x}_{A,g}}=\frac{{{\left. {{p}_{A}} \right}_{y={{0}^{}}}}}{p}</math>  
+  </center>  
+  {{EquationRef(44)}}  
}  }  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\omega }_{A,g}}=\frac{{{x}_{A,g}}{{M}_{A}}}{{{x}_{A,g}}{{M}_{A}}+{{x}_{B,g}}{{M}_{B}}}</math>  
+  </center>  
+  {{EquationRef(45)}}  
}  }  
  where M<sub>A</sub> and M<sub>B</sub> are molecular weights of species A and B. A conventional example for the case in Fig.  +  where M<sub>A</sub> and M<sub>B</sub> are molecular weights of species A and B. A conventional example for the case in Fig. 1(a) is the evaporation of water to a waterair mixture; obviously, one can imagine that the water absorbed a small amount of air. Sublimation of iodine in air is an application for the case of Fig. 1(b). 
Mass transfer from a solid to a gaseous state sometimes requires the specification of diffusion molar flux, rather than concentration, at the solid surface.  Mass transfer from a solid to a gaseous state sometimes requires the specification of diffusion molar flux, rather than concentration, at the solid surface.  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>J_{A}^{*}=c{{D}_{AB}}\frac{\partial {{x}_{A}}}{\partial y}</math>  
+  </center>  
+  {{EquationRef(46)}}  
}  }  
  where <math>J_{A}^{*}</math> can be a function of concentration and not necessarily constant. One application is related to catalytic surface reaction. For this case, <math>{  +  where <math>J_{A}^{*}</math> can be a function of concentration and not necessarily constant. One application is related to catalytic surface reaction. For this case, <math>{\dot{m}}''=0</math> because the mass flux of reactants equals the mass flux of products. Therefore, from eq. (34), 
  <math>{{{  +  <math>{{{\dot{m}}''}_{i}}={{J}_{k,i}}\cdot n</math>. Catalytic surfaces are used to promote heterogeneous reactions (see Chapter 3), which occur at the surface; the appropriate boundary condition is 
  +  
  { class='wikitable' border='0' width='100%' <center>  +  [[Image:Species concentration and mass transfer from gas mixture to liquid or solid.jpgthumb400 pxalt=Figure 2 Species concentration and mass transfer from gas mixture to liquid or solid Figure 2 Species concentration and mass transfer from gas mixture to liquid or solid]] 
  +  
+  { class='wikitable' border='0'  
+    
+   width='100%' <center>  
+  <math>J_{A}^{*}={{{k}'}_{1}}{{\left( {{\left. {{c}_{A}} \right}_{y=0}} \right)}^{n}}</math>  
+  </center>  
+  {{EquationRef(47)}}  
}  }  
  where <math>{{{k}  +  where <math>{{{k}'}_{1}}</math> is the reaction rate constant and n is the order of reaction. 
+  
+  [[Image:Species concentration and mass transfer from solid to liquid and liquid to solid.jpgthumb400 pxalt=Figure 3 Species concentration and mass transfer from solid to liquid and liquid to solid Figure 3 Species concentration and mass transfer from solid to liquid and liquid to solid]]  
b. Mass transfer from gases to liquids or solids  b. Mass transfer from gases to liquids or solids  
  There are two major forms of mass transfer from gas to liquid, shown in Fig.  +  There are two major forms of mass transfer from gas to liquid, shown in Fig. 2. Case I is for condensable species in the liquid phase, or for species that can deposit vapor into the solid phase. In this case, a species of low concentration in the gas phase can condense at a higher concentration in the liquid phase. This phenomenon is important in distillation processes. Case II is for species that are weakly soluble in a liquid or a solid. In the liquid (<math>{{x}_{A,\ell }}</math> is small), Henry’s law will relate the mole fraction of species A in the liquid to the partial vapor pressure of A in the gas mixture at the interface by following relation [Fig. 2(a)] 
  +  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\left. {{x}_{A,\ell }} \right}_{y={{0}^{+}}}}=\frac{{{\left. {{p}_{A}} \right}_{y={{0}^{}}}}}{H}</math>  
+  </center>  
+  {{EquationRef(48)}}  
}  }  
Line 386:  Line 592:  
The concentration of gas in a solid at the interface is usually obtained by the use of a property known as the solubility, S, defined below  The concentration of gas in a solid at the interface is usually obtained by the use of a property known as the solubility, S, defined below  
  { class='wikitable' border='0' width='100%' <center>  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{\left. {{c}_{A,s}} \right}_{y={{0}^{+}}}}=S{{\left. {{p}_{A,g}} \right}_{y={{0}^{}}}}</math>  
+  </center>  
+  {{EquationRef(49)}}  
}  }  
  where <math>{{c}_{A,s}}</math> is the mole concentration of species A in the solid at the interface and <math>{{  +  where <math>{{c}_{A,s}}</math> is the mole concentration of species A in the solid at the interface and <math>{{\left. {{p}_{A,g}} \right}_{y={{0}^{}}}}</math> is the vapor partial pressure of species A in the gas at the interface. The values of S (<math>\text{kmol/Pa}{{\text{m}}^{\text{3}}}</math>) for vapor gassolid combinations are given in chemistry handbooks; some are in Appendix B (Tables B.76, B.78 and B.79). Special care should be exercised when using eq. (49) for the variation in form and units – two different versions of S are presented in Appendix B. An example of this case is the diffusion of helium in a glass. The dissolution of gas in metals is much more complex and depends on the type of metal used. The dissolution of gas in some metals can be reversed, such as the case of hydrogen into titanium. In contrast to hydrogen, oxygen dissolution in titanium is irreversible but is complicated because it forms a layer of scale T<sub>iO2</sub> on the surface. This topic is beyond the scope of this book, and interested readers should refer to a metallurgy or materials handbook to consult appropriate phase diagrams. 
c. Mass transit from solid to liquid or liquid to solid  c. Mass transit from solid to liquid or liquid to solid  
  Two cases are also presented for the case of mass transfer between a solid and a liquid (Fig.  +  Two cases are also presented for the case of mass transfer between a solid and a liquid (Fig. 3). Case I represents a pure solid dissolving into a liquid, or a pure liquid diffusing into a solid. Case II represents a solid mixture melting into a liquid mixture, or a liquid mixture solidifying into a solid mixture. For the first case, the mass transfer is related to the solubility, which can be found in chemistry handbooks (Lide, 2004) and some are reproduced in Appendix B. A conventional example for this case is the dissolution of salt into water. 
For a liquid/solid interface during melting and solidification, the ratio of the species concentration of the solid in liquid phases is called the partition ratio, K<sub>p</sub>.  For a liquid/solid interface during melting and solidification, the ratio of the species concentration of the solid in liquid phases is called the partition ratio, K<sub>p</sub>.  
  <math>{{K}_{p}}=  +  { class='wikitable' border='0' 
  +    
+   width='100%' <center>  
+  <math>{{K}_{p}}=\frac{{{c}_{A,s}}}{{{c}_{A,\ell }}}</math>  
+  </center>  
+  {{EquationRef(50)}}  
+  }  
+  
The partition coefficient is a function of temperature. Furthermore, when the liquid and solid lines are nearly straight, it is a constant. This information can be obtained from phase diagrams for different solid/liquid mixtures.  The partition coefficient is a function of temperature. Furthermore, when the liquid and solid lines are nearly straight, it is a constant. This information can be obtained from phase diagrams for different solid/liquid mixtures.  
+  
+  ==Additional Equations==  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%" <center>  
+  <math>{{p}_{cap}}=\Delta p=\sigma \left( \frac{1}{{{R}_{I}}}+\frac{1}{{{R}_{II}}} \right)=\sigma ({{K}_{1}}+{{K}_{2}})</math>  
+  </center>  
+  {{EquationRef(51)}}  
+  }  
+  
+  ==References==  
+  
+  Faghri, A., and Zhang, Y., 2006, ''Transport Phenomena in Multiphase Systems'', Elsevier, Burlington, MA.  
+  
+  ==Further Reading==  
+  
+  ==External Links== 
Current revision as of 20:23, 15 July 2010
The conservation laws for transport phenomena can be reduced to local partial differential equations if they are considered at a point that does not belong to a surface of discontinuity, such as an interface. When considering a discontinuous point, appropriate jump conditions relating the values of the fundamental quantities on both sides of the interface should be considered. Jump conditions at an interface were discussed in Section 3.3.6, but the effects of surface tension and disjoining pressure associated with a nonflat liquidvapor interface were not taken into account. It is the objective of this subsection to specify mass, momentum, and energy balance at a nonflat liquidvapor interface, as well as species balance in solidliquidvapor interfaces.
Contents 
Mass Balance
At a liquidvapor interface, the mass balance is

where is mass flux at the interface due to phase change and VI is the velocity of the interface. For a threedimensional interface, there are three components of velocity: the normal direction, and two tangential directions, denoted by n, t_{1}, and t_{2}, respectively. Therefore, the velocity components should be defined according to these directions, as follows:



The interfacial mass balance can be rewritten in these terms by:

Momentum Balance
For a situation where the effects of surface tension and disjoining pressure are negligible, the momentum balance at the liquidvapor interface can be described by eq. (3.173), i.e.,

where n is the normal direction of the interface and points toward the vapor. The lefthand side of eq. (6) represents the force per unit area acting on the interface, while the righthand side represents the change of momentum across the interface. For applications involving thin film evaporation and condensation, the effects of surface tension and disjoining pressure will create additional forces on the interface; in these cases, the lefthand side of eq. (6) should be modified. The force per unit area created by the surface tension as indicated by eq. (51) is . The force per unit area created by the disjoining pressure is For the situations where the interfacial temperature is not a constant, the contribution of the Marangoni effect, , should also be included. Therefore, the momentum balance at the interface becomes

On the lefthand side, the first term is the stress tensor, the second term is the capillary pressure, the third term is the disjoining pressure, and the fourth term is the Marangoni stress. The righthand side is the momentum transfer due to inertia. In this equation, the tangential direction, t, can either be t_{1} or t_{2}. The stress tensor is:

The deformation tensor can be written for a reference frame that is adjusted to the interface:

The normal direction of the interface is [1 0 0], the first tangential direction is [0 1 0] and the second tangential direction is [0 0 1]. Therefore,

This can be reduced to the three components to obtain:



The momentum equation balance at the interface is then broken into its three components, as follows:
Normal Direction

Tangential 1

Tangential 2

The nonslip condition at the liquidvapor interface requires that
and The momentum balance at the tangential directions becomes


For most applications, the evaporation or condensation rate –
– is not very high; therefore, it can be assumed that . If the liquid and vapor phases are further assumed to be inviscid (), the momentum equation at the interface can be reduced to

Energy
The energy balance at the interface can be obtained from eq. (3.180), i.e.,

where the work done by surface tension and disjoining pressure is neglected. If the velocity of the reference frame is taken as the interfacial velocity VI, eq. (20) can be rewritten as

where the heat fluxes in the liquid and vapor phases have been determined by Fourier’s law of conduction.
Equation (21) can be rewritten in terms of enthalpy (h = e + p / ρ), i.e.,

where is the difference between the enthalpy of vapor and liquid at saturation, i.e., the latent heat of vaporization. The stress tensor in eq. (22) can be expressed as

Since the relative velocities at the interface satisfy and , we have

Substituting eqs. (23) and (24) into eq. (22), the energy balance can be written as

To simplify the energy equation, the kinetic energy terms are considered negligible and noslip conditions are assumed at the interface,

Therefore, the energy equation can be rewritten as

where are kinematic viscosities of liquid and vapor phases, respectively.
The energy balance at the interface can be simplified by assuming that the change in the kinetic energy across the interface is negligible, i.e.,

Equations (1), (19), and (28) are widely used in the analysis of evaporation and condensation.
Species Balance
In general, to obtain a detailed solution one needs to solve the continuity, species, momentum, and energy equations in each phase (presented in Chapters 3 and 4) and use the above interfacial balances to couple the conditions in each phase. This provides information for dependent variables, such as pressure, density, mass fraction, velocity, and temperature, as functions of time and space. Such an approach requires detailed numerical simulation in all phases. One may often be interested in a restricted solution for making an order of magnitude analysis or doing an analytical solution for a limiting case, in contrast to a detailed numerical solution. In some applications, dependent variables are also coupled, preventing the solution of one dependent variable independent of the others. In practice, one often must make assumptions to simplify the interface conditions or complexity of the multiphase problem. In other cases one can neglect, based on physical significance, the resistance to mass transfer and/or heat transfer in one of the phases; this permits simplification of the solution techniques and the solution of conservation equations in only one phase.
Some dependent variables, such as temperature, are continuous across the phases; other variables, such as concentration, are discontinuous. For a general interface between phases k and j in a multicomponent system, a local balance in mass flux of species i must be upheld. The total species mass flux, , at an interface is:

The velocity of species i in phase k and phase j is
and , respectively. These velocities are defined as:


Using the interfacial species mass balance, and substituting the definition of the species velocity, the interfacial species mass flux is:

Remembering the overall mass conservation at the interface,

The interfacial species mass flux is:

In some problems, the species mass flux will be specified, and the total mass flux is simply a sum of all the species mass fluxes.

If the species mass flux is not specified, the total mass flux at an interface can be calculated from the interfacial species mass flux equation, eq. (34):

Since the species equation is second order in space, and there are two phases, there needs to be a secondary condition relating the species mass fraction in phases k and j. This can be done by expressing the mass fraction of species i in phase k as a function of the species mass fraction i in phase j, or vice versa.


Note that, in general, the species mass fraction is not a continuous function, and there is almost always a jump condition at the interface when two or more species are present. To make this point clear, the simplest case of a binary mixture is considered. For a binary mixture, the species diffusion flux, J, can be calculated by Fick’s law.


There are several examples of different phenomena of a binary mixture in the presence of a liquid/vapor, solid/vapor, or liquid/solid interface. In other examples, the species flux in one of the phases is unimportant; therefore, only one phase must be considered.
A classical example of species flux is the condition of zero mass flux due to an impermeable surface where species A does not diffuse to the stationary media or . The case of constant surface concentration is more typical; however, all types of mass transfer boundary conditions are demonstrated in three categories below.
a. The mass transfer from solid or liquid to a gas stream Evaporation and sublimation are typical examples of mass transfer from liquid or solid to a gas mixture, as shown in Fig. 1, where , x_{A,s}, x_{A,g} are molar fractions in liquid, solid, and gas, respectively, and is mass flux for species A. There are two cases shown in Fig. 1. In case I, a pure liquid/solid is evaporating/sublimating into a gaseous mixture and the evaporation /sublimation rate is controlled purely by the species gradient in the gas phase. In case II, the liquid/solid is not pure, therefore concentration gradients also exist in the liquid/solid phases. The mass transfer rate is limited by both phases. Also note that the concentration gradient in the solid [Fig. 1(b)] is steeper and does not penetrate as deeply when compared to the liquid concentration gradient [Fig. 1(a)]. Also, the liquid concentration gradient is steeper and does not penetrate as far compared to the gas. These trends are due to the ratio of the mass diffusivities in each of the phases.
A simplified relationship for the mass concentration condition at the interface can be obtained assuming the gas mixture is approximated by the ideal gas and the solid or liquid has a high concentration of A. The condition within the gas stream is of interest, because the main resistance to mass transfer is within that region. With the preceding assumptions, the partial vapor pressure of A in the gas mixture at the interface can be approximated from Raoult’s Law:


where pA is the partial vapor pressure of A in the gas mixture, is the mole fraction of species A in liquid, is the mole fraction of species in solid, and pA,sat is the normal saturation pressure of species A at the surface interface. Clearly, if the solid or liquid is made of pure species A, then
and eqs. (41) and (42) reduce to

This means that the partial vapor pressure of species A in the gas mixture at the interface is equal to the normal saturation pressure for species A, which is a function of interfacial temperature and can be obtained from thermodynamic tables in Appendix B. At the interface, species A and B in the gas phase must be in equilibrium with species A and B in the liquid or solid phase, except in extreme circumstances. Knowing and total pressure p, one can easily calculate the mole fraction x_{A,g} and mass fraction ω_{A,g} of species A at the interface on the gas side by the following relation


where M_{A} and M_{B} are molecular weights of species A and B. A conventional example for the case in Fig. 1(a) is the evaporation of water to a waterair mixture; obviously, one can imagine that the water absorbed a small amount of air. Sublimation of iodine in air is an application for the case of Fig. 1(b).
Mass transfer from a solid to a gaseous state sometimes requires the specification of diffusion molar flux, rather than concentration, at the solid surface.

where can be a function of concentration and not necessarily constant. One application is related to catalytic surface reaction. For this case, because the mass flux of reactants equals the mass flux of products. Therefore, from eq. (34),
. Catalytic surfaces are used to promote heterogeneous reactions (see Chapter 3), which occur at the surface; the appropriate boundary condition is

where k'_{1} is the reaction rate constant and n is the order of reaction.
b. Mass transfer from gases to liquids or solids
There are two major forms of mass transfer from gas to liquid, shown in Fig. 2. Case I is for condensable species in the liquid phase, or for species that can deposit vapor into the solid phase. In this case, a species of low concentration in the gas phase can condense at a higher concentration in the liquid phase. This phenomenon is important in distillation processes. Case II is for species that are weakly soluble in a liquid or a solid. In the liquid ( is small), Henry’s law will relate the mole fraction of species A in the liquid to the partial vapor pressure of A in the gas mixture at the interface by following relation [Fig. 2(a)]

Henry’s constant, H, is dependent primarily on the temperature of the aqueous solution and the absorption of two species. The pressure dependence is usually small and in fact, negligible up to pressure equal to +5 bar. Table B.77 provides H for selected aqueous solutions. For binary soluble gases, Henry’s law is not appropriate, and solubility data are usually presented in term of gasphase partial pressure to liquidphase mole fraction. Such data are given in Appendix B for NH3water and SO2water systems.
In chemical engineering applications, there are many cases in which gas is absorbed into a liquid; two examples are the absorption of hydrogen sulfide from an H2Sair mixture into liquid water and the absorption of O2 or chlorine into liquid water.
The concentration of gas in a solid at the interface is usually obtained by the use of a property known as the solubility, S, defined below

where c_{A,s} is the mole concentration of species A in the solid at the interface and is the vapor partial pressure of species A in the gas at the interface. The values of S (kmol/Pam^{3}) for vapor gassolid combinations are given in chemistry handbooks; some are in Appendix B (Tables B.76, B.78 and B.79). Special care should be exercised when using eq. (49) for the variation in form and units – two different versions of S are presented in Appendix B. An example of this case is the diffusion of helium in a glass. The dissolution of gas in metals is much more complex and depends on the type of metal used. The dissolution of gas in some metals can be reversed, such as the case of hydrogen into titanium. In contrast to hydrogen, oxygen dissolution in titanium is irreversible but is complicated because it forms a layer of scale T_{iO2} on the surface. This topic is beyond the scope of this book, and interested readers should refer to a metallurgy or materials handbook to consult appropriate phase diagrams.
c. Mass transit from solid to liquid or liquid to solid Two cases are also presented for the case of mass transfer between a solid and a liquid (Fig. 3). Case I represents a pure solid dissolving into a liquid, or a pure liquid diffusing into a solid. Case II represents a solid mixture melting into a liquid mixture, or a liquid mixture solidifying into a solid mixture. For the first case, the mass transfer is related to the solubility, which can be found in chemistry handbooks (Lide, 2004) and some are reproduced in Appendix B. A conventional example for this case is the dissolution of salt into water.
For a liquid/solid interface during melting and solidification, the ratio of the species concentration of the solid in liquid phases is called the partition ratio, K_{p}.

The partition coefficient is a function of temperature. Furthermore, when the liquid and solid lines are nearly straight, it is a constant. This information can be obtained from phase diagrams for different solid/liquid mixtures.
Additional Equations

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