VentiFlux logoVentiFluxLearning Center
Heat Transfer Foundations

Forced and Mixed Convection in Ducts

Internal-flow heat transfer and pressure drop for real HVAC ducts: entrance effects, Dittus-Boelter and Gnielinski, rectangular ducts, mixed convection, and a full duct-sizing worked example.

Free Intermediate 5 hours 7 sessions

Session 1

Internal forced convection basics

Objectives

  • Explain how velocity and thermal boundary layers develop from a duct entrance
  • Estimate hydrodynamic and thermal entrance lengths for laminar and turbulent flow
  • Use the Graetz number to judge whether entrance effects matter in a given duct
  • Select laminar mean Nusselt number correlations that account for entrance effects

Internal forced convection describes heat transfer between a duct wall and a fluid that is driven through the duct by a fan, blower, or pump. Unlike external flow over a plate or cylinder, the fluid is confined: boundary layers growing from the walls eventually meet, and from that point on the flow is said to be fully developed. The convective heat transfer coefficient is not a property of the fluid alone. It depends on the flow regime, the fluid properties, and crucially on how far the flow has developed along the duct. A practicing engineer who ignores development effects will misjudge heat transfer in short ducts, in heat exchanger passages, and in the many short straight runs between fittings that make up a real duct system.

Consider air entering a straight duct with a nearly uniform velocity profile. At the wall the no-slip condition forces the fluid velocity to zero, so a velocity boundary layer begins to grow on each wall. Because the same mass flow must pass through every cross section, the fluid in the untouched core accelerates as the boundary layers thicken. This developing region ends when the boundary layers merge at the duct centerline. Beyond that point the velocity profile no longer changes with axial position, and the flow is hydrodynamically fully developed. In laminar flow that profile is the familiar parabola; in turbulent flow it is much flatter, with steep gradients confined to a thin layer near the wall.

The distance required for the profile to settle is the hydrodynamic entrance length. For laminar flow in a circular tube it grows in direct proportion to the Reynolds number:

Lfd,h = 0.05 * ReD * D

Lfd,h
Hydrodynamic entrance length, the axial distance from the duct inlet to the point where the velocity profile becomes fully developed (m)
ReD
Reynolds number based on duct diameter, ReD = V * D / nu, dimensionless
D
Inside diameter of the circular duct (m)

This relationship carries a warning that surprises many newcomers. At a Reynolds number of 2000 in a 200 mm duct, the laminar entrance length is 20 m. Laminar duct flow is almost never fully developed in practical equipment, which is one reason laminar correlations must handle entrance effects explicitly. Turbulent flow behaves very differently: intense cross-stream mixing completes the profile development quickly, and the hydrodynamic entrance length is nearly independent of Reynolds number, on the order of ten diameters:

Lfd,h ≈ 10 * D

Lfd,h
Hydrodynamic entrance length for turbulent internal flow (m)
D
Inside diameter of the duct (m)

Since HVAC duct flow is turbulent in the great majority of cases, a straight duct run becomes hydrodynamically developed after roughly ten diameters, and the velocity profile can be treated as settled for pressure drop purposes beyond that point.

Thermal development follows the same logic but on its own length scale. When the wall temperature differs from the entering fluid temperature, a thermal boundary layer grows alongside the velocity boundary layer. Whether it develops faster or slower depends on the Prandtl number, which compares the diffusion of momentum to the diffusion of heat. For laminar flow the thermal entrance length is well approximated by scaling the hydrodynamic entrance length with the Prandtl number:

Lfd,t ≈ Lfd,h * Pr

Lfd,t
Thermal entrance length, the axial distance to the point where the dimensionless temperature profile becomes fully developed (m)
Lfd,h
Hydrodynamic entrance length (m)
Pr
Prandtl number of the fluid, dimensionless

For gases with Pr near 0.7, the thermal entrance is somewhat shorter than the hydrodynamic entrance. For viscous oils with Prandtl numbers in the hundreds or thousands, the thermal entrance can be enormously longer, and the flow may be hydrodynamically developed while still thermally developing over the entire length of the equipment. In turbulent flow the vigorous mixing collapses both entrance lengths to roughly ten diameters, so the distinction matters mainly for laminar applications such as oil coolers, viscous process flows, and small liquid passages.

The Graetz number compresses the entrance question into a single dimensionless group. It can be read as the ratio of the time needed for heat to conduct across the duct radius to the time the fluid spends traveling the duct length:

Gz = (D / L) * ReD * Pr

Gz
Graetz number, dimensionless axial position that measures the importance of entrance effects
D
Duct diameter (m)
L
Duct length (m)
ReD
Reynolds number based on duct diameter, dimensionless
Pr
Prandtl number, dimensionless

A large Graetz number means a short duct in thermal terms: entrance effects dominate and the mean heat transfer coefficient is well above the fully developed value. A small Graetz number means a long duct where most of the surface sees fully developed conditions. The classic Hausen correlation for laminar flow with constant wall temperature is built directly on this idea:

Num = 3.66 + (0.0668 * Gz) / (1 + 0.04 * Gz2/3)

Num
Mean Nusselt number averaged over the duct length, dimensionless
Gz
Graetz number, dimensionless

As the duct length grows and Gz approaches zero, this expression recovers 3.66, the fully developed laminar Nusselt number for constant wall temperature in a circular tube. For short ducts the correction term raises the mean Nusselt number substantially, which is exactly the physics: thin developing boundary layers transfer heat far more effectively than the settled profile.

The practical message for duct design is that entrance effects almost always raise the heat transfer coefficient above the fully developed value, so using a fully developed correlation over an entire short duct underestimates heat transfer. In HVAC systems this matters in two ways. First, straight duct runs between elbows, takeoffs, and transitions are frequently only a few diameters long, so the flow is perpetually re-developing and the effective heat transfer exceeds textbook fully developed values. Second, many published correlations, including the widely used Dittus-Boelter correlation covered in Session 2, carry an explicit requirement such as L/D of at least 10, precisely because they were fitted to developed flow. Applying them to a two-diameter duct segment without correction is a misuse of the correlation, not a conservative simplification.

Common mistakes

Three errors appear repeatedly in practice. The first is applying the laminar entrance formula to turbulent flow, which predicts absurdly long entrance lengths; use the ten-diameter rule for turbulent flow instead. The second is confusing hydrodynamic development with thermal development: a flow can be fully developed in velocity while still thermally developing, particularly with high Prandtl number fluids. The third is treating an entire duct as fully developed when most of its length sits within entrance regions, for example a coil or a short duct run, which systematically underestimates the heat transfer coefficient.

