強制対流のCFD解析

そのとおり。強制対流では外部から与えられた流れ（ファン、ポンプ、配管内圧力差など）が主体で、浮力は無視できるほど小さい。CFDでは速度場・圧力場を解いた上で、壁面の熱流束やNusselt数を正確に予測することが目標になる。

Nusselt数 $Nu$ は壁面の対流伝熱と流体の熱伝導の比で、伝熱性能の無次元指標だ。

管内流れでよく使われるのがDittus-Boelter式とGnielinski式だ。

平板上の強制対流層流では $Nu_x = 0.332 Re_x^{1/2} Pr^{1/3}$ 、乱流では $Nu_x = 0.0296 Re_x^{0.8} Pr^{1/3}$ が基本だ。円柱まわりの流れではChurchill-Bernstein式が使われる。

その通り。壁面の熱伝達を正確に予測するには、温度境界層を十分に解像する必要がある。2つの主要アプローチがある。

(1) Low-Reynolds number approach: 壁面第一セルを $y^+ \approx 1$ に配置し、粘性底層から直接解く。SST k-ωモデルと組合せるのが標準。

定量的なNu数予測が必要なら $y^+ \approx 1$ のLow-Re approachを推奨する。壁面関数は流れの傾向把握には使えるけど、Nu数を10〜20%過大/過小評価するリスクがある。特に剥離を伴う流れ（後ろ向きステップ、ブラフボディ）では壁面関数の精度が悪化するよ。

完全発達した円管内乱流のNu数予測なら、標準k-εでもDittus-Boelter式と5%以内で一致する。ただし入口領域や急縮小・急拡大のある管路ではSST k-ωのほうが安全だ。遷移領域（$Re \approx 2300$〜$10000$）を扱うならTransition SSTモデルが必要になる。

そうだね。アルミ押出しフィン、ピンフィン、液冷コールドプレートなど、いずれも強制対流による放熱が主体だ。設計変数としてはフィンピッチ、フィン高さ、流路幅、流量などがあり、CFDで熱抵抗 $R_{th} = \Delta T / Q$ を評価してパラメトリックスタディを行う。

(1) フィンの対称性を利用して1ピッチ分だけモデル化する（周期境界条件を使用）。(2) 入口に一様流速、出口に圧力出口を設定。(3) フィン底面に一定熱流束を付与。(4) SST k-ωで定常計算。(5) フィン底面の平均温度から熱抵抗を算出。

短いヒートシンク（$L/D_h < 20$）では入口効果が無視できないので、全体モデルが必要だ。長いヒートシンクで完全発達流れが主体なら1ピッチ+周期条件で十分。Fluentではperiodic conditionにmass flow rateを指定できるし、STAR-CCM+でもperiodic interfaceで設定可能だ。

液冷コールドプレートは典型的なCHT問題だ。アルミや銅のプレートに微細流路（マイクロチャネル）を加工し、水やクーラントを流す。流路の水力直径 $D_h$ が1mm以下になるとRe数が低くなり層流になることも多い。その場合は層流モデルで十分だ。

非常に重要。ポンプの動力制約があるから、熱抵抗と圧力損失のトレードオフを最適化する必要がある。パレートフロント（多目的最適化）の考え方が有効で、STAR-CCM+のDesign ManagerやFluentのDesign ExplorationツールでDOEベースの最適化ができるよ。

直線流路ならDarcy-Weisbach式 $\Delta p = f \cdot (L/D_h) \cdot (\rho u^2/2)$ で検証できる。層流なら $f = 64/Re$ 、乱流ならMoody線図の値と比較する。CFDの圧損がこれらと10%以上ずれていたらメッシュか設定に問題がある可能性が高い。

定常・非圧縮性・乱流の標準ケースでの推奨を挙げよう。

最終的な解の精度は同等だけど、Coupledのほうが収束が速い。特にCHT問題や浮力が効く場合はCoupledが安定する。SIMPLEは反復ごとのメモリ使用量が少ないので、メモリ制約がある場合に選択する。

ソルバーはbuoyantSimpleFoam（定常、浮力項あり）またはsimpleFoam+エネルギー方程式追加が標準。乱流モデルはconstant/turbulencePropertiesでkOmegaSSTを指定する。壁面境界条件はkqRWallFunction、omegaWallFunction、nutUSpaldingWallFunctionの3点セットだ。

