Evaporation from adiabatic walls
From ThermalFluidsPedia
Yuwen Zhang (Talk  contribs) 

(5 intermediate revisions not shown)  
Line 1:  Line 1:  
  Evaporation from thin liquid films occurs in many industrial and natural processes, including drying, evaporative cooling, and sweating.  +  [[Image:NewChapter9 (9).jpgthumb400 pxalt= Evaporation from a thin liquid film on a horizontal surface.  '''Evaporation from a thin liquid film on a horizontal surface.''' ]] 
+  
+  Evaporation from thin liquid films occurs in many industrial and natural processes, including drying, evaporative cooling, and sweating. The figure on the right shows the physical model of horizontal thin film evaporation under consideration [[#References(Carey, 1992)]]. A film wets a horizontal surface over which flows a gas of ambient temperature <math>{{T}_{\infty }}</math> and mass fraction of water <math>{{\omega }_{\infty }}</math> at a velocity of <math>{{u}_{\infty }}</math>. Since the solid surface underneath the liquid is adiabatic, the latent heat of vaporization is provided by the vaporgas mixture flowing above the liquid film. The heat source (vaporgas mixture) and the evaporating liquid are in direct contact, making this scenario an example of direct contact evaporation. The liquid is evaporated and the latent heat is absorbed in the evaporation process. The resulting vapor injects into the boundary layer and is removed by the gas flow. The boundary layer becomes thicker, and a free stream of the gas flow is displaced from the surface being cooled. While phase change is the dominant mechanism of heat transfer at lower gas temperatures, vapor injection becomes more important at higher gas temperatures that correspond to higher evaporation rates.  
Due to diffusion and convection, the temperature rises and the concentration falls from their values at the film surface to their ambient values at the edge of the boundary layer. A nonslip boundary condition exists at the film surface, and the film itself is considered stationary with respect to the gas. Assuming a steady state, constant density, and incompressible flow, the continuity equation is  Due to diffusion and convection, the temperature rises and the concentration falls from their values at the film surface to their ambient values at the edge of the boundary layer. A nonslip boundary condition exists at the film surface, and the film itself is considered stationary with respect to the gas. Assuming a steady state, constant density, and incompressible flow, the continuity equation is  
  <center><math>\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0</math></center>  
+  {{EquationRef(1)}}  
+  }  
+  
With no pressure gradients and constant viscosity assumed, the boundary layer momentum equation is written as  With no pressure gradients and constant viscosity assumed, the boundary layer momentum equation is written as  
  <center><math>u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\tilde{\nu }\frac{{{\partial }^{2}}u}{\partial {{y}^{2}}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\tilde{\nu }\frac{{{\partial }^{2}}u}{\partial {{y}^{2}}}</math></center>  
+  {{EquationRef(2)}}  
+  }  
  
where ~ on top of viscosity signifies massaveraged properties of the mixture.  where ~ on top of viscosity signifies massaveraged properties of the mixture.  
Neglecting viscous dissipation and assuming constant thermal diffusivity and specific heats, the boundary layer energy balance is  Neglecting viscous dissipation and assuming constant thermal diffusivity and specific heats, the boundary layer energy balance is  
  <center><math>u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\tilde{\alpha }\frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\tilde{\alpha }\frac{{{\partial }^{2}}T}{\partial {{y}^{2}}}</math></center>  
+  {{EquationRef(3)}}  
+  }  
+  
where the axial conduction on the xdirection has been neglected because heat transfer occurred mainly in the ydirection.  where the axial conduction on the xdirection has been neglected because heat transfer occurred mainly in the ydirection.  
The mass fraction of water is also accounted for by assuming a constant mass diffusivity:  The mass fraction of water is also accounted for by assuming a constant mass diffusivity:  
  <center><math>u\frac{\partial \omega }{\partial x}+v\frac{\partial \omega }{\partial y}=D\frac{{{\partial }^{2}}\omega }{\partial {{y}^{2}}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>u\frac{\partial \omega }{\partial x}+v\frac{\partial \omega }{\partial y}=D\frac{{{\partial }^{2}}\omega }{\partial {{y}^{2}}}</math></center>  
+  {{EquationRef(4)}}  
+  }  
  As <math>y\to \infty </math>the boundary conditions can be taken directly from  +  As <math>y\to \infty </math>the boundary conditions can be taken directly from the above figure: 
  <math>u\to {{u}_{\infty }}</math>,  +  { class="wikitable" border="0" 
  <math>T\to {{T}_{\infty }}</math>,  +   
  +   width="100%"   
  +  <center><math>u\to {{u}_{\infty }}</math>,  
+  <math>T\to {{T}_{\infty }}</math>, <math>\omega \to {{\omega }_{\infty }}</math></center>  
+  {{EquationRef(5)}}  
+  }  
The boundary conditions are now laid out at the film surface ''y'' = 0. There is a nonslip condition, and the mass flux <math>{\dot{m}}''</math> evaporates normal to the surface.  The boundary conditions are now laid out at the film surface ''y'' = 0. There is a nonslip condition, and the mass flux <math>{\dot{m}}''</math> evaporates normal to the surface.  
  <math>u=0</math>,  +  { class="wikitable" border="0" 
+    
+   width="100%"   
+  <center><math>u=0</math>,  
<math>v=\frac{{{\dot{m}}''}}{{\tilde{\rho }}}</math>,  <math>v=\frac{{{\dot{m}}''}}{{\tilde{\rho }}}</math>,  
  <math>\omega ={{\omega }_{\delta }}</math>,  +  <math>\omega ={{\omega }_{\delta }}</math>, <math>T={{T}_{\delta }}={{T}_{sat}}\left( {{\omega }_{\delta }},{{p}_{\infty }} \right)</math></center> 
  +  {{EquationRef(6)}}  
  +  }  
  +  