Key takeaways

  • Internal flow develops from the entrance: boundary layers grow, merge, and then the profile stops changing, marking fully developed flow.
  • Laminar hydrodynamic entrance length scales with Reynolds number (0.05 * ReD * D); turbulent entrance length is roughly ten diameters regardless of Reynolds number.
  • Thermal entrance length follows the hydrodynamic length scaled by the Prandtl number in laminar flow, so high-Pr fluids stay thermally developing over long distances.
  • The Graetz number measures entrance importance: large Gz means entrance-dominated, small Gz means fully developed conditions prevail.
  • Developing flow transfers heat more effectively than fully developed flow, so entrance corrections raise, not lower, the predicted coefficient.

Session 2

Turbulent internal-flow correlations

Objectives

  • Apply the Dittus-Boelter correlation within its stated validity range
  • Apply the Gnielinski correlation using the Petukhov friction factor
  • Decide which correlation fits transition-region and rough-wall cases
  • Correct for property variation when wall and bulk temperatures differ strongly

Turbulent internal flow is the workhorse regime of HVAC engineering. The chaotic cross-stream motion of eddies carries heat from the wall deep into the fluid far more effectively than molecular conduction alone, so turbulent heat transfer coefficients are typically an order of magnitude above laminar values at the same mass flow. Because the governing equations cannot be solved in closed form for the mean heat transfer, engineers rely on empirical correlations built from the dimensionless groups that control the physics: the Reynolds number for the flow intensity, the Prandtl number for the fluid's thermal character, and the Nusselt number as the dimensionless heat transfer coefficient being predicted.

The most widely used of these is the Dittus-Boelter correlation, a simple power-law fit proposed in 1930 that has survived nearly a century of use because it is easy to evaluate and adequate for many design calculations:

NuD = 0.023 * ReD4/5 * Prn

NuD
Nusselt number based on duct diameter, NuD = h * D / k, dimensionless
ReD
Reynolds number based on duct diameter, dimensionless
Pr
Prandtl number of the fluid, dimensionless
n
Exponent: n = 0.4 when the wall heats the fluid (Ts greater than Tm), n = 0.3 when the wall cools the fluid (Ts less than Tm)

Fluid properties are evaluated at the bulk mean temperature, the mass-averaged fluid temperature at the cross section of interest. The correlation is stated to be valid for Reynolds numbers of 10,000 and above, Prandtl numbers from 0.6 to 160, and length-to-diameter ratios of at least 10. These limits are not decorative. Below Re of 10,000 the flow sits in the transition region where the correlation was never fitted, and predictions there can be substantially off. The L/D requirement restates the entrance lesson of Session 1: the correlation describes developed turbulent flow, and short ducts need entrance corrections.

The exponent n captures a real physical asymmetry. When the wall is hotter than the fluid, the fluid near the wall is warmer and, for liquids, less viscous than the bulk; the near-wall layer flows more easily and heat transfer improves slightly. When the wall is cooler, the effect reverses. For gases like air, viscosity increases with temperature, so the direction of the effect flips in detail, but the correlation's fitted exponents still give reasonable results across the stated range. In HVAC work with air at Prandtl number near 0.71, many engineers simply use n = 0.4 for heating and n = 0.3 for cooling without further thought, which is acceptable for ordinary design temperature differences.

When better accuracy is needed, or when the Reynolds number falls in the transition region between about 3000 and 10,000, the Gnielinski correlation from 1976 is the better tool. Rather than a pure power-law fit, it is built on the analogy between momentum transfer and heat transfer, which gives it a sounder physical footing across a wider range:

NuD = ((f / 8) * (ReD - 1000) * Pr) / (1 + 12.7 * (f / 8)1/2 * (Pr2/3 - 1))

NuD
Nusselt number based on duct diameter, dimensionless
f
Darcy friction factor for the duct, dimensionless
ReD
Reynolds number based on duct diameter, dimensionless
Pr
Prandtl number of the fluid, dimensionless

Gnielinski stated validity for Reynolds numbers from 3000 to 5,000,000 and Prandtl numbers from 0.5 to 2000, a range that comfortably covers everything from the transition region to highly turbulent flow, and from gases to viscous liquids. The correlation needs the Darcy friction factor, and for smooth tubes the Petukhov relation from 1970 is the standard companion:

f = (0.790 * ln(ReD) - 1.64)-2

f
Darcy friction factor for a smooth tube, dimensionless
ReD
Reynolds number based on duct diameter, dimensionless

Real HVAC ducts are not hydraulically smooth. Galvanized steel, duct liner, and flexible duct all carry surface roughness that raises the friction factor above the smooth-tube value. The accepted practice is to compute the friction factor for the actual roughness, for example with the Haaland equation covered in Session 5, and insert that friction factor into the Gnielinski correlation. This keeps the heat transfer prediction consistent with the pressure drop calculation, which matters because both feed into equipment selection and energy estimates.

A further refinement concerns property variation across the boundary layer. Both correlations evaluate properties at the bulk mean temperature, which is exact only when the wall and fluid temperatures are close. With large temperature differences, the fluid near the wall has markedly different viscosity and conductivity than the bulk, and the prediction drifts. Standard practice applies a property-ratio correction, most famously in the Sieder-Tate form that scales the Nusselt number by a viscosity ratio raised to a fitted power for laminar flow; analogous property-ratio corrections are used for turbulent flow. A simpler and common alternative for gases is to evaluate properties at the film temperature, the average of wall and bulk temperatures. Whichever approach is chosen, the key discipline is consistency: pick one property-evaluation scheme and apply it to both the heat transfer and the pressure drop calculations.

Choosing between the two correlations is straightforward in practice. For quick hand calculations in fully turbulent flow with moderate temperature differences, Dittus-Boelter is fast and adequate. For the transition region, for high-accuracy requirements, for rough ducts where the friction factor is known independently, or for fluids far from air, Gnielinski is the better choice. The numerical comparison below shows the two agreeing within a few percent in fully turbulent air flow, which is why Dittus-Boelter remains acceptable there, while only Gnielinski remains valid as the Reynolds number drops toward transition.

Worked example: comparing the two correlations

