Nusselt 相関 シミュレーター 戻る
熱伝達シミュレーター

Nusselt 相関 シミュレーター — 管内強制対流の Dittus-Boelter 式

管内乱流の Dittus-Boelter 相関 Nu = 0.023 Re^0.8 Pr^n を実時間に計算します。Reynolds 数・Prandtl 数・熱伝導率・管径から Nusselt 数、熱伝達率 h、熱流束、熱境界層厚さを求め、管内流れの可視化と Re-Nu 線図で強制対流の物理を直感的に理解できます。

パラメータ設定
流体プリセット
Reynolds 数 Re(流速)
-
Re<2300 は層流(Nu≈3.66 一定)、Re>4000 は乱流(Dittus-Boelter)。
Prandtl 数 Pr
-
熱伝導率 k
W/(m·K)
管径 D
mm
壁温 T_w
°C
熱伝達モード
アニメ速度
×
シナリオ

既定値 (Re=50000、Pr=7.0、k=0.60 W/(m·K)、D=25 mm、加熱) では Nu ≈ 287.7、h ≈ 6905 W/(m²·K)、熱流束 (ΔT=10 K) ≈ 69.0 kW/m²、熱境界層厚さ ≈ 0.087 mm。適用範囲は Re > 4000、0.6 < Pr < 160 で、Re < 4000 では層流域となり別の相関式が必要です。

計算結果(ライブ)
Nusselt 数 Nu
熱伝達率 h
熱流束 q
熱境界層厚さ δ_T
流れ領域
Reynolds 数 Re
Prandtl 数 Pr
モード / 流体
管内流れ・熱境界層・熱輸送(物理本体)

加熱された管壁(赤)から流体へ熱が伝わる様子をライブ表示。粒子は壁近傍で熱を受け取り(赤く)下流へ運びます。層流では流れが整い熱境界層が厚く Nu が低く、乱流では流れが混合し境界層が薄く Nu が高くなります。Re(流速)を上げると領域が遷移し、境界層が薄くなって Nu・h が上がる連鎖が一目で分かります。

Dittus-Boelter 線図 Re-Nu (両対数)

横軸:Reynolds 数 Re (4000〜200000、log10 表示) / 縦軸:Nusselt 数 Nu (log10 表示) / 青実線:Dittus-Boelter 直線 (傾き 0.8) / 黄色マーカー:現在の (Re, Nu) 動作点 / Pr を変えると線全体が上下にシフトし、Pr が大きいほど Nu が大きく対流伝熱が強くなります。

理論・主要公式

Dittus-Boelter 相関:管内乱流 (Re > 4000、0.6 < Pr < 160) における強制対流熱伝達の最も基本的な経験式です。

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

$n = 0.4$ は加熱中 (壁温 > 流体温度、流体が熱を受け取る)、$n = 0.3$ は冷却中 (壁温 < 流体温度)。Reynolds 数と Prandtl 数は次のように定義されます:

$$\mathrm{Re} = \frac{\rho U D}{\mu},\quad \mathrm{Pr} = \frac{\mu c_p}{k}$$

熱伝達率 $h$ は Nu の定義式 Nu = hD/k から:

$$h = \frac{\mathrm{Nu}\,k}{D}$$

壁面熱流束 $q$ は Newton の冷却法則から、熱境界層厚さは $\delta_T \approx D/\mathrm{Nu}$:

$$q = h\,(T_w - T_b),\quad \delta_T \approx \frac{D}{\mathrm{Nu}}$$

$\rho$ は密度、$U$ は平均流速、$D$ は管径、$\mu$ は粘性、$c_p$ は定圧比熱、$k$ は流体熱伝導率、$T_w$ は壁面温度、$T_b$ はバルク (混合平均) 温度です。

Nusselt 相関 シミュレーターとは