An energy balance at the interface accounts for the heat of vaporization and mass flux of the evaporating fluid. With the nonslip boundary condition, the only mode of heat transfer is conduction:  An energy balance at the interface accounts for the heat of vaporization and mass flux of the evaporating fluid. With the nonslip boundary condition, the only mode of heat transfer is conduction:  
  <  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{\dot{m}}''{{h}_{\ell v}}={q}''=\tilde{k}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}</math></center>  
+  {{EquationRef(7)}}  
+  }  
  Substituting eq. (  +  Substituting eq. (7) into eq. (6), the velocity component in the ydirection becomes 
  <center><math>{{\left. v \right}_{y=0}}=\frac{{\tilde{k}}}{\rho {{h}_{\ell v}}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{\left. v \right}_{y=0}}=\frac{{\tilde{k}}}{\rho {{h}_{\ell v}}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}</math></center>  
+  {{EquationRef(8)}}  
+  }  
The mass flux at the interface is the result of both diffusion and convection. Therefore, the mass balance at the interface is written as  The mass flux at the interface is the result of both diffusion and convection. Therefore, the mass balance at the interface is written as  
  <center><math>{\dot{m}}''=\rho D{{\left( \frac{\partial \omega }{\partial y} \right)}_{y=0}}+\rho {{\omega }_{\delta }}{{\left. v \right}_{y=0}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{\dot{m}}''=\rho D{{\left( \frac{\partial \omega }{\partial y} \right)}_{y=0}}+\rho {{\omega }_{\delta }}{{\left. v \right}_{y=0}}</math></center>  
+  {{EquationRef(9)}}  
+  }  
  Substituting eq. (  +  Substituting eq. (7) into eq. (9), conservation of mass at the interface becomes 
  <center><math>\tilde{\rho }D{{\left( \frac{\partial \omega }{\partial y} \right)}_{y=0}}+\tilde{\rho }{{\omega }_{\delta }}{{\left. v \right}_{y=0}}=\frac{{\tilde{k}}}{{{h}_{\ell v}}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\tilde{\rho }D{{\left( \frac{\partial \omega }{\partial y} \right)}_{y=0}}+\tilde{\rho }{{\omega }_{\delta }}{{\left. v \right}_{y=0}}=\frac{{\tilde{k}}}{{{h}_{\ell v}}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}</math></center>  
+  {{EquationRef(10)}}  
+  }  
  Equations (  +  Equations (1) – (2) are the same as the laminar boundary layer equations for forced convection over a flat plate with blowing on the liquid surface. By defining the stream function <math>\psi </math> as follows, 
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>u=\frac{\partial \psi }{\partial y}\begin{matrix}  <center><math>u=\frac{\partial \psi }{\partial y}\begin{matrix}  
{} & v=\frac{\partial \psi }{\partial x} \\  {} & v=\frac{\partial \psi }{\partial x} \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
  +  {{EquationRef(11)}}  
+  }  
  the continuity equation (  +  the continuity equation (1) is automatically satisfied. The momentum equation in terms of the stream function is then 
  <center><math>\frac{\partial \psi }{\partial y}\frac{{{\partial }^{2}}\psi }{\partial x\partial y}\frac{\partial \psi }{\partial x}\frac{{{\partial }^{2}}\psi }{\partial {{y}^{2}}}=\nu \frac{{{\partial }^{3}}\psi }{\partial {{y}^{3}}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{\partial \psi }{\partial y}\frac{{{\partial }^{2}}\psi }{\partial x\partial y}\frac{\partial \psi }{\partial x}\frac{{{\partial }^{2}}\psi }{\partial {{y}^{2}}}=\nu \frac{{{\partial }^{3}}\psi }{\partial {{y}^{3}}}</math></center>  
+  {{EquationRef(12)}}  
+  }  
Introducing the following similarity variable:  Introducing the following similarity variable:  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>\eta =y\sqrt{\frac{{{u}_{\infty }}}{\nu x}}\begin{matrix}  <center><math>\eta =y\sqrt{\frac{{{u}_{\infty }}}{\nu x}}\begin{matrix}  
, & f=\frac{\psi }{\sqrt{\nu {{u}_{\infty }}x}} \\  , & f=\frac{\psi }{\sqrt{\nu {{u}_{\infty }}x}} \\  
\end{matrix},\text{ }\theta =\frac{T{{T}_{\delta }}}{{{T}_{\infty }}{{T}_{\delta }}}\begin{matrix}  \end{matrix},\text{ }\theta =\frac{T{{T}_{\delta }}}{{{T}_{\infty }}{{T}_{\delta }}}\begin{matrix}  
, & \varphi =\frac{\omega {{\omega }_{\delta }}}{{{\omega }_{\infty }}{{\omega }_{\delta }}} \\  , & \varphi =\frac{\omega {{\omega }_{\delta }}}{{{\omega }_{\infty }}{{\omega }_{\delta }}} \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
  +  {{EquationRef(13)}}  
+  }  
  eq. (  +  eq. (12) and eqs. (3) – (4) can be reduced to a set of ordinary differential equations: 
  <center><math>{f}'''+\frac{1}{2}f{f}''=0</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{f}'''+\frac{1}{2}f{f}''=0</math></center>  
+  {{EquationRef(14)}}  
+  }  
  +  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{\theta }''+\frac{1}{2}\Pr f{\theta }'=0</math></center>  
+  {{EquationRef(15)}}  
+  }  
  +  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{\varphi }''+\frac{1}{2}\text{Sc}f{\varphi }'=0</math></center>  
