Sieder-Tate相関式

Sieder-Tate相関式の理論基礎

管内乱流熱伝達の相関式の系譜

円管内の乱流強制対流の熱伝達係数を無次元相関で与える試みは、Dittus-Boelter式（1930年代）に始まります。Dittus-Boelter式は物性を一定とみなすため、粘度が温度で大きく変わる流体（油・グリセリン・高粘性プロセス液）では、加熱時と冷却時で実測が系統的にずれる問題がありました。Sieder-Tate（1936）は、この物性変化の影響を「バルク粘度と壁面粘度の比」による補正項として取り込んだ古典で、今日でも熱交換器の初期設計・検算に現役で使われています。

式と各項の意味

乱流域のSieder-Tate相関式は次の形です。

$$ Nu = 0.027\, Re^{0.8}\, Pr^{1/3} \left( \frac{\mu_b}{\mu_w} \right)^{0.14} $$

ここで \( \mu_b \) はバルク（混合平均）温度での粘度、\( \mu_w \) は壁温での粘度です。物性はすべてバルク温度で評価し、粘度比の項だけが壁温の情報を持ちます。加熱時（壁が熱い）は壁近傍の粘度が下がって境界層が薄くなり熱伝達が増える——\( \mu_b/\mu_w > 1 \) で補正は増加側。冷却時はその逆です。適用範囲の目安は \( Re > 10^4 \)、\( 0.7 < Pr < 16700 \)、\( L/D > 10 \)（十分発達した流れ）です。

層流版と「どの式を選ぶか」

層流（\( Re < 2300 \)）で熱的助走区間が長い場合にはSieder-Tateの層流式 \( Nu = 1.86\,(Re\,Pr\,D/L)^{1/3} (\mu_b/\mu_w)^{0.14} \) があり、乱流域では後発のGnielinski式（摩擦係数ベース、遷移域 \( 3000 \lesssim Re \lesssim 5\times10^6 \) まで拡張、精度±10%程度）が現代の推奨です。実務の選択指針は、「精度重視・広範囲ならGnielinski、粘性液の粘度補正が主題ならSieder-Tate型補正を併用、手早い概算ならDittus-Boelter」という三択で覚えるのが実用的です。いずれも誤差帯（Sieder-Tateで±20〜25%程度）を明記して使います。

計算手順と数値例

壁温が未知のときの反復手順

数値例——油の加熱で補正はどれくらい効くか

機械油（バルク60℃: \( \mu_b \approx 0.072 \) Pa·s、壁温100℃: \( \mu_w \approx 0.017 \) Pa·s）を管内で加熱するケースを考えると、粘度補正は

$$ \left(\frac{\mu_b}{\mu_w}\right)^{0.14} = \left(\frac{0.072}{0.017}\right)^{0.14} \approx 4.2^{0.14} \approx 1.22 $$

つまり熱伝達率が22%増し——油では補正が測定誤差ではなく設計を左右する大きさになります。同じ条件を冷却（壁温20℃・\( \mu_w \approx 0.28 \) Pa·s）にすると補正は \( (0.072/0.28)^{0.14} \approx 0.83 \) で17%減。加熱と冷却で合計4割も違う——Dittus-Boelter式（補正なし）が粘性液で外す理由がこの数字に集約されています。一方、水（60℃→100℃で粘度比約1.7）では補正は \( 1.7^{0.14} \approx 1.08 \) と8%程度に留まります。

熱交換器計算への組み込み

Sieder-Tateで得た管内側 \( h_i \) は、汚れ係数・管壁伝導・外側熱伝達と直列に合成して総括伝熱係数 \( U \) を作り、LMTD法またはNTU-ε法で熱交換器の伝熱面積・出口温度を解きます。バルク温度は入口出口の算術平均（温度変化が大きい場合は区間分割）で評価します。相関式の±20%はUの不確かさに直結するため、設計余裕（面積マージン）の根拠として誤差帯を明記します。