🙋
熱交換器の設計で「Dittus-Boelter 式」って必ず出てくるんですけど、結局 Nu = 0.023 Re^0.8 Pr^0.4 の数字 0.023 とか 0.8 はどこから来てるんですか?理屈ですか実験ですか?
🎓
いい着眼点だ。Dittus-Boelter は純粋に経験式で、1930年に Dittus と Boelter が水・空気・油の管内乱流データを大量にフィッティングして作った。理屈で導出される指数 0.8 は「乱流の境界層理論」(普遍速度分布の1/7乗則) からほぼ予測でき、Pr^0.4 (加熱) や Pr^0.3 (冷却) は温度差による粘性変化の経験的補正だ。本ツール既定値 (Re=50000、Pr=7.0、水、D=25 mm、加熱) で実行すると Nu ≈ 287.7、熱伝達率 h ≈ 6905 W/(m²·K) と表示されるはずだよ。これが「水を流したパイプで壁から熱を取る」典型値だ。
🙋
h = 6905 W/(m²·K) ってピンと来ないんですけど、実用的にはどのくらいスゴい値なんですか?
🎓
ざっくり感覚を持つには対比が一番だ。自然対流の空気で h ≈ 5〜25、強制対流の空気で 25〜250、強制対流の水 (乱流) で 1000〜15000、沸騰水で 10^4〜10^5 という階層がある。だから h=6905 は「水を勢いよく流したパイプ」の典型値で、空気の数十倍だ。実務感覚としては、ΔT=10 K 取れれば q = h·ΔT = 69 kW/m² の熱を運べる — これは家庭用 IH コンロ (3 kW) を 1 m² の管壁で 20 台分逃がす規模だ。本ツールでスライダーを動かして、Re を 4000 に下げると h が約 915 W/(m²·K) まで落ちる。流速次第で 7〜8 倍の幅があると覚えておけ。
🙋
「熱境界層」のグラフで赤い帯がすごく薄くなったり厚くなったりするんですが、これは何を表してるんですか?
🎓
よく気づいた。熱境界層 δ_T は壁面から流体本体に向かう「温度が変化する薄い層」で、本ツールでは赤グラデーションで表示している。Nu の物理意味は実は δ_T ≈ D/Nu — つまり Nu が大きいほど境界層が薄くなる。既定値 Nu=287.7、D=25 mm なら δ_T ≈ 0.087 mm(87 μm)と髪の毛より細い。Re を上げる (流速 ↑) と乱流ミキシングが激しくなり境界層が薄くなる→Nu↑→h↑、という連鎖だ。試しに Re を 4000→200000 までスイープボタンで動かしてみて。境界層の赤帯が薄くなる様子と、右の Re-Nu 線図の黄色点が動く様子が同期するはずだ。
🙋
Prandtl 数を 7 から 1 に下げると Nu がガクンと落ちますね。Pr って結局何ですか?水と空気でそんなに違うんでしょうか?
🎓
大事な疑問だ。Prandtl 数 Pr = μc_p/k は「運動量拡散と熱拡散の比」で、流体ごとに固有の値だ。水で Pr ≈ 7 (20°C)、空気で Pr ≈ 0.71、エチレングリコールで Pr ≈ 200、液体ナトリウムで Pr ≈ 0.005、エンジンオイル (低温) で Pr ≈ 10000、と桁違いに違う。Pr が大きい液体ほど「熱境界層が速度境界層より薄い」ので、壁近傍の温度勾配が急で h が大きくなる。本ツールで Pr=0.71 (空気) にすると Nu が水の約 0.39 倍に落ちる — 同じ Re でも空気は水より遥かに伝熱が弱い、これが空冷より水冷が圧倒的に有利な理由だ。液体金属 (Pr<<1) は別の相関 (Lyon 式) が必要だが、本ツールの最小値 Pr=0.5 は Dittus-Boelter の適用範囲内ぎりぎりだ。

よくある質問

