Evaporation from wavy laminar falling film
From ThermalFluidsPedia
Line 1:  Line 1:  
Wavy flows of thin liquid films have higher heat transfer coefficients than smooth thin films. This effect is due to the former’s greater interfacial surface area and mixing action. [[#ReferencesFaghri and Seban (1985)]] analyzed a system in which a liquid film at an initial nondimensional temperature of zero flowed down a vertical wall and evaporated at the free liquidvapor interface. The laminar liquid film thickness varied sinusoidally and the Reynolds number ranged between 35 and 472. The energy equation for a twodimensional situation using the same vertical wall shown in Fig. 9.11 is as follows:  Wavy flows of thin liquid films have higher heat transfer coefficients than smooth thin films. This effect is due to the former’s greater interfacial surface area and mixing action. [[#ReferencesFaghri and Seban (1985)]] analyzed a system in which a liquid film at an initial nondimensional temperature of zero flowed down a vertical wall and evaporated at the free liquidvapor interface. The laminar liquid film thickness varied sinusoidally and the Reynolds number ranged between 35 and 472. The energy equation for a twodimensional situation using the same vertical wall shown in Fig. 9.11 is as follows:  
  +  
  <center><math>\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \left( \frac{{{\partial }^{2}}T}{\partial {{x}^{2}}}+\frac{{{\partial }^{2}}T}{\partial {{y}^{2}}} \right)</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \left( \frac{{{\partial }^{2}}T}{\partial {{x}^{2}}}+\frac{{{\partial }^{2}}T}{\partial {{y}^{2}}} \right)</math></center>  
+  {{EquationRef(1)}}  
+  }  
The continuity equation for this system is  The continuity equation for this system 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(2)}}  
+  }  
[[#ReferencesFaghri and Seban (1985)]] assumed a quasiparabolic velocity profile in the xdirection as follows:  [[#ReferencesFaghri and Seban (1985)]] assumed a quasiparabolic velocity profile in the xdirection as follows:  
  +  
  <center><math>u=3\bar{u}(x,t)\left( \eta \frac{{{\eta }^{2}}}{2} \right)</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>u=3\bar{u}(x,t)\left( \eta \frac{{{\eta }^{2}}}{2} \right)</math></center>  
+  {{EquationRef(3)}}  
+  }  
where the nondimensional coordinate ''η'' is defined as  where the nondimensional coordinate ''η'' is defined as  
  +  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\eta =\frac{y}{\delta }</math></center>  
+  {{EquationRef(4)}}  
+  }  
  and <math>\bar{u}</math> is the local mean velocity. Equation (  +  and <math>\bar{u}</math> is the local mean velocity. Equation (2) can be rearranged and integrated with respect to ''y'' to obtain the velocity component ''v'': 
  +  
  <center><math>v=\int\limits_{0}^{y}{\frac{\partial u}{\partial x}dy=3}\left[ \frac{\partial \bar{u}}{\partial x}\delta \left( \frac{{{\eta }^{2}}}{2}\frac{{{\eta }^{3}}}{6} \right)\bar{u}\frac{\partial \delta }{\partial x}\left( \frac{{{\eta }^{2}}}{2}\frac{{{\eta }^{3}}}{3} \right) \right]</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>v=\int\limits_{0}^{y}{\frac{\partial u}{\partial x}dy=3}\left[ \frac{\partial \bar{u}}{\partial x}\delta \left( \frac{{{\eta }^{2}}}{2}\frac{{{\eta }^{3}}}{6} \right)\bar{u}\frac{\partial \delta }{\partial x}\left( \frac{{{\eta }^{2}}}{2}\frac{{{\eta }^{3}}}{3} \right) \right]</math></center>  
+  {{EquationRef(5)}}  
+  }  
+  
+  It can be seen that eq. (3) has been substituted in order to find an expression for the <math>v</math> component of the velocity where the local mean velocity and film thickness δ are still unknown as functions of the wavy surface. To determine these two variables, the following mass balance can be written as a function of time and space:  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\frac{\partial \delta }{\partial t}=\frac{\partial }{\partial x}\int_{0}^{\delta }{udy=\frac{\partial }{\partial x}\left( \bar{u}\delta \right)}</math></center>  
+  {{EquationRef(6)}}  
+  }  
  
  
  
  
The film thickness can be written in terms of its average value,  The film thickness can be written in terms of its average value,  
<math>\overline{\delta }</math>  <math>\overline{\delta }</math>  
, and a local amplitude φ.  , and a local amplitude φ.  
  +  
  <center><math>\delta =\bar{\delta }(1+\varphi )</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\delta =\bar{\delta }(1+\varphi )</math></center>  
