# Evaporation from adiabatic walls

### From Thermal-FluidsPedia

(Created page with 'Evaporation from thin liquid films occurs in many industrial and natural processes, including drying, evaporative cooling, and sweating. Figure 9.6 shows the physical model of ho…') |
|||

Line 20: | Line 20: | ||

<center><math>u\frac{\partial \omega }{\partial x}+v\frac{\partial \omega }{\partial y}=D\frac{{{\partial }^{2}}\omega }{\partial {{y}^{2}}}</math> | <center><math>u\frac{\partial \omega }{\partial x}+v\frac{\partial \omega }{\partial y}=D\frac{{{\partial }^{2}}\omega }{\partial {{y}^{2}}}</math> | ||

(9.4)</center> | (9.4)</center> | ||

- | As | + | |

- | <math>y\to \infty </math> | + | As <math>y\to \infty </math>the boundary conditions can be taken directly from Fig. 9.6: |

- | + | ||

- | <math>u\to {{u}_{\infty }}</math> | + | <math>u\to {{u}_{\infty }}</math>, |

- | , | + | <math>T\to {{T}_{\infty }}</math>, |

- | <math>T\to {{T}_{\infty }}</math> | + | |

- | , | + | |

<center><math>\omega \to {{\omega }_{\infty }}</math> | <center><math>\omega \to {{\omega }_{\infty }}</math> | ||

(9.5)</center> | (9.5)</center> | ||

- | 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> | + | <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> | + | |

- | , | + | |

<center><math>T={{T}_{\delta }}={{T}_{sat}}\left( {{\omega }_{\delta }},{{p}_{\infty }} \right)</math> | <center><math>T={{T}_{\delta }}={{T}_{sat}}\left( {{\omega }_{\delta }},{{p}_{\infty }} \right)</math> | ||

(9.6)</center> | (9.6)</center> | ||

+ | |||

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: | ||

</center><math>{\dot{m}}''{{h}_{\ell v}}=-{q}''=\tilde{k}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}</math> | </center><math>{\dot{m}}''{{h}_{\ell v}}=-{q}''=\tilde{k}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}</math> | ||

(9.7)<center> | (9.7)<center> | ||

+ | |||

Substituting eq. (9.7) into eq. (9.6), the velocity component in the y-direction becomes | Substituting eq. (9.7) into eq. (9.6), the velocity component in the y-direction 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> | <center><math>{{\left. v \right|}_{y=0}}=\frac{{\tilde{k}}}{\rho {{h}_{\ell v}}}{{\left( \frac{\partial T}{\partial y} \right)}_{y=0}}</math> | ||

(9.8)</center> | (9.8)</center> | ||

+ | |||

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> | <center><math>{\dot{m}}''=-\rho D{{\left( \frac{\partial \omega }{\partial y} \right)}_{y=0}}+\rho {{\omega }_{\delta }}{{\left. v \right|}_{y=0}}</math> | ||

(9.9)</center> | (9.9)</center> | ||

+ | |||

