Basics of internal forced convection
From ThermalFluidsPedia
Line 30:  Line 30:  
[[Image:Fig5.2.pngthumb400 pxalt=Velocity profiles and friction factor in turbulent flow in a circular tube  Velocity profiles and friction factor in turbulent flow in a circular tube.]]  [[Image:Fig5.2.pngthumb400 pxalt=Velocity profiles and friction factor in turbulent flow in a circular tube  Velocity profiles and friction factor in turbulent flow in a circular tube.]]  
  The turbulent velocity profile and friction coefficient variation for a circular pipe are shown in figure to the right. Even for a very high inlet velocity, the boundary layer will be laminar over a part of the entrance. This transition from laminar to turbulent is clearly shown by the sudden increase in momentum boundary layer thickness as shown in  +  The turbulent velocity profile and friction coefficient variation for a circular pipe are shown in figure to the right. Even for a very high inlet velocity, the boundary layer will be laminar over a part of the entrance. This transition from laminar to turbulent is clearly shown by the sudden increase in momentum boundary layer thickness as shown in part (a) of the figure. The friction coefficient variation for turbulent flow in a pipe entrance is shown in part (b). 
+  
The hydrodynamic entry length required for fully developed flow should be obtained by a complete solution of the flow and thermal field in the entrance region. A rule of thumb to judge whether or not the flow is fully developed for circular pipes is  The hydrodynamic entry length required for fully developed flow should be obtained by a complete solution of the flow and thermal field in the entrance region. A rule of thumb to judge whether or not the flow is fully developed for circular pipes is  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 46:  Line 47:  
 {{EquationRef (3)}}   {{EquationRef (3)}}  
}  }  
  where  +  where ''L<sub>H</sub>'' is the hydrodynamic length and the Reynolds number is defined by<math>\operatorname{Re}=\frac{u_{m}D}{\nu }</math> 
  [[Image:Fig5.3.pngthumb400 pxalt=Temperature development along the flow in a circular tube   +  [[Image:Fig5.3.pngthumb400 pxalt=Temperature development along the flow in a circular tube Temperature development along the flow in a circular tube.]] 
+  
+  A similar behavior is expected for thermal cases with the thermal boundary layer growth at the entrance of a tube as shown in the figure to the right, which corresponds to a case where there may be an unheated length in which the velocity is fully developed before heating starts.  
  
One expects that the thermal boundary layer increases in the thermal entry region before the heat transfer coefficient becomes constant.  One expects that the thermal boundary layer increases in the thermal entry region before the heat transfer coefficient becomes constant.  
  It should be noted that the requirement for a fully developed thermal region is that the dimensionless temperature, <math>\theta </math>,  +  It should be noted that the requirement for a fully developed thermal region is that the dimensionless temperature, <math>\theta </math>, <math>\theta =\frac{T_{w}T}{T_{w}T_{m}}\quad \ \text{or}\quad \ \frac{T_{w}T}{T_{w}T_{c}}</math> does not change with distance along the flow direction. The mean and centerline temperatures are ''T<sub>m</sub>'' and ''T<sub>c</sub>'', respectively. 
  +  
  <math>\theta =\frac{T_{w}T}{T_{w}T_{m}}\quad \ \text{or}\quad \ \frac{T_{w}T}{T_{w}T_{c}}</math>  +  
  +  Similar requirements exist for the fully developed concentration profile where ''θ'' is replaced with <math>\varphi </math>,  
  Similar requirements exist for the fully developed concentration profile where θ is replaced with <math>\varphi </math>,  +  
  <math>\varphi =\frac{c_{w}c}{c_{w}c_{m}}\quad \ \text{or}\quad \ \frac{c_{w}c}{c_{w}c_{c}}</math>  +  <center><math>\varphi =\frac{c_{w}c}{c_{w}c_{m}}\quad \ \text{or}\quad \ \frac{c_{w}c}{c_{w}c_{c}}</math></center> 
where the mean and centerline concentrations (or mass fractions) are cm and cc, respectively.  where the mean and centerline concentrations (or mass fractions) are cm and cc, respectively.  
+  
In the subsequent sections, we use the following definitions to mathematically define the fully developed flow, temperature, and concentration profiles:  In the subsequent sections, we use the following definitions to mathematically define the fully developed flow, temperature, and concentration profiles:  
Fully developed flow profile  Fully developed flow profile  
Line 87:  Line 87:  
 {{EquationRef (6)}}   {{EquationRef (6)}}  
}  }  
  We can now define the local heat and mass transfer coefficients (h and  +  We can now define the local heat and mass transfer coefficients (''h'' and ''h<sub>m</sub>'') based on the mean temperature or concentration. 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
    
Line 102:  Line 102:  
 {{EquationRef (8)}}   {{EquationRef (8)}}  
}  }  
  where D is mass diffusivity.  +  where ''D'' is mass diffusivity. 