+  {{EquationRef(16)}}  
+  }  
  The velocity components in the x and ydirections in terms of the dimensionless stream function f can be obtained by eqs. (  +  The velocity components in the x and ydirections in terms of the dimensionless stream function f can be obtained by eqs. (11) and (13), i.e., 
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>u={{u}_{\infty }}{f}'(\eta )\begin{matrix}  <center><math>u={{u}_{\infty }}{f}'(\eta )\begin{matrix}  
, & v=\frac{1}{2}\sqrt{\frac{\nu {{u}_{\infty }}}{x}} \\  , & v=\frac{1}{2}\sqrt{\frac{\nu {{u}_{\infty }}}{x}} \\  
  \end{matrix}\left( \eta {f}'f \right)</math>  +  \end{matrix}\left( \eta {f}'f \right)</math></center> 
  +  {{EquationRef(17)}}  
+  }  
  The boundary conditions at y→∞ represented by eq. (  +  The boundary conditions at y→∞ represented by eq. (5) become 
  +  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>f'(\infty )=1\begin{matrix}  
, & \theta (\infty )=1, & \varphi (\infty )=1 \\  , & \theta (\infty )=1, & \varphi (\infty )=1 \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
+  {{EquationRef(18)}}  
+  }  
  The boundary conditions represented by eq. (  +  The boundary conditions represented by eq. (6), except the velocity component in the ydirection, can be rewritten to 
  +  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>f'(0)=0\begin{matrix}  <center><math>f'(0)=0\begin{matrix}  
, & \theta (0)=0, & \varphi (0)=0 \\  , & \theta (0)=0, & \varphi (0)=0 \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
+  {{EquationRef(19)}}  
+  }  
  Conservation of mass and energy at the interface, eqs. (  +  Conservation of mass and energy at the interface, eqs. (8) and (10), can be rewritten in terms of the similarity variables, i.e., 
  +  
  <center><math>f(0)=\frac{2\text{Ja}}{\Pr }{\theta }'(0)</math>  +  { class="wikitable" border="0" 
  +    
  +   width="100%"   
  <center><math>\frac{{{\omega }_{\delta }}{{\omega }_{\infty }}}{1{{\omega }_{\delta }}}{\phi }'(0)=\frac{\text{Sc}}{\Pr }Ja{\theta }'(0)</math>  +  <center><math>f(0)=\frac{2\text{Ja}}{\Pr }{\theta }'(0)</math></center> 
  +  {{EquationRef(20)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\frac{{{\omega }_{\delta }}{{\omega }_{\infty }}}{1{{\omega }_{\delta }}}{\phi }'(0)=\frac{\text{Sc}}{\Pr }Ja{\theta }'(0)</math></center>  
+  {{EquationRef(21)}}  
+  }  
where ''Ja'' is the Jakob number and is defined as  where ''Ja'' is the Jakob number and is defined as  
  +  
  <center><math>\text{Ja}=\frac{{{c}_{p}}\left[ {{T}_{\infty }}{{T}_{sat}}({{\omega }_{\delta }}) \right]}{{{h}_{\ell v}}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\text{Ja}=\frac{{{c}_{p}}\left[ {{T}_{\infty }}{{T}_{sat}}({{\omega }_{\delta }}) \right]}{{{h}_{\ell v}}}</math></center>  
+  {{EquationRef(22)}}  
+  }  
and Sc is the Schmidt number defined as  and Sc is the Schmidt number defined as  
  
  
  
  +  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\text{Sc}=\frac{{\tilde{\nu }}}{D}</math></center>  
+  {{EquationRef(23)}}  
+  }  
  The boundary value problem can be solved using the RungeKutta method in conjunction with a shooting method. The solution procedure begins with an assumed ω<sub>I</sub> and uses boundary conditions specified by eqs. (  +  The evaporation problem is now described by a set of ordinary differential equations (14) – (16) subjected to boundary conditions specified by eqs. (18) – (21). Since eq. (14) is a thirdorder ordinary differential equation, it requires three boundary conditions: in eq. (18), in eq. (19), and eq. (20). Equations (15) and (16) are both secondorder ordinary differential equations, each of them requiring two boundary conditions, which are specified in eqs. (18) and (19). Therefore, the problem is mathematically defined with boundary conditions stated by eqs. (18) – (20), which makes the boundary condition specified in eq. (21) an extra boundary condition – and makes the problem overstated. This happened because the heat and mass transfer are not independent of each other and the energy balance at the interface requires that eq. (21) be satisfied. 
+  
+  The boundary value problem can be solved using the RungeKutta method in conjunction with a shooting method. The solution procedure begins with an assumed ω<sub>I</sub> and uses boundary conditions specified by eqs. (18) – (20) to solve eqs. (14) – (16). Once a solution is obtained, eq. (21) is employed to find <math>{{\omega }_{\delta }}</math>. If <math>{{\omega }_{\delta }}</math> obtained from eq. (21) agrees with the assumed value, the solution is complete. Otherwise, the assumed value of <math>{{\omega }_{\delta }}</math> is corrected and the solution procedure is repeated until a converged solution is obtained.  
Once the converged solution is obtained, the local heat transfer coefficient can be evaluated by  Once the converged solution is obtained, the local heat transfer coefficient can be evaluated by  
  +  
  <center><math>{{h}_{x}}=\frac{{\tilde{k}}}{{{T}_{\infty }}{{T}_{\delta }}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}=\tilde{k}\sqrt{\frac{{{u}_{\infty }}}{\tilde{\nu }x}}{\theta }'(0)</math>  +  
  +  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{{h}_{x}}=\frac{{\tilde{k}}}{{{T}_{\infty }}{{T}_{\delta }}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}=\tilde{k}\sqrt{\frac{{{u}_{\infty }}}{\tilde{\nu }x}}{\theta }'(0)</math></center>  
+  {{EquationRef(24)}}  
+  }  
The local Nusselt number is then  The local Nusselt number is then  
  +  
  <center><math>N{{u}_{x}}=\frac{{{h}_{x}}x}{k}={\theta }'(0)\operatorname{Re}_{x}^{\frac{1}{2}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>N{{u}_{x}}=\frac{{{h}_{x}}x}{k}={\theta }'(0)\operatorname{Re}_{x}^{\frac{1}{2}}</math></center>  
+  {{EquationRef(25)}}  
+  }  
The local mass transfer coefficient can be evaluated by  The local mass transfer coefficient can be evaluated by  
  +  
  <center><math>{{h}_{mx}}=\frac{D}{{{\omega }_{\infty }}{{\omega }_{\delta }}}{{\left( \frac{\partial \omega }{\partial y} \right)}_{y=0}}=D\sqrt{\frac{{{u}_{\infty }}}{\nu x}}{\varphi }'(0)</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{h}_{mx}}=\frac{D}{{{\omega }_{\infty }}{{\omega }_{\delta }}}{{\left( \frac{\partial \omega }{\partial y} \right)}_{y=0}}=D\sqrt{\frac{{{u}_{\infty }}}{\nu x}}{\varphi }'(0)</math></center>  