Air flows in a duct at ReD = 20,000 and Pr = 0.7, with the wall heating the fluid. Dittus-Boelter with n = 0.4 gives NuD = 0.023 * (20,000)0.8 * (0.7)0.4. Evaluating: (20,000)0.8 is approximately 2759, and (0.7)0.4 is approximately 0.867, so NuD is approximately 0.023 * 2759 * 0.867 = 55.0. For Gnielinski, the Petukhov friction factor is f = (0.790 * ln(20,000) - 1.64)-2 = (0.790 * 9.903 - 1.64)-2 = (6.184)-2 = 0.0262. Then f/8 = 0.00327, and NuD = (0.00327 * 19,000 * 0.7) / (1 + 12.7 * (0.00327)0.5 * (0.70.667 - 1)) = 43.5 / 0.846 = 51.4. The two correlations differ by about 7 percent, which is well within the scatter expected of empirical heat transfer correlations. At ReD = 5000, Dittus-Boelter would be outside its stated range and unreliable, while Gnielinski continues to give a defensible prediction. The interpretation: use Dittus-Boelter for fast estimates in fully turbulent flow, and reach for Gnielinski whenever the Reynolds number is low, the accuracy requirement is high, or the wall roughness is significant.

Common mistakes

The most frequent error is applying Dittus-Boelter below Re of 10,000, deep in the transition region where it was never validated. Close behind is choosing the wrong exponent n, or forgetting that n depends on whether the wall heats or cools the fluid. Another common slip is evaluating fluid properties at the wall temperature instead of the bulk mean temperature, which biases the Reynolds and Prandtl numbers. Finally, engineers sometimes pair a smooth-tube friction factor with a rough duct, or mix Fanning and Darcy friction factors (they differ by a factor of four), silently corrupting both the Gnielinski prediction and the pressure drop.

Key takeaways

  • Dittus-Boelter (NuD = 0.023 * ReD4/5 * Prn) is a fast, adequate correlation for fully turbulent flow with ReD of at least 10,000, Pr from 0.6 to 160, and L/D of at least 10.
  • The exponent is n = 0.4 for wall heating and n = 0.3 for wall cooling, with properties evaluated at the bulk mean temperature.
  • Gnielinski covers ReD from 3000 to 5,000,000 and Pr from 0.5 to 2000, making it the right choice for transition flow and higher-accuracy work.
  • Gnielinski needs the Darcy friction factor: use Petukhov for smooth tubes and a roughness-based value such as Haaland for real ducts.
  • Large wall-to-fluid temperature differences call for a property-ratio correction or film-temperature property evaluation, applied consistently.

Session 3

Non-circular ducts

Objectives

  • Define the hydraulic diameter and explain why it is the correct length scale for non-circular ducts
  • Apply circular-tube heat transfer correlations to rectangular HVAC ducts
  • Describe how aspect ratio affects the laminar fully developed Nusselt number
  • Recognize the limits of the hydraulic-diameter approximation

Nearly every correlation in Sessions 1 and 2 was developed for circular tubes, yet the ducts on a real job site are overwhelmingly rectangular, with round duct reserved mainly for branches, risers, and high-pressure systems. The bridge between the two geometries is the hydraulic diameter, a single length scale that lets circular-tube correlations serve rectangular, triangular, and other cross sections with surprising accuracy. Understanding where it comes from, and where it stops working, keeps duct heat transfer calculations both simple and honest.

The hydraulic diameter is defined as four times the flow cross-sectional area divided by the wetted perimeter:

Dh = 4 * Ac / P

Dh
Hydraulic diameter, the equivalent length scale for a non-circular duct (m)
Ac
Cross-sectional area of the flow passage (m2)
P
Wetted perimeter, the length of duct wall in contact with the fluid at a cross section (m)

The factor of four is not arbitrary. For a circular tube, Ac = pi * D2 / 4 and P = pi * D, so 4 * Ac / P reduces exactly to D. The hydraulic diameter therefore reproduces the ordinary diameter for the geometry the correlations were built on, and extends the same scaling to anything else. Physically, the ratio Ac / P is proportional to the fluid volume per unit of wall surface area along the duct, which is the quantity that governs how wall friction and wall heat transfer scale with size. For a rectangular duct with sides a and b, the area is a * b and the perimeter is 2 * (a + b), giving Dh = 2 * a * b / (a + b). A 400 mm by 200 mm supply duct, the example worked in detail in Session 6, has a hydraulic diameter of about 267 mm.

In laminar flow, the cross-sectional shape genuinely matters. The fully developed Nusselt number depends on the geometry and on the thermal boundary condition, and its values for rectangular ducts of various aspect ratios are tabulated in standard references, notably the compilation by Shah and London. For constant wall temperature, the circular-tube value of 3.66 and the parallel-plate value of 7.54 bracket the rectangular-duct results: as a rectangular duct becomes flatter, its fully developed Nusselt number migrates from the circular value toward the parallel-plate limit. The constant-heat-flux values follow the same pattern at slightly higher levels. This shape sensitivity is a laminar-flow phenomenon, rooted in how the velocity profile fills the corners and how the thermal boundary condition couples to it.

In turbulent flow, the picture simplifies enormously, which is fortunate because HVAC duct flow is turbulent in essentially all practical cases. Experiments have repeatedly shown that circular-tube correlations predict turbulent heat transfer in non-circular ducts quite well when the Reynolds and Nusselt numbers are both based on the hydraulic diameter. The intense eddy mixing that dominates turbulent transport washes out the geometric details that matter so much in laminar flow. So the working procedure for a rectangular supply duct is direct: compute Dh, form Re based on Dh, check that it exceeds 10,000, and apply Dittus-Boelter or Gnielinski exactly as for a round tube. Session 6 demonstrates this end to end. The same approach serves the friction factor in pressure drop calculations, where the Reynolds number in the Haaland equation and the diameter in the Darcy-Weisbach equation are both replaced by their hydraulic-diameter counterparts.

The approximation does have limits worth knowing. For ducts of extreme aspect ratio, very flat wide ducts behave thermally like flow between parallel plates, and the proper length scale trends toward twice the plate spacing rather than the hydraulic diameter; in practice such geometries appear as duct liners, ceiling plenums, and double-skin facades rather than as supply ducts. Triangular and other exotic sections follow the hydraulic-diameter approach reasonably in turbulent flow but deserve caution in laminar flow, where tabulated shape-specific values should be consulted. Entrance effects, covered in Session 1, apply in the same way: the entrance length scales with the hydraulic diameter, so the ten-diameter rule for turbulent flow becomes ten hydraulic diameters.

