Dittus-Boelter相関式 — 管内乱流熱伝達の基本式
理論と物理
概要 — Dittus-Boelter式とは何か
先生、Dittus-Boelter式ってどこで使うんですか？ 教科書のどの章にも出てくるんですけど、正直ピンときてなくて…
管内乱流の熱伝達係数 $h$ を最も簡便に推算できる式だ。式の形はこう：
流体が加熱されるときは $n = 0.4$、冷却されるときは $n = 0.3$ を使う。冷却水配管、ヒートシンクのフィン間流路、化学プラントの熱交換器——こういった「管の中を流体が乱流で流れる」場面の初期設計で真っ先に使う式だよ。
え、たったこれだけの式で熱伝達係数が分かるんですか？ Re と Pr を入れるだけ？
そう、だからこそ「初期概算の王様」なんだ。例えば自動車のラジエーター配管を設計するとき、冷却水の流速と管の内径が決まれば、Re を計算して、水の Pr（常温で約7）を入れて、すぐに $h$ の概算値が得られる。これで「この配管径で十分冷やせるか？」の最初の判断ができるわけだ。
でも「概算」って言ってるってことは、当てにならない場面もあるってことですよね？
鋭いね。この式には前提条件がある。Re > 10,000 の完全発達乱流で、管が十分長い（$L/D > 60$）場合に限って使える。入口効果が強い短い管や、Re = 2,300〜10,000 の遷移域では Gnielinski式の方がずっと正確だ。この「使える範囲」を知らないまま適用すると、2〜3割の誤差を平気で抱え込むことになる。
核心の数式
Dittus-Boelter式を構成する2つの式を整理する。
ヌセルト数の相関式：
ここで各無次元数の定義は以下の通り：
熱伝達係数への変換：
| 記号 | 物理量 | SI単位 |
|---|---|---|
| $D$ | 管の内径 | m |
| $u$ | 管内平均流速 | m/s |
| $\rho$ | 流体密度 | kg/m³ |
| $\mu$ | 動粘度（粘性係数） | Pa·s |
| $c_p$ | 定圧比熱 | J/(kg·K) |
| $k$ | 流体の熱伝導率 | W/(m·K) |
| $h$ | 対流熱伝達係数 | W/(m²·K) |
| $n$ | Pr指数（加熱0.4／冷却0.3） | — |
適用条件と限界
適用条件、もうちょっと具体的に教えてもらえますか？ 「Re > 10,000」以外にも制限があるんですよね？
うん、まとめるとこうなる：
| 条件 | 範囲 | 違反した場合 |
|---|---|---|
| レイノルズ数 | $\mathrm{Re}_D > 10{,}000$ | 遷移域・層流では過大評価。Gnielinski式を使う |
| プラントル数 | $0.7 < \mathrm{Pr} < 160$ | 液体金属（Pr≪1）や高粘度油（Pr≫160）では大きな誤差 |
| 管長/直径比 | $L/D > 60$ | 入口効果で実際のhが高くなる。入口補正係数が必要 |
| 壁面温度と流体温度の差 | 中程度（$\Delta T$ が小） | 大温度差では粘度変化が大→Sieder-Tate式を使う |
$L/D > 60$ ってことは、内径20mmの管なら1.2m以上ないとダメってことですか？
そういうこと。電子機器のコールドプレートや小型ヒートシンクのフィン間流路だと $L/D = 10〜30$ 程度のことが多い。そういう短管では入口付近で熱伝達が局所的に非常に高くなるから、Dittus-Boelter式で均一な $h$ を仮定すると過小評価になる。実務では入口補正係数 $(1 + (D/L)^{0.7})$ を掛けるか、Gnielinski式に切り替えるのが安全だ。
各項の物理的意味
Re0.8 — なぜ0.8乗なのか
レイノルズ数は「慣性力/粘性力」の比であり、乱流の強さの指標。指数が1ではなく0.8なのは、乱流域でのNu数がReに対して線形ではなく、やや緩やかに増加するため。これは乱流境界層の発達パターン——壁面近傍の粘性底層が流速増加とともに薄くなるが完全には消えない——に由来する。実験データとの最良一致が0.8付近になることが多くの実験で確認されている。
Prn — 加熱と冷却でなぜ指数が違うのか
プラントル数は「運動量拡散/熱拡散」の比。流体が加熱されると壁面近傍の粘度が低下し速度境界層が薄くなる→熱伝達促進→大きめの指数（0.4）。冷却されると壁面近傍の粘度が増加し速度境界層が厚くなる→熱伝達抑制→小さめの指数（0.3）。この非対称性は、温度に依存する粘度変化の影響を簡易的に取り込んだものだ。
係数0.023の意味
この係数は理論的に導出されたものではなく、大量の円管流れの実験データへの回帰によって決定された経験的定数。Winterton (1998) の文献によれば、原著論文の実験データは±25%程度のばらつきがあり、0.023はその中心値である。
歴史的背景
式の真の発案者問題
「Dittus-Boelter式」として知られる $\mathrm{Nu} = 0.023\,\mathrm{Re}^{0.8}\,\mathrm{Pr}^n$ は、1930年にF.W. DittusとL.M.K. Boelterが発表した論文に由来する。しかし近年の文献調査（Winterton, 1998）により、指数 $n$ の加熱・冷却使い分け（0.4/0.3）はBoelterが1942年に単独で修正したものであり、Dittusの本来の論文には存在しないことが判明している。さらに元論文は正式な査読付きジャーナルではなくカリフォルニア大学の内部報告書であり、原本の入手自体が困難。にもかかわらず、あまりにも広く引用されているために「Dittus-Boelter式」の名称が定着した——工学史における面白いエピソードだ。
数値解法と実装
熱伝達係数の算出手順
実際に $h$ を求めるまでの手順を、ステップバイステップで教えてください。
手順はこうだ：
- 膜温度の決定：$T_f = (T_w + T_b)/2$（壁面温度 $T_w$ とバルク温度 $T_b$ の平均）
- 物性値の取得：$T_f$ における $\rho$, $\mu$, $c_p$, $k$ を物性表から読む
- Re の計算：$\mathrm{Re}_D = \rho u D / \mu$
- Pr の計算：$\mathrm{Pr} = \mu c_p / k$
- 適用条件の確認：Re > 10,000 かつ 0.7 < Pr < 160 かつ L/D > 60
- Nu の計算：$\mathrm{Nu}_D = 0.023\,\mathrm{Re}_D^{0.8}\,\mathrm{Pr}^n$
- $h$ の算出：$h = \mathrm{Nu}_D \cdot k / D$
膜温度がわからないとき——例えば設計の最初で壁面温度がまだ決まってない場合はどうするんですか？
いい質問だ。最初は壁面温度を仮定して計算し、得られた $h$ を使って壁面温度を再計算→もう一度物性値を更新→$h$ を再計算…という反復を回す。実務的には2〜3回で十分収束する。水のように温度による物性変化が小さい流体なら、バルク温度で物性値を取っても誤差は数%以内だ。
計算例 — 冷却水配管
具体的な数字で見てみたいです。例えばエンジンの冷却水系統とか。
じゃあ典型的なエンジン冷却水の条件でやってみよう。
条件：内径 $D = 0.02$ m、流速 $u = 1.5$ m/s、水温 80°C（加熱状態）
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 $$→ 完全乱流。OK。
Step 2 — Pr：
$$ \mathrm{Pr} = \frac{3.55 \times 10^{-4} \times 4{,}197}{0.670} = 2.22 $$→ 範囲内。OK。
Step 3 — Nu（加熱なので $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)} $$つまり約 $h \approx 10{,}700$ W/(m²·K) だ。自動車の冷却水系統で典型的な値だね。
おお、かなり大きな値ですね。空気の強制対流だとどのくらいですか？
空気だとPr ≈ 0.71、$k$ ≈ 0.03 W/(m·K)と水より桁違いに小さいから、同じ流速・管径でも $h = 50〜200$ W/(m²·K) 程度になる。水の $h$ が空気の50〜100倍あるのは、水の方が熱伝導率と熱容量がはるかに大きいからだ。だからエンジン冷却には水が使われる。
物性値の評価温度
物性値をどの温度で取るか——実務で最も見落とされるポイント
Dittus-Boelter式はもともと壁面温度とバルク温度の差が小さい実験データから導かれた。そのため、物性値は膜温度 $T_f = (T_w + T_b)/2$ で評価するのが標準。ただし壁温が未知の段階ではバルク温度で代用し、$h$ を求めた後に壁温を再計算して物性値を更新する反復法が一般的だ。
壁面と流体の温度差が大きい場合（例：高温の油がパイプ壁面で急冷される場合）、壁面近傍の粘度が大きく変化する。このとき Dittus-Boelter式では対応しきれないため、粘度補正項 $(\mu_b/\mu_w)^{0.14}$ を含む Sieder-Tate式 に切り替えるべきだ。
実践ガイド
CAEでの境界条件設定
CAEソフトで熱解析するとき、Dittus-Boelter式ってどう使うんですか？ 式をソフトに入力するんですか？
使い方は大きく2つある：
- 対流境界条件として $h$ を手動設定：流体側をモデリングせず、固体だけ解析する場合。Dittus-Boelter式で計算した $h$ と流体バルク温度 $T_\infty$ を壁面に指定する。Ansys MechanicalやAbaqusの熱解析で最も一般的なアプローチだ。
- CFDソルバーが内部で自動計算：FlunetやSTAR-CCM+で管内流れを直接解く場合、ソルバーは壁面関数を通じて自動的に $h$ を計算する。この場合はDittus-Boelter式を手動で使う必要はない。むしろ、CFD結果の妥当性検証としてDittus-Boelter式の概算値と比較する使い方になる。
つまり「流体を解かない熱解析」のときにDittus-Boelter式の出番があるってことですね。逆にCFDで直接解くなら検証用？
そのとおり。実務で多いのは「設計の初期段階でDittus-Boelter式による概算→詳細設計段階でCFD」という流れだ。初期概算の $h$ がCFD結果と大きく食い違ったら、適用条件を外れている可能性が高い。
Gnielinski式との使い分け
Gnielinski式との違いを、もう少し詳しく教えてください。どっちを使うか迷うことが多くて…
比較表を見てみよう：
| 項目 | Dittus-Boelter式 | Gnielinski式 |
|---|---|---|
| 式の形 | $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 適用範囲 | > 10,000 | 2,300 〜 5×10⁶ |
| Pr 適用範囲 | 0.7 〜 160 | 0.5 〜 2,000 |
| 入口効果 | 考慮なし（L/D > 60 前提） | 補正項で対応可 |
| 精度（実験との一致） | ±25% | ±10% |
| 計算の手間 | 電卓で即計算 | 摩擦係数 $f$ の計算が追加で必要 |
| 最適な用途 | 初期概算、感度分析 | 詳細設計、報告書の根拠 |
じゃあ精度が必要ならGnielinski一択で、Dittus-Boelterは「ざっくり見たいとき」に使う感じですか？
基本的にはそう。ただしRe > 10,000の完全乱流域では両者の差は数%以内に収まることが多い。つまり完全発達乱流で長い管ならDittus-Boelterで十分。逆にGnielinski式が真価を発揮するのは、遷移域（Re = 2,300〜10,000）や入口効果が大きい短管だ。
Sieder-Tate式との使い分け
Sieder-Tate式との違いも知りたいです。
Sieder-Tate式は $\mathrm{Nu} = 0.027\,\mathrm{Re}^{0.8}\,\mathrm{Pr}^{1/3}\,(\mu_b/\mu_w)^{0.14}$ という形で、壁面と流体のバルクの粘度比 $(\mu_b/\mu_w)^{0.14}$ が入っている。
つまり壁面温度と流体温度の差が大きく、粘度が大きく変わる場合に使う式だ。例えば原油パイプラインで壁面がヒーターで加熱されている場合、壁面近傍では粘度が激減する。こういう場面ではDittus-Boelter式では不十分で、Sieder-Tate式の粘度補正が効いてくる。
逆に水や空気のように温度変化による粘度変化が小さい流体では、Dittus-Boelter式とSieder-Tate式の結果はほぼ一致する。
産業応用事例
| 分野 | 典型的な条件 | 推算される $h$ [W/(m²·K)] | 備考 |
|---|---|---|---|
| 自動車ラジエーター配管 | 水、$u$ = 1〜2 m/s、$D$ = 15〜25 mm | 5,000〜15,000 | 80°C近傍で使用 |
| 化学プラント熱交換器 | 有機溶剤、$u$ = 0.5〜3 m/s | 500〜5,000 | Prが大きいため相対的にhが高い |
| 空調ダクト | 空気、$u$ = 5〜15 m/s、$D$ = 100〜300 mm | 30〜100 | 空気のkが小さいためhも小さい |
| データセンター液冷 | 不活性液体、$u$ = 0.5〜2 m/s | 2,000〜8,000 | マイクロチャネルではL/D要注意 |
| ガスタービン冷却通路 | 空気、$u$ = 30〜100 m/s | 200〜1,000 | 高Re・回転効果で補正が必要 |
給湯配管設計での常用
日本の給湯システム設計（JIS B 8417準拠）では管内強制対流の概算にDittus-Boelter式が広く使われ、流速 0.5〜2 m/s の温水配管で $h$ = 3,000〜8,000 W/(m²·K) 程度が得られる。配管の保温設計において「保温材を巻くことでどの程度放熱が減るか」を計算する際に、管内側 $h$ がまず必要になる。
ソフトウェアでの扱い
CFD直接解析 vs 相関式
CFDがあるなら、もう相関式なんて使わなくていいんじゃないですか？
それは「GPS があるなら地図は不要か？」みたいな問いだね。答えはNoだ。理由は3つ：
- 初期設計の速度：CFDは1ケースに数時間〜数日。Dittus-Boelter式なら数秒。設計パラメータを100通り振る感度分析はCFDでは非現実的
- V&V（検証と妥当性確認）：CFD結果が正しいことを確認するために、独立した手法（＝相関式）との照合が不可欠
- システムレベル解析：1D配管ネットワーク解析（FloMasterなど）では、各管セグメントのhをDittus-Boelter式で自動計算している
主要ツールでの実装
| ツール | Dittus-Boelter式の扱い |
|---|---|
| Ansys Fluent / CFX | 壁面関数モデルの内部でNu相関式を参照。ユーザーは直接入力しない。結果のNu値とDittus-Boelter値を比較して妥当性検証に使用 |
| STAR-CCM+ | 同上。Report機能でNu, hを出力し、相関式と比較する運用が推奨 |
| Ansys Mechanical / Abaqus | 熱解析の対流境界条件として $h$ を手動入力。Dittus-Boelter式で算出した値を使うのが一般的 |
| COMSOL Multiphysics | Pipe Flowモジュールに内蔵。管路解析で自動計算される |
| FloMASTER / Flowmaster | 1D配管ネットワーク解析で内部的にDittus-Boelterまたはgnielinski式を使用 |
| OpenFOAM | nutWallFunction系の壁面関数で類似の相関を使用。カスタマイズ可能 |
| AspenTech EDR | 熱交換器設計でDittus-Boelter/Gnielinski/Sieder-Tateの3式を自動選択 |
先端トピック
機械学習による相関式の改良
最近は機械学習で相関式を作り直す研究があるって聞いたんですけど、Dittus-Boelter式は時代遅れになるんですか？
確かにニューラルネットワークやGaussian Process回帰を使って、数千点の実験・CFDデータから新しい相関式を導く研究は増えている。例えばナノ流体（ナノ粒子を混ぜた冷却液）の熱伝達では、Dittus-Boelter式では予測できない粒子濃度の効果をMLモデルで取り込む試みがある。
ただし、Dittus-Boelter式が「時代遅れ」になることはないと思う。シンプルで物理的意味が明確な式は、初期設計・教育・コミュニケーションの場で不可欠だ。「MLモデルが $h$ = 12,345 と言ってます」では上司は納得しないが、「Dittus-Boelter式で約10,000、MLで12,000、差異はナノ粒子の効果」と言えば筋が通る。
マイクロチャネルへの拡張
半導体冷却でマイクロチャネルが注目されてますよね。Dittus-Boelter式は使えるんですか？
マイクロチャネル（水力直径 $D_h < 1$ mm）では、いくつかの追加効果が出てくる：
- 入口効果が支配的：$L/D$ が数十程度しかなく、完全発達流を仮定するDittus-Boelter式は不適切
- 軸方向熱伝導：管壁の軸方向熱伝導がバルク温度分布に影響する
- 表面粗さ：マクロ管では無視できた表面粗さが水力直径に対して相対的に大きくなる
これらの理由から、マイクロチャネルでは修正Gnielinski式や、Shah-Londonのエントリーレングス相関を使うのが標準だ。ただし「オーダー感をつかむ」目的でDittus-Boelter式を使うのは今でも有用で、$h$ が数万〜数十万 W/(m²·K) という桁の大きさを把握するのに役立つ。
トラブルシューティング
よくある間違いと対策
先生、Dittus-Boelter式を使うとき、初心者がやりがちな失敗パターンってありますか？
| 間違い | 結果への影響 | 対策 |
|---|---|---|
| 遷移域（Re = 2,300〜10,000）で使う | Nuを過大評価 → hが実際の1.5〜3倍 | Gnielinski式に切り替える |
| 加熱/冷却の $n$ を間違える | 5〜15%の系統誤差 | 流体が熱を受ける＝加熱（0.4）、放出＝冷却（0.3） |
| 物性値を常温で取る | 高温域では粘度が大幅に異なる → Re が数倍ずれる | 膜温度 $T_f$ で物性値を取得 |
| 短管（$L/D < 60$）に適用 | 入口効果を無視 → hを過小評価 | 入口補正係数を付加、またはGnielinski式 |
| 非円管断面にそのまま適用 | 矩形・環状断面では形状係数が異なる | 水力直径 $D_h = 4A/P$ で置き換え、さらに形状補正が必要 |
| 単位系の不整合（cm/m混在） | Re が2桁ずれる → 完全に誤った $h$ | すべてSI単位に統一してから計算 |
結果の妥当性チェック
計算した $h$ が合ってるかどうか、どうやって確認すればいいですか？
3つのチェックポイントがある：
- オーダーの確認：水の強制対流なら $h$ = 1,000〜20,000 W/(m²·K)。空気なら 20〜200。これを大幅に外れたら入力データを見直す
- 別の式との照合：同じ条件でGnielinski式でも計算してみる。完全乱流域なら両者の差は10%以内のはず。30%以上ずれたら何かがおかしい
- エネルギーバランス：得られた $h$ から壁面の放熱量 $q = h \cdot A \cdot (T_w - T_b)$ を計算し、流体の温度上昇 $\Delta T = q / (\dot{m} c_p)$ が妥当な範囲か確認
「$h$ を大きめに見積もった方が安全」は危険な誤解
冷却設計で $h$ を過大評価すると、「十分冷やせる」と判断して小さすぎる冷却系を設計してしまう。実際に運用するとオーバーヒートする——これは安全側ではなく危険側の誤りだ。逆に $h$ を過小評価すると過剰な冷却設備を設計することになり、コストは増えるが少なくとも安全。不確実性があるときは $h$ を小さめに見積もる（安全側に倒す）のが原則だ。
なるほど…Dittus-Boelter式って一見シンプルですけど、使いこなすには適用条件と物性値の扱いをちゃんと理解してないとダメなんですね。
そのとおり。「式を知っている」のと「正しく使える」のは全然違う。特にCAEの境界条件として $h$ を設定するとき、その根拠となる相関式の適用条件を理解していないと、どんなに精密なFEMモデルを作っても入力が間違っていれば結果はゴミだ。"Garbage in, garbage out"の典型例だね。
なった
詳しく
報告