Substituting eq. (9.7) into eq. (9.9), conservation of mass at the interface becomes | Substituting eq. (9.7) into eq. (9.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> | <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> | ||

(9.10)</center> | (9.10)</center> | ||

+ | |||

Equations (9.1) – (9.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, | Equations (9.1) – (9.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, | ||

Line 62: | Line 63: | ||

\end{matrix}</math> | \end{matrix}</math> | ||

(9.11)</center> | (9.11)</center> | ||

+ | |||

the continuity equation (9.1) is automatically satisfied. The momentum equation in terms of the stream function is then | the continuity equation (9.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> | <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> | ||

(9.12)</center> | (9.12)</center> | ||

+ | |||

Introducing the following similarity variable: | Introducing the following similarity variable: | ||

Line 74: | Line 77: | ||

\end{matrix}</math> | \end{matrix}</math> | ||

(9.13)</center> | (9.13)</center> | ||

+ | |||

eq. (9.12) and eqs. (9.3) – (9.4) can be reduced to a set of ordinary differential equations: | eq. (9.12) and eqs. (9.3) – (9.4) can be reduced to a set of ordinary differential equations: | ||

<center><math>{f}'''+\frac{1}{2}f{f}''=0</math> | <center><math>{f}'''+\frac{1}{2}f{f}''=0</math> | ||

(9.14)</center> | (9.14)</center> | ||

+ | |||

<center><math>{\theta }''+\frac{1}{2}\Pr f{\theta }'=0</math> (9.15)</center> | <center><math>{\theta }''+\frac{1}{2}\Pr f{\theta }'=0</math> (9.15)</center> | ||

+ | |||

<center><math>{\varphi }''+\frac{1}{2}\text{Sc}f{\varphi }'=0</math> (9.16)</center> | <center><math>{\varphi }''+\frac{1}{2}\text{Sc}f{\varphi }'=0</math> (9.16)</center> | ||

+ | |||

The velocity components in the x- and y-directions in terms of the dimensionless stream function f can be obtained by eqs. (9.11) and (9.13), i.e., | The velocity components in the x- and y-directions in terms of the dimensionless stream function f can be obtained by eqs. (9.11) and (9.13), i.e., | ||

Line 86: | Line 93: | ||

\end{matrix}\left( \eta {f}'-f \right)</math> | \end{matrix}\left( \eta {f}'-f \right)</math> | ||

(9.17)</center> | (9.17)</center> | ||

- | The boundary conditions at | + | |

+ | The boundary conditions at y→∞ represented by eq. (9.5) become | ||

<center><math>f'(\infty )=1\begin{matrix} | <center><math>f'(\infty )=1\begin{matrix} | ||

, & \theta (\infty )=1, & \varphi (\infty )=1 \\ | , & \theta (\infty )=1, & \varphi (\infty )=1 \\ | ||

\end{matrix}</math> (9.18)</center> | \end{matrix}</math> (9.18)</center> | ||

+ | |||

The boundary conditions represented by eq. (9.6), except the velocity component in the y-direction, can be rewritten to | The boundary conditions represented by eq. (9.6), except the velocity component in the y-direction, can be rewritten to | ||

- | <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> (9.19)</center> | \end{matrix}</math> (9.19)</center> | ||

+ | |||

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

Line 101: | Line 112: | ||

<center><math>\frac{{{\omega }_{\delta }}-{{\omega }_{\infty }}}{1-{{\omega }_{\delta }}}{\phi }'(0)=\frac{\text{Sc}}{\Pr }Ja{\theta }'(0)</math> | <center><math>\frac{{{\omega }_{\delta }}-{{\omega }_{\infty }}}{1-{{\omega }_{\delta }}}{\phi }'(0)=\frac{\text{Sc}}{\Pr }Ja{\theta }'(0)</math> | ||

(9.21)</center> | (9.21)</center> | ||

- | 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> | <center><math>\text{Ja}=\frac{{{c}_{p}}\left[ {{T}_{\infty }}-{{T}_{sat}}({{\omega }_{\delta }}) \right]}{{{h}_{\ell v}}}</math> | ||

(9.22)</center> | (9.22)</center> | ||

+ | |||

and Sc is the Schmidt number defined as | and Sc is the Schmidt number defined as | ||

<center><math>\text{Sc}=\frac{{\tilde{\nu }}}{D}</math> | <center><math>\text{Sc}=\frac{{\tilde{\nu }}}{D}</math> | ||

(9.23)</center> | (9.23)</center> | ||

- | The evaporation problem is now described by a set of ordinary differential equations (9.14) – (9.16) subjected to boundary conditions specified by eqs. (9.18) – (9.21). Since eq. (9.14) is a third-order ordinary differential equation, it requires three boundary conditions: in eq. (9.18), | + | |

- | The boundary value problem can be solved using the Runge-Kutta method in conjunction with a shooting method. The solution procedure begins with an assumed | + | The evaporation problem is now described by a set of ordinary differential equations (9.14) – (9.16) subjected to boundary conditions specified by eqs. (9.18) – (9.21). Since eq. (9.14) is a third-order ordinary differential equation, it requires three boundary conditions: in eq. (9.18), in eq. (9.19), and eq. (9.20). Equations (9.15) and (9.16) are both second-order ordinary differential equations, each of them requiring two boundary conditions, which are specified in eqs. (9.18) and (9.19). Therefore, the problem is mathematically defined with boundary conditions stated by eqs. (9.18) – (9.20), which makes the boundary condition specified in eq. (9.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. (9.21) be satisfied. |

+ | |||

+ | The boundary value problem can be solved using the Runge-Kutta method in conjunction with a shooting method. The solution procedure begins with an assumed ω<sub>I</sub> and uses boundary conditions specified by eqs. (9.18) – (9.20) to solve eqs. (9.14) – (9.16). Once a solution is obtained, eq. (9.21) is employed to find <math>{{\omega }_{\delta }}</math>. If <math>{{\omega }_{\delta }}</math> obtained from eq. (9.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> | <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> | ||

(9.24)</center> | (9.24)</center> | ||

+ | |||

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> | <center><math>N{{u}_{x}}=\frac{{{h}_{x}}x}{k}={\theta }'(0)\operatorname{Re}_{x}^{\frac{1}{2}}</math> | ||

(9.25)</center> | (9.25)</center> | ||

+ | |||

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> | <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> | ||

(9.26)</center> | (9.26)</center> | ||

+ | |||

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> | <center><math>S{{h}_{x}}=\frac{{{h}_{m}}x}{D}={\varphi }'(0)\operatorname{Re}_{x}^{\frac{1}{2}}</math> | ||

(9.27)</center> | (9.27)</center> | ||

+ | |||

where | where | ||

<center><math>{{\operatorname{Re}}_{x}}=\frac{{{u}_{\infty }}x}{\nu }</math> | <center><math>{{\operatorname{Re}}_{x}}=\frac{{{u}_{\infty }}x}{\nu }</math> | ||

(9.28)</center> | (9.28)</center> | ||

- | Although eqs. (9.25) and (9.27) have the same form as the case of single-phase 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 | + | |

+ | Although eqs. (9.25) and (9.27) have the same form as the case of single-phase 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 single-phase 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. (9.20) – (9.21). Since the blowing velocity at the surface of the liquid film, determined by eq. (9.8), is usually much lower than the incoming vapor-gas mixture velocity, <math>{{u}_{\infty }}</math>, one can assume that the blowing velocity at the liquid surface is negligible. Then eq. (9.20) can be replaced with | The evaporation problem presented here is a coupled problem because the solutions of the momentum, energy, and species equations are coupled with eqs. (9.20) – (9.21). Since the blowing velocity at the surface of the liquid film, determined by eq. (9.8), is usually much lower than the incoming vapor-gas mixture velocity, <math>{{u}_{\infty }}</math>, one can assume that the blowing velocity at the liquid surface is negligible. Then eq. (9.20) can be replaced with | ||

<center><math>f(0)=0</math> | <center><math>f(0)=0</math> | ||

(9.29)</center> | (9.29)</center> | ||

+ | |||

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

Line 143: | Line 166: | ||

<center><math>{\varphi }'(0)=0.332S{{c}^{\frac{1}{3}}}</math> | <center><math>{\varphi }'(0)=0.332S{{c}^{\frac{1}{3}}}</math> | ||

(9.31)</center> | (9.31)</center> | ||

+ | |||

Substituting eqs. (9.30) and (9.31) into eqs. (9.25) and (9.27) yields | Substituting eqs. (9.30) and (9.31) into eqs. (9.25) and (9.27) yields | ||

Line 150: | Line 174: | ||

<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>S{{h}_{x}}=\frac{{{h}_{m}}x}{k}=0.332\operatorname{Re}_{x}^{\frac{1}{2}}S{{c}^{\frac{1}{3}}}</math> | ||

(9.33)</center> | (9.33)</center> | ||

- | 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 | ||

<center><math>\overline{Nu}=\frac{\bar{h}L}{k}=0.664\operatorname{Re}_{L}^{\frac{1}{2}}{{\Pr }^{\frac{1}{3}}}</math> | <center><math>\overline{Nu}=\frac{\bar{h}L}{k}=0.664\operatorname{Re}_{L}^{\frac{1}{2}}{{\Pr }^{\frac{1}{3}}}</math> | ||

Line 157: | Line 182: | ||

<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> | ||

(9.35)</center> | (9.35)</center> | ||

+ | |||

where | where | ||

<center><math>{{\operatorname{Re}}_{L}}=\frac{{{u}_{\infty }}L}{\nu }</math> | <center><math>{{\operatorname{Re}}_{L}}=\frac{{{u}_{\infty }}L}{\nu }</math> | ||

(9.36)</center> | (9.36)</center> | ||

+ | |||

The mass fraction of water at the liquid surface can be obtained by substituting eqs. (9.30) and (9.31) into eq. (9.21), i.e., | The mass fraction of water at the liquid surface can be obtained by substituting eqs. (9.30) and (9.31) into eq. (9.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> | <center><math>{{\omega }_{\delta }}=\frac{{{\omega }_{\infty }}+{{(\text{Sc}/\Pr )}^{2/3}}\text{Ja}}{1+{{(\text{Sc}/\Pr )}^{2/3}}\text{Ja}}</math> | ||

(9.37)</center> | (9.37)</center> | ||

+ | |||

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. | ||

+ | |||

+ | ==References== |

## Revision as of 13:33, 1 June 2010

Evaporation from thin liquid films occurs in many industrial and natural processes, including drying, evaporative cooling, and sweating. Figure 9.6 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 vapor-gas mixture flowing above the liquid film. The heat source (vapor-gas 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 mass-averaged 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 x-direction has been neglected because heat transfer occurred mainly in the y-direction. The mass fraction of water is also accounted for by assuming a constant mass diffusivity:

As the boundary conditions can be taken directly from Fig. 9.6:

, ,

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.

*u* = 0,
,
ω = ω_{δ},

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:

</center>

(9.7)Substituting eq. (9.7) into eq. (9.6), the velocity component in the y-direction becomes

<center>

(9.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

Substituting eq. (9.7) into eq. (9.9), conservation of mass at the interface becomes

Equations (9.1) – (9.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 (9.1) is automatically satisfied. The momentum equation in terms of the stream function is then

Introducing the following similarity variable:

eq. (9.12) and eqs. (9.3) – (9.4) can be reduced to a set of ordinary differential equations:

The velocity components in the x- and y-directions in terms of the dimensionless stream function f can be obtained by eqs. (9.11) and (9.13), i.e.,

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

The boundary conditions represented by eq. (9.6), except the velocity component in the y-direction, can be rewritten to

Conservation of mass and energy at the interface, eqs. (9.8) and (9.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 (9.14) – (9.16) subjected to boundary conditions specified by eqs. (9.18) – (9.21). Since eq. (9.14) is a third-order ordinary differential equation, it requires three boundary conditions: in eq. (9.18), in eq. (9.19), and eq. (9.20). Equations (9.15) and (9.16) are both second-order ordinary differential equations, each of them requiring two boundary conditions, which are specified in eqs. (9.18) and (9.19). Therefore, the problem is mathematically defined with boundary conditions stated by eqs. (9.18) – (9.20), which makes the boundary condition specified in eq. (9.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. (9.21) be satisfied.

The boundary value problem can be solved using the Runge-Kutta method in conjunction with a shooting method. The solution procedure begins with an assumed ω_{I} and uses boundary conditions specified by eqs. (9.18) – (9.20) to solve eqs. (9.14) – (9.16). Once a solution is obtained, eq. (9.21) is employed to find ω_{δ}. If ω_{δ} obtained from eq. (9.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. (9.25) and (9.27) have the same form as the case of single-phase 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 single-phase 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. (9.20) – (9.21). Since the blowing velocity at the surface of the liquid film, determined by eq. (9.8), is usually much lower than the incoming vapor-gas mixture velocity, , one can assume that the blowing velocity at the liquid surface is negligible. Then eq. (9.20) can be replaced with

*f*(0) = 0 (9.29)

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

Substituting eqs. (9.30) and (9.31) into eqs. (9.25) and (9.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. (9.30) and (9.31) into eq. (9.21), i.e.,

The applications of the above empirical correlations can be demonstrated by the following example.