Evaporation from Cylindrical Pore under High Heat Flux
From ThermalFluidsPedia
(Created page with 'At high heat flux, high velocity vapor flow may create an extended thick liquid film attached to the evaporating meniscus in spite of the capillary pressure drop between the hemi…') 

(9 intermediate revisions not shown)  
Line 1:  Line 1:  
+  [[Image:NewChapter9 (5).jpgthumb200 pxalt= Schematic of extended thick liquid film in a circular microchannel.  Figure 9.25 Schematic of extended thick liquid film in a circular microchannel. ]]  
+  
At high heat flux, high velocity vapor flow may create an extended thick liquid film attached to the evaporating meniscus in spite of the capillary pressure drop between the hemispherical meniscus and the nearly flat thick film. [[#ReferencesKhrustalev and Faghri (1997)]] analyzed the evaporation of a pure liquid in the vicinity of a hemispherical liquidvapor meniscus formed within a circular micropore, as shown in Fig. 9.18(b). Since the radius of the pore and the liquid film thickness are of the same order of magnitude, it is necessary to describe the problem in a cylindrical coordinate system, as shown in Fig. 9.25. Both the vapor and the liquid flow along the zdirection. The origin of the coordinate system is fixed with respect to the minimum film thickness. It is assumed that the vapor is at saturation condition. Heat transfer in the liquid film is in the radial direction by conduction only, and convective effects are neglected. The velocity profile in the fluid region is assumed to be fully developed and to behave as a modified laminartype flow. This accounts for surface tension and mass transfer effects through the boundary conditions, but otherwise a single, axial direction momentum equation completely describes the flow field in the fluid.  At high heat flux, high velocity vapor flow may create an extended thick liquid film attached to the evaporating meniscus in spite of the capillary pressure drop between the hemispherical meniscus and the nearly flat thick film. [[#ReferencesKhrustalev and Faghri (1997)]] analyzed the evaporation of a pure liquid in the vicinity of a hemispherical liquidvapor meniscus formed within a circular micropore, as shown in Fig. 9.18(b). Since the radius of the pore and the liquid film thickness are of the same order of magnitude, it is necessary to describe the problem in a cylindrical coordinate system, as shown in Fig. 9.25. Both the vapor and the liquid flow along the zdirection. The origin of the coordinate system is fixed with respect to the minimum film thickness. It is assumed that the vapor is at saturation condition. Heat transfer in the liquid film is in the radial direction by conduction only, and convective effects are neglected. The velocity profile in the fluid region is assumed to be fully developed and to behave as a modified laminartype flow. This accounts for surface tension and mass transfer effects through the boundary conditions, but otherwise a single, axial direction momentum equation completely describes the flow field in the fluid.  
+  
Combining expressions for the conservation of mass and the conservation of energy in the liquid film, one obtains  Combining expressions for the conservation of mass and the conservation of energy in the liquid film, one obtains  
  +  
  <center><math>\int_{R\delta }^{R}{r{{w}_{\ell }}\left( r \right)dr}=\frac{1}{2\pi {{\rho }_{\ell }}}\left( {{{\dot{m}}}_{\ell ,in}}\frac{q}{{{h}_{\ell v}}} \right)</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\int_{R\delta }^{R}{r{{w}_{\ell }}\left( r \right)dr}=\frac{1}{2\pi {{\rho }_{\ell }}}\left( {{{\dot{m}}}_{\ell ,in}}\frac{q}{{{h}_{\ell v}}} \right)</math></center>  
+  {{EquationRef(1)}}  
+  }  
where <math>{{\dot{m}}_{\ell ,in}}</math> is the incoming liquid mass flow rate at z = 0 in the liquid film layer and q(z) is the heat flow rate through a given crosssection due to evaporation in the film from z = 0 to z , i.e.,  where <math>{{\dot{m}}_{\ell ,in}}</math> is the incoming liquid mass flow rate at z = 0 in the liquid film layer and q(z) is the heat flow rate through a given crosssection due to evaporation in the film from z = 0 to z , i.e.,  
  +  
  <center><math>q=2\pi R\int_{0}^{z}{{{{{q}''}}_{R}}dz}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>q=2\pi R\int_{0}^{z}{{{{{q}''}}_{R}}dz}</math></center>  
+  {{EquationRef(2)}}  
+  }  
For conduction heat transfer through the liquid film,  For conduction heat transfer through the liquid film,  
  +  
  <center><math>{{{q}''}_{R}}={{k}_{\ell }}\frac{{{T}_{w}}{{T}_{\delta }}}{R\ln \left[ R/(R\delta ) \right]}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{{q}''}_{R}}={{k}_{\ell }}\frac{{{T}_{w}}{{T}_{\delta }}}{R\ln \left[ R/(R\delta ) \right]}</math></center>  
