平板上の強制対流
平板強制対流の理論基礎
境界層という主役
一様流 \( U_\infty \) が平板に沿って流れるとき、壁面近傍には速度が0から \( U_\infty \) まで変化する速度境界層と、壁温 \( T_w \) から主流温度 \( T_\infty \) まで変化する温度境界層が発達します。熱伝達の抵抗はほぼこの薄い層の中に集中しており、両層の厚さ比を決めるのがプラントル数 \( Pr = \nu/\alpha \) です（空気 \( Pr \approx 0.7 \)：両層ほぼ同厚、油 \( Pr \gg 1 \)：温度境界層が薄い、液体金属 \( Pr \ll 1 \)：温度境界層が厚い）。平板はあらゆる対流熱伝達の「基準問題」であり、CFD・実験・相関式の三者を突き合わせられる数少ない構成として、V&Vの検証ベンチマークにも使われます。
層流域——Blasius解と局所Nu数
先端からの距離 \( x \) のレイノルズ数 \( Re_x = U_\infty x / \nu \) が臨界値（平板の目安 \( Re_{x,c} \approx 5\times10^5 \)、主流乱れで大きく変わる）以下では境界層は層流で、相似解（Blasius解）から局所ヌセルト数が導かれます（等温壁・\( Pr \gtrsim 0.6 \)）。
$$ Nu_x = 0.332\, Re_x^{1/2} Pr^{1/3}, \qquad \overline{Nu}_L = 0.664\, Re_L^{1/2} Pr^{1/3} $$
平均は局所の2倍——熱伝達率は先端で最大（境界層が薄い）、下流ほど \( x^{-1/2} \) で低下する分布です。等熱流束壁では係数が0.332→0.453に置き換わります。「壁の条件（等温か等熱流束か）で係数が変わる」ことは、CFD検証で相関式を突き合わせる際に最初に確認すべき前提です。
乱流域と混合平均
遷移後の乱流境界層では熱伝達が大きく増強され、代表的な局所式と、層流域・乱流域を通しで平均した混合平均式は次の形です（\( Re_{x,c} = 5\times10^5 \) 前提、\( 0.6 < Pr < 60 \)）。
$$ Nu_x = 0.0296\, Re_x^{4/5} Pr^{1/3}, \qquad \overline{Nu}_L = (0.037\, Re_L^{4/5} - 871)\, Pr^{1/3} $$
物性はすべて膜温度 \( T_f = (T_w + T_\infty)/2 \) で評価するのが規約です。相関式の誤差は理想条件でも±10〜20%程度あり、「相関式＝真値」ではなく「実験に裏打ちされた基準帯」として扱います。
手計算とCFDの計算手法
手計算の実例——空気・0.5m平板
式は分かったんですが、実際の数字でどれくらいのhになるのか感覚がなくて…。
一度手で通すと一生の財産になるよ。空気（膜温度50℃想定: \( \nu \approx 1.8\times10^{-5} \) m²/s、\( k \approx 0.028 \) W/m·K、\( Pr = 0.70 \)）、\( U_\infty = 10 \) m/s、板長 \( L = 0.5 \) mとしよう。まず \( Re_L = 10 \times 0.5 / 1.8\times10^{-5} \approx 2.8\times10^5 \)——臨界値以下だから全面層流だ。平均Nu数は \( \overline{Nu}_L = 0.664 \times \sqrt{2.8\times10^5} \times 0.70^{1/3} \approx 0.664 \times 529 \times 0.888 \approx 312 \)。平均熱伝達率は \( \bar{h} = \overline{Nu}_L\, k / L \approx 312 \times 0.028 / 0.5 \approx 17 \) W/m²K。強制対流の空気でhは数十W/m²Kのオーダー——この相場観があれば、CFDで h=500 が出た瞬間に「何かおかしい」と気づける。
CFDで解くときのメッシュと壁処理
- 低Re型（壁積分） — 壁第1セルを \( y^+ \le 1 \)、境界層内に15〜25層。熱伝達を精度良く出す基本形
- 壁関数型 — \( 30 \lesssim y^+ \lesssim 300 \) に第1セルを置く。粗いが大規模モデル向け。中間の \( y^+ \approx 5\sim30 \) が最悪で、どちらの仮定も破れて熱伝達が大きくずれる
- 流入条件 — 主流乱れ度は遷移位置を強く動かす。相関式と比べるなら低乱れ度（<1%）で層流域を作るか、遷移モデル（γ-Reθ等）を使う
- 先端の解像 — 境界層が最も薄い先端近傍は流れ方向にも細分。先端の数値的な鈍りは板全体の平均hに響く
検証ベンチマークとしての使い方
平板は「CFD設定の健全性テスト」として最も費用対効果が高い問題です。局所 \( Nu_x \) 分布を相関式と重ね、①層流域で \( x^{-1/2} \) の傾きが再現されるか、②遷移位置が想定と合うか、③乱流域の水準が±10%程度に入るか、を確認します。ここで合わない設定（メッシュ・壁処理・乱流モデル）のまま実形状に進むと、実形状での不一致の原因が切り分けられなくなります。
実務適用の指針
相関式の3つの使い道
- 概算設計 — 放熱面積・冷却風量の初期見積り。CFD前に答えの桁を掴む
- CFDの検算 — 解析結果の平均hを相関式と比較し、2倍以上ずれたら設定を疑う
- 境界条件の供給 — 詳細に解かない面の熱伝達係数 \( h \) として熱伝導解析・1D熱回路網に与える
実物とのずれの主因を知っておく
| 要因 | 影響 | 対処 |
|---|---|---|
| 主流乱れ度が高い（実機は数%） | 遷移が早まり平均hが増加 | 混合平均式の臨界Re数を下げて再評価、または安全側で全面乱流式 |
| 未加熱助走区間がある | 温度境界層が遅れて発達しhが上がる | 助走区間補正付きの式を使用 |
| 表面粗さ | 乱流域のhを増加（滑面式は下限側） | 粗面用の摩擦係数系相関へ |
| 有限幅・端効果 | 3次元流れで局所的にhが変動 | 幅/長さ比が小さい場合はCFDで確認 |
| 物性の評価温度 | 温度差が大きいと物性誤差が直撃 | 膜温度評価を徹底、必要なら物性補正係数 |
「等温でも等熱流束でもない」実物への態度
実際の放熱面は、等温壁（係数0.332系）と等熱流束壁（0.453系）の中間にあります。両式の差は約36%——つまり壁条件の理想化だけでその程度の幅が本質的にあるということです。精度が必要な場合は、固体側の熱伝導と連成した共役熱伝達解析に進むのが正道で、相関式は連成結果のサニティチェックとして残します。
ツールでの取り扱い
CFDツールの設定要点
| ツール | 要点 |
|---|---|
| Ansys Fluent | Enhanced Wall Treatment（y+適応壁処理）。遷移はTransition SST（γ-Reθ）モデル。壁面熱流束レポートで局所Nuを抽出 |
| OpenFOAM | 壁関数はnutkWallFunction系＋alphatWallFunction。低Re型はy+<1メッシュで壁関数を外す。wallHeatFluxユーティリティで後処理 |
| STAR-CCM+ | All y+壁処理。遷移モデル（Gamma ReTheta）内蔵。派生パート機能で局所h分布を抽出 |
| COMSOL | 共役伝熱インターフェースで流体固体一体解析。層流域なら直接解ける |
| 1D/熱回路網ツール | 平板相関式が組込み（h自動計算）。適用範囲（Re・Pr・壁条件）の内蔵仮定を確認して使う |
ツール間比較は「hの定義」を揃えてから
熱伝達率は \( h = q_w / (T_{ref} - T_w) \) で定義されますが、基準温度 \( T_{ref} \) の取り方（主流温度・断熱壁温度・入口温度）がツールと相関式で異なることがあります。特に高速流（粘性加熱が効く場合）や内部流との混在モデルでは定義差が数十%の見かけの不一致を生みます。ツール間・相関式との比較は、定義に依存しない壁面熱流束 \( q_w \) そのもので行うのが安全です。
先端研究の動向
遷移予測の高度化
層流から乱流への遷移は平板熱伝達の最大の不確定要素で、局所変数だけで遷移を判定するγ-Reθ型モデルの改良、感度の高い環境（乱れ度・圧力勾配・粗さ）への拡張が続いています。LES・DNSによる遷移機構（バイパス遷移・Tollmien-Schlichting波）のデータベースが整い、RANS遷移モデルの校正基盤として使われています。実務的には「遷移位置の予測誤差＝平均hの誤差の主因」という構図は変わっておらず、遷移が絡む案件では遷移位置の感度解析を標準に組み込むのが堅実です。
非理想条件への拡張
実機に近い条件——高乱れ度主流・粗面・脈動流・衝突噴流との複合——での平板熱伝達データの蓄積と相関式の拡張が続いています。また共役効果（板内の熱伝導が壁温分布を作り、それが境界層に影響する）を最初から含めた共役Blasius型問題の解析解・数値解も整備され、「等温/等熱流束の2択」からの脱却が理論側でも進んでいます。
データ駆動の相関式構築
実験・CFDデータから記号回帰やニューラルネットワークで熱伝達相関を構築する研究が増えています。古典相関式より適用範囲を広げられる可能性の一方、外挿の危険性と物理的整合性（次元整合・極限挙動）の担保が課題で、「無次元数の枠組みは保ち、係数・指数をデータ駆動で最適化する」中間的なアプローチが実用的な落とし所として有力です。
トラブル対応
症状別の原因と対策
| 症状 | 考えられる原因 | 対策 |
|---|---|---|
| CFDのNuが層流相関式より数十%高い | 意図せぬ遷移（乱れ度過大）、y+不適合 | 流入乱れ度を下げる/遷移モデル使用。y+分布を確認 |
| CFDのNuが相関式より大幅に低い | 境界層の解像不足（層数・先端）、壁関数の適用域違反 | 境界層メッシュ再設計。低Re型へ切替 |
| 実験と合わない（CFDは相関式と合う） | 実験側の乱れ度・粗さ・端効果 | 実験条件の無次元数を再確認。不確かさ幅内かを評価 |
| ツール間でhが2〜3割違う | 基準温度の定義差 | 熱流束で比較。定義をドキュメントで確認 |
| 局所Nu分布が先端で発散的に大きい | 先端特異性（境界層厚さ→0）の離散化 | 先端近傍の値は評価から外す（理論上もNu_x→∞） |
| 手計算と桁が合わない | 物性の評価温度ミス、特性長の取り違え（x vs L） | 膜温度で物性再評価。局所式と平均式の混同を確認 |
迷ったら最初に遷移を疑う
CFDと相関式と実験がバラバラの値を出したとき、どこから手を付ければいいですか？
平板の不一致の筆頭容疑者は遷移位置だよ。層流と乱流でhは数倍違うから、遷移がどこで起きているかの認識が三者で揃っていなければ、値が合わないのは当然なんだ。診断は局所分布で行う——CFDなら壁面熱流束の \( x \) 方向分布、実験なら壁温分布を見て、hが急増する位置（遷移）を特定する。三者の遷移位置を揃えてから比べ直すと、たいてい残差は±10〜20%の「正常な不確かさ」に収まる。値の比較の前に、まず分布の形の比較——これが対流熱伝達の検証の鉄則だね。
なった
詳しく
報告