+  {{EquationRef(7)}}  
+  }  
Assuming that the wave formation has a periodic characteristic with a velocity ''c'', the following relations can be written for the fluctuation of local mean velocity and film thickness:  Assuming that the wave formation has a periodic characteristic with a velocity ''c'', the following relations can be written for the fluctuation of local mean velocity and film thickness:  
  
  
  
  
  
  
  +  { class="wikitable" border="0"  
  +    
  <center><math>\  +   width="100%"  
  +  <center><math>\frac{\partial \bar{u}}{\partial t}=c\frac{\partial \bar{u}}{\partial x}</math></center>  
+  {{EquationRef(8)}}  
+  }  
  +  { class="wikitable" border="0"  
  +    
  <center><math>\  +   width="100%"  
  +  <center><math>\frac{\partial \delta }{\partial t}=c\frac{\partial \delta }{\partial x}</math></center>  
+  {{EquationRef(9)}}  
+  }  
  +  Combining eqs. (6) – (9) and integrating, the following expression is obtained:  
  +  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\left( c\bar{u} \right)\left( 1+\varphi \right)=\left[ c{{{\bar{u}}}_{0}} \right]</math></center>  
+  {{EquationRef(10)}}  
+  }  
+  
+  where <math>{{\bar{u}}_{0}}</math> is the average velocity for the average film thickness <math>\bar{\delta }</math>. For very small local amplitudes, φ < 1, <math>\bar{u}</math> may be approximated by expanding eq. (10) and neglecting the thirdorder terms.  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\bar{u}={{\bar{u}}_{0}}+(c{{\bar{u}}_{0}})\varphi (c{{\bar{u}}_{0}}){{\varphi }^{2}}</math></center>  
+  {{EquationRef(11)}}  
+  }  
+  
+  Substituting eqs. (7) and (11) into eq. (5), one obtains the velocity profile in terms of the wavy surface parameters:  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>\begin{align}  <center><math>\begin{align}  
& v=3{{{\bar{u}}}_{0}}\bar{\delta }\frac{\partial \phi }{\partial x}\left\{ \left( \frac{c}{{{{\bar{u}}}_{0}}}1 \right)(12\phi )\frac{\delta }{{\bar{\delta }}}\left( \frac{{{\eta }^{2}}}{2}\frac{{{\eta }^{3}}}{6} \right) \right. \\  & v=3{{{\bar{u}}}_{0}}\bar{\delta }\frac{\partial \phi }{\partial x}\left\{ \left( \frac{c}{{{{\bar{u}}}_{0}}}1 \right)(12\phi )\frac{\delta }{{\bar{\delta }}}\left( \frac{{{\eta }^{2}}}{2}\frac{{{\eta }^{3}}}{6} \right) \right. \\  
Line 59:  Line 108:  
{} & {} \\  {} & {} \\  
\end{matrix}\left[ 1+\left( \frac{c}{{{{\bar{u}}}_{0}}}1 \right)\phi \left( \frac{c}{{{{\bar{u}}}_{0}}}1 \right){{\phi }^{2}} \right]\left( \frac{{{\eta }^{2}}}{2}\frac{{{\eta }^{3}}}{3} \right) \right\} \\  \end{matrix}\left[ 1+\left( \frac{c}{{{{\bar{u}}}_{0}}}1 \right)\phi \left( \frac{c}{{{{\bar{u}}}_{0}}}1 \right){{\phi }^{2}} \right]\left( \frac{{{\eta }^{2}}}{2}\frac{{{\eta }^{3}}}{3} \right) \right\} \\  
  \end{align}</math>  +  \end{align}</math></center> 
  +  {{EquationRef(12)}}  
+  }  
Assuming the wave is sinusoidal,  Assuming the wave is sinusoidal,  
  +  
  <center><math>\varphi =A\sin \left[ \left( \frac{2\pi }{\lambda } \right)\left( xct \right) \right]</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\varphi =A\sin \left[ \left( \frac{2\pi }{\lambda } \right)\left( xct \right) \right]</math></center>  