+  {{EquationRef(3)}}  
+  }  
Drawn from the definition of total heat flow rate, an expression for the gradient of heat flow rate in the zdirection is  Drawn from the definition of total heat flow rate, an expression for the gradient of heat flow rate in the zdirection is  
  +  
  <center><math>\frac{dq}{dz}=2\pi {{k}_{\ell }}\frac{{{T}_{w}}{{T}_{\delta }}}{\ln \left[ R/(R\delta ) \right]}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{dq}{dz}=2\pi {{k}_{\ell }}\frac{{{T}_{w}}{{T}_{\delta }}}{\ln \left[ R/(R\delta ) \right]}</math></center>  
+  {{EquationRef(4)}}  
+  }  
The differential equation of the velocity profile from the conservation of momentum in the film is  The differential equation of the velocity profile from the conservation of momentum in the film is  
  +  
  <center><math>\frac{1}{r}\frac{\partial }{\partial r}\left( r\frac{\partial {{w}_{\ell }}}{\partial r} \right)=\frac{1}{{{\mu }_{\ell }}}\frac{d{{p}_{\ell }}}{dz}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{1}{r}\frac{\partial }{\partial r}\left( r\frac{\partial {{w}_{\ell }}}{\partial r} \right)=\frac{1}{{{\mu }_{\ell }}}\frac{d{{p}_{\ell }}}{dz}</math></center>  
+  {{EquationRef(5)}}  
+  }  
The boundary conditions for this equation are the nonslip condition at the wall and equality of shear stress at the liquidvapor interface:  The boundary conditions for this equation are the nonslip condition at the wall and equality of shear stress at the liquidvapor interface:  
  +  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>r=R\begin{matrix}  <center><math>r=R\begin{matrix}  
, & {{w}_{\ell }}=0 \\  , & {{w}_{\ell }}=0 \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
  +  {{EquationRef(6)}}  
  +  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>r=R\delta \begin{matrix}  <center><math>r=R\delta \begin{matrix}  
, & \frac{\partial {{w}_{\ell }}}{\partial r}=\frac{1}{{{\mu }_{\ell }}}\left( \frac{{{f}_{v}}}{2}{{\rho }_{v}}\bar{w}_{v}^{2}\frac{d\sigma }{dT}\frac{d{{T}_{\delta }}}{dz} \right)=E \\  , & \frac{\partial {{w}_{\ell }}}{\partial r}=\frac{1}{{{\mu }_{\ell }}}\left( \frac{{{f}_{v}}}{2}{{\rho }_{v}}\bar{w}_{v}^{2}\frac{d\sigma }{dT}\frac{d{{T}_{\delta }}}{dz} \right)=E \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
  +  {{EquationRef(7)}}  
+  }  
where <math>{{f}_{v}}</math> is the vapor fraction factor.  where <math>{{f}_{v}}</math> is the vapor fraction factor.  
  The solution to eq. (  +  The solution to eq. (5) with eqs. (6) and (7) as boundary conditions, gives 
  <center><math>{{w}_{\ell }}=\frac{1}{{{\mu }_{\ell }}}\frac{d{{p}_{\ell }}}{dz}\left[ \frac{1}{4}\left( {{R}^{2}}{{r}^{2}} \right)+\frac{{{\left( R\delta \right)}^{2}}}{2}\ln \frac{r}{R} \right]+E\left( R\delta \right)\ln \frac{r}{R}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{w}_{\ell }}=\frac{1}{{{\mu }_{\ell }}}\frac{d{{p}_{\ell }}}{dz}\left[ \frac{1}{4}\left( {{R}^{2}}{{r}^{2}} \right)+\frac{{{\left( R\delta \right)}^{2}}}{2}\ln \frac{r}{R} \right]+E\left( R\delta \right)\ln \frac{r}{R}</math></center>  
