Dittus-Boelter Correlation — Turbulent Pipe Flow Heat Transfer

Category: 熱解析 › 強制対流 | 更新 2026-04-12
Dittus-Boelter correlation for turbulent pipe flow heat transfer - Nusselt vs Reynolds number
図:管内乱流における熱伝達の概念図。流速・管径・流体物性からNu数を推算するDittus-Boelter式の位置づけ。

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:

$$ \mathrm{Nu}_D = 0.023\;\mathrm{Re}_D^{\,0.8}\;\mathrm{Pr}^{\,n} $$

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:

$$ \mathrm{Nu}_D = 0.023\;\mathrm{Re}_D^{\,0.8}\;\mathrm{Pr}^{\,n} $$

The definitions of each dimensionless number are as follows:

$$ \mathrm{Re}_D = \frac{\rho\,u\,D}{\mu} = \frac{u\,D}{\nu}, \qquad \mathrm{Pr} = \frac{\mu\,c_p}{k} = \frac{\nu}{\alpha} $$

Conversion to heat transfer coefficient:

$$ h = \frac{\mathrm{Nu}_D \cdot k}{D} $$
SymbolPhysical QuantitySI Unit
$D$Pipe inner diameterm
$u$Average flow velocity in pipem/s
$\rho$Fluid densitykg/m³
$\mu$Dynamic viscosity (viscosity coefficient)Pa·s
$c_p$Specific heat at constant pressureJ/(kg·K)
$k$Fluid thermal conductivityW/(m·K)
$h$Convective heat transfer coefficientW/(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:

ConditionRangeConsequence 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 temperatureModerate (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

Coffee Break Trivia Corner

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:

  1. Determine film temperature: $T_f = (T_w + T_b)/2$ (average of wall temperature $T_w$ and bulk temperature $T_b$)
  2. Obtain property values: Read $\rho$, $\mu$, $c_p$, $k$ at $T_f$ from property tables
  3. Calculate Re: $\mathrm{Re}_D = \rho u D / \mu$
  4. Calculate Pr: $\mathrm{Pr} = \mu c_p / k$
  5. Check applicability: Re > 10,000 and 0.7 < Pr < 160 and L/D > 60
  6. Calculate Nu: $\mathrm{Nu}_D = 0.023\,\mathrm{Re}_D^{0.8}\,\mathrm{Pr}^n$
  7. 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:

  1. 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.
  2. 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:

ItemDittus-Boelter EquationGnielinski 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,0002,300 – 5×10⁶
Pr Applicable Range0.7 – 1600.5 – 2,000
Entrance EffectNot considered (assumes L/D > 60)Can be addressed with correction term
Accuracy (vs. experiment)±25%±10%
Calculation EffortImmediate calculation with a calculatorAdditional calculation of friction factor $f$ required
Optimal UseInitial estimation, sensitivity analysisDetailed 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

FieldTypical ConditionsEstimated $h$ [W/(m²·K)]Notes
Automotive Radiator PipesWater, $u$ = 1–2 m/s, $D$ = 15–25 mm5,000–15,000Used around 80°C
Chemical Plant Heat ExchangersOrganic solvents, $u$ = 0.5–3 m/s500–5,000Relatively high h due to large Pr
Air Conditioning DuctsAir, $u$ = 5–15 m/s, $D$ = 100–300 mm30–100Low h due to low k of air
Data Center Liquid CoolingInert liquids, $u$ = 0.5–2 m/s2,000–8,000Caution with L/D in microchannels
Gas Turbine Cooling PassagesAir, $u$ = 30–100 m/s200–1,000High Re, rotation effects require correction
Coffee Break Trivia Corner

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

関連シミュレーター

この分野のインタラクティブシミュレーターで理論を体感しよう

シミュレーター一覧
この記事の評価
ご回答ありがとうございます!
参考に
なった
もっと
詳しく
誤りを
報告
参考になった
0
もっと詳しく
0
誤りを報告
0
Written by NovaSolver Contributors
Anonymous Engineers & AI — サイトマップ