Turbulent falling film evaporation
From ThermalFluidsPedia
(Created page with 'Stephan (1992) presented an empirical relationship between Reynolds number and Prandtl number that could be used to determine whether a falling film is completely…') 
Yuwen Zhang (Talk  contribs) 

(7 intermediate revisions not shown)  
Line 1:  Line 1:  
[[#ReferencesStephan (1992)]] presented an empirical relationship between Reynolds number and Prandtl number that could be used to determine whether a falling film is completely turbulent:  [[#ReferencesStephan (1992)]] presented an empirical relationship between Reynolds number and Prandtl number that could be used to determine whether a falling film is completely turbulent:  
  +  
  <center><math>\operatorname{Re}\ge {{\operatorname{Re}}_{turb}}=5840{{\Pr }^{1.05}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\operatorname{Re}\ge {{\operatorname{Re}}_{turb}}=5840{{\Pr }^{1.05}}</math></center>  
+  {{EquationRef(1)}}  
+  }  
In order to model turbulent flow, the following development from [[#ReferencesSeban and Faghri (1976)]] proceeds with assumptions similar to those made in the classical Nusselt laminar analysis. The momentum equation is  In order to model turbulent flow, the following development from [[#ReferencesSeban and Faghri (1976)]] proceeds with assumptions similar to those made in the classical Nusselt laminar analysis. The momentum equation is  
  +  
  <center><math>\frac{d}{dy}\left( {{v}_{\ell }}+{{\varepsilon }_{M}} \right)\frac{du}{dy}+g=0</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{d}{dy}\left( {{v}_{\ell }}+{{\varepsilon }_{M}} \right)\frac{du}{dy}+g=0</math></center>  
+  {{EquationRef(2)}}  
+  }  
where <math>{{\varepsilon }_{M}}</math> is the turbulent eddy diffusivity.  where <math>{{\varepsilon }_{M}}</math> is the turbulent eddy diffusivity.  
Defining the following dimensionless variables,  Defining the following dimensionless variables,  
  +  
  <center><math>{{u}_{\tau }}=\sqrt{{{{\tau }_{w}}}/{{{\rho }_{\ell }}}\;},\ \text{ }{{u}^{+}}={u}/{{{u}_{\tau }},\ \ {{y}^{+}}={y{{u}_{\tau }}}/{{{v}_{\ell }}}\;}\;</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{u}_{\tau }}=\sqrt{{{{\tau }_{w}}}/{{{\rho }_{\ell }}}\;},\ \text{ }{{u}^{+}}={u}/{{{u}_{\tau }},\ \ {{y}^{+}}={y{{u}_{\tau }}}/{{{v}_{\ell }}}\;}\;</math></center>  
+  {{EquationRef(3)}}  
+  }  
the momentum equation becomes  the momentum equation becomes  
  +  
  <center><math>\frac{d}{d{{y}^{+}}}\left( \frac{{{\varepsilon }_{M}}}{{{v}_{\ell }}}+1 \right)\frac{d{{u}^{+}}}{d{{y}^{+}}}+\frac{g{{v}_{\ell }}}{u_{\tau }^{3}}=0</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>\frac{d}{d{{y}^{+}}}\left( \frac{{{\varepsilon }_{M}}}{{{v}_{\ell }}}+1 \right)\frac{d{{u}^{+}}}{d{{y}^{+}}}+\frac{g{{v}_{\ell }}}{u_{\tau }^{3}}=0</math></center>  
+  {{EquationRef(4)}}  
+  }  
which is subject to the following boundary conditions:  which is subject to the following boundary conditions:  
  +  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>{{u}^{+}}=0\begin{matrix}  <center><math>{{u}^{+}}=0\begin{matrix}  
, & {{y}^{+}}=0 \\  , & {{y}^{+}}=0 \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
  +  {{EquationRef(5)}}  
  +  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>\frac{d{{u}^{+}}}{d{{y}^{+}}}=0{{\begin{matrix}  <center><math>\frac{d{{u}^{+}}}{d{{y}^{+}}}=0{{\begin{matrix}  
, & {{y}^{+}}=\delta \\  , & {{y}^{+}}=\delta \\  
  \end{matrix}}^{+}}</math>  +  \end{matrix}}^{+}}</math></center> 
  +  {{EquationRef(6)}}  