+  {{EquationRef(8)}}  
+  }  
  Substituting eq. (  +  Substituting eq. (8) into eq. (1), the following equation for liquid pressure in terms of the mass flow rate and heat flow rate is obtained: 
  +  
  <center><math>\frac{d{{p}_{\ell }}}{dz}=\frac{{{\mu }_{\ell }}\left[ \frac{1}{2\pi {{\rho }_{\ell }}}\left( \frac{q}{{{h}_{\ell v}}}{{{\dot{m}}}_{\ell ,in}} \right)+E\left( R\delta \right)F \right]}{\left[ \frac{{{R}^{4}}}{16}+\frac{{{\left( R\delta \right)}^{2}}}{2}\left( F+\frac{{{\left( R\delta \right)}^{2}}}{8}\frac{{{R}^{2}}}{4} \right) \right]}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{d{{p}_{\ell }}}{dz}=\frac{{{\mu }_{\ell }}\left[ \frac{1}{2\pi {{\rho }_{\ell }}}\left( \frac{q}{{{h}_{\ell v}}}{{{\dot{m}}}_{\ell ,in}} \right)+E\left( R\delta \right)F \right]}{\left[ \frac{{{R}^{4}}}{16}+\frac{{{\left( R\delta \right)}^{2}}}{2}\left( F+\frac{{{\left( R\delta \right)}^{2}}}{8}\frac{{{R}^{2}}}{4} \right) \right]}</math></center>  
+  {{EquationRef(9)}}  
+  }  
where  where  
  +  
  <center><math>F=\frac{{{\left( R\delta \right)}^{2}}}{2}\left( \ln \frac{R}{R\delta }+\frac{1}{2} \right)\frac{{{R}^{2}}}{4}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>F=\frac{{{\left( R\delta \right)}^{2}}}{2}\left( \ln \frac{R}{R\delta }+\frac{1}{2} \right)\frac{{{R}^{2}}}{4}</math></center>  