+  
Since we define the fully developed temperature profile as when the nondimensional temperature profile <math>(T_{w}T)/(T_{w}T_{m})</math> is invariant in the flow direction (xdirection), we can write the following equation:  Since we define the fully developed temperature profile as when the nondimensional temperature profile <math>(T_{w}T)/(T_{w}T_{m})</math> is invariant in the flow direction (xdirection), we can write the following equation:  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 111:  Line 112:  
 {{EquationRef (9)}}   {{EquationRef (9)}}  
}  }  
+  
The above conclusion, that the local heat transfer coefficient is constant along the flow direction for a fully developed temperature profile, is only valid for constant wall heat flux or constant wall temperature conditions.  The above conclusion, that the local heat transfer coefficient is constant along the flow direction for a fully developed temperature profile, is only valid for constant wall heat flux or constant wall temperature conditions.  
+  
The requirement for the dimensionless temperature to be invariant for a fully developed temperature profile can also be presented as  The requirement for the dimensionless temperature to be invariant for a fully developed temperature profile can also be presented as  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 128:  Line 131:  
 {{EquationRef (11)}}   {{EquationRef (11)}}  
}  }  
  In external flow, the heat and mass transfer coefficients are usually defined by a driving differential (  +  
+  In external flow, the heat and mass transfer coefficients are usually defined by a driving differential (''T<sub>w</sub>''  ''T<sub>∞</sub>'') or (''ω<sub>w</sub>'' – ''ω<sub>∞</sub>'') where ''T<sub>∞</sub>'' and ''ω<sub>∞</sub>'' are the temperature and mass fraction of the fluid in the free stream (far away from the wall). In most cases, ''T<sub>∞</sub>'' and ''ω<sub>∞</sub>'' are known and constant for external flows. However, in internal flow configurations, there is not usually a welldefined temperature or concentration (mass fraction), except at the inlet and/or the boundaries. In internal flow, the temperature and concentrations may change both in the axial direction and perpendicular to the flow direction. Therefore, there are several choices available for the driving differential for temperature and concentration in internal flow.  
+  
The most common choice for defining the driving temperature or concentration is based on mean temperature or concentration (mass fraction or mass density). The mixed mean fluid temperature or concentration is defined at a given local axial location based on the convective thermal energy or mass balance, i.e.,  The most common choice for defining the driving temperature or concentration is based on mean temperature or concentration (mass fraction or mass density). The mixed mean fluid temperature or concentration is defined at a given local axial location based on the convective thermal energy or mass balance, i.e.,  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 144:  Line 149:  
 {{EquationRef (13)}}   {{EquationRef (13)}}  
}  }  
  where <math>\rho _{A,m}</math> is the mean mass density for a given component A, and <math>\rho _{m}</math> is the mean density for the fluid where the mean velocity is defined as  +  where <math>\rho _{A,m}</math> is the mean mass density for a given component ''A'', and <math>\rho _{m}</math> is the mean density for the fluid where the mean velocity is defined as 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
    