+  }  
  Integrating eq. (  +  Integrating eq. (4) twice and using eqs. (5) and (6) to determine the integral constant, the dimensionless velocity becomes 
  +  
  <center><math>{{u}^{+}}=\int_{0}^{{{y}^{+}}}{\frac{\left( 1{{{y}^{+}}}/{{{\delta }^{+}}}\; \right)}{\left( 1+{{{\varepsilon }_{M}}}/{{{v}_{\ell }}}\; \right)}d{{y}^{+}}}</math>  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>{{u}^{+}}=\int_{0}^{{{y}^{+}}}{\frac{\left( 1{{{y}^{+}}}/{{{\delta }^{+}}}\; \right)}{\left( 1+{{{\varepsilon }_{M}}}/{{{v}_{\ell }}}\; \right)}d{{y}^{+}}}</math></center>  
+  {{EquationRef(7)}}  
+  }  
The energy equation is  The energy equation is  
  
  
  
  where <math>{{\varepsilon }_{H}}</math> is turbulent eddy thermal diffusivity. Equation (  +  { class="wikitable" border="0" 
  +    
+   width="100%"   
+  <center><math>u\frac{\partial T}{\partial x}=\frac{\partial }{\partial y}\left( {{\varepsilon }_{H}}+\alpha \right)\frac{\partial T}{\partial y}</math></center>  
+  {{EquationRef(8)}}  
+  }  
+  
+  where <math>{{\varepsilon }_{H}}</math> is turbulent eddy thermal diffusivity. Equation (8) is subject to the following boundary conditions:  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>\begin{matrix}  <center><math>\begin{matrix}  
T={{T}_{sat}} & x=0 \\  T={{T}_{sat}} & x=0 \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
  +  {{EquationRef(9)}}  
  +  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>\frac{\partial T}{\partial y}=\frac{{{{{q}''}}_{w}}}{{{k}_{\ell }}}\begin{matrix}  <center><math>\frac{\partial T}{\partial y}=\frac{{{{{q}''}}_{w}}}{{{k}_{\ell }}}\begin{matrix}  
, & y=0 \\  , & y=0 \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
  +  {{EquationRef(10)}}  
  +  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
<center><math>T={{T}_{sat}}\begin{matrix}  <center><math>T={{T}_{sat}}\begin{matrix}  
, & y=\delta \\  , & y=\delta \\  
  \end{matrix}</math>  +  \end{matrix}</math></center> 
  +  {{EquationRef(11)}}  
+  }  
By defining dimensionless temperature as  By defining dimensionless temperature as  
  
  
  
  and using the dimensionless variables from eq. (  +  { class="wikitable" border="0" 
  +    
  <center><math>{{u}^{+}}\frac{\partial T}{\partial {{x}^{+}}}=\frac{\partial }{\partial {{y}^{+}}}\left( \frac{\alpha }{{{v}_{\ell }}}+\frac{{{\varepsilon }_{H}}}{{{v}_{\ell }}} \right)\frac{\partial {{T}^{+}}}{\partial {{y}^{+}}}</math>  +   width="100%"  
  +  <center><math>{{T}^{+}}=\frac{T{{\rho }_{\ell }}c{{u}_{\tau }}}{{{{{q}''}}_{w}}}</math></center>  
+  {{EquationRef(12)}}  
+  }  
+  
+  and using the dimensionless variables from eq. (3), the energy equation becomes  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>{{u}^{+}}\frac{\partial T}{\partial {{x}^{+}}}=\frac{\partial }{\partial {{y}^{+}}}\left( \frac{\alpha }{{{v}_{\ell }}}+\frac{{{\varepsilon }_{H}}}{{{v}_{\ell }}} \right)\frac{\partial {{T}^{+}}}{\partial {{y}^{+}}}</math></center>  
+  {{EquationRef(13)}}  
+  }  
The solution requires choosing a turbulence model for the specification of <math>{{\varepsilon }_{M}}/\nu </math> and of <math>{{\varepsilon }_{H}}/{{\varepsilon }_{M}}.</math> The detailed turbulence model and solution procedure can be found in [[#ReferencesSeban and Faghri (1976)]]. The predicted result agrees very well with the experimental result of [[#ReferencesChun and Seban (1971)]] for constant heat flux at wall. To predict local heat transfer in turbulent falling film flow on a wall heated under constant heat flux, the following empirical correlation, recommended by [[#ReferencesChun and Seban (1971)]], can be used:  The solution requires choosing a turbulence model for the specification of <math>{{\varepsilon }_{M}}/\nu </math> and of <math>{{\varepsilon }_{H}}/{{\varepsilon }_{M}}.</math> The detailed turbulence model and solution procedure can be found in [[#ReferencesSeban and Faghri (1976)]]. The predicted result agrees very well with the experimental result of [[#ReferencesChun and Seban (1971)]] for constant heat flux at wall. To predict local heat transfer in turbulent falling film flow on a wall heated under constant heat flux, the following empirical correlation, recommended by [[#ReferencesChun and Seban (1971)]], can be used:  
  
  
  
  +  { class="wikitable" border="0"  
  +    
  <center><math>\frac  +   width="100%"  
  +  <center><math>\frac{h}{{{k}_{\ell }}}{{\left[ \frac{\mu _{\ell }^{2}}{{{\rho }_{\ell }}\left( {{\rho }_{\ell }}{{\rho }_{v}} \right)g} \right]}^{{1}/{3}\;}}=3.8\times {{10}^{3}}{{\operatorname{Re}}^{0.4}}{{\Pr }^{0.65}}</math></center>  
+  {{EquationRef(14)}}  
+  }  
  As was demonstrated in Example 9.4, eq. (  +  The average heat transfer coefficient can be obtained by substituting eq. (14) into eq. (9.112), i.e., 
  +  
  <center><math>\operatorname{Re}_{L}^{0.6}=\operatorname{Re}_{0}^{0.6}9.12\times {{10}^{3}}\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}\;}}=2.28\times {{10}^{3}}\frac{\left( {{\operatorname{Re}}_{o}}{{\operatorname{Re}}_{L}} \right)}{\left( \operatorname{Re}_{o}^{0.6}\operatorname{Re}_{L}^{0.6} \right)}{{\Pr }^{0.65}}</math></center>  