One subtle point concerns the wetted perimeter in ducts with internal obstructions. Turning vanes, splitter dampers, and internal liner surfaces all add wetted perimeter without adding flow area, which reduces the hydraulic diameter and raises both the friction factor and the heat transfer coefficient relative to a clean duct of the same outer dimensions. For ordinary design work this refinement is usually folded into fitting loss coefficients rather than computed explicitly, but it explains why a duct crowded with internal devices behaves as a "smaller" duct than its sheet-metal dimensions suggest.

Common mistakes

The classic error is replacing the duct with a circle of equal area and using that diameter in the correlations. Equal area does not mean equal perimeter, and since wall friction and heat transfer scale with the wall surface, the equal-area diameter gives the wrong answer; the hydraulic diameter is the correct equivalent. A second mistake is carrying laminar shape factors, such as the tabulated rectangular-duct Nusselt numbers, into turbulent calculations where they do not apply. A third is forgetting that Dh must be used consistently everywhere: in the Reynolds number, in the Nusselt number definition when converting back to the heat transfer coefficient, and in the pressure drop equation.

Key takeaways

  • The hydraulic diameter Dh = 4 * Ac / P reduces to the ordinary diameter for a circular tube and is the standard length scale for non-circular ducts.
  • For a rectangular duct with sides a and b, Dh = 2 * a * b / (a + b).
  • In laminar flow the fully developed Nusselt number depends on aspect ratio, lying between the circular-tube value 3.66 and the parallel-plate value 7.54 for constant wall temperature, with exact values tabulated by Shah and London.
  • In turbulent flow, circular-tube correlations applied with the hydraulic diameter predict rectangular-duct heat transfer well, which covers nearly all HVAC applications.
  • Use the hydraulic diameter consistently in the Reynolds number, the Nusselt number, and the pressure drop equations, and scale turbulent entrance lengths as ten hydraulic diameters.

Session 4

Mixed convection

Objectives

  • Define the Richardson number and interpret its magnitude
  • Distinguish aiding and opposing mixed convection in vertical ducts
  • Predict buoyancy-driven secondary-flow effects in horizontal ducts
  • Judge when buoyancy can be safely ignored in duct heat transfer calculations

Everything so far has assumed that the flow is driven entirely by the fan and that density differences within the fluid play no role. That assumption holds when velocities are high and temperature differences are modest, which describes most supply ductwork. But when temperature differences grow or velocities fall, the fluid's own weight begins to matter: warm fluid is lighter and tends to rise, cool fluid is heavier and tends to sink. These buoyancy forces superimpose a natural-convection motion on the forced flow, and the combined regime is called mixed convection. An engineer who never checks for it will occasionally, and expensively, be wrong about heat transfer in large low-velocity ducts, in ducts exposed to strong solar heating, and in any passage where the fan barely moves the air.

The strength of buoyancy relative to viscous forces is measured by the Grashof number:

GrD = g * beta * (Ts - Tm) * D3 / nu2

GrD
Grashof number based on duct diameter, the ratio of buoyancy forces to viscous forces, dimensionless
g
Gravitational acceleration (m/s2)
beta
Volume expansion coefficient of the fluid (1/K); for an ideal gas, beta = 1 / T where T is the absolute temperature
Ts
Wall surface temperature (K or C, as a difference)
Tm
Bulk mean fluid temperature (K or C, as a difference)
D
Duct diameter, or hydraulic diameter for non-circular ducts (m)
nu
Kinematic viscosity of the fluid (m2/s)

The quantity that decides the flow regime, though, is the Richardson number, which compares buoyancy directly against the inertia of the forced flow:

Ri = GrD / ReD2

Ri
Richardson number, the ratio of buoyancy forces to inertial forces, dimensionless
GrD
Grashof number based on duct diameter, dimensionless
ReD
Reynolds number based on duct diameter, dimensionless

The interpretation is deliberately simple. When Ri is much less than 1, inertia dominates and buoyancy can be ignored: pure forced convection applies, and every correlation in Sessions 2 and 3 is valid as written. When Ri is much greater than 1, buoyancy dominates and the flow behaves as natural convection in a confined passage. When Ri is of order 1, the two mechanisms are comparable and the heat transfer can differ markedly, usually upward, from either pure regime. In practice, engineers often treat Ri below about 0.1 as safely forced and Ri above about 10 as safely natural, with the decade in between requiring mixed-convection treatment. These are judgment thresholds, not sharp physical boundaries.

Orientation decides how buoyancy interacts with the forced flow. In a vertical duct, buoyancy acts along the flow direction, and two distinct cases arise. In aiding flow, buoyancy pushes in the same direction as the forced flow: the classic case is upward flow with heated walls, where the warm near-wall fluid is light and accelerates upward. The boundary layer thins, the wall shear rises, and the heat transfer coefficient increases above the forced-convection value. In opposing flow, buoyancy pushes against the forced flow: downward flow with heated walls is the standard example. The warm near-wall fluid wants to rise while the bulk flow descends, which decelerates the near-wall region, thickens the thermal boundary layer, and impairs heat transfer. If the opposition is strong enough, the near-wall fluid can stagnate or reverse, producing unstable, oscillatory behavior that no steady correlation can capture. Cooling the walls reverses the roles: upward flow with cooled walls becomes the opposing case.

Horizontal ducts present a different and very common situation. Here buoyancy acts across the duct rather than along it. With heated walls, the warm fluid near the wall rises toward the top of the duct while cooler core fluid sinks, setting up a pair of counter-rotating longitudinal vortices superimposed on the axial flow. This secondary motion continuously sweeps fresh fluid against the wall and typically enhances heat transfer relative to pure forced convection, sometimes substantially. With strong heating and low flow rates the stratification can become pronounced: a hot layer pools at the top of the duct while the bottom stays cool, and the circumferentially averaged heat transfer coefficient becomes a blunt instrument for describing a highly nonuniform reality. The top-to-bottom temperature asymmetry also matters mechanically, since it drives differential thermal expansion in large ducts.

For mainstream HVAC design, a quick Richardson number estimate usually settles the matter. Take a supply duct at 6 m/s with a 10 K wall-to-air temperature difference: the Grashof number is modest, the Reynolds number is large, and Ri comes out on the order of 0.002, as computed in the Session 6 worked example. Buoyancy is utterly negligible, and the analysis can proceed as pure forced convection with full confidence. Mixed convection earns attention in specific situations: large cross-section ducts at low velocity, where Re is small and D is large (note the D cubed in the Grashof number); rooftop or exposed ducts under strong solar heating, where the wall temperature can soar far above the air temperature; cold supply ducts running through hot ceiling plenums at part-load flow; and systems that rely partly on stack effect or natural ventilation, where the "forced" velocity itself is buoyancy-driven and the regime is mixed by construction.

