Natural convection in annular space between concentric cylinders and spheres
From ThermalFluidsPedia
Yuwen Zhang (Talk  contribs) 
Yuwen Zhang (Talk  contribs) 

Line 1:  Line 1:  
  [[Image:  +  [[Image:Natural convection in horizontal annular space between concentric cylinders.jpgthumb400 pxalt=Natural convection in horizontal annular space between concentric cylinders '''Natural convection in horizontal annular space between concentric cylinders.''']] 
  Natural convection in spaces between long horizontal concentric cylinders  +  Natural convection in spaces between long horizontal concentric cylinders or between spheres is very complicated. The only practical approach to analyze the problem is via numerical solution. The physical model of natural convection in annular space between concentric cylinders is shown in figure on the right. Since the problem is axisymmetric, one only needs to study the right half of the domain. In the coordinate system shown in the figure, the governing equations are 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
    
 width="100%" <center>   width="100%" <center>  
  <math>\frac{{\partial  +  <math>\frac{u}{r}\frac{\partial u}{\partial \theta }+\frac{\partial v}{\partial r}+\frac{v}{r}=0</math> 
</center>  </center>  
{{EquationRef(1)}}  {{EquationRef(1)}}  
}  }  
  
  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
    
 width="100%" <center>   width="100%" <center>  
  <math>g\beta (T{{T}_{ref}})\sin \theta </math>  +  <math>\frac{1}{r}\frac{\partial u}{\partial \theta }+v\frac{\partial u}{\partial r}+\frac{uv}{r}=\frac{1}{\rho r}\frac{\partial p}{\partial \theta }+\nu \left( \frac{1}{{{r}^{2}}}\frac{{{\partial }^{2}}u}{\partial {{\theta }^{2}}}+\frac{1}{r}\frac{\partial u}{\partial r}+\frac{{{\partial }^{2}}u}{\partial {{r}^{2}}}\frac{u}{{{r}^{2}}}+\frac{2}{{{r}^{2}}}\frac{\partial v}{\partial \theta } \right)g\beta (T{{T}_{ref}})\sin \theta </math> 
</center>  </center>  
{{EquationRef(2)}}  {{EquationRef(2)}}  
}  }  
  
  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
    
 width="100%" <center>   width="100%" <center>  
  <math>\  +  <math>\frac{u}{r}\frac{\partial v}{\partial \theta }+v\frac{\partial v}{\partial r}\frac{{{u}^{2}}}{r}=\frac{1}{\rho }\frac{\partial p}{\partial r}+\nu \left( \frac{1}{{{r}^{2}}}\frac{{{\partial }^{2}}v}{\partial {{\theta }^{2}}}+\frac{1}{r}\frac{\partial v}{\partial r}+\frac{{{\partial }^{2}}v}{\partial {{r}^{2}}}\frac{v}{{{r}^{2}}}\frac{2}{{{r}^{2}}}\frac{\partial u}{\partial \theta } \right)+g\beta (T{{T}_{ref}})\sin \theta </math> 
</center>  </center>  
{{EquationRef(3)}}  {{EquationRef(3)}}  
}  }  
  
  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 46:  Line 40:  
    
 width="100%" <center>   width="100%" <center>  
  <math>{{p}_{eff}}=p+\rho g\sin \theta \rho g\cos \theta </math>  +  <big><big><math>{{p}_{eff}}=p+\rho g\sin \theta \rho g\cos \theta </math></big></big> 
</center>  </center>  
{{EquationRef(5)}}  {{EquationRef(5)}}  
Line 64:  Line 58:  
    
 width="100%" <center>   width="100%" <center>  
  <math>u=v=0,\text{ }T={{T}_{i}},\text{ at }r={{r}_{i}}</math>  +  <big><big><math>u=v=0,\text{ }T={{T}_{i}},\text{ at }r={{r}_{i}}</math></big></big> 