+  {{EquationRef(15)}}  
+  }  
+  
+  As was demonstrated in Example 9.4, eq. (15) is valid for cases with either constant wall temperature or constant heat flux. However, the methods for determining the Reynolds number at x = L will vary for different thermal boundary conditions on the wall. For cases with constant wall temperature, the Reynolds number at ''x'' = ''L'', <math>{{\operatorname{Re}}_{L}}</math>, is obtained by substituting eq. (9.109) into eq. (15), i.e.,  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%"   
+  <center><math>\operatorname{Re}_{L}^{0.6}=\operatorname{Re}_{0}^{0.6}9.12\times {{10}^{3}}\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(16)}}  
+  }  
  
While gravitational force drives the falling film evaporation discussed here, liquid flow in a liquid film can also be driven by centrifugal force. [[#ReferencesRahman and Faghri (1992)]] analyzed the processes of heating and evaporation in a thin liquid film adjacent to a horizontal disk rotating about a vertical axis at a constant angular velocity. The fluid emanates axisymmetrically from a source at the center of the disk and then is carried downstream by inertial and centrifugal forces. Closedform analytical solutions were derived for fully developed flow and heat transfer. Simplified analyses were also presented for developing heat transfer in a fullydeveloped flow field. [[#ReferencesRice et al. (2005)]] provided a detailed analysis for evaporation from a thin liquid film on a rotating disk including conjugate effects.  While gravitational force drives the falling film evaporation discussed here, liquid flow in a liquid film can also be driven by centrifugal force. [[#ReferencesRahman and Faghri (1992)]] analyzed the processes of heating and evaporation in a thin liquid film adjacent to a horizontal disk rotating about a vertical axis at a constant angular velocity. The fluid emanates axisymmetrically from a source at the center of the disk and then is carried downstream by inertial and centrifugal forces. Closedform analytical solutions were derived for fully developed flow and heat transfer. Simplified analyses were also presented for developing heat transfer in a fullydeveloped flow field. [[#ReferencesRice et al. (2005)]] provided a detailed analysis for evaporation from a thin liquid film on a rotating disk including conjugate effects.  
==References==  ==References==  
+  
+  Faghri, A., and Zhang, Y., 2006, ''Transport Phenomena in Multiphase Systems'', Elsevier, Burlington, MA  
+  
+  Faghri, A., Zhang, Y., and Howell, J. R., 2010, ''Advanced Heat and Mass Transfer'', Global Digital Press, Columbia, MO.  
+  
+  Chun, K. R. and Seban, R. A., 1971, “Heat transfer to evaporating liquid films,” ASME Journal of Heat Transfer, Vol. 93, pp. 391396.  
+  
+  Rahman, M.M., and Faghri, A., 1992, “Analysis of Heating and Evaporation from a Liquid Film Adjacent to a Horizontal Rotating Disk,” International Journal of Heat and Mass Transfer, Vol. 35, pp. 26442655.  
+  
+  Rice, J., Faghri, A., and Cetegen, B.M., 2005, “Analysis of a Free Surface Flow for a Controlled Liquid Impinging Jet over a Rotating Disk Including Conjugate Effects, with and without Evaporation,” International Journal of Heat and Mass Transfer, Vol. 48, pp. 51925204.  
+  
+  Seban, R. A., and Faghri, A., 1976, “Evaporation and Heating with Turbulent Falling Liquid Films,” ASME Journal of Heat Transfer, Vol. 98, pp. 315318.  
+  
+  Stephan, K., 1992, Heat Transfer in Condensation and Boiling, SpringerVerlag, New York. 
Current revision as of 08:50, 17 July 2010
Stephan (1992) presented an empirical relationship between Reynolds number and Prandtl number that could be used to determine whether a falling film is completely turbulent:

In order to model turbulent flow, the following development from Seban and Faghri (1976) proceeds with assumptions similar to those made in the classical Nusselt laminar analysis. The momentum equation is

where is the turbulent eddy diffusivity.
Defining the following dimensionless variables,

the momentum equation becomes

which is subject to the following boundary conditions:


Integrating eq. (4) twice and using eqs. (5) and (6) to determine the integral constant, the dimensionless velocity becomes

The energy equation is

where is turbulent eddy thermal diffusivity. Equation (8) is subject to the following boundary conditions:



By defining dimensionless temperature as

and using the dimensionless variables from eq. (3), the energy equation becomes

The solution requires choosing a turbulence model for the specification of and of The detailed turbulence model and solution procedure can be found in Seban and Faghri (1976). The predicted result agrees very well with the experimental result of Chun and Seban (1971) for constant heat flux at wall. To predict local heat transfer in turbulent falling film flow on a wall heated under constant heat flux, the following empirical correlation, recommended by Chun and Seban (1971), can be used:

The average heat transfer coefficient can be obtained by substituting eq. (14) into eq. (9.112), i.e.,

As was demonstrated in Example 9.4, eq. (15) is valid for cases with either constant wall temperature or constant heat flux. However, the methods for determining the Reynolds number at x = L will vary for different thermal boundary conditions on the wall. For cases with constant wall temperature, the Reynolds number at x = L, , is obtained by substituting eq. (9.109) into eq. (15), i.e.,

While gravitational force drives the falling film evaporation discussed here, liquid flow in a liquid film can also be driven by centrifugal force. Rahman and Faghri (1992) analyzed the processes of heating and evaporation in a thin liquid film adjacent to a horizontal disk rotating about a vertical axis at a constant angular velocity. The fluid emanates axisymmetrically from a source at the center of the disk and then is carried downstream by inertial and centrifugal forces. Closedform analytical solutions were derived for fully developed flow and heat transfer. Simplified analyses were also presented for developing heat transfer in a fullydeveloped flow field. Rice et al. (2005) provided a detailed analysis for evaporation from a thin liquid film on a rotating disk including conjugate effects.
References
Faghri, A., and Zhang, Y., 2006, Transport Phenomena in Multiphase Systems, Elsevier, Burlington, MA
Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.
Chun, K. R. and Seban, R. A., 1971, “Heat transfer to evaporating liquid films,” ASME Journal of Heat Transfer, Vol. 93, pp. 391396.
Rahman, M.M., and Faghri, A., 1992, “Analysis of Heating and Evaporation from a Liquid Film Adjacent to a Horizontal Rotating Disk,” International Journal of Heat and Mass Transfer, Vol. 35, pp. 26442655.
Rice, J., Faghri, A., and Cetegen, B.M., 2005, “Analysis of a Free Surface Flow for a Controlled Liquid Impinging Jet over a Rotating Disk Including Conjugate Effects, with and without Evaporation,” International Journal of Heat and Mass Transfer, Vol. 48, pp. 51925204.
Seban, R. A., and Faghri, A., 1976, “Evaporation and Heating with Turbulent Falling Liquid Films,” ASME Journal of Heat Transfer, Vol. 98, pp. 315318.
Stephan, K., 1992, Heat Transfer in Condensation and Boiling, SpringerVerlag, New York.