+  {{EquationRef(10)}}  
+  }  
The pressure difference in the vapor and liquid phases due to capillary and disjoining pressure is  The pressure difference in the vapor and liquid phases due to capillary and disjoining pressure is  
  
  
  
  where the two terms in the braces represent the two radii of curvature in the microchannel. Equation (  +  { class="wikitable" border="0" 
  +    
  <center><math>\frac{d\delta }{dz}={\delta }'</math>  +   width="100%"  
  +  <center><math>{{p}_{v}}{{p}_{\ell }}=\sigma \left\{ \frac{{{d}^{2}}\delta }{d{{z}^{2}}}{{\left[ 1+{{\left( \frac{d\delta }{dz} \right)}^{2}} \right]}^{\frac{3}{2}}}+\frac{1}{R\delta }\cos \left( \arctan \frac{d\delta }{dz} \right) \right\}{{p}_{d}}</math></center>  
  +  {{EquationRef(11)}}  
  <center><math>\frac{d{\delta }'}{dz}={{\left( 1+{{{{\delta }'}}^{2}} \right)}^{3/2}}\left[ \frac{{{p}_{v}}{{p}_{\ell }}+{{p}_{d}}}{\sigma }\frac{\cos (\arctan {\delta }')}{R\delta } \right]</math>  +  } 
  +  
+  where the two terms in the braces represent the two radii of curvature in the microchannel. Equation (11) can be rewritten as the following two ordinary differential equations:  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\frac{d\delta }{dz}={\delta }'</math></center>  
+  {{EquationRef(12)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\frac{d{\delta }'}{dz}={{\left( 1+{{{{\delta }'}}^{2}} \right)}^{3/2}}\left[ \frac{{{p}_{v}}{{p}_{\ell }}+{{p}_{d}}}{\sigma }\frac{\cos (\arctan {\delta }')}{R\delta } \right]</math></center>  
+  {{EquationRef(13)}}  
+  }  
The vapor flow is assumed to be compressible and quasionedimensional [[#References(Faghri, 1995)]]:  The vapor flow is assumed to be compressible and quasionedimensional [[#References(Faghri, 1995)]]:  
  +  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>\begin{align}  <center><math>\begin{align}  
& \frac{d{{p}_{v}}}{dz}=\frac{1}{{{A}_{v}}}\left[ \frac{d}{dz}\left( {{\beta }_{v}}{{\rho }_{v}}\bar{w}_{v}^{2}{{A}_{v}} \right){{f}_{v}}\rho \bar{w}_{v}^{2}\left( R\delta \right) \right. \\  & \frac{d{{p}_{v}}}{dz}=\frac{1}{{{A}_{v}}}\left[ \frac{d}{dz}\left( {{\beta }_{v}}{{\rho }_{v}}\bar{w}_{v}^{2}{{A}_{v}} \right){{f}_{v}}\rho \bar{w}_{v}^{2}\left( R\delta \right) \right. \\  
Line 73:  Line 131:  
{} & {} & \left. 2\pi \left( R\delta \right){{\rho }_{v}}v_{v,\delta }^{2}\sin \left( \arctan \frac{d\delta }{dz} \right) \right] \\  {} & {} & \left. 2\pi \left( R\delta \right){{\rho }_{v}}v_{v,\delta }^{2}\sin \left( \arctan \frac{d\delta }{dz} \right) \right] \\  
\end{matrix} \\  \end{matrix} \\  
  \end{align}</math>  +  \end{align}</math></center> 
  +  {{EquationRef(14)}}  
+  }  
where <math>{{\beta }_{v}}</math> = 1.33 for small Reynolds numbers. The blowing velocity due to evaporation is  where <math>{{\beta }_{v}}</math> = 1.33 for small Reynolds numbers. The blowing velocity due to evaporation is  
  +  
  <center><math>{{v}_{v,\delta }}=\frac{dq}{dz}\frac{1}{2\pi \left( R\delta \right){{\rho }_{v}}{{h}_{\ell v}}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{v}_{v,\delta }}=\frac{dq}{dz}\frac{1}{2\pi \left( R\delta \right){{\rho }_{v}}{{h}_{\ell v}}}</math></center>  
+  {{EquationRef(15)}}  
+  }  
The vapor is assumed to behave as an ideal gas.  The vapor is assumed to behave as an ideal gas.  
  +  
  <center><math>{{\rho }_{v}}=\frac{{{p}_{v}}}{{{R}_{g}}{{T}_{v}}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{\rho }_{v}}=\frac{{{p}_{v}}}{{{R}_{g}}{{T}_{v}}}</math></center>  
+  {{EquationRef(16)}}  
+  }  
The density gradient in the zdirection is then  The density gradient in the zdirection is then  
  +  
  <center><math>\frac{d{{\rho }_{v}}}{dz}=\frac{1}{{{R}_{g}}}\left( \frac{d{{p}_{v}}}{dz}\frac{1}{{{T}_{v}}}\frac{{{p}_{v}}}{T_{v}^{2}}\frac{d{{T}_{v}}}{dz} \right)</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{d{{\rho }_{v}}}{dz}=\frac{1}{{{R}_{g}}}\left( \frac{d{{p}_{v}}}{dz}\frac{1}{{{T}_{v}}}\frac{{{p}_{v}}}{T_{v}^{2}}\frac{d{{T}_{v}}}{dz} \right)</math></center>  
+  {{EquationRef(17)}}  
+  }  
The temperature and pressure at saturation are related by the ClausiusClapeyron equation –  The temperature and pressure at saturation are related by the ClausiusClapeyron equation –  
  
  
  
  At this point, seven expressions – eqs. (  +  { class="wikitable" border="0" 
  +    
  +   width="100%"   
  +  <center><math>\frac{d{{T}_{v}}}{dz}=\frac{d{{p}_{v}}}{dz}\frac{{{R}_{g}}T_{v}^{2}}{{{p}_{v}}{{h}_{\ell v}}}</math></center>  
  +  {{EquationRef(18)}}  
  +  }  
  +  
  +  At this point, seven expressions – eqs. (4), (9), (13), (14), (15), (17), and (18) – involving seven unknown variables – <math>q</math>, <math>\delta </math>, <math>{\delta }'</math>, <math>{{p}_{\ell }}</math>, <math>{{p}_{v}}</math>, <math>{{\rho }_{v}}</math>, and Tv – completely describe the stated assumptions. Boundary conditions are required for each of the seven firstorder differential equations at z = 0:  
  +  
  +  
  +  
  +  
  +  
  <center><math>{{p}_{v}}={{p}_{v,in}}={{p}_{sat}}({{T}_{v,in}})</math>  +  { class="wikitable" border="0" 
  +    
  +   width="100%"   
  <center><math>{{\rho }_{v,in}}=\frac{{{p}_{v,in}}}{{{R}_{g}}{{T}_{v,in}}}</math>  +  <center><math>q=0</math></center> 
  +  {{EquationRef(19)}}  
  +  }  
  <center><math>{{T}_{v}}={{T}_{v,in}}</math>  +  
  +  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{{p}_{\ell }}={{p}_{v,in}}\frac{2\sigma }{R{{\delta }_{in}}}+{{p}_{d}}</math></center>  
+  {{EquationRef(20)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\delta ={{\delta }_{in}}</math></center>  
+  {{EquationRef(21)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{\delta }'=0</math></center>  
+  {{EquationRef(22)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{{p}_{v}}={{p}_{v,in}}={{p}_{sat}}({{T}_{v,in}})</math></center>  
+  {{EquationRef(23)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{{\rho }_{v,in}}=\frac{{{p}_{v,in}}}{{{R}_{g}}{{T}_{v,in}}}</math></center>  
+  {{EquationRef(24)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{{T}_{v}}={{T}_{v,in}}</math></center>  
+  {{EquationRef(25)}}  
+  }  
As can be seen from Fig. 9.25, the liquid film ends with a microfilm; therefore, its thickness at <math>z={{L}_{\delta }}</math> equals the nonevaporating film thickness δ0, determined by eq. (9.244). The interfacial temperature <math>{{T}_{\delta }}</math> is required to obtain the gradient of the heat transfer rate along the zdirection, as indicated by eq. (9.266). It can be determined by solving eqs. (9.241) and (9.242) simultaneously. The mass flow rate of the liquid at z = 0, <math>{{\dot{m}}_{\ell ,in}}</math>, can be determined by the following boundary condition at the end of the microfilm:  As can be seen from Fig. 9.25, the liquid film ends with a microfilm; therefore, its thickness at <math>z={{L}_{\delta }}</math> equals the nonevaporating film thickness δ0, determined by eq. (9.244). The interfacial temperature <math>{{T}_{\delta }}</math> is required to obtain the gradient of the heat transfer rate along the zdirection, as indicated by eq. (9.266). It can be determined by solving eqs. (9.241) and (9.242) simultaneously. The mass flow rate of the liquid at z = 0, <math>{{\dot{m}}_{\ell ,in}}</math>, can be determined by the following boundary condition at the end of the microfilm:  
  
  
  
  The problem described by seven ordinary differential equations – eqs. (  +  { class="wikitable" border="0" 
+    
+   width="100%"   
+  <center><math>{{\dot{m}}_{\ell ,in}}=\frac{q({{L}_{\delta }})}{{{h}_{\ell v}}}</math></center>  
+  {{EquationRef(26)}}  
+  }  
+  
+  [[Image:NewChapter9 (6).jpgthumb200 pxalt= Variation of characteristics of the liquidvapor interface along the cylindrical channel Figure 9.26 Variation of characteristics of the liquidvapor interface along the cylindrical channel: (a) thickness of the liquid film, and (b) pressure in the vapor and liquid phases. ]]  
+  [[Image:NewChapter9 (15).jpgthumb200 pxalt= Figure 9.27 Variation of characteristics of the liquidvapor interface along the cylindrical channel  Figure 9.27 Variation of characteristics of the liquidvapor interface along the cylindrical channel (a) heat flux at the wall, and (b) temperature of the vapor and liquidvapor interface. ]]  
+  
+  The problem described by seven ordinary differential equations – eqs. (4), (9), (13), (14), (15), (17) and (18) – with corresponding boundary conditions – eqs. (19) – (25) – can be solved using a RungeKutta method in conjunction with the shooting method to satisfy the constitutive equation, eq. (23). More detailed information about the numerical procedure can be found in [[#ReferencesKhrustalev and Faghri (1997)]].  
Numerical solution is performed for water evaporation from a cylindrical pore with an inner diameter of 20 <math>\text{ }\!\!\mu\!\!\text{ m}</math> and wall temperature of 388 K. The accommodation coefficient α is set equal to unity. Fig. 9.26(a) shows the variation in liquid film thickness along the zdirection. The liquid film thickness initially increases slowly along the zcoordinate, and then rapidly decreases approaching the microfilm region at the end of the film. The film is concave at z=0, where it attaches to the hemispherical meniscus; however, it is convex over most of its length. This is due to the pressure gradient in the liquid and vapor along the film. Figure 9.26(b) demonstrates that the vapor pressure is decreasing and the liquid pressure is increasing along the zcoordinate. The liquid flows in the film due mainly to frictional vaporliquid interaction at the interface.  Numerical solution is performed for water evaporation from a cylindrical pore with an inner diameter of 20 <math>\text{ }\!\!\mu\!\!\text{ m}</math> and wall temperature of 388 K. The accommodation coefficient α is set equal to unity. Fig. 9.26(a) shows the variation in liquid film thickness along the zdirection. The liquid film thickness initially increases slowly along the zcoordinate, and then rapidly decreases approaching the microfilm region at the end of the film. The film is concave at z=0, where it attaches to the hemispherical meniscus; however, it is convex over most of its length. This is due to the pressure gradient in the liquid and vapor along the film. Figure 9.26(b) demonstrates that the vapor pressure is decreasing and the liquid pressure is increasing along the zcoordinate. The liquid flows in the film due mainly to frictional vaporliquid interaction at the interface.  
  +  
+  Figure 9.27(a) demonstrates that the rate of evaporation decreases towards the thickest part of the liquid film and reaches its maximum in the microfilm region. However, the microfilm region at the end of the liquid film does not make any significant contribution to the total mass flow rate at the pore outlet, because (a) the length of the microfilm region is much smaller than in the case where a hemispherical meniscus ends directly with a microfilm; and (b) disjoining pressure causes the liquidvapor interfacial temperature <math>{{T}_{\delta }}</math> to sharply increase at the end of the film, as shown in Fig. 9.27 (b).  
+  
+  [[#ReferencesKhrustalev and Faghri’s (1997)]] work demonstrated, for the first time, that highvapor velocities during the evaporation of pure liquids in micropores could allow a thick liquid film to attach to a hemispherical meniscus. In addition to the evaporation of pure liquid in a circular channel discussed in this section, [[#ReferencesCoquard et al. (2005)]] investigated evaporation in a capillary tube with a square crosssection. The evaporation rate is much higher than in a circular tube because of the liquid flow along the corner induced by capillary force. They focused on slow evaporation controlled by mass transfer and identified three regimes: the capillary regime, capillaryviscous regime, and capillarygravity regime.  
+  
+  ==References==  
+  
+  Coquard, T., Camassel, B., and Prat, M., 2005, “Evaporation in Capillary Tubers of Squared Cross Section,” Proceedings of 2005 ASME Summer Heat Transfer Conference, San Francisco, CA.  
+  
+  Faghri, A., 1995, Heat Pipe Science and Technology, Taylor and Francis, Washington, D.C.  
+  
+  Khrustalev, D., and Faghri, A., 1997, “Thick Film Phenomenon in High Heat Flux Evaporation from Cylindrical Pores,” ASME Journal of Heat Transfer, Vol. 119, No. 2, pp. 272278. 
Current revision as of 20:07, 3 June 2010
At high heat flux, high velocity vapor flow may create an extended thick liquid film attached to the evaporating meniscus in spite of the capillary pressure drop between the hemispherical meniscus and the nearly flat thick film. Khrustalev and Faghri (1997) analyzed the evaporation of a pure liquid in the vicinity of a hemispherical liquidvapor meniscus formed within a circular micropore, as shown in Fig. 9.18(b). Since the radius of the pore and the liquid film thickness are of the same order of magnitude, it is necessary to describe the problem in a cylindrical coordinate system, as shown in Fig. 9.25. Both the vapor and the liquid flow along the zdirection. The origin of the coordinate system is fixed with respect to the minimum film thickness. It is assumed that the vapor is at saturation condition. Heat transfer in the liquid film is in the radial direction by conduction only, and convective effects are neglected. The velocity profile in the fluid region is assumed to be fully developed and to behave as a modified laminartype flow. This accounts for surface tension and mass transfer effects through the boundary conditions, but otherwise a single, axial direction momentum equation completely describes the flow field in the fluid.
Combining expressions for the conservation of mass and the conservation of energy in the liquid film, one obtains

where is the incoming liquid mass flow rate at z = 0 in the liquid film layer and q(z) is the heat flow rate through a given crosssection due to evaporation in the film from z = 0 to z , i.e.,

For conduction heat transfer through the liquid film,

Drawn from the definition of total heat flow rate, an expression for the gradient of heat flow rate in the zdirection is

The differential equation of the velocity profile from the conservation of momentum in the film is

The boundary conditions for this equation are the nonslip condition at the wall and equality of shear stress at the liquidvapor interface:


where f_{v} is the vapor fraction factor. The solution to eq. (5) with eqs. (6) and (7) as boundary conditions, gives

Substituting eq. (8) into eq. (1), the following equation for liquid pressure in terms of the mass flow rate and heat flow rate is obtained:

where

The pressure difference in the vapor and liquid phases due to capillary and disjoining pressure is

where the two terms in the braces represent the two radii of curvature in the microchannel. Equation (11) can be rewritten as the following two ordinary differential equations:


The vapor flow is assumed to be compressible and quasionedimensional (Faghri, 1995):

where β_{v} = 1.33 for small Reynolds numbers. The blowing velocity due to evaporation is

The vapor is assumed to behave as an ideal gas.

The density gradient in the zdirection is then

The temperature and pressure at saturation are related by the ClausiusClapeyron equation –

At this point, seven expressions – eqs. (4), (9), (13), (14), (15), (17), and (18) – involving seven unknown variables – q, δ, δ', , p_{v}, ρ_{v}, and Tv – completely describe the stated assumptions. Boundary conditions are required for each of the seven firstorder differential equations at z = 0:







As can be seen from Fig. 9.25, the liquid film ends with a microfilm; therefore, its thickness at z = L_{δ} equals the nonevaporating film thickness δ0, determined by eq. (9.244). The interfacial temperature T_{δ} is required to obtain the gradient of the heat transfer rate along the zdirection, as indicated by eq. (9.266). It can be determined by solving eqs. (9.241) and (9.242) simultaneously. The mass flow rate of the liquid at z = 0, , can be determined by the following boundary condition at the end of the microfilm:

The problem described by seven ordinary differential equations – eqs. (4), (9), (13), (14), (15), (17) and (18) – with corresponding boundary conditions – eqs. (19) – (25) – can be solved using a RungeKutta method in conjunction with the shooting method to satisfy the constitutive equation, eq. (23). More detailed information about the numerical procedure can be found in Khrustalev and Faghri (1997).
Numerical solution is performed for water evaporation from a cylindrical pore with an inner diameter of 20 and wall temperature of 388 K. The accommodation coefficient α is set equal to unity. Fig. 9.26(a) shows the variation in liquid film thickness along the zdirection. The liquid film thickness initially increases slowly along the zcoordinate, and then rapidly decreases approaching the microfilm region at the end of the film. The film is concave at z=0, where it attaches to the hemispherical meniscus; however, it is convex over most of its length. This is due to the pressure gradient in the liquid and vapor along the film. Figure 9.26(b) demonstrates that the vapor pressure is decreasing and the liquid pressure is increasing along the zcoordinate. The liquid flows in the film due mainly to frictional vaporliquid interaction at the interface.
Figure 9.27(a) demonstrates that the rate of evaporation decreases towards the thickest part of the liquid film and reaches its maximum in the microfilm region. However, the microfilm region at the end of the liquid film does not make any significant contribution to the total mass flow rate at the pore outlet, because (a) the length of the microfilm region is much smaller than in the case where a hemispherical meniscus ends directly with a microfilm; and (b) disjoining pressure causes the liquidvapor interfacial temperature T_{δ} to sharply increase at the end of the film, as shown in Fig. 9.27 (b).
Khrustalev and Faghri’s (1997) work demonstrated, for the first time, that highvapor velocities during the evaporation of pure liquids in micropores could allow a thick liquid film to attach to a hemispherical meniscus. In addition to the evaporation of pure liquid in a circular channel discussed in this section, Coquard et al. (2005) investigated evaporation in a capillary tube with a square crosssection. The evaporation rate is much higher than in a circular tube because of the liquid flow along the corner induced by capillary force. They focused on slow evaporation controlled by mass transfer and identified three regimes: the capillary regime, capillaryviscous regime, and capillarygravity regime.
References
Coquard, T., Camassel, B., and Prat, M., 2005, “Evaporation in Capillary Tubers of Squared Cross Section,” Proceedings of 2005 ASME Summer Heat Transfer Conference, San Francisco, CA.
Faghri, A., 1995, Heat Pipe Science and Technology, Taylor and Francis, Washington, D.C.
Khrustalev, D., and Faghri, A., 1997, “Thick Film Phenomenon in High Heat Flux Evaporation from Cylindrical Pores,” ASME Journal of Heat Transfer, Vol. 119, No. 2, pp. 272278.