Gnielinski相関式
Gnielinski相関式の理論基礎
概要
先生、Gnielinski式ってDittus-Boelter式より精度がいいって聞いたんですけど、本当ですか？
相関式の適用範囲は式の版と実験条件によって異なります。ここでは後述の研究例と滑らかな管の摩擦係数を組み合わせ、Re = 3,000〜10⁶、Pr = 0.5〜200を参考範囲とします。遷移域では特に流れの状態に注意が必要です。
精度はどれくらい違うんですか？
Gnielinski式とDittus–Boelter式の誤差は、実験条件、物性、摩擦係数や温度条件に依存します。公称誤差を全条件の保証とせず、対象条件で実測や適切な参考データと比較してください。 Gnielinski correlation: experimental study
え、そんなに違うんですか！ じゃあDittus-Boelter式ってもう使わなくていいんですか？
Dittus-Boelter式は概算や手計算の照合に役立ちます。設計では流動状態、管の形状、物性と境界条件を確認し、適用可能な相関式を選んで測定値とも比較してください。
Gnielinski式の導出背景
Gnielinski式ってどうやってできたんですか？ いきなり天才が思いついた？
いい質問だ。歴史的な流れがある。まず1970年にPetukhov（ペトゥホフ）が、運動量輸送と熱輸送のアナロジーに基づいた乱流熱伝達モデルを提案した。ただしPetukhov式はRe > 104 の完全乱流域限定だった。
1976年にGnielinski（グニエリンスキ）がPetukhov式を改良して、Reの項を (Re − 1000) に置き換えるという修正を加えた。たったこれだけの修正で遷移域まで精度よく予測できるようになったんだ。シンプルだけど物理的に深い修正で、「乱流が完全には発達していない領域では実効的なReが小さくなる」という現象をうまく捉えている。
Re を (Re − 1000) に変えるだけで遷移域もカバーできるって、なんか意外にシンプルですね。
そう、優れた工学モデルはシンプルなことが多い。Gnielinski本人はドイツのカールスルーエ工科大学で膨大な実験データと比較検証して、この修正が最も広範囲に高精度を保つことを示した。今では「内部流の乱流熱伝達と言えばGnielinski式」というくらいの定番になっている。
支配方程式
じゃあ実際の式を教えてください！
Gnielinski相関式は以下の通りだ:
ここで $\mathrm{Nu}$ はヌセルト数、$\mathrm{Re}$ はレイノルズ数、$\mathrm{Pr}$ はプラントル数、$f$ はDarcy摩擦係数だ。
分母の $12.7\sqrt{f/8}(\mathrm{Pr}^{2/3} - 1)$ って何を表してるんですか？
分母はPr数と摩擦係数の影響を含みます。Prが異なる流体のNuを比較するにはReや境界条件などを揃える必要があります。熱伝達係数hも熱伝導率と管径に依存します。
なるほど、Dittus-Boelter式のNu = 0.023 Re0.8 Prn と比べると、分母でPrの影響を補正してるから精度が良いわけですね！
両式はReとPrの扱いが異なります。計算結果に差があっても、それだけでどちらが高精度かは決まりません。流体と流動条件に対する適用範囲を確認し、実測値と比べてください。
Petukhov摩擦係数
式の中に出てくる摩擦係数 $f$ はどうやって求めるんですか？
Gnielinski式とセットで使うのが Petukhov（ペトゥホフ）の摩擦係数 だ:
これは滑らかな円管における乱流のDarcy摩擦係数を与える式で、Moody線図の滑面（smooth pipe）曲線とほぼ一致する。Re = 3,000〜5×106 で有効だ。
Moody線図って配管設計で使うやつですよね。あれの滑面カーブを式にしたものってことですか？
そう。ただし粗面（relative roughness ε/D > 0）の場合はPetukhov式のままでは使えない。粗面ではColebrook-White式やChurchill式から $f$ を求めて、それをGnielinski式に入れることになる。ただしその場合は「Gnielinski式の精度保証範囲外」だということは理解しておく必要がある。
なるほど、摩擦係数の選び方がNuの精度に直結するんですね。
適用条件と限界
Gnielinski式はどんな条件でも使えるんですか？
いや、しっかり適用条件がある。これを外すと精度が保証されない:
| パラメータ | 適用範囲 | 備考 |
|---|---|---|
| レイノルズ数 $\mathrm{Re}$ | 3,000 〜 106 | 遷移域〜乱流域 |
| プラントル数 $\mathrm{Pr}$ | 0.5 〜 200 | 気体〜高粘度液体 |
| 長さ/直径比 $L/D$ | ≧ 10 | 参考目安。十分な発達を保証しない |
| 管壁条件 | 滑面（smooth） | 粗面は別途補正が必要 |
Pr = 0.5〜2,000ってことは、液体金属（Na、Pr ≈ 0.005）には使えないんですね？
液体金属のような低Pr数流体や、表の範囲外の条件には、この例の相関式をそのまま使えません。対象流体・境界条件に適した相関式を文献で確認してください。
Dittus-Boelterとの精度比較
具体的にどれくらい差が出るか、数字で見たいです！
例えば水が円管内を流れるケースで比較してみよう。条件: D = 25 mm、Tbulk = 50°C（Pr ≈ 3.56）。
計算例では Pr = 3.56、滑らかな円管のDarcy摩擦係数 f = (0.790 ln Re − 1.64)⁻² を使用します。Dittus–Boelter式は流体加熱の指数 n = 0.4 とし、差分は (D-B − Gnielinski) / Gnielinski × 100% です。Re = 10,000は適用下限付近の参考値です。
| Re | Nu (Gnielinski) | Nu (Dittus-Boelter) | 差分 |
|---|---|---|---|
| 3,000（遷移域） | 17.8 | 適用外 | D-B不可 |
| 10,000 | 61.2 | 60.6 | -1.0% |
| 50,000 | 245.0 | 219.5 | -10.4% |
| 100,000 | 439.9 | 382.2 | -13.1% |
| 500,000 | 1,730.4 | 1,385.1 | -20.0% |
この条件ではDittus–Boelter式はGnielinski式より小さい値になります。式同士の差だけから、どちらが実測に近いかや安全側かは判断できません。遷移域は流れの状態に依存するため、適用範囲と実測条件を確認してください。
条件と式を明記すれば、比較結果を再計算できますね。
そのとおり。摩擦係数、物性、加熱・冷却条件で結果は変わります。機器設計では適用条件と実測との比較が必要です。
Gnielinski教授と800件の実験データ
Volker Gnielinski教授（カールスルーエ工科大学、ドイツ）は1976年の論文 "New Equations for Heat and Mass Transfer in Turbulent Pipe and Channel Flow" で、過去の研究者が蓄積した 約800組の実験データ を総括して相関式を提案した。それまでの相関式が特定の実験条件に合わせて提案されていたのに対し、Gnielinskiは「1つの式で全域をカバーする」ことにこだわった。この論文は被引用数が万を超え、熱流体工学分野で最も引用される論文の一つとなっている。
Gnielinski相関式の数値計算手法
CFDにおける壁関数検証
CFDでGnielinski式はどう使うんですか？ 直接CFDのソルバーに組み込まれてるわけじゃないですよね？
いい着眼点だ。CFDソルバーはNavier-Stokes方程式を直接解くから、Gnielinski式のような相関式は「ソルバーの内部計算」には使わない。Gnielinski式の出番は 検証（Verification & Validation） だ。具体的にはこういう使い方をする:
- 壁関数の妥当性検証：CFDで得られた壁面熱伝達係数 $h$ やNu数を、Gnielinski式の解析解と比較。2割以上ずれていたら、y+設定やメッシュ解像度を見直す
- メッシュ収束性確認：メッシュを段階的に細かくしたとき、CFDのNuがGnielinski式の値に収束していくかを確認
- 乱流モデルの選定：k-ε、k-ω SST、RSMなど異なるモデルでのNu予測値をGnielinski式と比較して、対象問題に適したモデルを選ぶ
つまり「答え合わせの正解」として使うわけですね。ところで y+ ってここで重要なんですか？
CFDと相関式の比較は妥当性を調べる一つの参考です。一致だけではモデルが正しいと証明できません。壁面処理に適したy+、格子・時間刻みの収束、物性と境界条件を確認し、可能なら実験データと比較します。
1D系統解析での利用
CFD以外でもGnielinski式って使いますか？
むしろ実務で最も多い使い方は 1D系統解析（system-level analysis）だ。例えば:
- 熱交換器設計：管内の熱伝達係数の概算や実験・CFDとの比較に相関式を利用します。製品の既定式と適用条件はバージョンとモデルの設定を確認してください。
- 配管系熱解析：FlownexやFLOWMASTERなどの1D熱流体コードで、管内強制対流のhを算出する標準モデル
- 原子力安全解析：RELAP5やTRACEなどの原子炉熱水力コードでGnielinski式ベースのモデルが使用されている
- 自動車冷却系：エンジン冷却ジャケット内の冷却水（水＋エチレングリコール混合液）の熱伝達予測
へぇ、原子力でも使われてるんですか。相関式1つでそこまで広く使われてるのは驚きです。
物性値の温度補正
物性値はどの温度で評価するんですか？ 壁温？ バルク温度？
Gnielinski式では原則 バルク平均温度 $T_b$（流体の断面平均温度）で物性値を評価する。ただし壁温と流体温度の差が大きいケースでは、Sieder-Tateの粘度補正を追加する:
ここで $\mu_b$ はバルク温度での粘度、$\mu_w$ は壁温での粘度だ。加熱の場合 $\mu_b > \mu_w$ なのでNuが増加し、冷却の場合はNuが減少する。特にオイル系流体で壁温差が50°C以上あるときはこの補正が効いてくる。
粘度って温度で大きく変わりますもんね。エンジンオイルなんか、冷間時と暖気後で全然違う。
そう。水なら20°C→80°Cで粘度は3分の1になる程度だが、エンジンオイルだと10分の1以上に変化する。この効果を無視するとNuが数十%ずれることもある。だから壁温差が大きいオイル冷却器の設計では粘度補正は必須だ。
入口効果の補正
管が短い場合（L/Dが小さい場合）はどうするんですか？
入口領域では速度・温度境界層が発達途上のため、局所Nuが軸方向位置に依存して大きくなる。Gnielinski自身が提案した入口効果補正はこうだ:
$L/D = 10$ でも補正は約 +21% と結構大きい。例えば自動車の冷却ジャケットのように短い流路が多い場合、入口効果を無視するとhを過小評価してしまう。
L/D = 10で+21%も変わるのは無視できないですね…。実務ではどっちを使うことが多いですか？
熱交換器のチューブは通常L/D > 50なので入口効果は小さいが、電子機器冷却のマイクロチャネル（L/D = 10〜30が多い）やガスタービン翼内部冷却流路では入口効果補正が重要だ。補正を入れるかどうかで設計マージンが変わるから、設計者は条件を見て判断する必要がある。
Gnielinski相関式の実務適用
熱伝達係数の算出手順
実際にGnielinski式で熱伝達係数hを計算するときのステップを教えてください！
手順は以下の5ステップだ:
- 流体物性値を取得：バルク温度 $T_b$ で $\rho$、$\mu$、$c_p$、$k$ を物性表から読む
- レイノルズ数を計算：$\mathrm{Re} = \rho u D / \mu$ （$u$ は断面平均流速、$D$ は管内径）
- プラントル数を計算：$\mathrm{Pr} = \mu c_p / k$
- 摩擦係数を計算：$f = (0.790 \ln \mathrm{Re} - 1.64)^{-2}$
- Nu → h を計算：Gnielinski式でNuを求め、$h = \mathrm{Nu} \cdot k / D$ で熱伝達係数を得る
ステップ4と5は、関数電卓かExcelがあれば簡単にできそうですね。
そう、計算自体は単純だ。一番の注意点は ステップ1の物性値。温度で大きく変わる粘度の値を正確に取ることが精度の鍵になる。NIST WebBookやRefpropが信頼できるデータソースだ。
設計計算の実例
具体的な数値で計算してみたいです！
じゃあ熱交換器のチューブ設計の例でやってみよう。
条件: 水、$T_b = 60°\text{C}$、管内径 $D = 20\,\text{mm}$、平均流速 $u = 2\,\text{m/s}$
60°Cの水の物性値:
- $\rho = 983\,\text{kg/m}^3$、$\mu = 4.67 \times 10^{-4}\,\text{Pa·s}$
- $c_p = 4,185\,\text{J/(kg·K)}$、$k = 0.654\,\text{W/(m·K)}$
Step 1: $\mathrm{Re} = 983 \times 2 \times 0.02 / (4.67 \times 10^{-4}) \approx 84,197$
Step 2: $\mathrm{Pr} = 4.67 \times 10^{-4} \times 4185 / 0.654 \approx 2.988$
Step 3: $f = (0.790 \ln 84197 - 1.64)^{-2} \approx 0.018666$（Darcy摩擦係数）
Step 4: $\mathrm{Nu} = \frac{(f/8)(\mathrm{Re}-1000)\mathrm{Pr}}{1+12.7\sqrt{f/8}(\mathrm{Pr}^{2/3}-1)} \approx 349.6$
Step 5: $h = 349.6 \times 0.654 / 0.02 \approx 11,430\,\text{W/(m}^2\text{·K)}$
h ≈ 11,430 W/(m²·K) となるのですね！
同じ条件でDittus-Boelter式の加熱用指数0.4を使うとNu ≈ 310.6、h ≈ 10,155 W/(m²·K)です。Gnielinski式の値より約11.2%小さくなりますが、式同士の比較だけでは、どちらが実測に近いかや安全側かは判断できません。
よくある失敗と対策
Gnielinski式の計算でやりがちなミスってありますか？
| 失敗パターン | 原因 | 対策 |
|---|---|---|
| Nu値が明らかにおかしい | Fanning摩擦係数とDarcy摩擦係数を混同（4倍の差） | Petukhov式はDarcy（Moody）摩擦係数。Fanningの場合は$f_\text{Darcy} = 4 f_\text{Fanning}$ |
| 遷移域でNuが負になる | Re < 1000で式に代入してしまった | Re < 2,300は適用外。層流式（Nu = 3.66 or 4.36）を使う |
| hが実験値と大きくずれる | 物性値の評価温度の誤り | バルク温度で評価。壁温差が大きい場合は粘度補正を追加 |
| 粗面管なのに滑面式を使用 | 管壁粗さの影響を無視 | 粗面ではColebrook-White式からfを求め、乖離を評価する |
| 流量から流速への換算ミス | 流路断面積の計算誤り | $u = Q / (\pi D^2/4)$。特に二重管やシェルアンドチューブでは内管径・外管径の区別に注意 |
Fanning摩擦係数とDarcy摩擦係数の混同って4倍もずれるんですか！ それは致命的ですね。
これは実務でも非常に多い間違いだ。アメリカの教科書はFanning（$f_F$）を使い、ヨーロッパの教科書はDarcy（$f_D = 4f_F$）を使うことが多い。Gnielinski式の原論文はDarcy摩擦係数なので、Petukhov式をそのまま使えばいい。ただしFluentの出力やHTRIのレポートがどちらの定義を使っているか、必ず確認すること。
ガスタービン翼冷却とGnielinski式
ガスタービンの冷却流路では形状、粗さ、回転、物性変化などが熱伝達に影響します。単純な円管相関式を実機設計に適用する前に、その適用性を実験や対象形状の解析で確認する必要があります。
Gnielinski相関式のソフトウェア比較
主要ソフトでの実装状況
Gnielinski式はCAEソフトにどう組み込まれてるんですか？
ソフトによって使い方が異なる。CFDソルバーでは「検証用の基準値」として、1D系統コードでは「熱伝達係数の計算モデル」として使われる:
| 製品・用途 | 確認できる内容 |
|---|---|
| COMSOL Pipe Flow | 公式の6.3資料では管内乱流の膜抵抗にGnielinski式を使用。適用条件と物性評価を確認。 |
| Flownex SE 8.10.0 | 2019年の公式資料では特定の熱交換器コンポーネントの管側で既定。対象機能とバージョンを確認。 |
| その他のCFD・熱交換器ソフト | 出力値の比較には利用できますが、内蔵式・既定値・設定方法は製品と版ごとの資料を確認してください。 |
各ソフトで、式の適用や既定値をどのように確認すればよいですか？
製品、バージョン、コンポーネント、流動条件ごとに公式資料と設定画面を確認します。相関式との一致だけでCFDや設計の精度は検証できません。
OpenFOAMでの実装
OpenFOAMでGnielinski式を使いたいときはどうするんですか？
OpenFOAMでは後処理で相関式と比較する方法や、条件が適切な場合に熱伝達係数を境界条件に与える方法があります。利用可能な機能は配布版とバージョンを確認してください:
- ポスト処理で比較：CFD結果の壁面熱流束と温度差からNuを算出し、Gnielinski式の解析値と比較するスクリプト（PythonやC++）を書く
- 境界条件として利用：conjugate heat transfer（共役熱伝達）の代わりに、外部流体側のhをGnielinski式で見積もり、
externalWallHeatFluxTemperature境界条件の $h$ 値として入力する
なるほど。利用する版の境界条件と適用範囲を確認する必要があるのですね。
Gnielinski相関式の先端研究
Gnielinski式の拡張と改良
Gnielinski式って1976年以降、改良されてるんですか？
Gnielinski自身が2013年に改訂版を発表している（VDI Heat Atlas, 2nd Edition）。主な改良点は:
- 遷移域の滑らかな接続：Re = 2,300〜4,000の範囲でNuを線形補間する手法の推奨。層流式のNulamと乱流式のNuturbを加重平均: $\mathrm{Nu} = (1-\gamma)\,\mathrm{Nu}_\text{lam} + \gamma\,\mathrm{Nu}_\text{turb}$ ここで$\gamma = (\mathrm{Re} - 2300)/(4000 - 2300)$
- 非円形断面への拡張：矩形、環状、三角形断面に対して水力直径 $D_h = 4A/P$ での適用条件を明示
- 回転流路への修正：ガスタービン翼内部のように回転する流路では、回転バースト数 $\mathrm{Bo} = \mathrm{Re}_\Omega^2 / \mathrm{Re}$ による補正が追加される
Gnielinski先生、自分で式をアップデートし続けてるんですね。すごい。
VDI Heat Atlas（ドイツ技術者協会の伝熱ハンドブック）の管内対流の章はGnielinskiが書いており、データの蓄積に応じて版を重ねるたびに推奨式を更新している。2013年版の遷移域補間法は、2300 < Re < 4000でNuが「飛ぶ」問題を解消するもので、数値シミュレーションとの整合性が向上した。
機械学習との融合
最近の機械学習ブームとGnielinski式の関係ってありますか？
興味深い研究が増えている。主な方向性は:
- PINN（Physics-Informed Neural Networks）：Gnielinski式を「物理的制約」としてニューラルネットに組み込み、実験データが少ない条件域でも信頼性の高い予測を行う
- サロゲートモデル：複雑な形状（フィン付き管、螺旋管）の伝熱性能をCFDデータベースから機械学習で構築。Gnielinski式の適用外条件（粗面、回転流路など）の補正係数を学習させる
- 逆問題：実測のNuデータからGnielinski式の修正係数を最適推定し、特定のプラントに特化した高精度モデルを構築する
古典的な相関式と最先端のAIが融合する時代なんですね。
重要なのは、機械学習は「ブラックボックス」だから、物理的な制約（Gnielinski式のような相関式）を組み込まないと非物理的な予測をしてしまうことがある。例えばReが2倍になったのにNuが半分になる、というような結果をMLモデルが出すことがある。Gnielinski式を制約として入れておけば、少なくとも傾向は物理的に妥当な範囲に収まる。これがPINNの強みだ。
Gnielinski相関式のトラブル対応
よくあるエラーと対策
Gnielinski式を使ってて「あれ？おかしいな」ってなる場面ってありますか？
あるある。代表的なものを挙げよう:
| 症状 | 考えられる原因 | 対策 |
|---|---|---|
| Nuが教科書の例題と合わない | Fanning/Darcy摩擦係数の混同 | Petukhov式はDarcy定義。$f_D = 4 f_F$ を確認 |
| Nuが負値になる | Re < 1000を代入 | Re < 2,300は適用外。適用範囲を必ずチェック |
| 計算値と実験値が30%以上乖離 | 入口効果・粗面効果・物性値誤り | 入口補正、粗面補正を適用。物性値の評価温度を確認 |
| オイル系でNuが過大 | 壁温補正（Sieder-Tate）の未適用 | $(\mu_b/\mu_w)^{0.14}$ 補正を追加 |
| Excelで計算するとln(Re)でエラー | Excelのln関数名がLN()であることの見落とし | LOG()は常用対数。自然対数はLN()を使う |
CFD結果との乖離
CFDのNuとGnielinski式のNuが20%以上ずれるときはどうすれば？
以下のチェックリストを上から順に確認していく:
- y+ の確認：壁関数使用時は y+ = 30〜300、Low-Reモデル使用時は y+ < 1 を満たしているか
- メッシュ収束性：メッシュを2倍に細かくしてNuが5%以上変化するなら、まだ収束していない
- 流れの発達区間：CFDの入口からどれだけ下流でNuを評価しているか。L/D < 10なら入口効果が残っている
- 乱流モデルの選択：k-εは壁面近傍の温度分布が不正確になりがち。k-ω SST の方がGnielinski式との一致度が高い傾向がある
- 物性値の一貫性：CFDに入力した物性値とGnielinski式に使った物性値が同じ温度で評価されているか
- 温度境界条件：Gnielinski式は等壁温（constant wall temperature）と等熱流束（constant heat flux）の両方に近い精度を持つが、CFDの境界条件がどちらかを確認する
y+ が30〜300の間に収まってなかったら壁関数の結果は信用できないってことですね。
y+の目安は使用する壁面処理や乱流モデルにより異なります。流速や粘度が変わるとy+も変化するので、ソルバーの推奨範囲を確認し、格子収束・物性・境界条件も併せて調べます。
関連トピック
なった
詳しく
報告