Line 174:  Line 179:  
 {{EquationRef (17)}}   {{EquationRef (17)}}  
}  }  
  We will now focus our attention on two special conventional boundary conditions; constant wall heat flux and constant surface temperature. First, consider the constant heat flux or heat rate at the wall, which occurs in many applications such as electronic cooling, electric resistance heating, and radiant heating. From eq. (  +  
+  We will now focus our attention on two special conventional boundary conditions; constant wall heat flux and constant surface temperature. First, consider the constant heat flux or heat rate at the wall, which occurs in many applications such as electronic cooling, electric resistance heating, and radiant heating. From eq. (7), since h and <math>{q}''_{w}</math> are constant, we can conclude ''T<sub>w</sub>'' – ''T<sub>m</sub>'' = constant. Differentiating leads to  
<math>\frac{dT_{w}}{dx}=\frac{dT_{m}}{dx}</math>  <math>\frac{dT_{w}}{dx}=\frac{dT_{m}}{dx}</math>  
  Substituting into eq. (  +  Substituting into eq. (11) gives us 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 187:  Line 193:  
{{EquationRef(18)}}  {{EquationRef(18)}}  
}  }  
  Now consider the case of constant surface or wall temperature, which also occurs in many applications including condensers, evaporators and any heat exchange surface where the heat transfer coefficient is extremely high. Using eq. (  +  
+  Now consider the case of constant surface or wall temperature, which also occurs in many applications including condensers, evaporators and any heat exchange surface where the heat transfer coefficient is extremely high. Using eq. (11) and the fact that ''dT''<sub>w</sub> /''dx'' = 0 for constant surface temperature, we get  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 196:  Line 203:  
{{EquationRef(19)}}  {{EquationRef(19)}}  
}  }  
  
+  [[Image:Fig5.4.pngthumb400 pxalt=Wall and mean temperature variation along the flow in a circular tube for fully developed flow and temperature profile  Wall and mean temperature variation along the flow in a circular tube for fully developed flow and temperature profile.]]  
+  
+  
+  It should be emphasized that eqs. (18) and (19) are only applicable when the temperature profile is fully developed.  
+  The variations of wall and mean temperature for the fully developed temperature profile along the flow for constant heat rate or surface temperature are shown in figure to the right.  
  
  
Finally, to obtain the convective heat and/or mass transfer coefficients, one needs to solve the continuity, mass, momentum, energy and appropriate species equations. In convective heat and mass transfer problems, it is important to obtain information about the flow by solving the continuity and momentum equations, in addition to the energy and species equations. These conservation equations are mostly decoupled, except for circumstances such as a variable property, or coupled governing equations or boundary conditions due to physical circumstances (which happens in applications such as natural convection, absorption, sublimation, evaporation and condensation problems).  Finally, to obtain the convective heat and/or mass transfer coefficients, one needs to solve the continuity, mass, momentum, energy and appropriate species equations. In convective heat and mass transfer problems, it is important to obtain information about the flow by solving the continuity and momentum equations, in addition to the energy and species equations. These conservation equations are mostly decoupled, except for circumstances such as a variable property, or coupled governing equations or boundary conditions due to physical circumstances (which happens in applications such as natural convection, absorption, sublimation, evaporation and condensation problems).  
  It is obviously more accurate to solve the complete transport conservation equations (elliptic form) for internal flow without making boundary layer assumptions (parabolic form) as discussed in  +  It is obviously more accurate to solve the complete transport conservation equations (elliptic form) for internal flow without making boundary layer assumptions (parabolic form) as discussed in [[external forced convection]]. However, in most cases it is not practical due to complexity of the geometry and/or solution techniquesm, as well as the requirement of additional boundary conditions in both analytical or numerical methods. 
  For the case of twodimensional fully developed steady laminar flow with constant properties, the momentum equation in a circular tube, including boundary conditions,  +  