When the Richardson number indicates that mixed convection matters, the honest approach is to consult mixed-convection correlations for the appropriate orientation and aiding or opposing condition, several of which appear in the standard heat transfer references. These correlations typically present the mixed Nusselt number as a correction to the forced-convection value, rising above it for aiding flow and horizontal ducts, and dipping below it for opposing flow. The correction is worth applying because the effect is not small: aiding and secondary-flow enhancement can raise heat transfer by tens of percent, while strong opposing flow can impair it by a similar margin. What the correlations cannot do is describe the unstable reversed-flow regime, which must be avoided by design rather than predicted precisely.

Common mistakes

The most common mistake is never computing the Richardson number at all, assuming forced convection always applies. It almost always does in supply ductwork, but "almost" is doing real work in that sentence. A second mistake is assuming buoyancy always helps: in opposing vertical flow it impairs heat transfer and can destabilize the flow. A third is applying vertical-duct aiding or opposing logic to horizontal ducts, where the physics is secondary vortices and stratification, not streamwise acceleration. Finally, using the duct diameter instead of the hydraulic diameter in the Grashof number for rectangular ducts silently corrupts the estimate, since Grashof scales with the cube of the length scale.

Key takeaways

  • Mixed convection arises when buoyancy forces are comparable to the inertia of the forced flow, measured by the Richardson number Ri = GrD / ReD2.
  • Ri much less than 1 means pure forced convection; Ri much greater than 1 means buoyancy-dominated flow; Ri near 1 requires mixed-convection treatment.
  • In vertical ducts, buoyancy aiding the flow enhances heat transfer, while buoyancy opposing the flow impairs it and can cause instability or reversal.
  • In horizontal ducts, buoyancy drives longitudinal secondary vortices and possible thermal stratification, usually enhancing the average heat transfer.
  • Typical HVAC supply ducts have Ri on the order of 0.001, so buoyancy is negligible; check Ri explicitly for large low-velocity ducts, solar-heated ducts, and part-load conditions.

Session 5

Pressure drop and fan power

Objectives

  • Apply the Darcy-Weisbach equation to straight-duct friction losses
  • Determine the Darcy friction factor from the Haaland equation
  • Account for fitting losses with loss coefficients
  • Compute fan power from volume flow rate and total pressure rise

Heat transfer determines what a duct does to the air thermally; pressure drop determines what the fan must do to the air mechanically. The two are coupled through the same velocity field and the same friction factor, which is why a course on duct convection must treat pressure drop as a first-class topic. Undersize the pressure analysis and the fan cannot deliver design airflow; oversize it and the client pays for fan energy and equipment for the life of the building. The calculation has three parts: friction along straight duct, losses in fittings, and the conversion of total pressure into fan power.

Friction loss in a straight duct section follows the Darcy-Weisbach equation, the central relation of duct and pipe flow:

Delta p = f * (L / D) * (rho * V2 / 2)

Delta p
Pressure drop due to wall friction over the duct length (Pa)
f
Darcy friction factor, dimensionless
L
Duct length (m)
D
Duct diameter, or hydraulic diameter for non-circular ducts (m)
rho
Fluid density (kg/m3)
V
Mean flow velocity in the duct (m/s)

The group rho * V2 / 2 is the velocity pressure, the kinetic energy per unit volume of the moving air, and the friction factor f is the dimensionless measure of how effectively the wall converts that kinetic energy into heat. In laminar flow the friction factor is known exactly, f = 64 / ReD, a result that follows directly from the parabolic velocity profile. In turbulent flow the friction factor depends on both the Reynolds number and the relative roughness of the wall, and it is read from the Moody chart or, more conveniently today, computed from an explicit approximation.

The Haaland equation, published in 1983, gives the Darcy friction factor explicitly without the iteration that the original Colebrook relation requires:

f-1/2 = -1.8 * log10[((epsilon / D) / 3.7)1.37 + 6.9 / ReD]

f
Darcy friction factor, dimensionless
epsilon
Equivalent sand-grain roughness height of the duct wall (m)
D
Duct diameter, or hydraulic diameter for non-circular ducts (m)
ReD
Reynolds number based on duct diameter, dimensionless

For new galvanized steel duct, the roughness height epsilon is commonly taken as 0.15 mm in design practice, a value that represents commercial steel with its seams, joints, and slight surface irregularities. Relative roughness epsilon / D then sets how far the friction factor sits above the smooth-tube line: in large ducts the relative roughness is small and the wall behaves as nearly smooth, while in small ducts the same absolute roughness is proportionally larger and friction rises. Duct liner, internally insulated duct, and especially flexible duct carry much higher effective roughness, and flexible duct additionally suffers when it is compressed, sagging, or bent, since every bend in a flexible run acts partly as a fitting. Manufacturer data or handbook values should be used for these materials rather than the bare-steel roughness.

Straight-duct friction is only half the story, and in many real systems it is the smaller half. Every elbow, transition, takeoff, damper, and grille disturbs the flow and dissipates mechanical energy, and these fitting losses are expressed through dimensionless loss coefficients:

Delta p = K * rho * V2 / 2

Delta p
Pressure drop across the fitting (Pa)
K
Minor loss coefficient for the fitting, dimensionless, referenced to the stated velocity
rho
Fluid density (kg/m3)
V
Reference mean velocity, usually the velocity in the duct at the fitting (m/s)

Loss coefficients are tabulated in standard design references such as ASHRAE Fundamentals and SMACNA duct design guidance, which devote extensive tables to them because fittings dominate so many systems. A few landmarks orient the intuition. A smooth radius elbow carries a small fraction of a velocity head, while a sharp mitered elbow without turning vanes can dissipate more than a full velocity head. A duct discharging directly into a room loses essentially one full velocity head, since the kinetic energy of the jet is dissipated in the space; this exit loss of K = 1.0 is a standard result. A sharp entrance from a plenum into a duct costs roughly half a velocity head. Dampers, fire dampers, and grilles each carry their own coefficients, usually obtained from manufacturer test data at the design flow. The critical discipline is consistency of the reference velocity: each K is defined against a particular velocity, normally the mean velocity in the duct at that fitting, and mixing reference velocities between fittings corrupts the total.