+  {{EquationRef(26)}}  
+  }  
and the local Sherwood number is  and the local Sherwood number is  
  <center><math>S{{h}_{x}}=\frac{{{h}_{m}}x}{D}={\varphi }'(0)\operatorname{Re}_{x}^{\frac{1}{2}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>S{{h}_{x}}=\frac{{{h}_{m}}x}{D}={\varphi }'(0)\operatorname{Re}_{x}^{\frac{1}{2}}</math></center>  
+  {{EquationRef(27)}}  
+  }  
where  where  
  <center><math>{{\operatorname{Re}}_{x}}=\frac{{{u}_{\infty }}x}{\nu }</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{\operatorname{Re}}_{x}}=\frac{{{u}_{\infty }}x}{\nu }</math></center>  
+  {{EquationRef(28)}}  
+  }  
  Although eqs. (  +  Although eqs. (25) and (27) have the same form as the case of singlephase forced convective heat or mass transfer over a flat plate, the solutions of <math>{\theta }'(0)</math> and <math>{\varphi }'(0)</math> depend on ''Pr'', ''Sc'', ''Ja'', and ''ω<sub>∞</sub>''. Consequently, simple empirical correlations similar to the case of singlephase forced convective heat or mass transfer over a flat plate are very difficult to obtain. 
  The evaporation problem presented here is a coupled problem because the solutions of the momentum, energy, and species equations are coupled with eqs. (  +  The evaporation problem presented here is a coupled problem because the solutions of the momentum, energy, and species equations are coupled with eqs. (20) – (21). Since the blowing velocity at the surface of the liquid film, determined by eq. (8), is usually much lower than the incoming vaporgas mixture velocity, <math>{{u}_{\infty }}</math>, one can assume that the blowing velocity at the liquid surface is negligible. Then eq. (20) can be replaced with 
  +  
  +  
  +  
  +  { class="wikitable" border="0"  
  +    
  <center><math>  +   width="100%"  
  +  <center><math>f(0)=0</math></center>  
  +  {{EquationRef(29)}}  
  +  }  
  +  
  +  The heat and mass transfer problem described by eqs. (14) – (16), (18) – (19) and (29) then cease to be a conjugated problem, because each equation in eq. (14) – (16) can be solved independently. Numerical solutions of these ordinary differential equations yield the following results (Kays et al., 2004):  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{\theta }'(0)=0.332{{\Pr }^{\frac{1}{3}}}</math></center>  
+  {{EquationRef(30)}}  
+  }  
  <center><math>N{{u}_{x}}=\frac{{{h}_{x}}x}{k}=0.332\operatorname{Re}_{x}^{\frac{1}{2}}{{\Pr }^{\frac{1}{3}}}</math>  +  { class="wikitable" border="0" 
  +    
  +   width="100%"   
  <center><math>S{{h}_{x}}=\frac{{{h}_{m}}x}{k}=0.332\operatorname{Re}_{x}^{\frac{1}{2}}S{{c}^{\frac{1}{3}}}</math>  +  <center><math>{\varphi }'(0)=0.332S{{c}^{\frac{1}{3}}}</math></center> 
  +  {{EquationRef(31)}}  
+  }  
+  
+  Substituting eqs. (30) and (31) into eqs. (25) and (27) yields  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>N{{u}_{x}}=\frac{{{h}_{x}}x}{k}=0.332\operatorname{Re}_{x}^{\frac{1}{2}}{{\Pr }^{\frac{1}{3}}}</math></center>  
+  {{EquationRef(32)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>S{{h}_{x}}=\frac{{{h}_{m}}x}{k}=0.332\operatorname{Re}_{x}^{\frac{1}{2}}S{{c}^{\frac{1}{3}}}</math></center>  
+  {{EquationRef(33)}}  
+  }  
The average heat and mass transfer characteristics are often relevant for practical applications. The average Nusselt and Sherwood numbers based on average heat and mass transfer coefficients for a liquid film with a length of ''L'' are obtained by  The average heat and mass transfer characteristics are often relevant for practical applications. The average Nusselt and Sherwood numbers based on average heat and mass transfer coefficients for a liquid film with a length of ''L'' are obtained by  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\overline{Nu}=\frac{\bar{h}L}{k}=0.664\operatorname{Re}_{L}^{\frac{1}{2}}{{\Pr }^{\frac{1}{3}}}</math></center>  
+  {{EquationRef(34)}}  
+  }  
  +  { class="wikitable" border="0"  
  +    
  +   width="100%"   
  <center><math>\overline{Sh}=\frac{{{{\bar{h}}}_{m}}L}{k}=0.664\operatorname{Re}_{L}^{\frac{1}{2}}S{{c}^{\frac{1}{3}}}</math>  +  <center><math>\overline{Sh}=\frac{{{{\bar{h}}}_{m}}L}{k}=0.664\operatorname{Re}_{L}^{\frac{1}{2}}S{{c}^{\frac{1}{3}}}</math></center> 
  +  {{EquationRef(35)}}  
+  }  
where  where  
  <center><math>{{\operatorname{Re}}_{L}}=\frac{{{u}_{\infty }}L}{\nu }</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{\operatorname{Re}}_{L}}=\frac{{{u}_{\infty }}L}{\nu }</math></center>  
