In heat transfer, the Nusselt number is a nondimensionalization of the convective heat transfer coefficient. The Nusselt number is named in honor of Wilhelm Nusselt.

The heat transfer coefficient, <math>h</math>, at surface cooled by a fluid is the ratio of the heat flux <math>q</math> from the surface to the difference <math>\Delta T</math>

between the surface temperature and the temperature of the fluid far from the surface: <math>h = q/\Delta T</math>. The Nusselt is formed by nondimensionalizing <math>h</math> with a length <math>L</math> and the thermal conductivity of the fluid <math>k</math>:

:<math>\mathrm{Nu}_L = \frac{qL}{k\Delta T} = \frac{hL}{k}</math>

The Nusselt number may be defined locally, at a single position on a surface, or as an average value that represents the heat flow from the entire surface.

A Nusselt number between one and 10 is characteristic of slug flow or laminar flow. A larger Nusselt number corresponds to more active convection, with turbulent flow typically in the 100–1000 range.

  • Selection of the characteristic length should be in the direction of growth (or thickness) of the boundary layer; some examples of characteristic length are: the outer diameter of a cylinder in (external) cross flow (perpendicular to the cylinder axis), the height of a vertical plate undergoing natural convection, or the diameter of a sphere. For complex shapes, the length may be defined as the volume of the fluid body divided by the surface area.
  • The thermal conductivity of the fluid is typically (but not always) evaluated at the film temperature, which for engineering purposes may be calculated as the mean-average of the bulk fluid temperature and wall surface temperature.

In contrast to the definition given above, known as average Nusselt number, the local Nusselt number is scaled using the distance from the leading edge of a boundary layer to the point of interest on the surface.

