助走区間効果
助走区間の理論基礎
2種類の助走区間——流力的と熱的
管に流体が入ると、壁の影響が中心へ広がりきるまでの区間で速度分布が変化し続けます。これが流力的助走区間 \( L_h \)。同様に、加熱・冷却が始まってから温度分布の形が定まるまでが熱的助走区間 \( L_t \) です。層流の目安は
$$ L_h \approx 0.05\, Re\, D, \qquad L_t \approx 0.05\, Re\, Pr\, D $$
で、熱的助走はPr倍長いのがポイントです。乱流では乱れの混合が速く、どちらも \( 10\sim60\,D \) 程度で発達します。「発達した流れ」の相関式(Sieder-Tate・Gnielinski等)はこの助走が終わった後の世界の式——適用の前提を破っていないかを判定するのが、この長さの実務的な意味です。
助走区間では熱伝達が「高い」
助走区間って「まだ整っていない区間」ですよね。熱伝達的には悪い区間なんですか?
逆なんだ、そこが面白いところ。加熱開始直後は温度境界層が薄いから、局所Nu数は入口で理論上無限大から始まって、下流へ向かって発達値(層流の等壁温なら \( Nu = 3.66 \))へ単調に減衰していく。つまり助走区間は熱伝達のボーナス区間。短い熱交換器ほど平均Nuが発達値より高くなる——この効果を無視すると性能を過小評価するし、逆に発達値の適用範囲外なのに気づかず短管を設計すると、実測が「なぜか良すぎる」ことになる。良い方向のズレでも、原因不明のズレは検証としては failed だからね。
Graetz問題という古典
層流の熱的助走を厳密に扱ったのがGraetz問題(速度分布は発達済み・温度が発達していく設定)で、無次元距離 \( x^* = x/(Re\,Pr\,D) \)(Graetz数の逆数)が支配パラメータです。\( x^* \gtrsim 0.05 \) で発達域とみなせる——先の \( L_t \) の式はこの厳密解の要約です。平均Nu数への助走効果を織り込む実用式としてHausen式などの助走補正付き相関が整備されており、短管では発達域相関ではなくこちらを使うのが正解です。
計算での取り扱い
数値例で「長さの感覚」を持つ
水(\( Pr \approx 7 \))、\( Re = 1000 \)(層流)、管径 \( D = 10 \) mmとすると:流力的助走 \( L_h \approx 0.05 \times 1000 \times 0.01 = 0.5 \) m、熱的助走 \( L_t \approx 0.05 \times 1000 \times 7 \times 0.01 = \) 3.5 m。つまり長さ1 mの熱交換器なら全域が熱的助走区間で、発達域の \( Nu = 3.66 \) を使うと熱伝達を大幅に過小評価します。油(\( Pr = 100 \)超)ではさらに極端で、層流の油クーラーは実質的に「常に助走区間」です。逆に空気(\( Pr = 0.7 \))や乱流では助走は短く、\( L/D > 60 \) なら無視して構いません。設計計算の最初に \( L_t/L \) を見積もる——このひと手間が式の誤用を防ぎます。
CFDでの入口条件設計
CFDでは助走の扱いは「入口境界条件の設計問題」になります。
| 目的 | 入口条件 | 注意 |
|---|---|---|
| 実機の入口形状の影響も含めて評価 | 上流の実形状(ヘッダ・曲がり)ごとモデル化 | 最も忠実だが計算域が大きい |
| 発達流の性能だけ欲しい | 発達プロファイルを入口に直接指定、または周期境界(fully-developed periodic) | 助走効果は評価されない |
| 助走効果そのものを評価 | 一様流入+十分な助走距離をモデル化 | 「一様流入直後に加熱開始」は実機に無い理想化=hが過大に出る |
頻出ミスは、一様流入のすぐ下流から加熱壁を始める設定です。これは流力・熱の助走が同時に始まる特殊条件で、実機(上流配管で速度は発達済み、加熱部から温度助走が始まる)と異なります。「速度は発達・温度は入口から」というGraetz型の状況を再現したいなら、断熱の助走管を前置するのが正しいモデル化です。
局所Nu分布の抽出と検証
CFDの検証は平均値でなく局所 \( Nu_x \) の軸方向分布で行います。層流ならGraetz解と重ね、入口近傍の減衰カーブと発達値への漸近の両方が再現されているかを確認します。局所Nuの定義(基準温度に混合平均温度を使う)に注意——壁温と入口温度の差で定義すると発達域でNuが一定にならず、「収束しない」ように見えます(平板で述べたhの定義問題の管内版です)。
実務適用の指針
設計計算のチェックリスト
- \( L_h, L_t \) を見積もり、使う相関式の前提(発達域か助走込みか)と整合させたか
- 短管・高Pr流体では助走補正付き相関(Hausen型・Sieder-Tate層流式)を使ったか
- 乱流でも \( L/D < 10 \) の極短管では入口効果(縮流・剥離)を別途考慮したか
- 直列する曲がり・バルブの下流では助走がリセットされる(再発達)ことを考慮したか
- CFDの入口条件が想定する物理状況(実機のどこを切り出したか)を文書化したか
助走効果は両刃——性能とリスク
助走区間の高い熱伝達は、コンパクト熱交換器では積極的に活用されます(流路を短く分割し、何度も助走を「やり直させる」——オフセットフィンやルーバーフィンの原理はまさにこれです)。一方でリスク側もあります:入口近傍は熱流束・温度勾配が最大になるため、熱応力・スケール析出・局所沸騰が入口に集中します。「平均で設計し、入口で壊れる」を避けるため、局所値の評価を設計プロセスに残すことが重要です。
実験と比較するときの注意
実験装置の「入口」は理想的な一様流でも完全発達流でもありません。整流格子・縮流ノズル・上流配管長で助走状態が変わり、これが実験室間でNuデータがばらつく古典的な原因の一つです。CFD・相関式と実験を突き合わせる際は、装置の入口構成(上流直管長 \( L/D \))を必ず記録・確認します。不一致の調停では、まず三者の「入口条件の想定」を揃えるのが先決です。
ツールでの取り扱い
ツール別の機能
| ツール | 助走関連の機能 |
|---|---|
| Ansys Fluent | 周期境界(translational periodic)で発達流を直接計算。profile書き出し→別解析の入口へマップ |
| OpenFOAM | mappedPatch(下流断面を入口へ写像して発達流を自給)、codedFixedValueで解析プロファイル指定 |
| STAR-CCM+ | fully-developed interface・テーブルによるプロファイル指定 |
| COMSOL | 「完全発達流」入口オプション(内部的に補助問題を解く) |
| 1D熱流体ツール/熱交換器設計ソフト | 助走補正付き相関の内蔵。どの相関か(発達域式に補正係数か、Graetz系か)を確認して使う |
発達プロファイルの作り方の定石
「発達済みの速度・温度を入口に与えたい」場合の定石は3通り:①短い周期境界解析を先に回してプロファイルを抽出、②同じ断面の長い直管を1回解いて出口プロファイルを再利用、③層流なら解析解(ポアズイユ分布)を直接式で与える。③が使えるのは速度だけで、温度の発達プロファイルは壁条件に依存するため①②で作るのが安全です。作ったプロファイルは質量流量・エネルギー流量が目標値と一致することを検算してから使います。
先端動向
マイクロチャネルへの展開
電子機器冷却のマイクロチャネルは高Re化が難しく層流運用が基本のため、助走区間の設計的活用が性能の主要因になります。流路分割・再合流による境界層リスタート、発達抑制形状(波状流路)の最適化が活発で、従来「補正項」だった助走効果が設計の主役になっています。気体マイクロ流路では希薄化(スリップ境界)と助走の複合も研究対象です。
非定常・脈動流の助走
往復動機器・生体流(血流)のような脈動流では、助走長が時間変化し、位相によって逆流も生じます。時間平均Nuが定常より増減どちらにも振れることが知られており、脈動条件(Womersley数)ごとの整理が進んでいます。定常の助走理論をそのまま適用できない領域として、CFD検証の定番題材になっています。
データ駆動の高速予測
複雑断面・分岐を含むネットワーク流路で、断面形状→助走含みの局所Nu分布を学習するサロゲート(GNN・演算子学習)が研究されています。1D熱流体ネットワーク解析の「管路要素」を、助走・再発達を含む学習済みモデルで置き換え、系全体の設計探索を高速化する方向です。
トラブル対応
症状別の原因と対策
| 症状 | 考えられる原因 | 対策 |
|---|---|---|
| CFDのhが相関式より大幅に高い | 助走区間を発達域相関と比較している、一様流入直後の加熱開始 | L_t/Lを確認。助走補正付き相関と比較、入口条件を実機相当に |
| 実験が相関式より高い(短管) | 助走効果(正常) | Hausen型の助走込み相関で再評価 |
| 発達域のはずなのに局所Nuが一定にならない | Nuの基準温度の定義誤り(混合平均温度を使っていない) | 各断面の混合平均温度で再定義 |
| 周期境界解析が収束しない・非物理的 | 周期区間に正味の加熱があるのに温度周期条件が不整合 | 温度はスケール化周期条件(または線形成分除去)を使用 |
| 曲がり下流の実測hが計算より高い | 二次流れと再発達(助走リセット)の未考慮 | 曲がりを含めてモデル化、または再発達込みの局所相関を使用 |
| 入口近傍で壁温・熱応力が想定超過 | 入口の高熱流束集中(助走の裏面) | 局所値で設計評価。入口の熱流束緩和(段階加熱)を検討 |
実務の一言まとめ
助走区間について、設計者として覚えておくべきことを一言でまとめるなら?
「相関式を使う前に \( L_t = 0.05\,Re\,Pr\,D \) を暗算せよ」——これに尽きるね。この長さが装置より短ければ発達域の式でよし、同等以上なら助走込みの式かCFDへ。たった1行の見積りで、層流×高Pr流体という「発達域相関が最も派手に外れる組合せ」を確実に検知できる。熱設計のレビューで最初に聞くべき質問のひとつだよ。ついでに言えば、乱流化(Re>4000)すれば助走もPr依存もほぼ消えて設計が簡単になる——「迷ったら乱流域で設計」も、先人の知恵として覚えておいて損はない。
関連記事:Sieder-Tate相関式、平板上の強制対流、強制対流の記事一覧。
なった
詳しく
報告