+  {{EquationRef(36)}}  
+  }  
  The mass fraction of water at the liquid surface can be obtained by substituting eqs. (  +  The mass fraction of water at the liquid surface can be obtained by substituting eqs. (30) and (31) into eq. (21), i.e., 
  <center><math>{{\omega }_{\delta }}=\frac{{{\omega }_{\infty }}+{{(\text{Sc}/\Pr )}^{2/3}}\text{Ja}}{1+{{(\text{Sc}/\Pr )}^{2/3}}\text{Ja}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{\omega }_{\delta }}=\frac{{{\omega }_{\infty }}+{{(\text{Sc}/\Pr )}^{2/3}}\text{Ja}}{1+{{(\text{Sc}/\Pr )}^{2/3}}\text{Ja}}</math></center>  
+  {{EquationRef(37)}}  
+  }  
The applications of the above empirical correlations can be demonstrated by the following example.  The applications of the above empirical correlations can be demonstrated by the following example.  
Line 199:  Line 353:  
Carey, V.P., 1992, LiquidVapor PhaseChange Phenomena: An Introduction to the Thermophysics of Vaporization and Condensation Processes in Heat Transfer Equipment, Hemisphere Publishing Corp., Washington, D. C.  Carey, V.P., 1992, LiquidVapor PhaseChange Phenomena: An Introduction to the Thermophysics of Vaporization and Condensation Processes in Heat Transfer Equipment, Hemisphere Publishing Corp., Washington, D. C.  
+  
+  Faghri, A., and Zhang, Y., 2006, ''Transport Phenomena in Multiphase Systems'', Elsevier, Burlington, MA  
+  
+  Faghri, A., Zhang, Y., and Howell, J. R., 2010, ''Advanced Heat and Mass Transfer'', Global Digital Press, Columbia, MO. 
Current revision as of 08:21, 7 July 2010
Evaporation from thin liquid films occurs in many industrial and natural processes, including drying, evaporative cooling, and sweating. The figure on the right shows the physical model of horizontal thin film evaporation under consideration (Carey, 1992). A film wets a horizontal surface over which flows a gas of ambient temperature and mass fraction of water at a velocity of . Since the solid surface underneath the liquid is adiabatic, the latent heat of vaporization is provided by the vaporgas mixture flowing above the liquid film. The heat source (vaporgas mixture) and the evaporating liquid are in direct contact, making this scenario an example of direct contact evaporation. The liquid is evaporated and the latent heat is absorbed in the evaporation process. The resulting vapor injects into the boundary layer and is removed by the gas flow. The boundary layer becomes thicker, and a free stream of the gas flow is displaced from the surface being cooled. While phase change is the dominant mechanism of heat transfer at lower gas temperatures, vapor injection becomes more important at higher gas temperatures that correspond to higher evaporation rates.
Due to diffusion and convection, the temperature rises and the concentration falls from their values at the film surface to their ambient values at the edge of the boundary layer. A nonslip boundary condition exists at the film surface, and the film itself is considered stationary with respect to the gas. Assuming a steady state, constant density, and incompressible flow, the continuity equation is

