Thermally developing laminar flow
From ThermalFluidsPedia
(Created page with 'In the previous section, we considered problems where the velocity and temperature profile were fully developed, so that the heat transfer coefficient was constant with distance …') 
Yuwen Zhang (Talk  contribs) 

(11 intermediate revisions not shown)  
Line 1:  Line 1:  
  +  In this article, we consider problems in which only velocity is fully developed at the point where the heat transfer starts. Furthermore, we consider two cases of constant wall temperature and wall heat flux, both by assuming uniform temperature at the inlet. Under these conditions, the heat transfer coefficient is not constant but varies along the tube.  
  Whiteman and Drake  +  
  The following assumptions are made in order to obtain a closed form solution for heat transfer analysis for fully developed flow and developing temperature profile in a circular tube:  +  Whiteman and Drake <ref name="WD1980">Whiteman, I. R., and Drake, W. B., 1980, Trans. ASME, Vol. 80, pp. 728732.</ref>, Lyche and Bird <ref name="LB1956">Lyche, B. C., and Bird, R. B., 1956, “GraetzNusselt Problem for a PowerLaw NonNewtonian Fluid,” Chem. Eng. Sci., Vol. 6, pp. 3541.</ref> and Blackwell <ref name="B1985">Blackwell, B. F., 1985, “Numerical Solution of the Graetz Problem for a Bingham Plastic in Laminar Tube Flow with Constant Wall temperature,” ASME J. Heat Transfer, Vol. 107, pp. 466468.</ref> studied the case of fully developed flow with thermal entry effects for nonNewtonian fluids. Sellars et al. <ref name="S1956">Sellars, J.R., Tribus, M., and Klein, J.S., 1956, “Heat Transfer to Laminar Flow in a Flat Conduit –The Graetz Problem Extended,” Trans. ASME, Vol. 78, pp. 441448.</ref> obtained thermal entry length solutions for the case of a Newtonian fluid with constant wall temperature and fully developed flow, which are presented below. 
  1. Incompressible Newtonian fluid  +  
  2. Laminar flow  +  The following assumptions are made in order to obtain a closed form solution for heat transfer analysis for fully developed flow and developing temperature profile in a circular tube<ref>Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.</ref>:<br> 
  3. Twodimensional steady state  +  1. Incompressible Newtonian fluid<br> 
  4. Axial heat conduction and viscous dissipation are neglected  +  2. Laminar flow<br> 
  5. Constant properties  +  3. Twodimensional steady state<br> 
+  4. Axial heat conduction and viscous dissipation are neglected<br>  
+  5. Constant properties<br>  
+  
This does not mean that one cannot obtain analytical solutions when one or more of the above assumptions is valid, but the solution will be much easier by making the above assumptions.  This does not mean that one cannot obtain analytical solutions when one or more of the above assumptions is valid, but the solution will be much easier by making the above assumptions.  
  Since the fully developed velocity was already obtained in  +  Since the fully developed velocity was already obtained in [[Basics of Internal Forced Convection]], we will focus on the solution of the energy equation and boundary conditions for a developing temperature profile. 
==Constant Wall Temperature==  ==Constant Wall Temperature==  
  The dimensionless energy  +  The dimensionless energy equation and boundary conditions using the above assumptions for the case of constant wall temperature are reduced to 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 26:  Line 29:  
<math>\begin{align} & \theta \left( {{r}^{+}},0 \right)=1 \\ & \theta \left( 1,{{x}^{+}} \right)=0 \\ & \theta \left( 0,{{x}^{+}} \right)=\quad \text{finite}\quad \text{or}\quad \frac{\partial \theta }{\partial {{r}^{+}}}\left( 0,{{x}^{+}} \right)=0 \\ \end{align}</math>  <math>\begin{align} & \theta \left( {{r}^{+}},0 \right)=1 \\ & \theta \left( 1,{{x}^{+}} \right)=0 \\ & \theta \left( 0,{{x}^{+}} \right)=\quad \text{finite}\quad \text{or}\quad \frac{\partial \theta }{\partial {{r}^{+}}}\left( 0,{{x}^{+}} \right)=0 \\ \end{align}</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(2)}} 
}  }  
where  where  
  <math>{{r}^{+}}=\frac{r}{{{r}_{o}}},\quad \theta =\frac{T{{T}_{w}}}{{{T}_{in}}{{T}_{w}}},\quad {{u}^{+}}=\frac{u}{{{u}_{m}}},\quad {{x}^{+}}=\frac{x/{{r}_{0}}}{\operatorname{Re}\Pr }</math>  +  <center><math>{{r}^{+}}=\frac{r}{{{r}_{o}}},\quad \theta =\frac{T{{T}_{w}}}{{{T}_{in}}{{T}_{w}}},\quad {{u}^{+}}=\frac{u}{{{u}_{m}}},\quad {{x}^{+}}=\frac{x/{{r}_{0}}}{\operatorname{Re}\Pr }</math></center> 
