衝突噴流熱伝達

衝突噴流（jet impingement）はノズルから噴出したジェットが壁面に衝突する流れで、衝突点（stagnation point）近傍で非常に高い熱伝達係数が得られる。通常の管内強制対流の2〜3倍の熱伝達率を達成できるため、ガスタービン翼内部冷却、鋼板焼入れ、電子部品冷却などに幅広く使われているよ。

流れは3つの領域に分かれる。(1) 自由噴流領域（free jet region）：ノズルから壁面に向かう領域でポテンシャルコアが存在。(2) 衝突域（impingement zone）：壁面近傍で流れが方向転換する領域。(3) 壁面噴流領域（wall jet region）：壁面に沿って放射状に広がる領域。

ある。Martin（1977）の相関式が広く使われている。単一円形ノズルの停滞点Nu数は

ポテンシャルコアの長さ（通常4〜6D）を超えると噴流が拡散して衝突時の速度が低下し、Nu数が減少する。逆に $H/D < 4$ でもconfined geometry（閉じた空間）ではcross-flow effect（使用済みの噴流が新しい噴流に干渉）が出て性能低下することがある。

そうだ。配列ジェット（array impingement）ではジェット間の干渉とcross-flowの影響が重要になる。孔間隔 $S/D$が小さいほど面積平均のNu数は上がるけど、cross-flowが強くなって下流側のジェットの冷却性能が劣化する。典型的には $S/D = 4$〜$8$ の範囲で設計する。Florschuetz et al.（1981）の相関式が配列噴流の標準的な参照データだよ。

そのとおり。衝突噴流は乱流モデルのベンチマーク問題として有名で、多くのモデルが衝突域のNu数を過大予測する。最大の原因は衝突域での乱流エネルギー生成項の過大評価だ。

具体的には、標準k-εモデルは停滞点Nu数を実験値の2倍以上過大予測することがある。v2f モデルや SST k-ω はかなり改善するが、それでも20〜30%の過大予測が残ることがある。$\omega$方程式ベースのモデルのほうがk-ε系より一般的に良好だ。

文献調査の結果、v2f モデル（Fluent: $\overline{v^2}$-$f$ モデル）が衝突噴流で最も実験値に近い予測を出す傾向がある。次いで SST k-ω。ただし v2f はFluentとOpenFOAM（v2f turbulence model）で利用可能だがSTAR-CCM+には標準実装されていない。

衝突域では壁面垂直方向に $y^+ < 1$ を確保し、壁面平行方向にもノズル直径 $D$ に対して $\Delta r / D \approx 0.02$〜$0.05$ の解像度が必要だ。ノズル出口から壁面までの間にも十分なセルを配置して噴流の発達を解像する。典型的な2D軸対称計算でも5万〜20万セル、3D配列噴流なら数百万セルが必要になる。

単一円形ノズルの場合はaxisymmetric（軸対称）モデルが非常に効率的だ。FluentでもSTAR-CCM+でもaxisymmetric solverが用意されている。ただしノズルが矩形の場合やcross-flowがある場合は3D必須だよ。

研究レベルでは盛んに使われている。噴流のシアー層に発生するKelvin-Helmholtz不安定性と、衝突域での大規模渦構造がNu数の非定常変動を引き起こす。LESはこの渦構造を直接解像できるので、時間平均Nu数分布がRANSより実験に近い。DES（Detached Eddy Simulation）やSBESも中間的な選択肢として有用だ。

高密度のデータセンターや次世代パワー半導体では、空冷の限界を超えて液体マイクロジェットが検討されている。直径0.5mmのノズル配列でチップ表面に直接噴射し、$Nu \sim 100$〜$500$ を達成する研究が進んでいるよ。

まず単一ノズルの軸対称問題で乱流モデルの検証を行う。Baughn & Shimizu（1989）や Cooper et al.（1993）の実験データが標準的な検証データセットだ。停滞点Nu数と径方向Nu数分布の両方を比較する。

配列ジェットの場合はFlorschuetz et al.の実験データと比較する。CFD結果は面積平均Nu数に加えて、ジェット中心線上のNu数分布、ジェット間の最小Nu数も評価ポイントだ。

(1) まずメッシュ感度分析で数値誤差を排除する。(2) 異なる乱流モデル（v2f、SST k-ω、RSM）を試す。(3) それでも合わなければ壁面処理（壁面関数 vs Low-Re）を変更する。(4) 最終手段としてLESやDESで非定常計算を行う。実務では「全ての条件で20%以内」を目標にすることが多いよ。

定常・軸対称の単一ノズルケースを例にしよう。(1) AxisymmetricモデルをMesh Setupで選択。(2) 乱流モデルはViscous Models > k-omega > SSTを選択。可能ならv2fモデル（Viscous Models > Transition > v2f）を試す。(3) Enhanced Wall Treatmentを有効化。(4) エネルギー方程式をONにする。

境界条件は、ノズル入口にvelocity inlet（噴流速度と温度を指定）、壁面にconstant heat fluxまたはconstant temperatureを設定。出口はpressure outlet。拘束壁がない場合（自由噴流衝突）は、上方と側方にpressure outletを配置する。

ゼロ速度場で開始すると収束に時間がかかることがある。噴流方向に一様速度を初期値として与えると収束が速い。Hybrid Initializationで自動設定してもよい。