With no pressure gradients and constant viscosity assumed, the boundary layer momentum equation is written as

where ~ on top of viscosity signifies massaveraged properties of the mixture. Neglecting viscous dissipation and assuming constant thermal diffusivity and specific heats, the boundary layer energy balance is

where the axial conduction on the xdirection has been neglected because heat transfer occurred mainly in the ydirection. The mass fraction of water is also accounted for by assuming a constant mass diffusivity:

As the boundary conditions can be taken directly from the above figure:

The boundary conditions are now laid out at the film surface y = 0. There is a nonslip condition, and the mass flux evaporates normal to the surface.
, ω = ω_{δ}, 
An energy balance at the interface accounts for the heat of vaporization and mass flux of the evaporating fluid. With the nonslip boundary condition, the only mode of heat transfer is conduction:

Substituting eq. (7) into eq. (6), the velocity component in the ydirection becomes

The mass flux at the interface is the result of both diffusion and convection. Therefore, the mass balance at the interface is written as

Substituting eq. (7) into eq. (9), conservation of mass at the interface becomes

Equations (1) – (2) are the same as the laminar boundary layer equations for forced convection over a flat plate with blowing on the liquid surface. By defining the stream function ψ as follows,

the continuity equation (1) is automatically satisfied. The momentum equation in terms of the stream function is then