For a fully developed laminar flow, the parabolic velocity profile previously developed is applicable, i.e.,  For a fully developed laminar flow, the parabolic velocity profile previously developed is applicable, i.e.,  
  <math>u=2{{u}_{m}}\left( 1\frac{{{r}^{2}}}{r_{o}^{2}} \right)\quad \text{or}\quad {{u}^{+}}=2\left( 1{{r}^{{{+}^{2}}}} \right)</math>  +  <center><math>u=2{{u}_{m}}\left( 1\frac{{{r}^{2}}}{r_{o}^{2}} \right)\quad \text{or}\quad {{u}^{+}}=2\left( 1{{r}^{{{+}^{2}}}} \right)</math></center> 
  Substituting the above equation into the energy eq. (  +  Substituting the above equation into the energy eq. (1), we get 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 43:  Line 46:  
<math>\left( 1{{r}^{{{+}^{2}}}} \right)\frac{\partial \theta }{\partial {{x}^{+}}}=\frac{{{\partial }^{2}}\theta }{\partial {{r}^{{{+}^{2}}}}}+\frac{1}{{{r}^{+}}}\frac{\partial \theta }{\partial {{r}^{+}}}</math>  <math>\left( 1{{r}^{{{+}^{2}}}} \right)\frac{\partial \theta }{\partial {{x}^{+}}}=\frac{{{\partial }^{2}}\theta }{\partial {{r}^{{{+}^{2}}}}}+\frac{1}{{{r}^{+}}}\frac{\partial \theta }{\partial {{r}^{+}}}</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(3)}} 
}  }  
Since the above partial differential equation is linear and homogeneous, one can apply the method of separation of variables. The separation of variables solution is assumed of the form  Since the above partial differential equation is linear and homogeneous, one can apply the method of separation of variables. The separation of variables solution is assumed of the form  
Line 52:  Line 55:  
<math>\theta \left( {{r}^{+}},{{x}^{+}} \right)=R\left( {{r}^{+}} \right)X\left( {{x}^{+}} \right)</math>  <math>\theta \left( {{r}^{+}},{{x}^{+}} \right)=R\left( {{r}^{+}} \right)X\left( {{x}^{+}} \right)</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(4)}} 
}  }  
  The substitution of the above equation into eq. (  +  The substitution of the above equation into eq. (3) yields two ordinary differential equations 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 61:  Line 64:  
<math>{X}'+{{\lambda }^{2}}X=0</math>  <math>{X}'+{{\lambda }^{2}}X=0</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(5)}} 
}  }  
Line 69:  Line 72:  
<math>{R}''+\frac{1}{{{r}^{+}}}{R}'+{{\lambda }^{2}}R\left( 1{{r}^{{{+}^{2}}}} \right)=0</math>  <math>{R}''+\frac{1}{{{r}^{+}}}{R}'+{{\lambda }^{2}}R\left( 1{{r}^{{{+}^{2}}}} \right)=0</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(6)}} 
}  }  
where  where  
  <math>{X}'=\frac{dX}{d{{x}^{+}}},\quad {X}''=\frac{{{d}^{2}}X}{d{{x}^{{{+}^{2}}}}},\quad {R}'=\frac{dR}{d{{r}^{+}}},\quad {R}''=\frac{{{d}^{2}}R}{d{{r}^{{{+}^{2}}}}}</math>  +  <center><math>{X}'=\frac{dX}{d{{x}^{+}}},\quad {X}''=\frac{{{d}^{2}}X}{d{{x}^{{{+}^{2}}}}},\quad {R}'=\frac{dR}{d{{r}^{+}}},\quad {R}''=\frac{{{d}^{2}}R}{d{{r}^{{{+}^{2}}}}}</math></center> 
  and –  +  and –''λ''<sup>2</sup> is the separation constant or eigenvalue. 
  The solution for eq. (5  +  
  <math>{{e}^{{{\lambda }^{2}}{{x}^{+}}}}</math> while the solution of eq. (  +  The solution for eq. (5) is a simple exponential function of the form <math>{{e}^{{{\lambda }^{2}}{{x}^{+}}}}</math> while the solution of eq. (6) is of infinite series referred to by the SturmLiouville theory. The solution is of the form 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 84:  Line 87:  
<math>\theta \left( {{r}^{+}},{{x}^{+}} \right)=\sum\limits_{n=0}^{\infty }{{{c}_{n}}{{R}_{n}}\left( {{r}^{+}} \right)}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)</math>  <math>\theta \left( {{r}^{+}},{{x}^{+}} \right)=\sum\limits_{n=0}^{\infty }{{{c}_{n}}{{R}_{n}}\left( {{r}^{+}} \right)}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(7)}} 
}  }  
  where  +  where ''λ<sub>n</sub>'' are the eigenvalues, ''R<sub>n</sub>'' are the eigenfunctions corresponding to eq. (6), and ''c<sub>n</sub>'' are constants. 
