Gnielinski相関式

ざっくり言うと、そのとおりだ。Gnielinski式は Re = 2,300〜5×10⁶ という非常に広いレンジで使え、しかも遷移域（Re = 2,300〜10,000あたり）もカバーする。Dittus-Boelter式はRe > 10,000の完全乱流域しか想定していないから、遷移域では使えないんだ。

Gnielinski式は実験データに対して ±10% の精度。対してDittus-Boelter式は ±25% だ。つまり精度が2倍以上良い。だから現在の教科書や設計基準では、特別な理由がない限りGnielinski式が推奨されている。CFDの壁関数検証でもGnielinski式をベンチマーク基準にすることが多いよ。

いや、Dittus-Boelter式は形がシンプルだから、概算や手計算のクロスチェックには今でも便利だ。ただし設計の最終値を出すときはGnielinski式を使うのが現場のスタンダードだね。

いい質問だ。歴史的な流れがある。まず1970年にPetukhov（ペトゥホフ）が、運動量輸送と熱輸送のアナロジーに基づいた乱流熱伝達モデルを提案した。ただしPetukhov式はRe > 10⁴ の完全乱流域限定だった。

1976年にGnielinski（グニエリンスキ）がPetukhov式を改良して、Reの項を (Re − 1000) に置き換えるという修正を加えた。たったこれだけの修正で遷移域まで精度よく予測できるようになったんだ。シンプルだけど物理的に深い修正で、「乱流が完全には発達していない領域では実効的なReが小さくなる」という現象をうまく捉えている。

そう、優れた工学モデルはシンプルなことが多い。Gnielinski本人はドイツのカールスルーエ工科大学で膨大な実験データと比較検証して、この修正が最も広範囲に高精度を保つことを示した。今では「内部流の乱流熱伝達と言えばGnielinski式」というくらいの定番になっている。

Gnielinski相関式は以下の通りだ:

この項は 粘性底層の熱抵抗 を表している。壁面のごく近くには粘性が支配する薄い層（粘性底層）があって、そこでは乱流による混合がほぼ効かない。Prが1から大きく離れる流体ほどこの層の影響が大きくなる。例えばオイル（Pr ≈ 100〜1000）では粘性底層の熱抵抗がかなり大きく、Nuは水（Pr ≈ 7）に比べて下がる。逆に気体（Pr ≈ 0.7）では粘性底層が薄いから影響は小さい。

そのとおり。Dittus-Boelter式はPrをべき乗で近似しているだけだが、Gnielinski式は粘性底層と乱流コアの熱抵抗を分離してモデル化している。だから特にPrが1から離れた流体（オイル、高粘度冷媒など）で精度が格段に良くなるんだ。

Gnielinski式とセットで使うのが Petukhov（ペトゥホフ）の摩擦係数 だ:

そう。ただし粗面（relative roughness ε/D > 0）の場合はPetukhov式のままでは使えない。粗面ではColebrook-White式やChurchill式から $f$ を求めて、それをGnielinski式に入れることになる。ただしその場合は「Gnielinski式の精度保証範囲外」だということは理解しておく必要がある。

いや、しっかり適用条件がある。これを外すと精度が保証されない:

そう。液体金属は粘性底層がほとんど存在しないので、熱伝導がバルク流れの中まで支配的になる。Gnielinski式の分母にある粘性底層モデルが成り立たないんだ。液体金属にはLyon-Martinelli式やSeban-Shimazaki式を使う。逆に超高粘度の流体（Pr > 2,000）ではSieder-Tate式の壁温補正付きモデルの方が信頼できることが多い。

例えば水が円管内を流れるケースで比較してみよう。条件: D = 25 mm、T_bulk = 50°C（Pr ≈ 3.56）。

そのとおり。特にオイル（Pr ≈ 100以上）になるとその差はさらに広がる。Dittus-Boelter式のPr^0.4 近似が粗くなるからだ。機器メーカーの設計基準書でも、HTRIやASPENのような熱交換器設計ツールでも、デフォルトはGnielinski式だよ。

いい着眼点だ。CFDソルバーはNavier-Stokes方程式を直接解くから、Gnielinski式のような相関式は「ソルバーの内部計算」には使わない。Gnielinski式の出番は 検証（Verification & Validation） だ。具体的にはこういう使い方をする:

非常に重要だ。壁関数アプローチ（y+ = 30〜300）を使う場合、Gnielinski式のNuとCFDのNuが±10%以内で一致していれば、壁関数が適切に機能していると判断できる。壁面解像（Low-Re）メッシュ（y+ < 1）の場合は、Gnielinski式とほぼ完全に一致するはずだ。一致しなければ乱流モデルか境界条件に問題がある。

むしろ実務で最も多い使い方は 1D系統解析（system-level analysis）だ。例えば:

Gnielinski式では原則 バルク平均温度 $T_b$（流体の断面平均温度）で物性値を評価する。ただし壁温と流体温度の差が大きいケースでは、Sieder-Tateの粘度補正を追加する:

そう。水なら20°C→80°Cで粘度は3分の1になる程度だが、エンジンオイルだと10分の1以上に変化する。この効果を無視するとNuが数十%ずれることもある。だから壁温差が大きいオイル冷却器の設計では粘度補正は必須だ。

入口領域では速度・温度境界層が発達途上のため、局所Nuが軸方向位置に依存して大きくなる。Gnielinski自身が提案した入口効果補正はこうだ:

熱交換器のチューブは通常L/D > 50なので入口効果は小さいが、電子機器冷却のマイクロチャネル（L/D = 10〜30が多い）やガスタービン翼内部冷却流路では入口効果補正が重要だ。補正を入れるかどうかで設計マージンが変わるから、設計者は条件を見て判断する必要がある。

手順は以下の5ステップだ:

そう、計算自体は単純だ。一番の注意点は ステップ1の物性値。温度で大きく変わる粘度の値を正確に取ることが精度の鍵になる。NIST WebBookやRefpropが信頼できるデータソースだ。

じゃあ熱交換器のチューブ設計の例でやってみよう。

水の乱流強制対流では典型的な値だ。ちなみにDittus-Boelter式で同条件を計算すると Nu ≈ 390、h ≈ 12,750 W/(m²·K) で 約11%過大評価 になる。この差は熱交換器の伝熱面積の設計に直結する。11%の過大評価は面積の過小設計を招き、結果として性能未達のリスクがある。

これは実務でも非常に多い間違いだ。アメリカの教科書はFanning（$f_F$）を使い、ヨーロッパの教科書はDarcy（$f_D = 4f_F$）を使うことが多い。Gnielinski式の原論文はDarcy摩擦係数なので、Petukhov式をそのまま使えばいい。ただしFluentの出力やHTRIのレポートがどちらの定義を使っているか、必ず確認すること。

ソフトによって使い方が異なる。CFDソルバーでは「検証用の基準値」として、1D系統コードでは「熱伝達係数の計算モデル」として使われる:

そう。HTRIは会員企業の実験データに基づいて独自の係数修正を加えている場合もあるが、ベースはGnielinski式だ。ただしフィン付き管やコルゲート管など特殊な伝熱促進管では、Gnielinski式に追加の補正因子を掛ける独自モデルを使っている。

OpenFOAMはN-S方程式を直接解くCFDコードなので、内部にGnielinski式は組み込まれていない。使い方は2つ:

Gnielinski自身が2013年に改訂版を発表している（VDI Heat Atlas, 2nd Edition）。主な改良点は:

VDI Heat Atlas（ドイツ技術者協会の伝熱ハンドブック）の管内対流の章はGnielinskiが書いており、データの蓄積に応じて版を重ねるたびに推奨式を更新している。2013年版の遷移域補間法は、2300 < Re < 4000でNuが「飛ぶ」問題を解消するもので、数値シミュレーションとの整合性が向上した。

興味深い研究が増えている。主な方向性は:

重要なのは、機械学習は「ブラックボックス」だから、物理的な制約（Gnielinski式のような相関式）を組み込まないと非物理的な予測をしてしまうことがある。例えばReが2倍になったのにNuが半分になる、というような結果をMLモデルが出すことがある。Gnielinski式を制約として入れておけば、少なくとも傾向は物理的に妥当な範囲に収まる。これがPINNの強みだ。

あるある。代表的なものを挙げよう:

以下のチェックリストを上から順に確認していく:

そう。y+ = 5〜30の「バッファ層」に最初のセルが位置していると、壁関数もLow-Reモデルも精度が出ない「デッドゾーン」になる。Gnielinski式との乖離の原因で最も多いのがこのy+問題だ。流速や粘度が変わればy+も変わるから、設計変更のたびに確認が必要だよ。