Introducing the following similarity variable:

eq. (12) and eqs. (3) – (4) can be reduced to a set of ordinary differential equations:



The velocity components in the x and ydirections in terms of the dimensionless stream function f can be obtained by eqs. (11) and (13), i.e.,

The boundary conditions at y→∞ represented by eq. (5) become

The boundary conditions represented by eq. (6), except the velocity component in the ydirection, can be rewritten to

Conservation of mass and energy at the interface, eqs. (8) and (10), can be rewritten in terms of the similarity variables, i.e.,


where Ja is the Jakob number and is defined as

and Sc is the Schmidt number defined as

The evaporation problem is now described by a set of ordinary differential equations (14) – (16) subjected to boundary conditions specified by eqs. (18) – (21). Since eq. (14) is a thirdorder ordinary differential equation, it requires three boundary conditions: in eq. (18), in eq. (19), and eq. (20). Equations (15) and (16) are both secondorder ordinary differential equations, each of them requiring two boundary conditions, which are specified in eqs. (18) and (19). Therefore, the problem is mathematically defined with boundary conditions stated by eqs. (18) – (20), which makes the boundary condition specified in eq. (21) an extra boundary condition – and makes the problem overstated. This happened because the heat and mass transfer are not independent of each other and the energy balance at the interface requires that eq. (21) be satisfied.
The boundary value problem can be solved using the RungeKutta method in conjunction with a shooting method. The solution procedure begins with an assumed ω_{I} and uses boundary conditions specified by eqs. (18) – (20) to solve eqs. (14) – (16). Once a solution is obtained, eq. (21) is employed to find ω_{δ}. If ω_{δ} obtained from eq. (21) agrees with the assumed value, the solution is complete. Otherwise, the assumed value of ω_{δ} is corrected and the solution procedure is repeated until a converged solution is obtained.
Once the converged solution is obtained, the local heat transfer coefficient can be evaluated by