divスキームはbounded Gauss linearUpwindが安定と精度のバランスが良い。Gauss linearは精度が高いけど非定常流れでは振動しやすい。gradスキームはGauss linearで、limiterを入れる場合はcellLimitedを使おう。

Physics ModelでSegregated Flowを選択し、エネルギーモデルをSegregated FluidTemperatureで有効化する。乱流はSST k-ω + All y+ Wall Treatment（$y^+$ によらず自動対応）が推奨。Mesh ModelではPolyhedral Mesher + Prism Layer Mesherの組合せが品質と効率の両立に優れているよ。

$Re \approx 2300$〜$10000$ の遷移領域はRANS乱流モデルの弱点だ。通常のk-εやk-ωでは遷移の位置を予測できない。Ansys FluentのTransition SSTモデル（$\gamma$-$Re_{\theta}$モデル）がこの問題に対応しており、遷移の開始位置と長さを実験的相関式に基づいて予測する。

STAR-CCM+にもGamma Transition Modelがある。OpenFOAMではkOmegaSSTLM（Langtry-Menter遷移モデル）が利用可能だ。いずれも自由流れ乱流強度 $Tu$ と自由流れ粘性比 $\mu_t/\mu$ の入口条件が遷移位置に大きく影響するので、これらを実験条件に合わせて正確に設定することが重要だよ。

もちろん。内燃機関の冷却通路やガスタービンの内部冷却ではリブ付き流路が標準的だ。リブの高さ $e$、ピッチ $P$、角度 $\alpha$ が設計パラメータで、Nu数の増大比（smooth管比）とfriction factorの増大比の関係を評価する。

リブ付き流路のCFDで注意すべきなのは、リブ間の再循環域を正確に捉えること。SST k-ωは再付着点の位置を比較的よく予測するが、標準k-εは再付着長さを過小評価する傾向がある。メッシュはリブ高さ方向に20層以上配置し、リブ間を流れ方向に少なくとも40分割するのが目安だ。

トポロジー最適化がホットなトピックだ。Ansys Fluent 2023R2以降ではAdjoint Topology Optimizationが使えるようになり、流路形状を目的関数（熱抵抗最小化+圧損制約など）に基づいて自動的に進化させられる。COMSOL MultiphysicsのTopology Optimization Moduleも同様の機能を持っている。得られた形状は有機的で従来設計では思いつかないような構造になることが多いよ。

3Dプリンティング（SLMやDMLS）なら製造可能だ。実際にAdditive Manufacturingを前提としたトポロジー最適化ヒートシンクの研究事例が増えている。製造制約（最小壁厚、オーバーハング角度など）を最適化に組み込むのが実用上のポイントだよ。

まずメッシュと流れの状態を確認しよう。完全発達しているか？ 入口から十分な助走距離（乱流なら $L > 10D$ 程度）をとっているか？ 助走区間ではNu数が発達値より高くなるのは正常だ。

次に壁面処理を確認しよう。壁面関数を使っていて $y^+$ が5〜30の「バッファ層」に入っていないか？ バッファ層は壁面関数の適用範囲外で、熱伝達の予測精度が悪化する。Fluentなら$y^+$のコンター図を確認して、30以上（壁面関数用）か1以下（Low-Re用）のどちらかに統一すべきだ。

メッシュの非対称性が原因であることが多い。テトラメッシュでは面の向きがランダムなので、粗いメッシュだと非対称な数値拡散が生じる。ヘキサメッシュに変更するか、メッシュを十分に細かくすれば改善する。また、Second Order Upwindを使っているか確認しよう。First Orderだと数値拡散が大きく非対称になりやすい。

温度差が大きい場合（$\Delta T > 20$度C程度）は入れるべきだ。特に粘度 $\mu$ は温度依存性が大きく、壁面近傍の速度勾配ひいてはNu数に直結する。Fluentではmaterial propertiesでpolynomialやpiecewise-linearで温度依存性を入力できる。OpenFOAMではtransportProperties内のtransportModelをSutherlandやpolynomialに設定する。

空気や窒素を冷媒とする場合、Ma数が0.3以下（圧縮性効果が5%以下）なら非圧縮性近似で十分。ただし大きな温度差がある場合はBoussinesq近似ではなく理想気体として密度の温度依存性を扱ったほうがよい。FluentのIncompressible Ideal Gas設定がこの用途に適している。