+  
The local heat flux, dimensionless mean temperature, local Nusselt number and mean Nusselt number can be obtained from the following equations, using the above temperature distribution  The local heat flux, dimensionless mean temperature, local Nusselt number and mean Nusselt number can be obtained from the following equations, using the above temperature distribution  
Line 94:  Line 98:  
<math>\begin{align} & {{q}_{w}}^{\prime \prime }=k{{\left. \frac{\partial T}{\partial r} \right}_{r={{r}_{o}}}}=k\frac{\left( {{T}_{w}}{{T}_{in}} \right)}{{{r}_{o}}}{{\left. \frac{\partial \theta }{\partial {{r}^{+}}} \right}_{{{r}^{+}}=1}} \\ & \text{ }=\frac{2k}{{{r}_{o}}}\left( {{T}_{w}}{{T}_{in}} \right)\sum\limits_{n=0}^{\infty }{{{G}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)} \\ \end{align}</math>  <math>\begin{align} & {{q}_{w}}^{\prime \prime }=k{{\left. \frac{\partial T}{\partial r} \right}_{r={{r}_{o}}}}=k\frac{\left( {{T}_{w}}{{T}_{in}} \right)}{{{r}_{o}}}{{\left. \frac{\partial \theta }{\partial {{r}^{+}}} \right}_{{{r}^{+}}=1}} \\ & \text{ }=\frac{2k}{{{r}_{o}}}\left( {{T}_{w}}{{T}_{in}} \right)\sum\limits_{n=0}^{\infty }{{{G}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)} \\ \end{align}</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(8)}} 
}  }  
Line 102:  Line 106:  
<math>{{\theta }_{m}}=\frac{{{T}_{m}}{{T}_{w}}}{{{T}_{in}}{{T}_{w}}}=8\sum\limits_{n=0}^{\infty }{{{G}_{n}}\left[ \frac{\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)}{{{\lambda }_{n}}^{2}} \right]}</math>  <math>{{\theta }_{m}}=\frac{{{T}_{m}}{{T}_{w}}}{{{T}_{in}}{{T}_{w}}}=8\sum\limits_{n=0}^{\infty }{{{G}_{n}}\left[ \frac{\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)}{{{\lambda }_{n}}^{2}} \right]}</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(9)}} 
}  }  
Line 110:  Line 114:  
<math>N{{u}_{x}}=\frac{{{h}_{x}}\left( 2{{r}_{o}} \right)}{k}=\frac{{{q}_{w}}^{\prime \prime }\left( 2{{r}_{o}} \right)}{\left( {{T}_{w}}{{T}_{in}} \right)k{{\theta }_{m}}}=\frac{2}{{{\theta }_{m}}}{{\left. \frac{\partial \theta }{\partial {{r}^{+}}} \right}_{{{r}^{+}}=1}}=\frac{\sum\limits_{n=0}^{\infty }{{{G}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)}}{2\sum\limits_{n=0}^{\infty }{{{G}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)/{{\lambda }_{n}}^{2}}}</math>  <math>N{{u}_{x}}=\frac{{{h}_{x}}\left( 2{{r}_{o}} \right)}{k}=\frac{{{q}_{w}}^{\prime \prime }\left( 2{{r}_{o}} \right)}{\left( {{T}_{w}}{{T}_{in}} \right)k{{\theta }_{m}}}=\frac{2}{{{\theta }_{m}}}{{\left. \frac{\partial \theta }{\partial {{r}^{+}}} \right}_{{{r}^{+}}=1}}=\frac{\sum\limits_{n=0}^{\infty }{{{G}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)}}{2\sum\limits_{n=0}^{\infty }{{{G}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)/{{\lambda }_{n}}^{2}}}</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(10)}} 
}  }  
Line 118:  Line 122:  
<math>\text{N}{{\text{u}}_{\text{m}}}=\frac{{{{\bar{h}}}_{x}}\left( 2{{r}_{o}} \right)}{k}=\frac{1}{{{x}^{+}}}\int_{0}^{{{x}^{+}}}{N{{u}_{x}}d{{x}^{+}}=}\frac{1}{2{{x}^{+}}}\ln \left[ 8\sum\limits_{n=0}^{\infty }{\frac{{{G}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)}{{{\lambda }_{n}}^{2}}} \right]</math>  <math>\text{N}{{\text{u}}_{\text{m}}}=\frac{{{{\bar{h}}}_{x}}\left( 2{{r}_{o}} \right)}{k}=\frac{1}{{{x}^{+}}}\int_{0}^{{{x}^{+}}}{N{{u}_{x}}d{{x}^{+}}=}\frac{1}{2{{x}^{+}}}\ln \left[ 8\sum\limits_{n=0}^{\infty }{\frac{{{G}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)}{{{\lambda }_{n}}^{2}}} \right]</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(11)}} 
}  }  
where  where  
+  <center><math>{{G}_{n}}=\frac{1}{2}{{c}_{n}}{{{R}'}_{n}}\left( 1 \right)</math></center>  
  +  The first five terms in eqs. (8) – (11) are sufficient to provide accurate solutions to the above infinite series. The eigenvalues, ''λ<sub>n</sub>'' and ''G<sub>n</sub>'', used to calculate <math>{{q}_{w}}^{\prime \prime }</math>, ''θ<sub>m</sub>'', ''Nu<sub>x</sub>'' and ''Nu<sub>m</sub>'' for the above problem are presented in the following table.  
  +  
  The first five terms in eqs. (  +  
  <math>{{q}_{w}}^{\prime \prime }</math>  +  
  ,  +  
<center>  <center>  
<div style="display:inline;">  <div style="display:inline;">  
  '''Table 1''' Eigenvalues and Eigenfunctions of a Circular Duct; Thermal Entry Effect with Fully Developed Laminar Flow and Constant Wall Temperature  +  '''Table 1''' Eigenvalues and Eigenfunctions of a Circular Duct; Thermal Entry Effect with Fully Developed Laminar Flow and Constant Wall Temperature 
{ class="wikitable" border="1"  { class="wikitable" border="1"  
   align="center" style="background:#f0f0f0;" width="  +   align="center" style="background:#f0f0f0;" width="10%"  n 
   align="center" style="background:#f0f0f0;" width="  +   align="center" style="background:#f0f0f0;" width="10%"  <math>{{\lambda }_{n}}^{2}/2</math> 
   align="center" style="background:#f0f0f0;" width="  +   align="center" style="background:#f0f0f0;" width="10%"  ''G<sub>n</sub>'' 
    
0  0  
Line 160:  Line 161:  
'''Table 2''' Nusselt Solution for Thermal Entry Effect of a Circular Tube for Fully Developed Laminar Flow and Constant Wall Temperature  '''Table 2''' Nusselt Solution for Thermal Entry Effect of a Circular Tube for Fully Developed Laminar Flow and Constant Wall Temperature  
{ class="wikitable" border="1"  { class="wikitable" border="1"  
   align="center" style="background:#f0f0f0;" width="  +   align="center" style="background:#f0f0f0;" width="15%"  ''x''<sup>+</sup> 
   align="center" style="background:#f0f0f0;" width="  +   align="center" style="background:#f0f0f0;" width="15%"  Nu<sub>x</sub> 
   align="center" style="background:#f0f0f0;" width="  +   align="center" style="background:#f0f0f0;" width="15%"  Nu<sub>m</sub> 
   align="center" style="background:#f0f0f0;" width="  +   align="center" style="background:#f0f0f0;" width="15%"  ''θ<sub>m</sub>'' 
    
0  0  
Line 211:  Line 212:  
}  }  
</div></center>  </div></center>  
  +  The above table provides the variations of ''Nu<sub>x</sub>'', ''Nu<sub>m</sub>'' and ''θ<sub>m</sub>'' with distance along the tube. It can be easily observed from the table that the fully developed temperature profile starts at approximately:  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 218:  Line 219:  
<math>{{x}^{+}}=\frac{x/{{r}_{0}}}{\operatorname{Re}\Pr }=0.1</math>  <math>{{x}^{+}}=\frac{x/{{r}_{0}}}{\operatorname{Re}\Pr }=0.1</math>  
</center>  </center>  
  {{EquationRef(  +  {{EquationRef(12)}} 
}  }  
  Therefore,  +  Therefore, 
  (  +  <center><math>\left( {{L}_{T,T}}/D \right)=0.05\operatorname{Re}\Pr </math></center> 
  where  +  
+  where ''L<sub>T,T</sub>'' is the thermal entrance length for constant wall temperature.  
The thermal entry length increases as both the Reynolds number and Prandtl number increase. A very long thermal entry length is needed for fluids with a high Prandtl number, such as oil. Therefore, care should be taken to make a fully developed temperature profile assumption for fluids with a high Prandtl number.  The thermal entry length increases as both the Reynolds number and Prandtl number increase. A very long thermal entry length is needed for fluids with a high Prandtl number, such as oil. Therefore, care should be taken to make a fully developed temperature profile assumption for fluids with a high Prandtl number.  
==Constant Heat Flux at the Wall==  ==Constant Heat Flux at the Wall==  
+  The laminar fully developed flow with thermal entry length effects (developing temperature profile) for constant wall heat flux is very similar to that in the case of constant wall temperature, except that the dimensionless temperature and boundary conditions are defined as  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%" <center>  
+  <math>\theta =\frac{{{T}_{in}}T}{{{{{q}''}}_{w}}D/k}</math>  
+  </center>  
+  {{EquationRef(13)}}  
+  }  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%" <center>  
+  <math>{{{q}''}_{w}}=k{{\left. \frac{\partial T}{\partial r} \right}_{r={{r}_{o}}}}=\text{constant}</math>  
+  </center>  
+  {{EquationRef(14)}}  
+  }  
+  Siegel et al. <ref name="S1958">Siegel, R., Sparrow, E.M., and Hallman, T.M., 1958, “Steady Laminar Heat Transfer in a Circular Tube with Prescribed Wall Heat Flux,” Appl. Sci. Res., Ser. A, Vol. 7, pp. 386392.</ref> solved the above problem for laminar fully developed flow using separation of variables and the StrumLiouville theory, of which the result is presented below:  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%" <center>  
+  <math>\theta ={{\theta }^{*}}\left( {{r}^{+}},{{x}^{+}} \right)+4{{x}^{+}}+{{r}^{{{+}^{2}}}}\frac{{{r}^{{{+}^{2}}}}}{4}\frac{7}{24}</math>  
+  </center>  
+  {{EquationRef(15)}}  
+  }  
+  where θ* is given below  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%" <center>  
+  <math>{{\theta }^{*}}\left( {{r}^{+}},{{x}^{+}} \right)=\sum\limits_{n=1}^{\infty }{{{c}_{n}}{{R}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right)}</math>  
+  </center>  
+  {{EquationRef(16)}}  
+  }  
+  The eigenvalues, ''λ<sub>n</sub>'', eigenfunctions, ''R<sub>n</sub>'', and constants, ''c<sub>n</sub>'', are presented in the following table.  
+  <center>  
+  <div style="display:inline;">  
+  '''Table 3''' The local friction coefficient for laminar boundary layer flow over a wedge, with U = cx<sup>m</sup> and an impermeable wall  
+  { class="wikitable" border="1"  
+   align="center" style="background:#f0f0f0;" width="15%"  n  
+   align="center" style="background:#f0f0f0;" width="15%"  λn<sub>2</sub>/2  
+   align="center" style="background:#f0f0f0;" width="15%"  Rn(1)  
+   align="center" style="background:#f0f0f0;" width="15%"  c<sub>n</sub>  
+    
+  1  
+  25.6796  
+   0.492517  
+  0.403483  
+    
+  2  
+  83.8618  
+  0.395508  
+   0.175111  
+    
+  3  
+  174.167  
+   0.345872  
+  0.105594  
+    
+  4  
+  296.536  
+  0.314047  
+   0.0732804  
+    
+  5  
+  450.947  
+   0.291252  
+  0.0550357  
+    
+  6  
+  637.387  
+  0.273808  
+   0.043483  
+    
+  7  
+  855.850  
+   0.259852  
+  0.035597  
+  }  
+  </div></center>  
+  The local Nusselt number, based on the above solution, is given below and numerical values are presented in the following table  
+  <center>  
+  <div style="display:inline;">  
+  '''Table 4''' Local Nusselt Number for Thermal Entry Effect with Fully Developed Flow of a Circular Tube with Constant Wall Heat Flux  
+  { class="wikitable" border="1"  
+   align="center" style="background:#f0f0f0;" width="20%" ''x''<sup>+</sup>  
+   align="center" style="background:#f0f0f0;" width="20%" Nu<sub>x</sub>  
+    
+  0  
+  ∞  
+    
+  0.0025  
+  11.5  
+    
+  0.005  
+  9.0  
+    
+  0.01  
+  7.5  
+    
+  0.02  
+  6.1  
+    
+  0.05  
+  5.0  
+    
+  0.1  
+  4.5  
+    
+  0.2  
+  4.364  
+    
+  ∞  
+  4.364  
+  }  
+  </div></center>  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%" <center>  
+  <math>N{{u}_{x}}=\frac{\left( 48/11 \right)}{1+\left( 24/11 \right)\sum\limits_{n=0}^{\infty }{{{c}_{n}}\exp \left( {{\lambda }_{n}}^{2}{{x}^{+}} \right){{R}_{n}}\left( 1 \right)}}</math>  
+  </center>  
+  {{EquationRef(17)}}  
+  }  
+  The thermal entrance for constant wall heat flux based on the numerical results presented in the above table is <math>{{x}^{+}}\approx 0.05</math> or  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%" <center>  
+  <math>\left( {{L}_{T,H}} \right)=0.05\left( \operatorname{Re}\Pr \right)\left( D \right)</math>  
+  </center>  
+  {{EquationRef(18)}}  
+  }  
+  where ''L<sub>T,H</sub>'' is the thermal entry length for fully developed flow with constant wall heat flux.  
+  
+  The mean temperature variation can be obtained from the Nusselt number, eq. (17), using the following equation:  
+  
+  { class="wikitable" border="0"  
+    
+   width="100%" <center>  
+  <math>{{T}_{w}}{{T}_{m}}=\frac{{{q}_{w}}^{\prime \prime }}{{{h}_{x}}}=\frac{{{q}_{w}}^{\prime \prime }D}{N{{u}_{x}}k}</math>  
+  </center>  
+  {{EquationRef(19)}}  
+  }  
+  The thermal entry length solutions presented above for hydrodynamically fully developed flows are based only on either constant wall temperature or constant heat flux. Although the wall temperature or heat flux may be assumed constant along the tube length in the entrance region for certain internal convection problems, there are cases where the wall temperature or heat flux varies considerably with tube length. In these cases, a decision must be made whether to assume constant wall temperature or constant heat flux, use a mean overall wall temperature or heat flux with respect to length, or attempt to take into account the effect of varying wall temperature or heat flux. If the latter is desirable, the solution can be obtained by superposing the thermal entry length solutions at infinitesimal and finite surfacetemperature steps<ref>Kays, W.M., Crawford, M.E., and Weigand, B., 2005, Convective Heat Transfer, 4th ed., McGrawHill, New York, NY</ref>.  
+  
+  ==References==  
+  {{Reflist}} 
Current revision as of 07:27, 27 July 2010
In this article, we consider problems in which only velocity is fully developed at the point where the heat transfer starts. Furthermore, we consider two cases of constant wall temperature and wall heat flux, both by assuming uniform temperature at the inlet. Under these conditions, the heat transfer coefficient is not constant but varies along the tube.
Whiteman and Drake ^{[1]}, Lyche and Bird ^{[2]} and Blackwell ^{[3]} studied the case of fully developed flow with thermal entry effects for nonNewtonian fluids. Sellars et al. ^{[4]} obtained thermal entry length solutions for the case of a Newtonian fluid with constant wall temperature and fully developed flow, which are presented below.
The following assumptions are made in order to obtain a closed form solution for heat transfer analysis for fully developed flow and developing temperature profile in a circular tube^{[5]}:
1. Incompressible Newtonian fluid
2. Laminar flow
3. Twodimensional steady state
4. Axial heat conduction and viscous dissipation are neglected
5. Constant properties
This does not mean that one cannot obtain analytical solutions when one or more of the above assumptions is valid, but the solution will be much easier by making the above assumptions. Since the fully developed velocity was already obtained in Basics of Internal Forced Convection, we will focus on the solution of the energy equation and boundary conditions for a developing temperature profile.
Constant Wall Temperature
The dimensionless energy equation and boundary conditions using the above assumptions for the case of constant wall temperature are reduced to


where
For a fully developed laminar flow, the parabolic velocity profile previously developed is applicable, i.e.,
Substituting the above equation into the energy eq. (1), we get

Since the above partial differential equation is linear and homogeneous, one can apply the method of separation of variables. The separation of variables solution is assumed of the form

The substitution of the above equation into eq. (3) yields two ordinary differential equations
X' + λ^{2}X = 0 

where
and –λ^{2} is the separation constant or eigenvalue.
The solution for eq. (5) is a simple exponential function of the form while the solution of eq. (6) is of infinite series referred to by the SturmLiouville theory. The solution is of the form

where λ_{n} are the eigenvalues, R_{n} are the eigenfunctions corresponding to eq. (6), and c_{n} are constants.
The local heat flux, dimensionless mean temperature, local Nusselt number and mean Nusselt number can be obtained from the following equations, using the above temperature distribution




where
The first five terms in eqs. (8) – (11) are sufficient to provide accurate solutions to the above infinite series. The eigenvalues, λ_{n} and G_{n}, used to calculate , θ_{m}, Nu_{x} and Nu_{m} for the above problem are presented in the following table.
Table 1 Eigenvalues and Eigenfunctions of a Circular Duct; Thermal Entry Effect with Fully Developed Laminar Flow and Constant Wall Temperature
n  G_{n}  
0  3.656  0.749 
1  22.31  0.544 
2  56.9  0.463 
3  107.6  0.414 
4  174.25  0.383 
Table 2 Nusselt Solution for Thermal Entry Effect of a Circular Tube for Fully Developed Laminar Flow and Constant Wall Temperature
x^{+}  Nu_{x}  Nu_{m}  θ_{m} 
0  ∞  ∞  1 
0.001  10.1  15.4  0.940 
0.004  8.06  12.2  0.907 
0.01  6.00  8.94  0.836 
0.04  4.17  5.82  0.628 
0.08  3.79  4.89  0.457 
0.1  3.71  4.64  0.395 
0.2  3.658  4.16  0.190 
∞  3.657  3.657  0 
The above table provides the variations of Nu_{x}, Nu_{m} and θ_{m} with distance along the tube. It can be easily observed from the table that the fully developed temperature profile starts at approximately:

Therefore,
where L_{T,T} is the thermal entrance length for constant wall temperature. The thermal entry length increases as both the Reynolds number and Prandtl number increase. A very long thermal entry length is needed for fluids with a high Prandtl number, such as oil. Therefore, care should be taken to make a fully developed temperature profile assumption for fluids with a high Prandtl number.
Constant Heat Flux at the Wall
The laminar fully developed flow with thermal entry length effects (developing temperature profile) for constant wall heat flux is very similar to that in the case of constant wall temperature, except that the dimensionless temperature and boundary conditions are defined as


Siegel et al. ^{[6]} solved the above problem for laminar fully developed flow using separation of variables and the StrumLiouville theory, of which the result is presented below:

where θ* is given below

The eigenvalues, λ_{n}, eigenfunctions, R_{n}, and constants, c_{n}, are presented in the following table.
Table 3 The local friction coefficient for laminar boundary layer flow over a wedge, with U = cx^{m} and an impermeable wall
n  λn_{2}/2  Rn(1)  c_{n} 
1  25.6796  0.492517  0.403483 
2  83.8618  0.395508  0.175111 
3  174.167  0.345872  0.105594 
4  296.536  0.314047  0.0732804 
5  450.947  0.291252  0.0550357 
6  637.387  0.273808  0.043483 
7  855.850  0.259852  0.035597 
The local Nusselt number, based on the above solution, is given below and numerical values are presented in the following table
Table 4 Local Nusselt Number for Thermal Entry Effect with Fully Developed Flow of a Circular Tube with Constant Wall Heat Flux
x^{+}  Nu_{x} 
0  ∞ 
0.0025  11.5 
0.005  9.0 
0.01  7.5 
0.02  6.1 
0.05  5.0 
0.1  4.5 
0.2  4.364 
∞  4.364 

The thermal entrance for constant wall heat flux based on the numerical results presented in the above table is or

where L_{T,H} is the thermal entry length for fully developed flow with constant wall heat flux.
The mean temperature variation can be obtained from the Nusselt number, eq. (17), using the following equation:

The thermal entry length solutions presented above for hydrodynamically fully developed flows are based only on either constant wall temperature or constant heat flux. Although the wall temperature or heat flux may be assumed constant along the tube length in the entrance region for certain internal convection problems, there are cases where the wall temperature or heat flux varies considerably with tube length. In these cases, a decision must be made whether to assume constant wall temperature or constant heat flux, use a mean overall wall temperature or heat flux with respect to length, or attempt to take into account the effect of varying wall temperature or heat flux. If the latter is desirable, the solution can be obtained by superposing the thermal entry length solutions at infinitesimal and finite surfacetemperature steps^{[7]}.
References
 ↑ Whiteman, I. R., and Drake, W. B., 1980, Trans. ASME, Vol. 80, pp. 728732.
 ↑ Lyche, B. C., and Bird, R. B., 1956, “GraetzNusselt Problem for a PowerLaw NonNewtonian Fluid,” Chem. Eng. Sci., Vol. 6, pp. 3541.
 ↑ Blackwell, B. F., 1985, “Numerical Solution of the Graetz Problem for a Bingham Plastic in Laminar Tube Flow with Constant Wall temperature,” ASME J. Heat Transfer, Vol. 107, pp. 466468.
 ↑ Sellars, J.R., Tribus, M., and Klein, J.S., 1956, “Heat Transfer to Laminar Flow in a Flat Conduit –The Graetz Problem Extended,” Trans. ASME, Vol. 78, pp. 441448.
 ↑ Faghri, A., Zhang, Y., and Howell, J. R., 2010, Advanced Heat and Mass Transfer, Global Digital Press, Columbia, MO.
 ↑ Siegel, R., Sparrow, E.M., and Hallman, T.M., 1958, “Steady Laminar Heat Transfer in a Circular Tube with Prescribed Wall Heat Flux,” Appl. Sci. Res., Ser. A, Vol. 7, pp. 386392.
 ↑ Kays, W.M., Crawford, M.E., and Weigand, B., 2005, Convective Heat Transfer, 4th ed., McGrawHill, New York, NY