+  For the case of twodimensional fully developed steady laminar flow with constant properties, the momentum equation in a circular tube, including boundary conditions, is:  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 240:  Line 250:  
{{EquationRef(22)}}  {{EquationRef(22)}}  
}  }  
  Using the definition of mean velocity,  +  Using the definition of mean velocity, <math>u_{m}</math>, for constant properties and the above equation, we obtain 
  <math>u_{m}</math>  +  
  , for constant properties and the above equation, we obtain  +  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 251:  Line 259:  
{{EquationRef(23)}}  {{EquationRef(23)}}  
}  }  
  Equation (  +  
+  Equation (22) in terms of mean velocity is  
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 269:  Line 278:  
{{EquationRef(25)}}  {{EquationRef(25)}}  
}  }  
  The above result can be presented in terms of the friction coefficient,  +  The above result can be presented in terms of the friction coefficient, ''c<sub>f</sub>''. 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 287:  Line 296:  
{{EquationRef(27)}}  {{EquationRef(27)}}  
}  }  
  It follows from eq. (  +  It follows from eq. (25) that 
{ class="wikitable" border="0"  { class="wikitable" border="0"  
Line 296:  Line 305:  
{{EquationRef(28)}}  {{EquationRef(28)}}  
}  }  
  Let’s make an analysis of the energy equation to get a feeling about the importance of various terms, since determination of the temperature field in the fluid is required for the heat transfer coefficients. To simplify the analysis, lets consider a twodimensional cylindrical geometry with the following assumptions:  +  
  1. Steady laminar flow  +  Let’s make an analysis of the energy equation to get a feeling about the importance of various terms, since determination of the temperature field in the fluid is required for the heat transfer coefficients. To simplify the analysis, lets consider a twodimensional cylindrical geometry with the following assumptions:<br> 
  2. Constant properties  +  1. Steady laminar flow<br> 
  3. Fully developed flow  +  2. Constant properties<br> 
  4. Newtonian incompressible fluid  +  3. Fully developed flow<br> 
+  4. Newtonian incompressible fluid<br>  
+  
The energy equation under the above assumptions is  The energy equation under the above assumptions is  
Line 321:  Line 332:  
 {{EquationRef(30)}}   {{EquationRef(30)}}  
}  }  
  where  +  where ''T<sub>r</sub>'' is a reference temperature and E and Pe are Eckert and Peclet numbers, respectively. 
+  
The resulting dimensionless energy equation is  The resulting dimensionless energy equation is  
Line 331:  Line 343:  
{{EquationRef(31)}}  {{EquationRef(31)}}  
}  }  
  The second term on the right hand side of the above equation is due to axial heat conduction, and the last term is due to the viscous dissipation effect. If E Pr is small, viscous dissipation can be neglected. This is true for flow with a low velocity and low Prandtl number. The second term on the right hand side (axial heat conduction) is neglected when the Peclet number, Pe, is greater than 100. Axial heat conduction should be accounted for when the Peclet number is small, in the case of liquid metals.  +  The second term on the right hand side of the above equation is due to axial heat conduction, and the last term is due to the viscous dissipation effect. If ''E'' Pr is small, viscous dissipation can be neglected. This is true for flow with a low velocity and low Prandtl number. The second term on the right hand side (axial heat conduction) is neglected when the Peclet number, Pe, is greater than 100. Axial heat conduction should be accounted for when the Peclet number is small, in the case of liquid metals. 
  +  
  +  ==References==  