実務適用の指針

適用前チェックリスト

• Reの確認 — \( Re < 2300 \)（層流）なら層流式へ、\( 2300 < Re < 10^4 \)（遷移域）ならSieder-Tate乱流式の適用外——Gnielinski式を使う
• 発達長 — 入口から \( L/D < 10 \) の区間は助走効果でhが高い。短管では助走補正を検討
• 非円形断面 — 水力直径 \( D_h = 4A/P \) で代用可（乱流のみ。層流は断面形状依存の厳密解を使う）
• 自然対流の混入 — 低流速・大温度差ではリチャードソン数を確認。混合対流域は強制対流相関の適用外
• 物性評価温度 — 主物性はバルク温度、\( \mu_w \) のみ壁温。この規約の混同が最頻出ミス

CFDとの役割分担

管内流はCFDでも解けますが、実務では「単純な円管はSieder-Tate/Gnielinskiの方が速くて信頼できる」が原則です。CFDに進む価値があるのは、①複雑断面・曲がり・フィン付きで相関式が存在しない、②入口効果・局所ホットスポットの分布が欲しい、③物性の温度依存が極端で相関式の外挿になる、の場合です。その際も単純円管でCFD設定を相関式に対して検証してから実形状へ進むのが、平板検証と同じ作法です。

誤差帯込みの報告

相関式ベースの設計報告には、①使用した式と適用範囲の充足確認、②物性の出典と評価温度、③相関式の公称誤差帯（Sieder-Tate ±20〜25%）、④その誤差が設計量（伝熱面積・出口温度）に与える影響、を含めます。「h = 1234 W/m²K」と点値だけ書く報告は、相関式の性質上、精度を装った過剰表示になります。

ツールでの取り扱い

ツール別の実装状況

CFDで粘度補正効果を再現する条件

CFDでSieder-Tateの粘度補正に相当する物理を捉えるには、粘度の温度依存性を物性モデルに入れることが必須です（定物性のCFDはDittus-Boelter相当の答えしか出しません）。加えて壁近傍の粘度変化は境界層内で起きるため、低Re型壁処理（y+≤1）でないと効果が壁関数の仮定に埋もれます。油系のCHT解析で「CFDが相関式より系統的に低い」ときは、まずこの2点（温度依存粘度・壁解像）を確認します。

先端動向

適用限界の拡張——超臨界・ナノ流体・非ニュートン

超臨界圧流体（sCO2サイクル・超臨界水）では擬臨界点近傍の物性激変により古典相関が大きく外れ、専用相関と可変物性CFDの突き合わせが活発な研究領域です。ナノ流体・非ニュートン流体（食品・ポリマー溶液）でも、Sieder-Tate型の「粘度比補正」という発想を拡張した相関の提案が続いています。共通する教訓は、相関式は開発時のデータベースの範囲でのみ信頼できるということ——適用範囲確認の重要性は新分野ほど高まります。

データ駆動相関とハイブリッド化

大量の実験・DNSデータから機械学習で管内熱伝達を回帰する研究では、無次元数（Re・Pr・粘度比）を入力特徴量に保つ設計が主流です。古典相関の関数形を骨格に係数をデータで最適化するアプローチは、外挿安全性と精度改善のバランスが良く、産業ツールへの実装が始まっています。一方で学習データ外の流体・条件での劣化は避けられず、古典相関を「サニティチェックの下限・上限」として併走させる使い方が推奨されます。

局所計測との融合

赤外線サーモグラフィや分布型光ファイバ温度計測により、管軸方向の局所熱伝達分布が実測できるようになり、「平均hの相関式」から「局所分布の検証」へ検証の解像度が上がっています。助走区間・曲がり下流・脈動の影響など、平均値に埋もれていた局所現象の相関化・モデル化が進行中です。

トラブル対応

症状別の原因と対策

式の選択に迷ったら

関連記事：平板上の強制対流、強制対流の記事一覧、共役熱伝達の反復カップリング。