The Nusselt number is a nondimensionalization of the convective heat transfer coefficient. Like the heat transfer coefficient, 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.
The Nusselt number is named in honor of Wilhelm Nusselt, who first identified this dimensionless group in 1915.
Analytical results and empirical correlations allow the Nusselt number to be estimated in many situations. For forced convection, these expressions depend on the Reynolds number and the Prandtl number. For natural convection (or "free" convection) the predictions use the Grashof number or Rayleigh number along with the Prandtl number.
The mass transfer analog of the Nusselt number is the Sherwood number.
The group now known as the Nusselt number was first identified by Wilhelm Nusselt in a 1915 paper. In that paper, he also identified the groups that became the Grashof number, the Prandtl number, and the Reynolds number. Nusselt used these groups effectively to correlate data for natural convection and, in a 1917 paper, forced convection.
In 1931, the Association of German Engineers formally named the dimensionless group the Nusselt number. The decision was announced in paper by Max Jakob. The term was quickly adopted internationally.
The heat transfer coefficient, h {\displaystyle h} , at a surface cooled by a fluid is the ratio of the heat flux q {\displaystyle q} [W/m2] from the surface to the temperature difference Δ T {\displaystyle \Delta T} [K] between the surface and the fluid far from the surface: h = q / Δ T {\displaystyle h=q/\Delta T} [W/(m2·K)]. The Nusselt number, Nu {\displaystyle {\textrm {Nu}}} , is formed by nondimensionalizing h {\displaystyle h} with a characteristic length L {\displaystyle L} [m] and the thermal conductivity of the fluid k {\displaystyle k} [W/(m·K)]:
The Nusselt number is a dimensionless group. The subscript, L, indicates which characteristic length is used.
The characteristic length should be related to the thickness of the thermal boundary layer. In a tube, the thermal boundary layer will be limited by the tube diameter, D {\displaystyle D} , so that the diameter represents the scale of the boundary layer. On a flat plate, the length of the plate, L {\displaystyle L} , determines the average thickness of the boundary layer on the plate. On the other hand, the local thickness of the flat plate boundary layer, at a position x {\displaystyle x} , is determined by the value of x {\displaystyle x} , so that x {\displaystyle x} is the characteristic length for the local heat transfer coefficient on the plate.
Some additional examples of characteristic length are: the outer diameter of a cylinder in (external) crossflow (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.
At a location x {\displaystyle x} on a surface, the local heat transfer coefficient, h x {\displaystyle h_{x}} , relates the local heat flux, q {\displaystyle q} , to the local temperature difference Δ T x {\displaystyle \Delta T_{x}} :
The local heat transfer coefficient varies along the surface of a body. It becomes lower as the thermal boundary layer becomes thicker. It increases suddenly where the boundary layer transitions from laminar to turbulent flow.
The local Nusselt number is a nondimensionalization of the local heat transfer coefficient, h x {\displaystyle h_{x}} , at a position x {\displaystyle x} on a surface:
The average heat transfer coefficient relates the average temperature difference to the average heat flux over the body. For an isothermal body the temperature difference is constant, and if the body has length L {\displaystyle L} the average heat flux is
The average Nusselt number is then defined as
For a uniform flux surface, the temperature difference is not constant, and so the temperature difference must be averaged to find the average heat transfer coefficient:
The average temperature difference is another integral:
The average Nusselt number is defined just as for an isothermal surface:
The term "average Nusselt number" can cause confusion because the Nusselt number is not what has been averaged. For an isothermal surface, the heat flux has been averaged, while for a uniform flux surface, the temperature difference has been averaged.
A closely related dimensionless group is the Stanton number used in some studies of forced convection:
where ρ , c p , and u {\displaystyle \rho ,c_{p},{\text{and }}u} are the fluid density, specific heat capacity, and free stream speed.
A similar nondimensional group is the Biot number, which compares the thermal resistance of heat conduction in the body to the thermal resistance of convection outside the body. The Biot number uses the thermal conductivity of the solid body rather than the fluid. The Biot number should not be confused with the Nusselt number.
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 is created by the energy exchange resulting from this temperature difference.
The local heat transfer rate can be written using Newton's law of cooling as
At the surface of the body, the fluid velocity is zero as a result of the no-slip condition, so heat transfer in the fluid at the surface is by conduction alone:
These two terms are equal; thus
For a laminar boundary layer (see figure), the temperature gradient at the surface has a temperature change T s − T 0 {\displaystyle T_{s}-T_{0}} over a distance proportional to the thermal boundary layer thickness, δ t {\displaystyle \delta _{t}} , so a crude approximation is:
Hence, the laminar Nusselt number is greater when the thermal boundary layer is thinner:
For forced convection, the Nusselt number is generally a function of the Reynolds number and the Prandtl number, or
For free, or natural, convection, the average Nusselt number is usually expressed as a function of the Rayleigh number or Grashof number and the Prandtl number:
or
since Ra = GrPr.
Analytical results and empirical correlations that express the Nusselt number in the aforementioned forms are available for a wide variety of geometries.
When calculating the Nusselt number, the thermal conductivity, the Prandtl number, and the properties in the Reynolds or Grashof numbers are usually evaluated at the film temperature. The film temperature is the average of the wall temperature and free stream or bulk temperature. When the bulk temperature changes significantly along the length of the tube, the average bulk temperature can be used.[citation needed]
Churchill and Chu correlated a wide range of data with the following expression, which includes both laminar and turbulent flow:
The transition from a laminar to a turbulent boundary occurs at R a L ≈ 10 9 {\displaystyle \mathrm {Ra} _{L}\approx 10^{9}} . For laminar flows, Churchill and Chu recommended the following, slightly more accurate correlation:
If the characteristic length is defined
where A s {\displaystyle \mathrm {A} _{s}} is the surface area of the plate and P {\displaystyle P} 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
And for the bottom surface of a hot object in a colder environment or top surface of a cold object in a hotter environment
In 1959, Globe and Dropkin reported the following correlation for wide fluid layers between a lower heated and upper cooled plate:
This equation "holds when the horizontal layer is sufficiently wide so that the effect of the short vertical sides is minimal."
The solution of the energy and momentum equations for laminar flow over an isothermal flat plate leads to the following equation for the local Nusselt number at a distance x {\displaystyle x} downstream from the leading edge of the plate.
The average Nusselt number for laminar flow over an isothermal flat plate, from the edge of the plate to a downstream distance L {\displaystyle L} , is given by
For a sphere in forced flow, such as an evaporating droplet, Faeth suggests:
The result limits to N u D = 2 {\displaystyle \mathrm {Nu} _{D}\ ={2}} for small Reynolds numbers, which corresponds to heat conduction into a sphere in an infinite medium.
Gnielinski's correlation (1975) for flow in smooth tubes is:
where f {\displaystyle f} is the Darcy friction factor that can either be obtained from the Moody chart or from the correlation of Filonenko:
The Gnielinski correlation is valid for:
This correlation is much more accurate than 1930s power-law correlations such as the Dittus–Boelter equation.
Gnielinski also developed power-law correlations for limited ranges of Prandtl number, which agree very closely with his full-range correlation. These equations are simpler to use, while remaining accurate:
These equations apply for smooth tubes with 3000 ≤ R e D ≤ 5 × 10 6 {\displaystyle 3000\leq \mathrm {Re} _{D}\leq 5\times 10^{6}} .
The Dittus–Boelter equation (for turbulent flow in smooth-walled tubes), as introduced by W.H. McAdams, is one of several power-law correlations that were published in the 1930s. These equations are based on smaller data sets than the Gnielinski correlations, represent smaller ranges of Reynolds and Prandtl number, and are generally less accurate. The Dittus–Boelter equation is:
where:
The Dittus–Boelter equation is valid for
The Sieder–Tate correlation for turbulent pipe flow (1936) takes account of the change in viscosity ( μ b {\displaystyle \mu _{b}} and μ w {\displaystyle \mu _{w}} ) due to temperature change between the average bulk temperature of the fluid and wall temperature. The Sieder–Tate correlation may need to be solved iteratively if the bulk or wall temperature is unknown.
where:
The Sieder–Tate correlation is valid for
Taking water with a bulk fluid average temperature of 20 °C (68 °F), viscosity 10.07×10−4 Pa⋅s and a heat transfer surface temperature of 40 °C (104 °F) (viscosity 6.96×10−4 Pa⋅s, a viscosity correction factor for ( μ / μ s ) {\displaystyle ({\mu }/{\mu _{s}})} can be obtained as 1.45. This increases to 3.57 with a heat transfer surface temperature of 100 °C (212 °F) (viscosity 2.82×10−4 Pa⋅s), making a significant difference to the Nusselt number and the heat transfer coefficient.
Solution of the governing equations for fully developed, laminar pipe flow leads to very simple results. The Nusselt numbers tend towards a constant value for long tubes:
where:
From Incropera & DeWitt,
OEIS sequence A282581 gives this value as N u D = 3.6567934577632923619... {\displaystyle \mathrm {Nu} _{D}=3.6567934577632923619...} .
For the case of constant surface heat flux,
Churchill–Bernstein equation
Convective heat transfer