Dittus-Boelter 相関は、管内乱流における強制対流熱伝達の最も基本的な経験式で、Nu = 0.023 Re^0.8 Pr^n (加熱で n=0.4、冷却で n=0.3) と表されます。Re は Reynolds 数 (ρUD/μ)、Pr は Prandtl 数 (μc_p/k) です。適用範囲は Re > 4000 (乱流域) かつ 0.6 < Pr < 160、L/D > 10 (発達流れ)。本ツール既定値 (Re=50000、Pr=7.0、水、D=25 mm) では Nu ≈ 287.7、熱伝達率 h ≈ 6905 W/(m²·K) と表示されます。
Nusselt 数 Nu = hD/k は対流熱伝達と流体内純粋熱伝導の比率を表す無次元数で、Nu=1 なら対流効果はなく純粋熱伝導と同じ、Nu>>1 なら強い対流ミキシングで伝熱が増強されていることを意味します。同時に Nu ≈ D/δ_T (D は管径、δ_T は熱境界層厚さ) と解釈でき、熱境界層が管径より十分薄いほど Nu は大きくなります。本ツールでは熱境界層厚さも併記しており、Nu=288 で約 0.087 mm と非常に薄いことが確認できます。
Dittus-Boelter 式の Pr^n の指数が加熱で n=0.4、冷却で n=0.3 と異なるのは、温度変化に伴う流体粘性の変化を経験的に補正するためです。液体は温度上昇で粘性が下がるため、加熱時は壁近傍の薄い高温層で粘性が低下し境界層が薄くなって伝熱が増強 (指数を大きく)、冷却時は逆に粘性が増加して境界層が厚くなり伝熱が低下 (指数を小さく) します。物性変化を別途補正する Sieder-Tate 式 (μ/μ_w)^0.14 と等価な簡易補正と理解できます。
Dittus-Boelter 式は (1) Re < 4000 (層流・遷移流)、(2) Pr < 0.6 (液体金属) または Pr > 160 (高粘性油)、(3) L/D < 10 (入口効果支配)、(4) 大きな温度差で物性が大きく変わる場合、(5) 環状管・コイル管などの非単純形状、には不適切です。これらでは Gnielinski 式 (Re=3000〜5×10^6 で精度高い)、Sieder-Tate 式 (粘性比補正)、Hausen 式 (層流入口流)、液体金属用の Lyon 式などを使い分けます。本ツールでは適用範囲を守るため Re スライダーの下限を 4000 に制限しています。

実世界での応用

ボイラ・蒸気発生器の伝熱面設計:火力発電所のボイラ管 (D=50 mm) では給水を 200°C→300°C に加熱しますが、Re≈100000、Pr≈1 程度になり、本ツールでこの値(k=0.68 W/(m·K)、D=50 mm)を入れると Nu ≈ 230、h ≈ 3100 W/(m²·K) と概算できます。実際の燃焼ガス側 h との総括 U を考えて A = Q/(U·LMTD) で必要管長を決めますが、Dittus-Boelter の精度 (±25%) を意識した余裕設計が業界標準です。発展形の Gnielinski 式 (誤差 ±10%) を使う場合もあります。

原子力発電プラント一次冷却管:PWR の一次冷却水 (300°C、加圧水) を高熱流束で炉心から取り出すため、Re ≈ 5×10^5、Pr ≈ 0.9 で h ≈ 30000 W/(m²·K) 級の強力な対流伝熱を実現します。Dittus-Boelter 式では Nu ≈ 800 となり、D=12 mm なら δ_T ≈ 15 μm という極薄境界層が形成されます(この Re は本ツールのスライダー上限 200,000 を超えるため式での確認値です)。CHF (限界熱流束) 設計では Dittus-Boelter の予測値を出発点に、Tong 式や Bowring 式で安全評価します。

自動車エンジン冷却水路:シリンダブロックのウォータージャケット (D≈8 mm) では Re≈20000、Pr≈4、k≈0.33 W/(m·K) (50% 水/グリコール混合液、80°C)、本ツールで h ≈ 4600 W/(m²·K) と推算できます。燃焼室壁面から 200 kW/m² の熱流束を逃がすため、ΔT ≈ 45 K の駆動力で運転されます。冷却水温が上がる (Pr が下がる) と h も低下するため、低速走行時の冷却劣化を考慮した設計が必要です。

HVAC 冷温水コイル:ビル空調機の銅管コイル (D=12 mm) を流れる冷温水 (Re≈10000、Pr≈7、k=0.6 W/(m·K)) で h ≈ 4000 W/(m²·K) と本ツールで推算でき、フィン側の空気対流 h≈50 W/(m²·K) と組み合わせて UA を最適化します。水側の流量を上げる (Re↑) と h は上がりますがポンプ動力も増えるため、伝熱と動力のバランス点を Dittus-Boelter で評価して、フィン側 (律速側) 強化との優先順位を決めます。