:<math>\overline{\mathrm{Nu=\frac{\frac{1}{L} \int_0^L h_x\ dx\ L}{k}=\frac{\overline{h} L}{k}</math>

For a constant flux surface, the temperature difference is not constant, and so the average process is more complicated.

Context

An understanding of convection boundary layers is necessary to understand convective heat transfer between a surface and a fluid flowing past it. A thermal boundary layer develops if the fluid free stream temperature and the surface temperatures differ. A temperature profile exists due to the energy exchange resulting from this temperature difference.

400px|thumb|Momentum boundary layer (left) and thermal boundary layer (right)

The heat transfer rate can be written using Newton's law of cooling as

:<math>Q_y=hA\left( T_s-T_\infty \right)</math>,

where h is the heat transfer coefficient and A is the heat transfer surface area. Because heat transfer at the surface is by conduction, the same quantity can be expressed in terms of the thermal conductivity k:

:<math>Q_y=-kA\frac{\partial }{\partial y}\left. \left( T-T_s \right) \right|_{y=0}</math>.

These two terms are equal; thus

:<math>-kA\frac{\partial }{\partial y}\left. \left( T-T_s \right) \right|_{y=0}=hA\left( T_s-T_\infty \right)</math>.

Rearranging,

:<math>\frac{h}{k}=\frac{\left. \frac{\partial \left( T_s-T \right)}{\partial y} \right|_{y=0{\left( T_s-T_\infty \right)}</math>.

Multiplying by a representative length L gives a dimensionless expression:

:<math>\frac{hL}{k}=\frac{\left. \frac{\partial \left( T_s-T \right)}{\partial y} \right|_{y=0{\frac{\left( T_s-T_\infty \right)}{L</math>.

The right-hand side is now the ratio of the temperature gradient at the surface to the reference temperature gradient, while the left-hand side is similar to the Biot number. This becomes the ratio of conductive thermal resistance to the convective thermal resistance of the fluid, otherwise known as the Nusselt number, Nu.

:<math>\mathrm{Nu} = \frac{h}{k/L} = \frac{hL}{k}</math>.

Empirical correlations

Typically, for free convection, the average Nusselt number is expressed as a function of the Rayleigh number and the Prandtl number, written as:

:<math>\mathrm{Nu} = f(\mathrm{Ra}, \mathrm{Pr})</math>

Otherwise, for forced convection, the Nusselt number is generally a function of the Reynolds number and the Prandtl number, or

:<math>\mathrm{Nu} = f(\mathrm{Re}, \mathrm{Pr})</math>

Empirical correlations for a wide variety of geometries are available that express the Nusselt number in the aforementioned forms.

Free convection

Free convection at a vertical wall

Cited as coming from Churchill and Chu:

:<math>\overline{\mathrm{Nu_L \ = \left({0.825 + \frac{0.387 \mathrm{Ra}_L^{1/6{\left(1 + (0.492/\mathrm{Pr})^{9/16} \right)^{8/27} \right)^2 \, \quad \mathrm{Ra}_L < 10^{12}</math>

For laminar flows, the following correlation is slightly more accurate. It is observed that a transition from a laminar to a turbulent boundary occurs when Ra<sub>L</sub> exceeds around 10<sup>9</sup>.

:<math>\overline{\mathrm{Nu_L \ = \left(0.68 + \frac{0.67 \mathrm{Ra}_L^{1/4{\left(1 + (0.492/\mathrm{Pr})^{9/16}\right)^{4/9\right) \, \quad \mathrm10^{-1} < \mathrm{Ra}_L < 10^9 </math>

Free convection from horizontal plates

If the characteristic length is defined

:<math>L \ = \frac{A_s}{P}</math>

where <math>\mathrm{A}_s</math> is the surface area of the plate and <math>P</math> is its perimeter.

Then for the top surface of a hot object in a colder environment or bottom surface of a cold object in a hotter environment

:<math>\overline{\mathrm{Nu_L \ = 0.54\, \mathrm{Ra}_L^{1/4} \, \quad 10^4 \le \mathrm{Ra}_L \le 10^7</math>

:<math>\overline{\mathrm{Nu_L \ = 0.15\, \mathrm{Ra}_L^{1/3} \, \quad 10^7 \le \mathrm{Ra}_L \le 10^{11}</math>

And for the bottom surface of a hot object in a colder environment or top surface of a cold object in a hotter environment

:<math>\overline{\mathrm{Nu_L \ = 0.27\, \mathrm{Ra}_L^{1/4} \, \quad 10^5 \le \mathrm{Ra}_L \le 10^{10}</math>

Free convection from enclosure heated from below

Cited as coming from Bejan:

:<math>\overline{\mathrm{Nu_L \ = 0.069\, \mathrm{Ra}_L^{1/3}Pr^{0.074} \, \quad 3 * 10^5 \le \mathrm{Ra}_L \le 7 * 10^{9}</math>

This equation <i>"holds when the

horizontal layer is sufficiently wide so that the effect of the short vertical sides

is minimal."</i>

It was empirically determined by Globe and Dropkin in 1959: <i>"Tests were made in cylindrical containers having copper tops and bottoms and insulating walls."</i> The containers used were around 5" in diameter and 2" high.

Flat plate in laminar flow

The local Nusselt number for laminar flow over a flat plate, at a distance <math>x</math> downstream from the edge of the plate, is given by

:<math>\mathrm{Nu}_x\ = 0.332\, \mathrm{Re}_x^{1/2}\, \mathrm{Pr}^{1/3}, (\mathrm{Pr} > 0.6) </math>

The average Nusselt number for laminar flow over a flat plate, from the edge of the plate to a downstream distance <math>x</math>, is given by

:<math>\overline{\mathrm{Nu_x \ = {2} \cdot 0.332\, \mathrm{Re}_x^{1/2}\, \mathrm{Pr}^{1/3}\ = 0.664\, \mathrm{Re}_x^{1/2}\, \mathrm{Pr}^{1/3}, (\mathrm{Pr} > 0.6) </math>

Sphere in convective flow

In some applications, such as the evaporation of spherical liquid droplets in air, the following correlation is used:

:<math>\mathrm{Nu}_D \ = {2} + 0.4\, \mathrm{Re}_D^{1/2}\, \mathrm{Pr}^{1/3}\, </math>

Forced convection in turbulent pipe flow

Gnielinski correlation

Gnielinski's correlation for flow in tubes with <math>3000 \le \mathrm{Re}_D \le 5 \times 10^{6}</math>:

:<math>\mathrm{Nu}_D = \frac{ \left( f/8 \right) \left( \mathrm{Re}_D - 1000 \right) \mathrm{Pr} } {1 + 12.7(f/8)^{1/2} \left( \mathrm{Pr}^{2/3} - 1 \right)}</math>

where f is the Darcy friction factor that can either be obtained from the Moody chart or for smooth tubes from correlation developed by Petukhov:

:<math>f= \left( 0.79 \ln \left(\mathrm{Re}_D \right)-1.64 \right)^{-2}</math>

The Gnielinski Correlation is valid for:

:<math>0.5 \le \mathrm{Pr} \le 2000</math>

:<math>3000 \le \mathrm{Re}_D \le 5 \times 10^{6}</math>

Dittus–Boelter equation

The Dittus–Boelter equation (for turbulent flow) as introduced by W.H. McAdams is an explicit function for calculating the Nusselt number. It is easy to solve but is less accurate when there is a large temperature difference across the fluid. It is tailored to smooth tubes, so use for rough tubes (most commercial applications) is cautioned. The Dittus–Boelter equation is:

:<math>\mathrm{Nu}_D = 0.023\, \mathrm{Re}_D^{4/5}\, \mathrm{Pr}^{n}</math>

where:

:<math>D</math> is the inside diameter of the circular duct

:<math>\mathrm{Pr}</math> is the Prandtl number

:<math>n = 0.4</math> for the fluid being heated, and <math>n = 0.3</math> for the fluid being cooled.

The Dittus–Boelter equation is valid for

:<math>0.6 \le \mathrm{Pr} \le 160</math>

:<math>\mathrm{Re}_D \gtrsim 10\,000</math>

:<math>\frac{L}{D} \gtrsim 10</math>

The Dittus–Boelter equation is a good approximation where temperature differences between bulk fluid and heat transfer surface are minimal, avoiding equation complexity and iterative solving. Taking water with a bulk fluid average temperature of , viscosity and a heat transfer surface temperature of (viscosity , a viscosity correction factor for <math>({\mu} / {\mu_s})</math> can be obtained as 1.45. This increases to 3.57 with a heat transfer surface temperature of (viscosity ), making a significant difference to the Nusselt number and the heat transfer coefficient.

Sieder–Tate correlation

The Sieder–Tate correlation for turbulent flow is an implicit function, as it analyzes the system as a nonlinear boundary value problem. The Sieder–Tate result can be more accurate as it takes into account the change in viscosity (<math>\mu</math> and <math>\mu_s</math>) due to temperature change between the bulk fluid average temperature and the heat transfer surface temperature, respectively. The Sieder–Tate correlation is normally solved by an iterative process, as the viscosity factor will change as the Nusselt number changes.

:<math>\mathrm{Nu}_D = 0.027\,\mathrm{Re}_D^{4/5}\, \mathrm{Pr}^{1/3}\left(\frac{\mu}{\mu_s}\right)^{0.14}</math>

where:

:<math>\mu</math> is the fluid viscosity at the bulk fluid temperature

:<math>\mu_s</math> is the fluid viscosity at the heat-transfer boundary surface temperature

The Sieder–Tate correlation is valid for

:<math>0.7 \le \mathrm{Pr} \le 16\,700</math>

:<math>\mathrm{Re}_D \ge 10\,000</math>

:<math>\frac{L}{D} \gtrsim 10</math>

Forced convection in fully developed laminar pipe flow

For fully developed internal laminar flow, the Nusselt numbers tend towards a constant value for long pipes.

For internal flow:

:<math>\mathrm{Nu} = \frac{h D_h}{k_f}</math>

where:

:D<sub>h</sub> = Hydraulic diameter

:k<sub>f</sub> = thermal conductivity of the fluid

:h = convective heat transfer coefficient

Convection with uniform temperature for circular tubes

From Incropera & DeWitt,

:<math>\mathrm{Nu}_D = 3.66</math>

OEIS sequence gives this value as <math>\mathrm{Nu}_D = 3.6567934577632923619...</math>.

Convection with uniform heat flux for circular tubes

For the case of constant surface heat flux,

:<math>\mathrm{Nu}_D = 4.36</math>

See also

  • Sherwood number (mass transfer Nusselt number)
  • Churchill–Bernstein equation
  • Biot number
  • Reynolds number
  • Convective heat transfer
  • Heat transfer coefficient
  • Thermal conductivity

References

  • Simple derivation of the Nusselt number from Newton's law of cooling (Accessed 23 September 2009)