The total pressure drop the fan must overcome is the sum of all straight-duct friction terms and all fitting losses along the critical path, the run from the fan to the terminal with the highest resistance. In systems with many elbows, takeoffs, and dampers, the fitting total frequently exceeds the straight-duct friction by a large factor, which is why experienced designers obsess over fitting selection and duct routing rather than shaving a few pascals off straight runs. Once the total pressure is known, the air power required is the product of volume flow rate and pressure rise, divided by the overall fan efficiency:

P = Q * Delta p / etafan

P
Electrical or shaft power required by the fan (W)
Q
Volume flow rate of air (m3/s)
Delta p
Total pressure rise the fan must produce, straight friction plus fittings (Pa)
etafan
Overall efficiency from air power to input power, accounting for the fan, motor, and drive, dimensionless

The efficiency term deserves respect: it folds together the aerodynamic efficiency of the fan wheel, the motor efficiency, and any belt or drive losses, and it is properly taken from manufacturer selection data at the operating point rather than guessed. Two scaling laws, the fan affinity laws, govern how power responds to changes in flow and are among the most economically consequential relations in HVAC. At fixed system geometry, pressure drop varies with the square of the flow rate and fan power varies with the cube. Reducing airflow by 20 percent cuts fan power nearly in half. This cubic relationship is the entire economic logic behind variable-air-volume systems and variable-speed drives: trimming flow to match the actual load, rather than running constant volume against throttled dampers, converts directly into large energy savings.

For preliminary sizing, many designers work toward a straight-duct friction rate on the order of 1 Pa/m, then add fitting losses explicitly and check noise. This friction rate is a practice target that balances duct size against fan energy, not a code requirement, and it is adjusted up or down with project priorities: lower for quiet or energy-critical projects, higher where space is tight. The equal-friction sizing method extends the idea across a whole system by sizing each duct section for roughly the same pressure drop per meter, which tends to produce self-balancing layouts, with balancing dampers trimming the residual differences between branches.

Common mistakes

The costliest mistake is omitting fitting losses or estimating them with a crude length allowance, then discovering the fan cannot make the design flow. Next is confusing the Darcy friction factor with the Fanning friction factor: the two differ by a factor of four, and grabbing the wrong one from a reference silently scales every pressure drop. Mixing reference velocities between loss coefficients is another quiet corruptor. Finally, many analyses forget that flexible duct is not rigid duct: its pressure drop per meter can be several times higher, and a compressed or kinked flexible run behaves like a partially closed damper.

Key takeaways

  • Straight-duct friction follows Darcy-Weisbach: Delta p = f * (L / D) * (rho * V2 / 2), with f = 64 / ReD in laminar flow.
  • The Haaland equation gives the turbulent Darcy friction factor explicitly from Reynolds number and relative roughness; commercial steel duct is commonly taken at epsilon = 0.15 mm.
  • Fitting losses use Delta p = K * rho * V2 / 2 with K values from ASHRAE, SMACNA, or manufacturer data, always referenced to a consistent velocity.
  • Fitting losses frequently dominate straight-duct friction in real systems, so fitting selection and routing deserve the most design attention.
  • Fan power is P = Q * Delta p / etafan, and the affinity laws (pressure with the square, power with the cube of flow) make flow reduction the most powerful energy lever available.

Session 6

Worked example: sizing a rectangular supply duct

Objectives

  • Carry out a complete hand calculation for a rectangular supply duct from geometry to heat gain
  • Compute the hydraulic diameter, Reynolds number, and friction factor for a real duct
  • Determine the pressure drop per meter and the internal convection coefficient
  • Interpret the results against HVAC design practice: velocity, friction rate, and insulation needs

Theory becomes engineering judgment only when the numbers are worked by hand at least once. This session does exactly that for a completely ordinary piece of HVAC hardware: a rectangular galvanized steel supply duct. Every step uses the tools built in Sessions 1 through 5, and every result is checked against the kind of practice rules introduced in Session 7. Follow the arithmetic closely, because the pattern here, geometry to Reynolds number to friction factor to pressure drop to heat transfer coefficient to thermal consequence, is the same pattern behind every duct calculation you will ever perform.

The scenario is a straight horizontal supply duct, 400 mm by 200 mm in cross section, carrying air at 20 C with a mean velocity of 6 m/s. The duct is commercial galvanized steel. We want four answers: the pressure drop per meter of straight duct, the internal convection coefficient, the wall heat gain per meter for a 10 K temperature difference between the air and the surroundings, and a judgment on whether this duct is sensibly sized.

Worked example: 400 mm x 200 mm supply duct at 6 m/s

Given data. Duct cross section: a = 0.40 m, b = 0.20 m. Air at 20 C: rho = 1.204 kg/m3, mu = 1.81e-5 kg/m-s, nu = 15.1e-6 m2/s, k = 0.0259 W/m-K, Pr = 0.713, cp = 1007 J/kg-K. Mean velocity V = 6 m/s. Wall roughness for commercial steel: epsilon = 0.15 mm = 1.5e-4 m, the commonly used roughness value for this material. Temperature difference between air and surroundings: Delta T = 10 K.

Step 1: Geometry and hydraulic diameter

Cross-sectional area Ac = 0.40 * 0.20 = 0.080 m2. Wetted perimeter P = 2 * (0.40 + 0.20) = 1.20 m. Hydraulic diameter Dh = 4 * Ac / P = 4 * 0.080 / 1.20 = 0.267 m (267 mm). The volume flow rate is Q = Ac * V = 0.080 * 6 = 0.48 m3/s, roughly 1020 cfm, a typical flow for a medium-size zone.

Step 2: Reynolds number

ReDh = V * Dh / nu = 6 * 0.267 / 15.1e-6 = 1.60 / 15.1e-6 = 106,000 (approximately). The flow is fully turbulent, well above the 10,000 threshold, so the Dittus-Boelter correlation is applicable provided the straight run is at least about ten hydraulic diameters, roughly 2.7 m, long.

Step 3: Friction factor from Haaland

Relative roughness epsilon / Dh = 1.5e-4 / 0.267 = 5.63e-4. Insert into the Haaland equation: f-1/2 = -1.8 * log10[((5.63e-4) / 3.7)1.37 + 6.9 / 106,000]. The roughness term evaluates to (1.52e-4)1.37 = 5.88e-6 and the Reynolds term to 6.51e-5, giving a bracket sum of 7.10e-5. Then f-1/2 = -1.8 * (-4.149) = 7.47, so f = (1 / 7.47)2 = 0.0179. For reference, the smooth-tube Petukhov value at this Reynolds number is 0.0178, confirming that this duct is nearly hydraulically smooth at the design flow.