よくある誤解と注意点

最も多い誤解が 「Dittus-Boelter は全ての管内流れに使える万能式」 というものです。実際の適用条件は厳しく、Re > 4000 (乱流)、0.6 < Pr < 160、L/D > 10 (発達流れ)、温度差による物性変化が小さい、の全てを満たす必要があります。たとえばエンジンオイル (Pr>1000) や液体ナトリウム (Pr<0.01) では誤差が 100% を超え、それぞれ Sieder-Tate 式や Lyon 式に切り替える必要があります。本ツールで Pr スライダーを 0.5 や 100 の両端に動かすと、適用範囲外であることを意識してください。

次に多いのが 「Nu が大きいほど熱量 Q が大きい」 という単純な理解です。実際には Q = h·A·ΔT = (Nu·k/D)·A·ΔT なので、Nu を上げる手段によって結果が変わります。流速を上げて Re↑ で Nu↑ させればポンプ動力も増え、管径を小さくして D↓ させると同じ Nu でも h は上がりますが圧力損失が二乗で増えます。実務では「熱伝達 vs 圧力損失」のトレードオフ (Reynolds 比較・Colburn の j 因子など) で最適点を探します。本ツールで D を 5 mm に下げると h が劇的に上がる様子と、それが実は危険である (圧損が深刻に増える) ことを忘れないでください。

最後に 「Dittus-Boelter は最新の標準式」 という誤解です。実は 1976 年に提案された Gnielinski 式 Nu = (f/8)(Re-1000)Pr / [1+12.7(f/8)^0.5 (Pr^(2/3)-1)] のほうが現代の標準で、誤差 ±10% (Dittus-Boelter は ±25%) と精度が高く、適用範囲も 3000 < Re < 5×10^6、0.5 < Pr < 2000 と広がっています。ただし計算が複雑なので、概算や教育用途では今でも Dittus-Boelter が広く使われます。本ツールは Dittus-Boelter に特化していますが、実際の設計では Gnielinski や設備固有の補正係数を加える慣習があることを覚えておいてください。

使い方ガイド

  1. Reynolds 数(Re)スライダーを 10,000~100,000 の範囲で調整し、管内流速に対応させます。工業用水冷却では Re = 50,000 が標準です
  2. Prandtl 数(Pr)を設定します。水の場合 Pr ≈ 6.5、エチレングリコールの場合 Pr ≈ 25 を目安に選択してください
  3. 熱伝導率 k(W/m·K)と管内径 D(mm)を入力し、Dittus-Boelter 式 Nu = 0.023 Re^0.8 Pr^0.4 を実行して Nusselt 数を計算します
  4. 熱伝達率 h = Nu·k / D から管壁での対流熱伝達係数を導出し、熱流束 q = h·ΔT(ΔT = 10 K)を求めます

具体的な計算例

鋼製熱交換器で水が流通する場合を想定します。Re = 50,000、Pr = 6.5、k = 0.65 W/m·K、D = 12 mm の条件で計算すると、Nu = 0.023 × 50,000^0.8 × 6.5^0.4 ≈ 279 となります。これより h = 279 × 0.65 / 0.012 ≈ 15,100 W/m²·K、熱流束 q = 15,100 × 10 ≈ 151 kW/m² を得ます。熱境界層厚さ δ_th ≈ D / Nu ≈ 0.043 mm となり、管径に対し約 0.4% の極めて薄い層が形成されることが確認できます

実務での注意点

  1. 厳密な Dittus-Boelter の適用範囲は Re ≳ 10,000 とされ、本ツール下限の Re = 4,000~10,000 の遷移域では精度が落ちます。層流領域(Re < 2,300)では Hausen 式や Shah 式を使用してください
  2. 管壁温度が 80℃ を超える場合、流体の物性値(特に Pr、k)が変動するため、スライダーで動的に修正して精度を確保します
  3. 塩水や高粘度油など特殊流体では Pr が大幅に変わるため、メーカー仕様書から正確な値を入力することが必須です
  4. 管内径が 6 mm 以下の微細管では圧力損失が急増するため、Re = 30,000 以下に制限する設計が推奨されます