</center>  </center>  
{{EquationRef(7)}}  {{EquationRef(7)}}  
Line 72:  Line 66:  
    
 width="100%" <center>   width="100%" <center>  
  <math>u=v=0,\text{ }T={{T}_{o}},\text{ at }r={{r}_{o}}</math>  +  <big><big><math>u=v=0,\text{ }T={{T}_{o}},\text{ at }r={{r}_{o}}</math></big></big> 
</center>  </center>  
{{EquationRef(8)}}  {{EquationRef(8)}}  
}  }  
  [[#ReferencesDate (1986)]] solved this problem numerically using a modified SIMPLE algorithm for <math>\delta /{{D}_{i}}=0.8\text{ and }0.15</math>.  +  [[#ReferencesDate (1986)]] solved this problem numerically using a modified SIMPLE algorithm for <math>\delta /{{D}_{i}}=0.8\text{ and }0.15</math>. The streamlines and isotherms obtained by [[#ReferencesDate (1986)]] agreed very well with results obtained by [[#ReferencesKuehn and Goldstein (1976)]] using vortexstream function method. 
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  +  
  [[  +  
  +  
  +  
The rate of heat transfer per unit length of the annulus can be calculated by the following correlation ([[#ReferencesRaithby and Hollands, 1975]]; [[#ReferencesBejan, 2004]]):  The rate of heat transfer per unit length of the annulus can be calculated by the following correlation ([[#ReferencesRaithby and Hollands, 1975]]; [[#ReferencesBejan, 2004]]):  
Line 108:  Line 80:  
<math>{q}'\cong \frac{2.425k({{T}_{i}}{{T}_{o}})}{{{[1+{{({{D}_{i}}/{{D}_{o}})}^{3/5}}]}^{5/4}}}{{\left( \frac{\Pr \text{R}{{\text{a}}_{{{D}_{i}}}}}{0.861+\Pr } \right)}^{1/4}}</math>  <math>{q}'\cong \frac{2.425k({{T}_{i}}{{T}_{o}})}{{{[1+{{({{D}_{i}}/{{D}_{o}})}^{3/5}}]}^{5/4}}}{{\left( \frac{\Pr \text{R}{{\text{a}}_{{{D}_{i}}}}}{0.861+\Pr } \right)}^{1/4}}</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(9)}} 
}  }  
Line 118:  Line 90:  
<math>\text{R}{{\text{a}}_{{{D}_{i}}}}=\frac{g\beta ({{T}_{i}}{{T}_{o}})D_{i}^{3}}{\nu \alpha }</math>  <math>\text{R}{{\text{a}}_{{{D}_{i}}}}=\frac{g\beta ({{T}_{i}}{{T}_{o}})D_{i}^{3}}{\nu \alpha }</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(10)}} 
}  }  
  Equation (  +  Equation (9) is valid for <math>0.7<\Pr <6000\text{ and Ra}<{{10}^{7}}</math>. The thermophysical properties of the fluid should be evaluated at the mean temperature <math>({{T}_{i}}+{{T}_{o}})/2</math>. The scales of the thermal boundary layer on the inner surface of the outer cylinder and on the outer surface of the inner cylinder are:<math>{{\delta }_{o}}\tilde{\ }{{D}_{o}}Ra_{{{D}_{o}}}^{1/4},\text{ }{{\delta }_{i}}\tilde{\ }{{D}_{i}}Ra_{{{D}_{i}}}^{1/4}</math> 
  It is obvious that <math>{{\delta }_{o}}>{{\delta }_{i}}</math> since Do > Di. Equation (  +  It is obvious that <math>{{\delta }_{o}}>{{\delta }_{i}}</math> since Do > Di. Equation (9) will be valid only if the boundary layer thickness is less than the gap between the two cylinders, i.e. only if δo < Do – Di. Under lower Rayleigh numbers, on the other hand, we have 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 130:  Line 102:  
<math>{{D}_{o}}\text{Ra}_{{{D}_{o}}}^{1/4}>{{D}_{o}}{{D}_{i}}</math>  <math>{{D}_{o}}\text{Ra}_{{{D}_{o}}}^{1/4}>{{D}_{o}}{{D}_{i}}</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(11)}} 
}  }  
  and the heat transfer mechanism between two cylinders will approach pure conduction. In this case, the heat transfer rate obtained from eq. (  +  and the heat transfer mechanism between two cylinders will approach pure conduction. In this case, the heat transfer rate obtained from eq. (9) may be less than that obtained from the pure conduction model, which does not make sense. Instead of using eq. (11) to check the validity of eq. (9), an alternative method is to calculate the heat transfer rate via eq. (9) and pure conduction model, and the larger of the two heat transfer rates should be used. 
For natural convection in the annulus between two concentric spheres, the trends for the evolution of the flow pattern and isotherms are similar to the concentric cylinder except the circulation between concentric spheres has the shape of a doughnut. The empirical correlation for the heat transfer rate is ([[#ReferencesRaithby and Hollands, 1975]]; [[#ReferencesBejan, 2004]]):  For natural convection in the annulus between two concentric spheres, the trends for the evolution of the flow pattern and isotherms are similar to the concentric cylinder except the circulation between concentric spheres has the shape of a doughnut. The empirical correlation for the heat transfer rate is ([[#ReferencesRaithby and Hollands, 1975]]; [[#ReferencesBejan, 2004]]):  
  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
    
Line 143:  Line 114:  
<math>q\cong \frac{2.325k{{D}_{i}}({{T}_{i}}{{T}_{o}})}{{{[1+{{({{D}_{i}}/{{D}_{o}})}^{7/5}}]}^{5/4}}}{{\left( \frac{\Pr \text{R}{{\text{a}}_{{{D}_{i}}}}}{0.861+\Pr } \right)}^{1/4}}</math>  <math>q\cong \frac{2.325k{{D}_{i}}({{T}_{i}}{{T}_{o}})}{{{[1+{{({{D}_{i}}/{{D}_{o}})}^{7/5}}]}^{5/4}}}{{\left( \frac{\Pr \text{R}{{\text{a}}_{{{D}_{i}}}}}{0.861+\Pr } \right)}^{1/4}}</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(12)}} 
}  }  
  where the definition of Rayleigh number is the same as for eq. (  +  where the definition of Rayleigh number is the same as for eq. (11). Equation (12) is valid for 0.7<Pr<4000 and Ra<10<sup>4</sup>. The thermophysical properties of the fluid should be evaluated at the mean temperature, <math>({{T}_{i}}+{{T}_{o}})/2</math>. Similar to the annulus between two cylinders, one should also make use of both eq. (12) and the pure conduction model and retain the larger of the two values obtained. 
==References==  ==References==  
Line 153:  Line 124:  
Date, A.W., 1986, “Numerical Prediction of Natural Convection Heat Transfer in Horizontal Annulus,” International Journal of Heat and Mass Transfer, Vol. 29, pp. 14571464.  Date, A.W., 1986, “Numerical Prediction of Natural Convection Heat Transfer in Horizontal Annulus,” International Journal of Heat and Mass Transfer, Vol. 29, pp. 14571464.  
+  
+  Faghri, A., Zhang, Y., and Howell, J. R., 2010, ''Advanced Heat and Mass Transfer'', Global Digital Press, Columbia, MO.  
Kuehn, T.H., and Goldstein, R.J., 1976, “An Experimental and Theoretical Study of Natural Convection in the Annulus between Horizontal Concentric Cylinders,” J. Fluid Mech., Vol. 74, pp. 695719.  Kuehn, T.H., and Goldstein, R.J., 1976, “An Experimental and Theoretical Study of Natural Convection in the Annulus between Horizontal Concentric Cylinders,” J. Fluid Mech., Vol. 74, pp. 695719.  
Raithby, G.D., and Hollands, K.G.T., 1975, “A General Method of Obtaining Approximate Solutions to Laminar and Turbulent Free Convection Problems,” Irvine, T.F., and Hartnett, J.P., Eds., Advances in Heat Transfer, Vol. 11, pp. 265315, Academic Press, New York, NY.  Raithby, G.D., and Hollands, K.G.T., 1975, “A General Method of Obtaining Approximate Solutions to Laminar and Turbulent Free Convection Problems,” Irvine, T.F., and Hartnett, J.P., Eds., Advances in Heat Transfer, Vol. 11, pp. 265315, Academic Press, New York, NY. 
Current revision as of 12:56, 12 July 2010
Natural convection in spaces between long horizontal concentric cylinders or between spheres is very complicated. The only practical approach to analyze the problem is via numerical solution. The physical model of natural convection in annular space between concentric cylinders is shown in figure on the right. Since the problem is axisymmetric, one only needs to study the right half of the domain. In the coordinate system shown in the figure, the governing equations are




where T_{ref} is a reference temperature and:
p_{eff} = p + ρgsinθ − ρgcosθ 
is the effective pressure. Equations (1) – (4) are subject to the following boundary conditions:

u = v = 0, T = T_{i}, at r = r_{i} 
u = v = 0, T = T_{o}, at r = r_{o} 
Date (1986) solved this problem numerically using a modified SIMPLE algorithm for δ / D_{i} = 0.8 and 0.15. The streamlines and isotherms obtained by Date (1986) agreed very well with results obtained by Kuehn and Goldstein (1976) using vortexstream function method.
The rate of heat transfer per unit length of the annulus can be calculated by the following correlation (Raithby and Hollands, 1975; Bejan, 2004):

where the Rayleigh number is defined as

Equation (9) is valid for . The thermophysical properties of the fluid should be evaluated at the mean temperature (T_{i} + T_{o}) / 2. The scales of the thermal boundary layer on the inner surface of the outer cylinder and on the outer surface of the inner cylinder are:
It is obvious that δ_{o} > δ_{i} since Do > Di. Equation (9) will be valid only if the boundary layer thickness is less than the gap between the two cylinders, i.e. only if δo < Do – Di. Under lower Rayleigh numbers, on the other hand, we have

and the heat transfer mechanism between two cylinders will approach pure conduction. In this case, the heat transfer rate obtained from eq. (9) may be less than that obtained from the pure conduction model, which does not make sense. Instead of using eq. (11) to check the validity of eq. (9), an alternative method is to calculate the heat transfer rate via eq. (9) and pure conduction model, and the larger of the two heat transfer rates should be used.
For natural convection in the annulus between two concentric spheres, the trends for the evolution of the flow pattern and isotherms are similar to the concentric cylinder except the circulation between concentric spheres has the shape of a doughnut. The empirical correlation for the heat transfer rate is (Raithby and Hollands, 1975; Bejan, 2004):

where the definition of Rayleigh number is the same as for eq. (11). Equation (12) is valid for 0.7<Pr<4000 and Ra<10^{4}. The thermophysical properties of the fluid should be evaluated at the mean temperature, (T_{i} + T_{o}) / 2. Similar to the annulus between two cylinders, one should also make use of both eq. (12) and the pure conduction model and retain the larger of the two values obtained.
References
Bejan, A., 2004, Convection Heat Transfer, 3rd ed., John Wiley & Sons, New York.
Date, A.W., 1986, “Numerical Prediction of Natural Convection Heat Transfer in Horizontal Annulus,” International Journal of Heat and Mass Transfer, Vol. 29, pp. 14571464.
Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.
Kuehn, T.H., and Goldstein, R.J., 1976, “An Experimental and Theoretical Study of Natural Convection in the Annulus between Horizontal Concentric Cylinders,” J. Fluid Mech., Vol. 74, pp. 695719.
Raithby, G.D., and Hollands, K.G.T., 1975, “A General Method of Obtaining Approximate Solutions to Laminar and Turbulent Free Convection Problems,” Irvine, T.F., and Hartnett, J.P., Eds., Advances in Heat Transfer, Vol. 11, pp. 265315, Academic Press, New York, NY.