Step 4: Pressure drop per meter

Velocity pressure rho * V2 / 2 = 1.204 * 36 / 2 = 21.7 Pa. Darcy-Weisbach per meter: Delta p / L = f * (1 / Dh) * (rho * V2 / 2) = 0.0179 * (1 / 0.267) * 21.7 = 0.0179 * 81.3 = 1.46 Pa/m. One meter of this straight duct costs about one and a half pascals of pressure.

Step 5: Convection coefficient from Dittus-Boelter

Take n = 0.4, corresponding to the duct wall warming the air. NuDh = 0.023 * ReDh0.8 * Pr0.4 = 0.023 * (106,000)0.8 * (0.713)0.4. Now (106,000)0.8 is approximately 10,474 and (0.713)0.4 is approximately 0.873, so NuDh = 0.023 * 10,474 * 0.873 = 210. The heat transfer coefficient is h = NuDh * k / Dh = 210 * 0.0259 / 0.267 = 20.4 W/m2-K.

Step 6: Wall heat gain per meter

With a 10 K difference between the air and the surroundings, the heat gain per meter of duct is Q' = h * P * Delta T = 20.4 * 1.20 * 10 = 245 W/m. This is stated as an upper-bound estimate: it neglects the conduction resistance of the steel wall, which is negligible anyway, and the external convection and radiation resistances, which are not. The true heat gain with external resistances included would be lower, but this bound shows the scale of the effect.

Step 7: Buoyancy check

GrDh = g * beta * Delta T * Dh3 / nu2 with beta = 1 / 293 = 0.00341 1/K: Gr = 9.81 * 0.00341 * 10 * (0.267)3 / (15.1e-6)2 = 2.78e7. Richardson number Ri = Gr / Re2 = 2.78e7 / (106,000)2 = 0.0025, far below 1. Buoyancy is negligible and the pure forced-convection analysis stands.

Now the interpretation, which is where the calculation earns its keep. First, the velocity: 6 m/s sits at the bottom of the typical 6 to 9 m/s range used for commercial main supply ducts, so it is entirely reasonable, quiet enough for most applications and not wasteful of duct material. Second, the pressure drop: 1.46 Pa/m for straight duct is acceptable and sits near the friction-rate targets designers commonly aim for, though slightly above a nominal 1 Pa/m target. Remember that this is straight duct only. In a real run, elbows, takeoffs, and dampers will add substantially, often doubling or tripling the total, so the system curve must include the fittings before a fan is selected.

Third, the heat gain: 245 W/m is a large number wearing an innocent expression. Over a 15 m run through an unconditioned ceiling plenum or attic, that is roughly 3.7 kW of unwanted heat picked up by the supply air, a meaningful fraction of the cooling capacity serving a typical zone, and it arrives as a direct penalty on both energy use and delivered air temperature. Because the estimate is an upper bound, the real figure is smaller, but not small enough to ignore. The conclusion is unambiguous: insulate supply ducts that run through unconditioned spaces, and where condensation is possible on cold duct surfaces, the insulation needs a proper vapor retarder as well. An uninsulated 400 by 200 duct at these conditions is not a minor oversight; it is a measurable load the mechanical system must carry for its entire life.

Finally, note what the calculation did not need: no CFD, no charts, no software. Six hand steps with standard correlations sized the duct thermally and hydraulically, checked the flow regime, and produced a defensible insulation decision. That is the working method this course is built to instill.

Common mistakes

In hand calculations like this one, the errors are almost always unit or definition slips rather than conceptual ones. Forgetting to convert millimeters to meters before forming the hydraulic diameter is the classic. Using the duct's longest side, or the diameter of an equal-area circle, instead of the hydraulic diameter corrupts the Reynolds number, the friction factor, and the Nusselt number all at once. Evaluating air properties at the wall temperature instead of the bulk mean temperature biases everything downstream. And applying Dittus-Boelter to a short duct run without checking the L/D requirement, or to a flow below Re of 10,000, uses the correlation outside the conditions it was built for.

Key takeaways

  • A 400 mm x 200 mm duct has Dh = 267 mm; at 6 m/s the Reynolds number is about 106,000, fully turbulent.
  • The Haaland friction factor is 0.0179, giving a straight-duct pressure drop of 1.46 Pa/m, acceptable but only part of the system total.
  • Dittus-Boelter gives Nu = 210 and h = 20.4 W/m2-K for air being warmed by the duct wall.
  • The upper-bound heat gain of 245 W/m shows why ducts in unconditioned spaces must be insulated: a 15 m run approaches 3.7 kW of unwanted load.
  • The Richardson number of 0.0025 confirms buoyancy is negligible, so the forced-convection analysis is complete and self-consistent.

Session 7

HVAC design rules of thumb

Objectives

  • Apply the duct velocity ranges commonly used in design practice
  • Balance pressure drop, noise, and fan energy when selecting duct sizes
  • Size duct systems to avoid the most frequent failure modes
  • Identify and avoid the classic mistakes of duct design

The correlations of the previous sessions give exact answers to well-posed questions. Real design adds a harder problem: choosing the questions. Duct velocity, friction rate, aspect ratio, and insulation level are not dictated by physics alone; they are compromises between first cost, operating cost, space, noise, and comfort. The rules of thumb in this session are the distilled compromises of the industry, the ranges within which experienced designers work before refining with calculation. They come from design references such as ASHRAE Fundamentals and SMACNA duct design guidance, and they should be understood as typical practice ranges, not as code mandates. Codes govern fire, smoke, and life safety; velocity selection remains an engineering decision.

Velocity is the master variable because nearly everything scales with it. The commonly taught practice ranges run roughly as follows: main supply ducts in commercial systems, about 6 to 9 m/s; branch ducts serving individual zones, lower, typically in the 4 to 6 m/s region; risers and shafts, moderate values between the two; residential systems, lower still, reflecting smaller fans, shorter runs, and stricter noise expectations. Terminal devices, the grilles and diffusers that actually deliver air to rooms, operate far below duct velocities, and their selection follows manufacturer throw and noise data rather than any duct rule. These ranges exist because velocity is a three-way bargain. Higher velocity means smaller ducts, which saves sheet metal, hangers, and precious ceiling space. But pressure drop grows with the square of velocity and fan power with the cube, so each extra meter per second is bought with permanently higher energy bills. Lower velocity reverses the trade: cheap energy, expensive and bulky ducts.

