温度境界層 — 壁面熱伝達予測のトラブルシューティング
より充実した内容を thermal-boundary-layer.html でご覧いただけます。
壁面熱伝達係数が理論値と合わない
層流平板で $Nu_x$ を計算したのに、Pohlhausen解と10%以上ずれています。
まず基本を確認しよう。(1) 流れは本当に層流か？ Re数を確認して $Re_x < 5 \times 10^5$ であること。(2) 壁面温度条件は等温か等熱流束か？ Pohlhausen解は等温壁の式なので、等熱流束なら $Nu_x = 0.453 Re_x^{1/2} Pr^{1/3}$ を使うべき。(3) 入口条件は一様流か？ 入口にフリーストリーム乱流があると遷移が早まってNu数が増加する。
条件は全部合っているのにずれます。
離散化スキームの精度を確認しよう。First Order Upwindは数値拡散が大きく、温度境界層を実際より厚く（=Nu数を低く）する。Second Order Upwindに変更すべきだ。また、流れ方向のメッシュが粗いと前縁付近の急峻なNu数変化を捉えられない。前縁近傍で $\Delta x$ を特に小さくすること。
壁面関数領域でのNu数の不連続
壁面に沿ってNu数を追跡すると、セルサイズが変わるところでNu数がジャンプします。
壁面第一セルの $y^+$ がセルサイズの変化に伴って変わり、壁面関数の切り替えが起きている可能性がある。FluentのEnhanced Wall Treatmentでもブレンディング領域で微小な不連続が出ることがある。対策は壁面近傍のメッシュを均一化すること。Inflation layer の成長率を壁面に沿って一定に保つのがポイントだ。
高Pr数流体でNu数が大幅に過大
エンジンオイル（$Pr \approx 500$）の管内流れでNu数がDittus-Boelterの2倍になります。
高Pr数では温度境界層が極めて薄いので、壁面第一セルが温度境界層外に出てしまっている可能性がある。$y^+ \cdot Pr < 5$ を目安にメッシュを設計しよう。$Pr = 500$ なら $y^+ < 0.01$ が理想だが、これは現実的に厳しい。
どうすればいいですか？
$y^+ < 1$ まで細かくした上で、壁面関数に高Pr数補正が含まれているモデルを使う。FluentのEnhanced Wall Treatment + Thermal Effectsの組合せが最も信頼性が高い。OpenFOAMならalphatJayatillekeWallFunctionを必ず使うこと。それでもGnielinski式との乖離が大きい場合は、物性値の温度依存性（特に粘度）を入力しているか確認しよう。
y+を間違えると「壁が断熱材」になる
温度境界層のトラブルシュートで最も多いのが「熱伝達係数が実験の半分以下になる」という事例だ。原因を調べると、ほぼ確実にy+の設定ミスが出てくる。標準k-εで壁関数を使う場合、y+は30〜300の範囲に入れないといけないのに、「メッシュを細かくするほど良い」という先入観でy+ ＜ 5の細かいメッシュを使ってしまうケースが多い。壁関数は「壁から離れた粗いメッシュ」を前提とした近似式なので、細かすぎると物理が正しく再現されず、実質的に壁が断熱材になる。対策は「使うモデルに合ったy+を徹底する」か「低Re数モデルに切り替えてy+＜1」のどちらかに統一すること。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——温度境界層の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告