+  {{Reflist}} 
Revision as of 05:02, 22 July 2010
It is important to clarify some basic definitions, terminologies and criteria that are often used in internal convective heat and mass transfer. These include:
1. Mean velocity, temperature, and concentration
2. Fully developed flow, temperature, and concentration profiles
3. Hydrodynamic, thermal, and concentration entrance lengths
The figure to the right shows the development of a velocity profile inside a duct or tube with uniform inlet velocity for laminar flow of an incompressible Newtonian fluid. The velocity profile at some distance away from the tube’s inlet no longer changes along the flow direction, where it is referred to as the fully developed flow condition.
The fully developed condition is often met at some distance away from the inlet. However, there are also applications in which fully developed flow is never reached. Momentum, thermal, and concentration boundary layers form on the inside surface of the tube. The thickness of the layers increases in a similar manner as boundary layer flow over a flat plate (which was presented in detail in external forced convection).
Part (a) of the figure to the right shows how the momentum boundary layer builds up in a pipe along the flow direction. At some distance away from the inlet, the boundary layer fills the flow area. The flow downstream from this point is referred to as fully developed flow since the velocity slope does not change after this point. The distance downstream from the inlet to where the flow becomes fully developed is called the hydrodynamic entrance length. If the flow is laminar (Re < 2300 for flow inside circular tubes), the fully developed velocity is a parabolic shape. It should be noted that the fluid velocity outside the boundary layer increases with x, which is required to satisfy the conservation of mass (or continuity) equation.
The center line velocity finally reaches a value two times the inlet velocity, u_{in}, for fully developed, steady, incompressible, laminar flow inside tubes. It should be noted that the hydrodynamic entrance length for fully developed flow does not start from the point where the friction coefficient,

does not change along the flow. The friction coefficient variation for laminar flow inside a circular tube with uniform inlet velocity is shown in part (b) of the above figure.
The friction coefficient is highest at the entrance and then decreases smoothly to a constant value, corresponding to fully developed flow. Two factors cause the friction coefficient to be higher in the entrance region of tubes than in the fully developed region. The first factor is the larger velocity gradient at the entrance on the wall. The gradient decreases along the pipe and becomes constant before the velocity becomes fully developed. The second factor is the velocity outside the boundary layer, which must increase to satisfy the conservation of mass or continuity equation. Accelerating velocity in the core produces an additional drag force when its effect is considered in the friction coefficient.
The turbulent velocity profile and friction coefficient variation for a circular pipe are shown in figure to the right. Even for a very high inlet velocity, the boundary layer will be laminar over a part of the entrance. This transition from laminar to turbulent is clearly shown by the sudden increase in momentum boundary layer thickness as shown in part (a) of the figure. The friction coefficient variation for turbulent flow in a pipe entrance is shown in part (b).
The hydrodynamic entry length required for fully developed flow should be obtained by a complete solution of the flow and thermal field in the entrance region. A rule of thumb to judge whether or not the flow is fully developed for circular pipes is


where L_{H} is the hydrodynamic length and the Reynolds number is defined by
A similar behavior is expected for thermal cases with the thermal boundary layer growth at the entrance of a tube as shown in the figure to the right, which corresponds to a case where there may be an unheated length in which the velocity is fully developed before heating starts.
One expects that the thermal boundary layer increases in the thermal entry region before the heat transfer coefficient becomes constant. It should be noted that the requirement for a fully developed thermal region is that the dimensionless temperature, θ, does not change with distance along the flow direction. The mean and centerline temperatures are T_{m} and T_{c}, respectively.
Similar requirements exist for the fully developed concentration profile where θ is replaced with ,
where the mean and centerline concentrations (or mass fractions) are cm and cc, respectively.
In the subsequent sections, we use the following definitions to mathematically define the fully developed flow, temperature, and concentration profiles: Fully developed flow profile

Fully developed temperature profile

Fully developed concentration profile

We can now define the local heat and mass transfer coefficients (h and h_{m}) based on the mean temperature or concentration.


where D is mass diffusivity.
Since we define the fully developed temperature profile as when the nondimensional temperature profile (T_{w} − T) / (T_{w} − T_{m}) is invariant in the flow direction (xdirection), we can write the following equation:

The above conclusion, that the local heat transfer coefficient is constant along the flow direction for a fully developed temperature profile, is only valid for constant wall heat flux or constant wall temperature conditions.
The requirement for the dimensionless temperature to be invariant for a fully developed temperature profile can also be presented as

Differentiating the above equation yields