Noise is the constraint that most often overrules the energy optimum. Air noise is generated wherever the flow is disturbed: at elbows, at dampers, at takeoffs, and especially at terminal devices discharging into occupied rooms. Room background noise is rated against NC, noise criteria, curves, and each space type carries an expectation: an open office tolerates more background sound than a conference room or a classroom. The design response is to keep velocities low near occupied spaces and at terminals, to prefer smooth radius fittings over abrupt ones in noise-sensitive paths, and to locate dampers and takeoffs away from diffusers so that the turbulence they generate has duct length in which to decay. A system that meets its pressure budget but whistles at every grille has failed just as surely as one whose fan cannot make the flow.

The pressure side of the design is managed through the friction rate, the pressure drop per meter of straight duct. Many designers begin sizing with a target on the order of 1 Pa/m for straight duct, then add fitting losses explicitly and check the noise implications. This target is a starting compromise, not a law: energy-critical or noise-critical projects size lower, space-constrained retrofits size higher. The equal-friction method systematizes the approach by sizing every section of the system for approximately the same friction rate, which produces layouts that are nearly self-balancing, with opposed-blade or splitter dampers trimming the residual differences between branches. A related classical method, static regain, sizes downstream sections so that the drop in velocity pressure at each expansion is converted back into static pressure, holding the available static pressure roughly constant along the run. Both methods then face the same reality check from Session 5: fittings often contribute more than the straight duct, so the fitting schedule deserves as much care as the duct sizes.

Fan energy ties the whole design together through the relation P = Q * Delta p / etafan

P = Q * Delta p / etafan

P
Fan input power (W)
Q
Volume flow rate of air (m3/s)
Delta p
Total pressure rise across the fan (Pa)
etafan
Overall fan, motor, and drive efficiency, dimensionless

and the affinity laws that make flow reduction the most powerful efficiency measure available. A variable-air-volume system that trims flow to match the actual load exploits the cubic relationship between flow and power, which is why modern commercial systems rarely run constant volume. Duct leakage deserves mention in the same breath: every cubic meter per second that escapes through unsealed joints is flow the fan paid to move that never reaches the room, so sealing ducts to a proper leakage class is among the cheapest energy measures in the building. Insulation closes the thermal side: supply ducts in unconditioned spaces gain or lose heat along their length, as Session 6 quantified, and cold duct surfaces in humid environments condense, so insulation with a continuous vapor retarder is both an energy and a moisture-control decision.

Aspect ratio is a quieter but real design variable. A square duct carries a given flow with less sheet metal per unit of delivered air and lower pressure drop than a flat wide duct of the same area, because the flat duct has more wall surface per unit of flow area. Deep ceiling plenums therefore favor near-square ducts, and flat ducts should be reserved for cases where the available depth truly forces them. When flat ducts are unavoidable, the higher pressure drop and material cost should be carried explicitly in the design rather than discovered during balancing.

Balancing deserves a final word because it is where good designs are confirmed or quietly ruined. Dampers at each branch allow the airflows to be set to design values during commissioning, but dampers can only redistribute the flow the fan provides; they cannot create airflow the system cannot deliver. A design that depends on heavily throttled dampers to force balance is a design whose duct sizes were wrong: the throttling wastes fan energy permanently to correct a sizing error made once. Better to size the branches so that their natural pressure drops are close, leaving the dampers for fine trimming.

Common mistakes

The classic mistakes form a familiar gallery. Oversized ducts waste money, steal ceiling space, and can run at velocities so low that air distribution at the terminals suffers. Undersized ducts scream: noise complaints, excessive fan power, and comfort problems follow. Ignoring fitting losses produces fans that cannot deliver design flow on day one. Forgetting insulation on ducts in unconditioned spaces, or omitting the vapor retarder on cold ducts, creates energy waste and condensation damage that persists for decades. Sizing an entire system at one velocity ignores the different noise and pressure economics of mains, branches, and terminals. Overusing flexible duct, especially with compression, sag, or sharp bends, introduces pressure drops far above the rigid-duct calculation. And skipping proper sealing turns the duct system into a leaky distribution network that no amount of fan power can fully correct.

Key takeaways

  • Typical practice velocities are about 6 to 9 m/s for commercial main supply ducts, lower for branches, and lower still for residential systems; these are design-practice ranges from references such as ASHRAE Fundamentals and SMACNA guidance, not code mandates.
  • Velocity is a three-way trade between duct size and first cost, fan energy through the square and cube laws, and noise against NC criteria expectations.
  • Size straight duct toward a friction rate on the order of 1 Pa/m as a starting target, then add fitting losses explicitly, since fittings often dominate the system total.
  • Keep ducts near square where plenum depth allows; flat high-aspect-ratio ducts cost more material and pressure drop per unit of airflow.
  • Seal ducts against leakage, insulate ducts in unconditioned spaces with a proper vapor retarder on cold surfaces, and size branches to balance naturally rather than by heavy damper throttling.

References

References

Every source below is an established engineering textbook, an ASHRAE or ISO standard or handbook, or a peer-reviewed paper from a Q1/Q2 journal. Inline citations in the sessions above refer to these entries.

  1. Dittus, F. W. and Boelter, L. M. K., "Heat transfer in automobile radiators of the tubular type," University of California Publications in Engineering, Vol. 2, p. 443, 1930.
  2. Gnielinski, V., "New equations for heat and mass transfer in turbulent pipe and channel flow," International Chemical Engineering, Vol. 16, pp. 359-368, 1976.
  3. Petukhov, B. S., "Heat transfer and friction in turbulent pipe flow with variable physical properties," Advances in Heat Transfer, Vol. 6, pp. 503-564, 1970.
  4. Haaland, S. E., "Simple and explicit formulas for the friction factor in turbulent pipe flow," Journal of Fluids Engineering (ASME), Vol. 105, pp. 89-90, 1983.
  5. Shah, R. K. and London, A. L., Laminar Flow Forced Convection in Ducts, Advances in Heat Transfer, Supplement 1, Academic Press, 1978.
  6. Incropera, F. P., DeWitt, D. P., Bergman, T. L. and Lavine, A. S., Fundamentals of Heat and Mass Transfer, 6th ed., Wiley, 2007.
  7. ASHRAE Handbook: Fundamentals, American Society of Heating, Refrigerating and Air-Conditioning Engineers (duct design chapter).