The local Nusselt number is then

The local mass transfer coefficient can be evaluated by

and the local Sherwood number is

where

Although eqs. (25) and (27) have the same form as the case of singlephase forced convective heat or mass transfer over a flat plate, the solutions of θ'(0) and depend on Pr, Sc, Ja, and ω_{∞}. Consequently, simple empirical correlations similar to the case of singlephase forced convective heat or mass transfer over a flat plate are very difficult to obtain.
The evaporation problem presented here is a coupled problem because the solutions of the momentum, energy, and species equations are coupled with eqs. (20) – (21). Since the blowing velocity at the surface of the liquid film, determined by eq. (8), is usually much lower than the incoming vaporgas mixture velocity, , one can assume that the blowing velocity at the liquid surface is negligible. Then eq. (20) can be replaced with

The heat and mass transfer problem described by eqs. (14) – (16), (18) – (19) and (29) then cease to be a conjugated problem, because each equation in eq. (14) – (16) can be solved independently. Numerical solutions of these ordinary differential equations yield the following results (Kays et al., 2004):


Substituting eqs. (30) and (31) into eqs. (25) and (27) yields


The average heat and mass transfer characteristics are often relevant for practical applications. The average Nusselt and Sherwood numbers based on average heat and mass transfer coefficients for a liquid film with a length of L are obtained by


where

The mass fraction of water at the liquid surface can be obtained by substituting eqs. (30) and (31) into eq. (21), i.e.,

The applications of the above empirical correlations can be demonstrated by the following example.
References
Carey, V.P., 1992, LiquidVapor PhaseChange Phenomena: An Introduction to the Thermophysics of Vaporization and Condensation Processes in Heat Transfer Equipment, Hemisphere Publishing Corp., Washington, D. C.
Faghri, A., and Zhang, Y., 2006, Transport Phenomena in Multiphase Systems, Elsevier, Burlington, MA
Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.