In external flow, the heat and mass transfer coefficients are usually defined by a driving differential (T_{w}  T_{∞}) or (ω_{w} – ω_{∞}) where T_{∞} and ω_{∞} are the temperature and mass fraction of the fluid in the free stream (far away from the wall). In most cases, T_{∞} and ω_{∞} are known and constant for external flows. However, in internal flow configurations, there is not usually a welldefined temperature or concentration (mass fraction), except at the inlet and/or the boundaries. In internal flow, the temperature and concentrations may change both in the axial direction and perpendicular to the flow direction. Therefore, there are several choices available for the driving differential for temperature and concentration in internal flow.
The most common choice for defining the driving temperature or concentration is based on mean temperature or concentration (mass fraction or mass density). The mixed mean fluid temperature or concentration is defined at a given local axial location based on the convective thermal energy or mass balance, i.e.,


where ρ_{A,m} is the mean mass density for a given component A, and ρ_{m} is the mean density for the fluid where the mean velocity is defined as

Assuming constant properties for mean velocity, temperature and mass density in the above equation, we obtain



We will now focus our attention on two special conventional boundary conditions; constant wall heat flux and constant surface temperature. First, consider the constant heat flux or heat rate at the wall, which occurs in many applications such as electronic cooling, electric resistance heating, and radiant heating. From eq. (7), since h and q''_{w} are constant, we can conclude T_{w} – T_{m} = constant. Differentiating leads to
Substituting into eq. (11) gives us

Now consider the case of constant surface or wall temperature, which also occurs in many applications including condensers, evaporators and any heat exchange surface where the heat transfer coefficient is extremely high. Using eq. (11) and the fact that dT_{w} /dx = 0 for constant surface temperature, we get

It should be emphasized that eqs. (18) and (19) are only applicable when the temperature profile is fully developed.
The variations of wall and mean temperature for the fully developed temperature profile along the flow for constant heat rate or surface temperature are shown in figure to the right.
Finally, to obtain the convective heat and/or mass transfer coefficients, one needs to solve the continuity, mass, momentum, energy and appropriate species equations. In convective heat and mass transfer problems, it is important to obtain information about the flow by solving the continuity and momentum equations, in addition to the energy and species equations. These conservation equations are mostly decoupled, except for circumstances such as a variable property, or coupled governing equations or boundary conditions due to physical circumstances (which happens in applications such as natural convection, absorption, sublimation, evaporation and condensation problems).
It is obviously more accurate to solve the complete transport conservation equations (elliptic form) for internal flow without making boundary layer assumptions (parabolic form) as discussed in external forced convection. However, in most cases it is not practical due to complexity of the geometry and/or solution techniquesm, as well as the requirement of additional boundary conditions in both analytical or numerical methods.
For the case of twodimensional fully developed steady laminar flow with constant properties, the momentum equation in a circular tube, including boundary conditions, is:


Integrating the above equation twice and using the boundary conditions yield a parabolic velocity profile:

Using the definition of mean velocity, u_{m}, for constant properties and the above equation, we obtain

Equation (22) in terms of mean velocity is

The shear stress at the wall can be calculated from the velocity gradient at the wall.

The above result can be presented in terms of the friction coefficient, c_{f}.

In addition to the above friction coefficient, the following friction factor is also widely used:

It follows from eq. (25) that

Let’s make an analysis of the energy equation to get a feeling about the importance of various terms, since determination of the temperature field in the fluid is required for the heat transfer coefficients. To simplify the analysis, lets consider a twodimensional cylindrical geometry with the following assumptions:
1. Steady laminar flow
2. Constant properties
3. Fully developed flow
4. Newtonian incompressible fluid
The energy equation under the above assumptions is

The above equation is nondimensionalized using the following variables to show the effect of axial conduction and viscous dissipation:

where T_{r} is a reference temperature and E and Pe are Eckert and Peclet numbers, respectively.
The resulting dimensionless energy equation is

The second term on the right hand side of the above equation is due to axial heat conduction, and the last term is due to the viscous dissipation effect. If E Pr is small, viscous dissipation can be neglected. This is true for flow with a low velocity and low Prandtl number. The second term on the right hand side (axial heat conduction) is neglected when the Peclet number, Pe, is greater than 100. Axial heat conduction should be accounted for when the Peclet number is small, in the case of liquid metals.