# Evaporation in an Inclined Microchannel

### From Thermal-FluidsPedia

The analysis in the preceding subsection dealt with evaporation of a liquid to pure vapor when the evaporation was driven by heat transfer. Chakraborty and Som (2005) analytically investigated heat transfer in an evaporating thin liquid film moving slowly along the walls of an inclined microchannel, and this work will be introduced in this subsection.

Figure 9.28 shows the lower half of the microchannel, in which quiescent air with a temperature and mass fraction is in the core of the channel while liquid film forms on the microchannel wall. It is assumed that the transport phenomena are two-dimensional and that the diffusion in the axial direction is negligible. The temperature distribution in the liquid film is assumed to be linear along the y-direction. The momentum equation that governs the liquid flow is

and is subject to the following boundary conditions

*u*= 0 at

*y*= 0 (9.291)

Integrating eq. (9.290) twice and considering eqs. (9.291) and (9.292) yields

The energy equation in the vapor phase is

where *V*_{s} is the Stefan velocity due to evaporation at the liquid surface. Equation (9.294) is subject to the following boundary conditions

*T*

_{v}=

*T*

_{δ}at

*y*= δ (9.295)

The solution of eq. (9.294) with boundary conditions eqs. (9.295) and (9.296) is

Similarly, the mass fraction of the vapor is

where *D*_{v} is the mass diffusivity of the vapor in the air.

The mass balance at the interface is

where is the mass flow rate per unit width. Considering the velocity profile in eq. (9.293), the mass balance becomes

Assuming the interface is impermeable to nonevaporating species, the mass balance at the interface can be written as

Substituting eq. (9.298) into eq. (9.301) yields

The energy balance at the interface is

Substituting eq. (9.297) into eq. (9.303) and considering the linear temperature profile in the liquid film, the interfacial temperature becomes

where

The pressure in the liquid film is

where is the vapor pressure, *R*(*x*) is the radius of curvature at *x*, and the last term on the right-hand side accounts for the effect of disjoining pressure.

Figure 9.29 Local Nusselt number (Chakraborty and Som, 2005; Reprinted with permission from Elsevier).
Therefore, the function *F*(*x*) in eq. (9.290) becomes

Substituting eqs. (9.302) and (9.304) into eq. (9.300) and considering the variation of , a fourth order ordinary differential equation of δ can be obtained. The boundary conditions of the ordinary differential equation include δ = δ_{0} at *x* = 0 and . The liquid film thickness can be obtained by a fourth order Runge-Kutta method.
The local Nusselt number is then obtained by

Figure 9.29 shows the variation of the local Nusselt number along the axial direction for different δ_{0} / *h*. It can be seen that the Nusselt number increases with increasing x and eventually reaches a very high value. This represents the location where dryout occurs and the liquid film becomes very small. The Nusselt number increases with decreasing initial film thickness δ_{0} / *h*, but its effect becomes insignificant for δ_{0} / *h* < 0.01.