衝突噴流熱伝達 — CFDの典型的な問題と対策

これは衝突噴流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^+$ の値に応じて自動調整するので比較的安定だ。