+  {{EquationRef(13)}}  
+  }  
and introducing the nondimensional variable ξ defined as  and introducing the nondimensional variable ξ defined as  
  +  
  <center><math>\xi =\left( \frac{2\pi }{\lambda } \right)\left( xct \right)</math>  +  { class="wikitable" border="0" 
  +    
  where ''A'' is magnitude of the wave, and λ is the wavelength, eq. (  +   width="100%"  
  +  <center><math>\xi =\left( \frac{2\pi }{\lambda } \right)\left( xct \right)</math></center>  
  <center><math>{{C}_{1}}\frac{\partial T}{\partial \xi }+{{C}_{2}}\frac{\partial T}{\partial \eta }={{C}_{3}}\frac{{{\partial }^{2}}T}{\partial {{\eta }^{2}}}+{{C}_{4}}\frac{{{\partial }^{2}}T}{\partial \eta \partial \xi }+{{C}_{5}}\frac{{{\partial }^{2}}T}{\partial {{\xi }^{2}}}</math>  +  {{EquationRef(14)}} 
  +  }  
+  
+  where ''A'' is magnitude of the wave, and λ is the wavelength, eq. (1) is transformed from the (t, x, y) to the (ξ, η) coordinate system,  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{{C}_{1}}\frac{\partial T}{\partial \xi }+{{C}_{2}}\frac{\partial T}{\partial \eta }={{C}_{3}}\frac{{{\partial }^{2}}T}{\partial {{\eta }^{2}}}+{{C}_{4}}\frac{{{\partial }^{2}}T}{\partial \eta \partial \xi }+{{C}_{5}}\frac{{{\partial }^{2}}T}{\partial {{\xi }^{2}}}</math></center>  
+  {{EquationRef(15)}}  
+  }  
where the constants are as follows:  where the constants are as follows:  
  +  
  <center><math>{{C}_{1}}=\frac{2\pi }{\lambda }(uc)</math>  +  { class="wikitable" border="0" 
  +    
  +   width="100%"   
+  <center><math>{{C}_{1}}=\frac{2\pi }{\lambda }(uc)</math></center>  
+  {{EquationRef(16)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>\begin{align}  <center><math>\begin{align}  
& {{C}_{2}}=\frac{2\pi \eta }{\lambda }\frac{{\bar{\delta }}}{\delta }(cu)A\cos \xi {{\left( \frac{2\pi }{\lambda } \right)}^{2}}\frac{{\bar{\delta }}}{\delta }\alpha \eta A\sin \xi \\  & {{C}_{2}}=\frac{2\pi \eta }{\lambda }\frac{{\bar{\delta }}}{\delta }(cu)A\cos \xi {{\left( \frac{2\pi }{\lambda } \right)}^{2}}\frac{{\bar{\delta }}}{\delta }\alpha \eta A\sin \xi \\  
Line 86:  Line 156:  
{} & {} \\  {} & {} \\  
\end{matrix}2{{\left( \frac{2\pi }{\lambda } \right)}^{2}}\alpha \left( \frac{{\bar{\delta }}}{\delta } \right)\eta {{A}^{2}}{{\cos }^{2}}\xi +\frac{v}{\delta } \\  \end{matrix}2{{\left( \frac{2\pi }{\lambda } \right)}^{2}}\alpha \left( \frac{{\bar{\delta }}}{\delta } \right)\eta {{A}^{2}}{{\cos }^{2}}\xi +\frac{v}{\delta } \\  
  \end{align}</math>  +  \end{align}</math></center> 
  +  {{EquationRef(17)}}  
  +  }  
  <center><math>{{C}_{3}}=\frac{\alpha }{{{\delta }^{2}}}+\alpha {{\left( \frac{2\pi }{\lambda } \right)}^{2}}{{\left( \frac{{\bar{\delta }}}{\delta } \right)}^{2}}{{\eta }^{2}}{{A}^{2}}{{\cos }^{2}}\xi </math>  +  
  +  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{{C}_{3}}=\frac{\alpha }{{{\delta }^{2}}}+\alpha {{\left( \frac{2\pi }{\lambda } \right)}^{2}}{{\left( \frac{{\bar{\delta }}}{\delta } \right)}^{2}}{{\eta }^{2}}{{A}^{2}}{{\cos }^{2}}\xi </math></center>  
+  {{EquationRef(18)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{{C}_{4}}=2\alpha {{\left( \frac{2\pi }{\lambda } \right)}^{2}}\eta \frac{{\bar{\delta }}}{\delta }A\cos \xi </math></center>  
+  {{EquationRef(19)}}  
+  }  
  <center><math>{{C}_{  +  { class="wikitable" border="0" 
  +    
  +   width="100%"   
  +  <center><math>{{C}_{5}}=\alpha {{\left( \frac{2\pi }{\lambda } \right)}^{2}}</math></center>  
  +  {{EquationRef(20)}}  
+  }  
  These coefficients are evaluated with <math>\bar{u}</math> from eqs. (  +  These coefficients are evaluated with <math>\bar{u}</math> from eqs. (3) and (11), and v from eq. (12). 
  In the evaporating film with a fluid with a high latent heat of vaporization, the average film thickness over a wave length, <math>\bar{\delta },</math> will not vary substantially with distance, x. Furthermore, the temporal average temperature profile will not vary with x far from where heating begins. Therefore, the boundary condition for eq. (  +  In the evaporating film with a fluid with a high latent heat of vaporization, the average film thickness over a wave length, <math>\bar{\delta },</math> will not vary substantially with distance, x. Furthermore, the temporal average temperature profile will not vary with x far from where heating begins. Therefore, the boundary condition for eq. (15) that describes evaporation on a wavy film is approximated as below. 
The equality of temperatures at corresponding points in the period gives  The equality of temperatures at corresponding points in the period gives  
  +  
  <center><math>T(0,\eta )=T(2\pi ,\eta )</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>T(0,\eta )=T(2\pi ,\eta )</math></center>  
+  {{EquationRef(21)}}  
+  }  
Boundary conditions for dimensionless temperature at the heated vertical wall and wavy surface are, respectively,  Boundary conditions for dimensionless temperature at the heated vertical wall and wavy surface are, respectively,  
  +  
  <center><math>T(\xi ,0)=1</math>  +  { class="wikitable" border="0" 
  +    
  +   width="100%"   
  <center><math>T(\xi ,1)=0</math>  +  <center><math>T(\xi ,0)=1</math></center> 
  +  {{EquationRef(22)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>T(\xi ,1)=0</math></center>  
+  {{EquationRef(23)}}  
+  }  
Therefore, the energy equation with the above boundary conditions can be solved to obtain the wavy condensate film’s temperature profile. With this profile, the heat flux at the wall can be found from  Therefore, the energy equation with the above boundary conditions can be solved to obtain the wavy condensate film’s temperature profile. With this profile, the heat flux at the wall can be found from  
  +  
  <center><math>{{{q}''}_{w}}=\frac{k}{\delta }{{\left. \frac{\partial T}{\partial \eta } \right}_{\eta =0}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{{q}''}_{w}}=\frac{k}{\delta }{{\left. \frac{\partial T}{\partial \eta } \right}_{\eta =0}}</math></center>  
+  {{EquationRef(24)}}  
+  }  
and the heat flux normal to the wavy surface (η = 1) is as follows:  and the heat flux normal to the wavy surface (η = 1) is as follows:  
  
  
  
  where <math>\mathbf{i}</math> and <math>\mathbf{j}</math> are unit vectors in the x and ydirections, respectively. The average heat flow over a period 0 < ξ < 2π is found by integrating eqs. (  +  { class="wikitable" border="0" 
  +    
  <center><math>\frac{{{{\bar{{q}''}}}_{w}}}{k}=\frac{1}{2\pi }\int_{0}^{2\pi }{\frac{1}{\delta }{{\left. \frac{\partial T}{\partial \eta } \right}_{0}}d\xi }</math>  +   width="100%"  
  +  <center><math>\frac{{{{\mathbf{{q}''}}}_{\delta }}}{k}=\left[ \frac{2\pi }{\lambda }\left( \frac{{\bar{\delta }}}{\delta }A\cos \xi \right)\mathbf{i}\frac{1}{\delta }\mathbf{j} \right]{{\left. \frac{\partial T}{\partial \eta } \right}_{\eta =1}}</math></center>  
+  {{EquationRef(25)}}  
+  }  
+  
+  where <math>\mathbf{i}</math> and <math>\mathbf{j}</math> are unit vectors in the x and ydirections, respectively. The average heat flow over a period 0 < ξ < 2π is found by integrating eqs. (24) and (25). For a sinusoidal wavy layer, such as the one considered here, this average is identical to the average obtained for the actual surface length. The average heat flux at the wall (η = 0) is approximated by  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\frac{{{{\bar{{q}''}}}_{w}}}{k}=\frac{1}{2\pi }\int_{0}^{2\pi }{\frac{1}{\delta }{{\left. \frac{\partial T}{\partial \eta } \right}_{0}}d\xi }</math></center>  
+  {{EquationRef(26)}}  
+  }  
and the average normal heat flux at the wavy surface is as follows:  and the average normal heat flux at the wavy surface is as follows:  
  +  
  <center><math>\frac{{{{\bar{{q}''}}}_{\delta }}}{k}=\frac{1}{2\pi }\int_{0}^{2\pi }{\frac{1}{\delta }{{\left. \frac{\partial T}{\partial \eta } \right}_{\eta =1}}{{\left[ 1+{{\left( \frac{\partial \delta }{\partial x} \right)}^{2}} \right]}^{1/2}}d\xi }</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{{{{\bar{{q}''}}}_{\delta }}}{k}=\frac{1}{2\pi }\int_{0}^{2\pi }{\frac{1}{\delta }{{\left. \frac{\partial T}{\partial \eta } \right}_{\eta =1}}{{\left[ 1+{{\left( \frac{\partial \delta }{\partial x} \right)}^{2}} \right]}^{1/2}}d\xi }</math></center>  
+  {{EquationRef(27)}}  
+  }  
Table 9.1 Comparison of wavy film analysis with experimental data [[#References(Faghri and Seban, 1985)]].  Table 9.1 Comparison of wavy film analysis with experimental data [[#References(Faghri and Seban, 1985)]].  
Line 238:  Line 349:  
* Conduction. Reprinted with permission from Elsevier.  * Conduction. Reprinted with permission from Elsevier.  
  The calculations were made using 40 increments in both <math>\xi </math> and <math>\eta </math> for the value of <math>2\pi /\lambda </math>, ''A'', and <math>c/{{\bar{u}}_{0}}</math>. Table 9.1 lists these experimental determinations made by the cited authors. Rows 1 through 7 give the experimental conditions and measurements. Row 8 gives the Prandtl numbers for which the calculations were made. One, on the order of 7, corresponds to the experimental conditions, and the other, 1.7, was selected as a comparison to show the effect of Prandtl number. Row 9 gives the calculated average Nusselt number at the wall, and Row 10 gives the average Nusselt number at the outer edge of the layer. These values should be the same, and the difference indicates the failure of the calculation to satisfy the energy balance. This difference is small; it increases with the Reynolds number, reflecting some increase in truncation error. Row 11 gives the average Nusselt number as evaluated for unidimensional conduction. On this basis, the local Nusselt number is <math>h\delta /k=1</math>, and then <math>\frac{{\bar{h}}}{k}=\frac{1}{2\pi }\int_{0}^{2\pi }{\frac{d\xi }{\delta }}</math>  +  The calculations were made using 40 increments in both <math>\xi </math> and <math>\eta </math> for the value of <math>2\pi /\lambda </math>, ''A'', and <math>c/{{\bar{u}}_{0}}</math>. Table 9.1 lists these experimental determinations made by the cited authors. Rows 1 through 7 give the experimental conditions and measurements. Row 8 gives the Prandtl numbers for which the calculations were made. One, on the order of 7, corresponds to the experimental conditions, and the other, 1.7, was selected as a comparison to show the effect of Prandtl number. Row 9 gives the calculated average Nusselt number at the wall, and Row 10 gives the average Nusselt number at the outer edge of the layer. These values should be the same, and the difference indicates the failure of the calculation to satisfy the energy balance. This difference is small; it increases with the Reynolds number, reflecting some increase in truncation error. Row 11 gives the average Nusselt number as evaluated for unidimensional conduction. On this basis, the local Nusselt number is <math>h\delta /k=1</math>, and then 
  +  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\frac{{\bar{h}}}{k}=\frac{1}{2\pi }\int_{0}^{2\pi }{\frac{d\xi }{\delta }}</math></center>  
+  {{EquationRef(28)}}  
+  }  
+  
+  For the sinusoidal wave,  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>\begin{align}  <center><math>\begin{align}  
& \frac{{\bar{h}}}{k}=\frac{1}{2\pi \bar{\delta }}\int_{0}^{2\pi }{\frac{d\xi }{1+A\sin \xi }} \\  & \frac{{\bar{h}}}{k}=\frac{1}{2\pi \bar{\delta }}\int_{0}^{2\pi }{\frac{d\xi }{1+A\sin \xi }} \\  
Line 245:  Line 368:  
{} & =\frac{1}{2\pi \bar{\delta }}\left[ \int_{0}^{\pi }{\frac{d\xi }{1+A\sin \xi }+\int_{0}^{\pi }{\frac{d\xi }{1A\sin \xi }}} \right]=\frac{1}{\bar{\delta }\sqrt{1{{A}^{2}}}} \\  {} & =\frac{1}{2\pi \bar{\delta }}\left[ \int_{0}^{\pi }{\frac{d\xi }{1+A\sin \xi }+\int_{0}^{\pi }{\frac{d\xi }{1A\sin \xi }}} \right]=\frac{1}{\bar{\delta }\sqrt{1{{A}^{2}}}} \\  
\end{matrix} \\  \end{matrix} \\  
  \end{align}</math>  +  \end{align}</math></center> 
+  {{EquationRef(29)}}  
+  }  
The difference between the calculated average Nusselt numbers is due to the contribution of convection and to the twodimensional nature of the conduction that exists because of the variation in layer thickness. The result for ''w'' = ''v'' = ''c'' = 0 corresponds to the twodimensional conduction solution for the wave, and the Nusselt numbers for such a calculation are shown in Rows 12 and 13. These Rows should be identical, and the truncation error in the calculation creates a discrepancy between them, as it does in the cases in which there is fluid motion. These values are essentially the same as Row 11, for unidimensional conduction, and this correspondence shows that twodimensional conduction effects are negligible. Therefore, the difference between the Nusselt numbers in Rows 9 and 11 reflect only the effect of convection.  The difference between the calculated average Nusselt numbers is due to the contribution of convection and to the twodimensional nature of the conduction that exists because of the variation in layer thickness. The result for ''w'' = ''v'' = ''c'' = 0 corresponds to the twodimensional conduction solution for the wave, and the Nusselt numbers for such a calculation are shown in Rows 12 and 13. These Rows should be identical, and the truncation error in the calculation creates a discrepancy between them, as it does in the cases in which there is fluid motion. These values are essentially the same as Row 11, for unidimensional conduction, and this correspondence shows that twodimensional conduction effects are negligible. Therefore, the difference between the Nusselt numbers in Rows 9 and 11 reflect only the effect of convection.  
Line 255:  Line 380:  
Figure 9.13 shows the variation of the local heat transfer coefficient as a function of <math>\xi </math>, normalized with respect to the average value over a period. The results are for the Prandtl number of 8.2 and the Reynolds numbers of 35 and 472 (the same portrayal for the Prandtl number of 1.7 is not much different.) The ratio <math>{{h}_{w}}/\bar{h}</math> for the wall varies considerably for the low Reynolds number, but not very much for the high Reynolds number. For the surface, the ratio <math>{{h}_{\delta }}/\bar{h}</math> varies substantially for both cases. The numerical results for the Nusselt numbers of [[#ReferencesFaghri and Seban (1985)]] are on the order of, but tend to be higher than, the correlation equations of [[#ReferencesKutateladze (1963, 1982)]]. Comparison of Rows 9 and 11 in Table 9.1 reveals the considerable effects of convection and twodimensional conduction.  Figure 9.13 shows the variation of the local heat transfer coefficient as a function of <math>\xi </math>, normalized with respect to the average value over a period. The results are for the Prandtl number of 8.2 and the Reynolds numbers of 35 and 472 (the same portrayal for the Prandtl number of 1.7 is not much different.) The ratio <math>{{h}_{w}}/\bar{h}</math> for the wall varies considerably for the low Reynolds number, but not very much for the high Reynolds number. For the surface, the ratio <math>{{h}_{\delta }}/\bar{h}</math> varies substantially for both cases. The numerical results for the Nusselt numbers of [[#ReferencesFaghri and Seban (1985)]] are on the order of, but tend to be higher than, the correlation equations of [[#ReferencesKutateladze (1963, 1982)]]. Comparison of Rows 9 and 11 in Table 9.1 reveals the considerable effects of convection and twodimensional conduction.  
+  
For practical purposes, [[#ReferencesChun and Seban (1971)]] suggested that the following correlation based on wavy laminar film condensation is valid for wavy laminar falling film evaporation:  For practical purposes, [[#ReferencesChun and Seban (1971)]] suggested that the following correlation based on wavy laminar film condensation is valid for wavy laminar falling film evaporation:  
  
  
  
  An empirical correlation for the local heat transfer coefficient for laminarwavy flow is obtained by combining eq. (  +  { class="wikitable" border="0" 
  +    
  <center><math>\frac{{{h}_{x}}}{{{k}_{\ell }}}{{\left[ \frac{\mu _{\ell }^{2}}{{{\rho }_{\ell }}\left( {{\rho }_{\ell }}{{\rho }_{v}} \right)g} \right]}^{{1}/{3}\;}}=0.828{{\operatorname{Re}}^{0.22}},\text{ }\operatorname{Re}>{{\operatorname{Re}}_{wavy}}</math>  +   width="100%"  
  +  <center><math>{{h}_{x}}=0.876{{\left( \frac{\operatorname{Re}}{4} \right)}^{0.11}}{{h}_{Nusselt}}</math></center>  
+  {{EquationRef(30)}}  
+  }  
+  
+  An empirical correlation for the local heat transfer coefficient for laminarwavy flow is obtained by combining eq. (30) with eq. (9.108),  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\frac{{{h}_{x}}}{{{k}_{\ell }}}{{\left[ \frac{\mu _{\ell }^{2}}{{{\rho }_{\ell }}\left( {{\rho }_{\ell }}{{\rho }_{v}} \right)g} \right]}^{{1}/{3}\;}}=0.828{{\operatorname{Re}}^{0.22}},\text{ }\operatorname{Re}>{{\operatorname{Re}}_{wavy}}</math></center>  
+  {{EquationRef(31)}}  
+  }  
[[#ReferencesChun and Seban (1971)]] suggested the following empirical correlation to predict the onset of wavy laminar flow:  [[#ReferencesChun and Seban (1971)]] suggested the following empirical correlation to predict the onset of wavy laminar flow:  
  +  
  <center><math>\operatorname{Re}\ge {{\operatorname{Re}}_{wavy}}=2.43K{{a}^{1/11}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\operatorname{Re}\ge {{\operatorname{Re}}_{wavy}}=2.43K{{a}^{1/11}}</math></center>  
+  {{EquationRef(32)}}  
+  }  
where Ka is Kapitza number, a ratio of surface tension to viscous force:  where Ka is Kapitza number, a ratio of surface tension to viscous force:  
  
  
  
  +  { class="wikitable" border="0"  
  +    
  <center><math>  +   width="100%"  
  +  <center><math>Ka=\frac{\mu _{\ell }^{4}g}{\left( {{\rho }_{\ell }}{{\rho }_{v}} \right){{\sigma }^{3}}}</math></center>  
+  {{EquationRef(33)}}  
+  }  
  +  The average heat transfer coefficient is often relevant to the practical application and can be obtained by using eq. (9.112), which is also valid for laminar flow with waves. Substituting eq. (31) into eq. (9.112), one obtains an empirical correlation of the average heat transfer coefficient as follows:  
  +  
  <center><math>\operatorname{Re}_{L}^{1.22}=\operatorname{Re}_{0}^{1.22}4\frac{{{k}_{\ell }}L\left( {{T}_{w}}{{T}_{v}} \right)}{{{\mu }_{\ell }}{{h}_{\ell v}}}{{\left[ \frac{{{\rho }_{\ell }}\left( {{\rho }_{\ell }}{{\rho }_{v}} \right)g}{\mu _{\ell }^{2}} \right]}^{{1}/{3}\;}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{\overset{\_}{\mathop{h}}\,}{{{k}_{\ell }}}{{\left[ \frac{\mu _{\ell }^{2}}{{{\rho }_{\ell }}\left( {{\rho }_{\ell }}{{\rho }_{v}} \right)g} \right]}^{{1}/{3}\;}}=\frac{\left( {{\operatorname{Re}}_{o}}{{\operatorname{Re}}_{L}} \right)}{\left( \operatorname{Re}_{o}^{1.22}\operatorname{Re}_{L}^{1.22} \right)}</math></center>  
+  {{EquationRef(34)}}  
+  }  
+  
+  Substituting eq. (9.109) into eq. (34), an equation correlating <math>{{\operatorname{Re}}_{L}}</math>, the Reynolds number at x = L, is obtained:  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\operatorname{Re}_{L}^{1.22}=\operatorname{Re}_{0}^{1.22}4\frac{{{k}_{\ell }}L\left( {{T}_{w}}{{T}_{v}} \right)}{{{\mu }_{\ell }}{{h}_{\ell v}}}{{\left[ \frac{{{\rho }_{\ell }}\left( {{\rho }_{\ell }}{{\rho }_{v}} \right)g}{\mu _{\ell }^{2}} \right]}^{{1}/{3}\;}}</math></center>  
+  {{EquationRef(35)}}  
+  }  
==References==  ==References== 
Revision as of 18:45, 3 June 2010
Wavy flows of thin liquid films have higher heat transfer coefficients than smooth thin films. This effect is due to the former’s greater interfacial surface area and mixing action. Faghri and Seban (1985) analyzed a system in which a liquid film at an initial nondimensional temperature of zero flowed down a vertical wall and evaporated at the free liquidvapor interface. The laminar liquid film thickness varied sinusoidally and the Reynolds number ranged between 35 and 472. The energy equation for a twodimensional situation using the same vertical wall shown in Fig. 9.11 is as follows:

The continuity equation for this system is

Faghri and Seban (1985) assumed a quasiparabolic velocity profile in the xdirection as follows:

where the nondimensional coordinate η is defined as

and is the local mean velocity. Equation (2) can be rearranged and integrated with respect to y to obtain the velocity component v:

It can be seen that eq. (3) has been substituted in order to find an expression for the v component of the velocity where the local mean velocity and film thickness δ are still unknown as functions of the wavy surface. To determine these two variables, the following mass balance can be written as a function of time and space:

The film thickness can be written in terms of its average value, , and a local amplitude φ.

Assuming that the wave formation has a periodic characteristic with a velocity c, the following relations can be written for the fluctuation of local mean velocity and film thickness:


Combining eqs. (6) – (9) and integrating, the following expression is obtained:

where is the average velocity for the average film thickness . For very small local amplitudes, φ < 1, may be approximated by expanding eq. (10) and neglecting the thirdorder terms.

Substituting eqs. (7) and (11) into eq. (5), one obtains the velocity profile in terms of the wavy surface parameters:

Assuming the wave is sinusoidal,

and introducing the nondimensional variable ξ defined as

where A is magnitude of the wave, and λ is the wavelength, eq. (1) is transformed from the (t, x, y) to the (ξ, η) coordinate system,

where the constants are as follows:





These coefficients are evaluated with from eqs. (3) and (11), and v from eq. (12).
In the evaporating film with a fluid with a high latent heat of vaporization, the average film thickness over a wave length, will not vary substantially with distance, x. Furthermore, the temporal average temperature profile will not vary with x far from where heating begins. Therefore, the boundary condition for eq. (15) that describes evaporation on a wavy film is approximated as below. The equality of temperatures at corresponding points in the period gives

Boundary conditions for dimensionless temperature at the heated vertical wall and wavy surface are, respectively,


Therefore, the energy equation with the above boundary conditions can be solved to obtain the wavy condensate film’s temperature profile. With this profile, the heat flux at the wall can be found from

and the heat flux normal to the wavy surface (η = 1) is as follows:

where and are unit vectors in the x and ydirections, respectively. The average heat flow over a period 0 < ξ < 2π is found by integrating eqs. (24) and (25). For a sinusoidal wavy layer, such as the one considered here, this average is identical to the average obtained for the actual surface length. The average heat flux at the wall (η = 0) is approximated by

and the average normal heat flux at the wavy surface is as follows:

Table 9.1 Comparison of wavy film analysis with experimental data (Faghri and Seban, 1985).
References  Rogovan et al. (1969)  Kapitza and Kapitza (1975)  Rogovan et al. (1969)  Rogovan et al. (1969)  

1  T(°C)  20  25  20  20 
2  4Γ / μ  35  92  268  472 
3  (104ft)  4.132  5.96  9.41  11.5 
4  (ft/s)  0.279  0.383  0.836  1.164 
5  2π / λ(ft1)  191  241  159  136 
6  2.12  1.93  1.65  1.55  
7  A  0.30  0.52  0.55  0.40 (0.55) 
8  ν / α  7.2, 1.7  6.2, 1.7  7.2, 1.7  7.2, 1.7, 7.0, 1.7 
9  1.059, 1.057  1.309, 1.296  1.567, 1.543  71.393, 1.372, 1.724, 1.694  
10  1.058, 1.056  1.314, 1.300  1.611, 1.580  1.433, 1.405, 1.816, 1.772  
11  1.05  1.17  1.19  1.09, 1.19  
12  1.049, 1.048  1.171, 1.169  1.199, 1.199  1.092, 1.092, 1.199, 1.199  
13  1.049, 1.048  1.171, 1.169  1.197, 1.197  1.091, 1.091, 1.197, 1.197 
 Conduction. Reprinted with permission from Elsevier.
The calculations were made using 40 increments in both ξ and η for the value of 2π / λ, A, and . Table 9.1 lists these experimental determinations made by the cited authors. Rows 1 through 7 give the experimental conditions and measurements. Row 8 gives the Prandtl numbers for which the calculations were made. One, on the order of 7, corresponds to the experimental conditions, and the other, 1.7, was selected as a comparison to show the effect of Prandtl number. Row 9 gives the calculated average Nusselt number at the wall, and Row 10 gives the average Nusselt number at the outer edge of the layer. These values should be the same, and the difference indicates the failure of the calculation to satisfy the energy balance. This difference is small; it increases with the Reynolds number, reflecting some increase in truncation error. Row 11 gives the average Nusselt number as evaluated for unidimensional conduction. On this basis, the local Nusselt number is hδ / k = 1, and then

For the sinusoidal wave,

The difference between the calculated average Nusselt numbers is due to the contribution of convection and to the twodimensional nature of the conduction that exists because of the variation in layer thickness. The result for w = v = c = 0 corresponds to the twodimensional conduction solution for the wave, and the Nusselt numbers for such a calculation are shown in Rows 12 and 13. These Rows should be identical, and the truncation error in the calculation creates a discrepancy between them, as it does in the cases in which there is fluid motion. These values are essentially the same as Row 11, for unidimensional conduction, and this correspondence shows that twodimensional conduction effects are negligible. Therefore, the difference between the Nusselt numbers in Rows 9 and 11 reflect only the effect of convection.
Figure 9.12 shows the values of the Nusselt number as given by Row 9 of Table 9.1 for a Prandtl number of about 7. The value of 472 for (4Γ / μ) is also shown for an arbitrarilyhigher value of A = 0.55, because the value of A probably should not decrease as the Reynolds number increases. Figure 9.11 also contains lines A – to show the specification of Zazuli given by Kutateladze (1963) – and B – to show that of Kutateladze (1982). Curve C is that of Hirschburg and Florschuetz (1982) for the intermediate wave solution designated by them as f+ = 0.65, which best fits the average of the data. Three data points from Chun and Seban (1971) for evaporation of water with a Prandtl number of about 5.6 are shown by plus (+) symbols, and four numerical results of Faghri and Seban (1985) are shown by circles (◦) in the same figure.
Figure 9.13 shows the variation of the local heat transfer coefficient as a function of ξ, normalized with respect to the average value over a period. The results are for the Prandtl number of 8.2 and the Reynolds numbers of 35 and 472 (the same portrayal for the Prandtl number of 1.7 is not much different.) The ratio for the wall varies considerably for the low Reynolds number, but not very much for the high Reynolds number. For the surface, the ratio varies substantially for both cases. The numerical results for the Nusselt numbers of Faghri and Seban (1985) are on the order of, but tend to be higher than, the correlation equations of Kutateladze (1963, 1982). Comparison of Rows 9 and 11 in Table 9.1 reveals the considerable effects of convection and twodimensional conduction.
For practical purposes, Chun and Seban (1971) suggested that the following correlation based on wavy laminar film condensation is valid for wavy laminar falling film evaporation:

An empirical correlation for the local heat transfer coefficient for laminarwavy flow is obtained by combining eq. (30) with eq. (9.108),

Chun and Seban (1971) suggested the following empirical correlation to predict the onset of wavy laminar flow:

where Ka is Kapitza number, a ratio of surface tension to viscous force:

The average heat transfer coefficient is often relevant to the practical application and can be obtained by using eq. (9.112), which is also valid for laminar flow with waves. Substituting eq. (31) into eq. (9.112), one obtains an empirical correlation of the average heat transfer coefficient as follows:

Substituting eq. (9.109) into eq. (34), an equation correlating , the Reynolds number at x = L, is obtained:

References
Chun, K. R. and Seban, R. A., 1971, “Heat transfer to evaporating liquid films,” ASME Journal of Heat Transfer, Vol. 93, pp. 391396.
Faghri, A. and Seban, R., 1985, “Heat Transfer in Wavy Liquid Films,” International Journal of Heat Mass Transfer, Vol. 28, pp. 506508.
Kutateladze, S.S., 1963, Fundamentals of Heat Transfer, Academic Press Inc., New York.
Kutateladze, S.S., 1982, “SemiEmpirical Theory of Film Condensation of Pure Vapors,” International Journal of Heat Mass Transfer, Vol. 25, pp. 653660.