Dittus-Boelter Correlation — Turbulent Pipe Flow Heat Transfer
Theory and Physics
Overview — What is the Dittus-Boelter Equation?
Professor, where is the Dittus-Boelter equation used? It appears in every chapter of the textbook, but honestly, I don't quite grasp it...
It's the simplest equation for estimating the convective heat transfer coefficient $h$ for turbulent flow inside pipes. The form of the equation is:
Use $n = 0.4$ when the fluid is being heated and $n = 0.3$ when it's being cooled. Cooling water pipes, flow channels between heat sink fins, heat exchangers in chemical plants—this is the first equation used in the initial design of any scenario where "fluid flows turbulently inside a pipe."
Wait, can you really find the heat transfer coefficient with just this simple equation? Just plug in Re and Pr?
Yes, that's precisely why it's the "king of initial estimates." For example, when designing an automotive radiator pipe, once the coolant flow velocity and the pipe's inner diameter are decided, you calculate Re, plug in the Pr of water (about 7 at room temperature), and immediately get an estimated value for $h$. This allows you to make the first judgment: "Is this pipe diameter sufficient for cooling?"
But you said "estimate," which means there are situations where it's not reliable, right?
Sharp observation. This equation has prerequisites. It's only applicable for fully developed turbulent flow with Re > 10,000 and for sufficiently long pipes ($L/D > 60$). For short pipes with strong entrance effects or for the transition region (Re = 2,300–10,000), the Gnielinski equation is much more accurate. If you apply this equation without knowing its "applicable range," you can easily end up with errors of 20–30%.
Core Equations
Let's organize the two equations that constitute the Dittus-Boelter equation.
Nusselt number correlation:
The definitions of each dimensionless number are as follows:
Conversion to heat transfer coefficient:
| Symbol | Physical Quantity | SI Unit |
|---|---|---|
| $D$ | Pipe inner diameter | m |
| $u$ | Average flow velocity in pipe | m/s |
| $\rho$ | Fluid density | kg/m³ |
| $\mu$ | Dynamic viscosity (viscosity coefficient) | Pa·s |
| $c_p$ | Specific heat at constant pressure | J/(kg·K) |
| $k$ | Fluid thermal conductivity | W/(m·K) |
| $h$ | Convective heat transfer coefficient | W/(m²·K) |
| $n$ | Pr exponent (heating 0.4 / cooling 0.3) | — |
Applicability and Limitations
Could you tell me more specifically about the applicability conditions? There are limitations besides "Re > 10,000," right?
Yes, to summarize:
| Condition | Range | Consequence of Violation |
|---|---|---|
| Reynolds number | $\mathrm{Re}_D > 10{,}000$ | Overestimation in transition/laminar flow. Use Gnielinski equation. |
| Prandtl number | $0.7 < \mathrm{Pr} < 160$ | Large error for liquid metals (Pr≪1) or high-viscosity oils (Pr≫160) |
| Pipe length/diameter ratio | $L/D > 60$ | Actual h is higher due to entrance effects. Entrance correction factor needed. |
| Difference between wall and fluid temperature | Moderate (small $\Delta T$) | Large temperature difference causes significant viscosity change → Use Sieder-Tate equation. |
So $L/D > 60$ means for a pipe with an inner diameter of 20mm, it needs to be longer than 1.2m?
Exactly. In cold plates for electronics or flow channels between fins of small heat sinks, $L/D$ is often around 10–30. In such short pipes, heat transfer becomes locally very high near the entrance, so assuming a uniform $h$ with the Dittus-Boelter equation leads to underestimation. In practice, it's safer to multiply by an entrance correction factor $(1 + (D/L)^{0.7})$ or switch to the Gnielinski equation.
Physical Meaning of Each Term
Re0.8 — Why the 0.8 exponent?
The Reynolds number is the ratio of "inertial force/viscous force" and is an indicator of turbulence intensity. The exponent is 0.8, not 1, because in the turbulent region, the Nusselt number increases somewhat more gently than linearly with Re. This originates from the development pattern of the turbulent boundary layer—the viscous sublayer near the wall becomes thinner with increasing velocity but never completely disappears. Many experiments have confirmed that the best fit with experimental data occurs around 0.8.
Prn — Why is the exponent different for heating and cooling?
The Prandtl number is the ratio of "momentum diffusivity/thermal diffusivity." When a fluid is heated, viscosity near the wall decreases, thinning the velocity boundary layer → enhanced heat transfer → larger exponent (0.4). When cooled, viscosity near the wall increases, thickening the velocity boundary layer → suppressed heat transfer → smaller exponent (0.3). This asymmetry is a simplified way to incorporate the effect of temperature-dependent viscosity changes.
Meaning of the coefficient 0.023
This coefficient is not derived theoretically; it is an empirical constant determined by regression against a large amount of experimental data for circular pipe flow. According to literature by Winterton (1998), the original paper's experimental data had a scatter of about ±25%, and 0.023 is the central value.
Historical Background
The True Originator Problem
The equation known as the "Dittus-Boelter equation," $\mathrm{Nu} = 0.023\,\mathrm{Re}^{0.8}\,\mathrm{Pr}^n$, originates from a paper published in 1930 by F.W. Dittus and L.M.K. Boelter. However, recent literature review (Winterton, 1998) revealed that the differentiation of the exponent $n$ for heating/cooling (0.4/0.3) was a modification made by Boelter alone in 1942 and was not present in Dittus's original paper. Furthermore, the original paper was not a formally peer-reviewed journal but an internal report from the University of California, making the original itself difficult to obtain. Nevertheless, because it has been cited so widely, the name "Dittus-Boelter equation" has become established—an interesting episode in engineering history.
Numerical Solution and Implementation
Procedure for Calculating the Heat Transfer Coefficient
Could you teach me the step-by-step procedure to actually find $h$?
The procedure is as follows:
- Determine film temperature: $T_f = (T_w + T_b)/2$ (average of wall temperature $T_w$ and bulk temperature $T_b$)
- Obtain property values: Read $\rho$, $\mu$, $c_p$, $k$ at $T_f$ from property tables
- Calculate Re: $\mathrm{Re}_D = \rho u D / \mu$
- Calculate Pr: $\mathrm{Pr} = \mu c_p / k$
- Check applicability: Re > 10,000 and 0.7 < Pr < 160 and L/D > 60
- Calculate Nu: $\mathrm{Nu}_D = 0.023\,\mathrm{Re}_D^{0.8}\,\mathrm{Pr}^n$
- Calculate $h$: $h = \mathrm{Nu}_D \cdot k / D$
What if the film temperature is unknown—for example, at the very beginning of design when the wall temperature isn't decided yet?
Good question. Initially, assume a wall temperature, calculate $h$, then recalculate the wall temperature using the obtained $h$ → update property values again → recalculate $h$... iterate this loop. In practice, 2–3 iterations are sufficient for convergence. For fluids like water where property changes with temperature are small, using property values at the bulk temperature keeps the error within a few percent.
Calculation Example — Cooling Water Pipe
I'd like to see it with concrete numbers. For example, an engine cooling water system.
Let's try with typical engine cooling water conditions.
Conditions: Inner diameter $D = 0.02$ m, flow velocity $u = 1.5$ m/s, water temperature 80°C (heating state)
Water properties at 80°C:
- $\rho = 972$ kg/m³
- $\mu = 3.55 \times 10^{-4}$ Pa·s
- $c_p = 4,197$ J/(kg·K)
- $k = 0.670$ W/(m·K)
Step 1 — Re:
$$ \mathrm{Re} = \frac{972 \times 1.5 \times 0.02}{3.55 \times 10^{-4}} = 82{,}028 $$→ Fully turbulent. OK.
Step 2 — Pr:
$$ \mathrm{Pr} = \frac{3.55 \times 10^{-4} \times 4{,}197}{0.670} = 2.22 $$→ Within range. OK.
Step 3 — Nu (heating, so $n = 0.4$):
$$ \mathrm{Nu} = 0.023 \times 82{,}028^{0.8} \times 2.22^{0.4} = 0.023 \times 10{,}173 \times 1.37 = 320.6 $$Step 4 — $h$:
$$ h = \frac{320.6 \times 0.670}{0.02} = 10{,}740 \;\text{W/(m²·K)} $$So approximately $h \approx 10{,}700$ W/(m²·K). A typical value for automotive cooling water systems.
Wow, that's quite a large value. What about forced convection with air?
For air, Pr ≈ 0.71 and $k$ ≈ 0.03 W/(m·K), which are orders of magnitude smaller than water, so even with the same velocity and pipe diameter, $h$ would be around 50–200 W/(m²·K). Water's $h$ is 50–100 times larger than air's because water has much higher thermal conductivity and heat capacity. That's why water is used for engine cooling.
Property Evaluation Temperature
At Which Temperature to Take Property Values—The Most Overlooked Point in Practice
The Dittus-Boelter equation was originally derived from experimental data where the difference between wall and bulk temperatures was small. Therefore, property values should be evaluated at the film temperature $T_f = (T_w + T_b)/2$ as standard. However, when the wall temperature is unknown in the initial stage, it's common to substitute the bulk temperature, then use an iterative method: calculate $h$, recalculate the wall temperature, and update the property values.
When the temperature difference between the wall and fluid is large (e.g., hot oil being rapidly cooled at a pipe wall), the viscosity near the wall changes significantly. In such cases, the Dittus-Boelter equation cannot adequately handle this, and one should switch to the Sieder-Tate equation, which includes a viscosity correction term $(\mu_b/\mu_w)^{0.14}$.
Practical Guide
Setting Boundary Conditions in CAE
When doing thermal analysis with CAE software, how is the Dittus-Boelter equation used? Do you input the equation into the software?
There are two main ways to use it:
- Manually set $h$ as a convection boundary condition: When modeling only the solid without modeling the fluid side. Specify the $h$ calculated by the Dittus-Boelter equation and the fluid bulk temperature $T_\infty$ on the wall. This is the most common approach in thermal analysis with Ansys Mechanical or Abaqus.
- Automatic internal calculation by CFD solver: When directly solving pipe flow in Fluent or STAR-CCM+, the solver automatically calculates $h$ through wall functions. In this case, there's no need to manually use the Dittus-Boelter equation. Rather, it's used for validating the reasonableness of CFD results by comparing them with the Dittus-Boelter equation's estimated values.
So the Dittus-Boelter equation comes into play for "thermal analysis without solving the fluid," right? Conversely, if solving directly with CFD, it's for verification?
Exactly. In practice, a common flow is "initial estimation using the Dittus-Boelter equation in the early design stage → CFD in the detailed design stage." If the initial estimated $h$ differs significantly from the CFD result, it's likely outside the applicable conditions.
Choosing Between the Gnielinski Equation
Could you explain the difference with the Gnielinski equation in more detail? I often get confused about which one to use...
Let's look at a comparison table:
| Item | Dittus-Boelter Equation | Gnielinski Equation |
|---|---|---|
| Equation Form | $0.023\,\mathrm{Re}^{0.8}\,\mathrm{Pr}^n$ | $\frac{(f/8)(\mathrm{Re}-1000)\mathrm{Pr}}{1+12.7\sqrt{f/8}(\mathrm{Pr}^{2/3}-1)}$ |
| Re Applicable Range | > 10,000 | 2,300 – 5×10⁶ |
| Pr Applicable Range | 0.7 – 160 | 0.5 – 2,000 |
| Entrance Effect | Not considered (assumes L/D > 60) | Can be addressed with correction term |
| Accuracy (vs. experiment) | ±25% | ±10% |
| Calculation Effort | Immediate calculation with a calculator | Additional calculation of friction factor $f$ required |
| Optimal Use | Initial estimation, sensitivity analysis | Detailed design, basis for reports |
So if accuracy is needed, Gnielinski is the only choice, and Dittus-Boelter is for "when you want a rough idea"?
Basically, yes. However, in the fully turbulent region with Re > 10,000, the difference between the two is often within a few percent. That is, for fully developed turbulent flow in long pipes, Dittus-Boelter is sufficient. Conversely, the Gnielinski equation truly shines in the transition region (Re = 2,300–10,000) or for short pipes with significant entrance effects.
Choosing Between the Sieder-Tate Equation
I'd also like to know the difference with the Sieder-Tate equation.
The Sieder-Tate equation has the form $\mathrm{Nu} = 0.027\,\mathrm{Re}^{0.8}\,\mathrm{Pr}^{1/3}\,(\mu_b/\mu_w)^{0.14}$, which includes the viscosity ratio $(\mu_b/\mu_w)^{0.14}$ between the wall and the fluid bulk.
In other words, it's an equation used for cases where the temperature difference between the wall and fluid is large, causing significant viscosity changes. For example, in a crude oil pipeline where the wall is heated by a heater, viscosity near the wall drastically decreases. In such scenarios, the Dittus-Boelter equation is insufficient, and the viscosity correction of the Sieder-Tate equation becomes effective.
Conversely, for fluids like water or air where viscosity changes with temperature are small, the results of the Dittus-Boelter and Sieder-Tate equations are almost identical.
Industrial Application Examples
| Field | Typical Conditions | Estimated $h$ [W/(m²·K)] | Notes |
|---|---|---|---|
| Automotive Radiator Pipes | Water, $u$ = 1–2 m/s, $D$ = 15–25 mm | 5,000–15,000 | Used around 80°C |
| Chemical Plant Heat Exchangers | Organic solvents, $u$ = 0.5–3 m/s | 500–5,000 | Relatively high h due to large Pr |
| Air Conditioning Ducts | Air, $u$ = 5–15 m/s, $D$ = 100–300 mm | 30–100 | Low h due to low k of air |
| Data Center Liquid Cooling | Inert liquids, $u$ = 0.5–2 m/s | 2,000–8,000 | Caution with L/D in microchannels |
| Gas Turbine Cooling Passages | Air, $u$ = 30–100 m/s | 200–1,000 | High Re, rotation effects require correction |
Common Use in Hot Water Pipe Design
In Japanese hot water system design (based on JIS B 8417), the Dittus-Boelter equation is widely used for estimating forced convection inside pipes, yielding $h$ = 3,000–8,000 W/(m²·K) for hot water pipes with flow velocities of 0.5–2 m/s. When calculating "how much heat loss is reduced by wrapping insulation" in pipe insulation design, the inside $h$ is first required.
Related Topics
ThermalGnielinski相関式ThermalSieder-Tate相関式Thermal入口助走区間Thermal乱流熱伝達用語集Dittus-Boelter式 — CAE用語解説Thermal混合対流(強制+自然)
この記事の評価ご回答ありがとうございます!参考に
なったもっと
詳しく誤りを
報告
なった
詳しく
報告