Physics ModelsでSteady + Segregated Flow + Segregated Fluid Temperature + Turbulent + SST k-ωを選択。All y+ Wall Treatmentが有効になっていることを確認する。衝突噴流では壁面近傍のメッシュが特に重要なので、Prism Layer Mesherで壁面総厚さと層数を十分に設定する。

STAR-CCM+にはv2fモデルの標準実装がないので、衝突域の過大予測が気になる場合はElliptic Blending k-epsilon（EB k-ε）モデルを検討しよう。これはv2fと同等の壁面近傍処理を持つ。

ソルバーはsimpleFoamまたはbuoyantSimpleFoam。乱流モデルはconstant/turbulencePropertiesでkOmegaSSTを指定。v2fモデルはOpenFOAM v2212以降のv2f classで利用可能だ。壁面境界条件はnutUSpaldingWallFunction、kqRWallFunction、omegaWallFunctionを使用する。

blockMeshDictでwedge形状（角度5度程度）のメッシュを作成し、前面と背面にwedge境界条件を指定する。あるいはextrudeMeshDictでaxisymmetricオプションを使う方法もある。中心軸にはempty（2D）ではなくaxisの境界条件を設定するのを忘れないように。

回転はコリオリ力と遠心浮力の2つの効果で衝突噴流の熱伝達に影響する。回転数 $Ro = \Omega D / u_j$ が0.1以上になるとNu数分布が非回転時と有意に異なってくる。圧力面側ではコリオリ力がジェットを壁面に押し付けてNu数が増加し、吸い込み面側では引き離してNu数が減少する。

回転座標系（Rotating Reference Frame / Moving Reference Frame）を使用する。FluentのFrame Motionまたはsliding mesh、STAR-CCM+のRotating Physics、OpenFOAMのSRFSimpleFoamなどで設定できる。回転座標系ではコリオリ力項と遠心力項がソース項として運動方程式に追加される。

沸騰衝突噴流はさらに複雑だ。核沸騰領域では気泡の生成と離脱が熱伝達を大幅に増強し、$h$ が単相の10倍以上になることもある。CHF（臨界熱流束）を超えると膜沸騰に遷移して伝熱が急激に悪化する。

VOF法やEulerian多相流モデルで気液界面を追跡し、蒸発・凝縮モデルを組合せることで近似的に解ける。FluentのEvaporation-Condensation ModelやSTAR-CCM+のBoiling ModelがLee modelベースの蒸発/凝縮を扱う。ただし気泡の核生成を直接解像するのは計算コストが膨大で、実務ではWall Boiling Model（RPI model）のような壁面沸騰モデルを使うのが現実的だよ。

IBMやIntelの研究グループが発表しているマイクロジェット冷却では、ノズル直径 $D = 50$〜$500$μm、$H/D = 2$〜$5$ の配列で1000 W/cm2超の除熱を目指している。CFDでは層流から乱流遷移が重要で、$Re_D = 500$〜$5000$ の範囲をカバーする。この領域ではTransition SSTモデルが有用だ。

シリコンチップ基板にTSV（Through-Silicon Via）を加工してジェットノズルとする「integrated microjet cooling」も研究されており、3D-IC（積層IC）の冷却ソリューションとして注目されている。このような微細構造のCFDでは、壁面粗さと滑り速度（slip velocity）の影響も検討が必要になるよ。

これは衝突噴流CFDで最もよくある問題だ。原因は乱流モデルが衝突域の乱流エネルギー生成を過大評価していること。対策を優先順に並べよう。

(1) 乱流モデルを変更する: 標準k-ε → SST k-ω → v2f（Fluent）/ EB k-ε（STAR-CCM+）の順に試す。(2) Kato-Launderの修正を適用する: Fluentでは乱流モデルのOptionsでProduction Limiterを有効化。STAR-CCM+ではRealizability修正をONにする。(3) RSM（Reynolds Stress Model）を試す: 計算コストは上がるが、衝突域の異方性を直接扱える。

標準的な乱流エネルギー生成項 $P_k = \mu_t S^2$ を $P_k = \mu_t S \Omega$ に置き換える修正だ（$S$ は歪み速度、$\Omega$ は渦度）。衝突域では $S$ は大きいが $\Omega$ は小さいので、生成項が抑制される。単純だが効果が大きいよ。

cross-flow（使用済み噴流が溜まって新しいジェットに干渉する効果）が正しくモデル化されていない可能性がある。計算領域の出口境界条件を確認しよう。出口が近すぎるとcross-flowの蓄積が不十分になる。また、RANSはcross-flowによるジェット偏向を過小評価する傾向がある。配列噴流の正確な予測にはLESが有効だ。

衝突噴流は壁面近傍の速度・温度勾配が急峻なので、メッシュ依存性が出やすい。3水準以上のメッシュでRichardson外挿を行い、GCIが5%以内になるまでメッシュを細かくすべきだ。特に衝突点近傍の壁面平行方向のメッシュが不足していることが多い。$\Delta r / D < 0.03$ を衝突域で確保するのが目安だよ。

$y^+$ が1を大幅に下回ると（$y^+ < 0.1$）、ω方程式の壁面境界値が非常に大きくなり数値的に不安定になることがある。Fluentではωの under-relaxation factor を0.6程度に下げてみよう。STAR-CCM+のAll y+ Wall Treatmentは $y^+$ の値に応じて自動調整するので比較的安定だ。