DSMC-Neutrals　解析事例 - 極超音速希薄流解析

ウェーブフロント＞ CAEソリューション＞ DSMC-Neutlrals ＞極超音速希薄流シミュレーション

DSMC-Neutralsはモンテカルロ直接法（Direct Simulation Monte Carlo Method；DSMC法）を用いた3次元希薄気体解析ソフトウェアです．非構造メッシュを用いることにより，複雑な形状のシミュレーションが可能です．また，化学反応の計算も可能であり，真空チャンバー内の希薄なガス流れはもちろんのこと，化学蒸着（Chemical Vapor Deposition；CVD）などの半導体製造における薄膜生成のシミュレーションにも適しています．本記事を読んで希薄気体解析・ガス流れシミュレーションに少しでも興味をお持ちになりましたら 資料請求や 無料セミナーなど承りますので，いつでもお気軽にお問い合わせください．

極超音速希薄流シミュレーションHypersonic Rarefied Flow Simulation

以下に極超音速希薄流解析のシミュレーション結果を示す。流速は 7.5 km/s と音速の20倍である。化学反応計算として窒素・酸素の解離、再結合、原子交換反応を考慮している。物体表面は温度が非常に高いために、ほとんどの酸素が解離していることがわかる。流体モデルの解析ソフトでは計算できない極超音速のガス流れを再現することがDSMC-Neutralsでは可能である。また、DSMC-Neutralsの計算結果はDogra等の計算結果 [1] をよく再現している。

▲ シミュレーション条件

▲ モデル形状とメッシュ

衝突モデルCollision Model

弾性衝突
非弾性衝突
- 回転
  Larsen-Bolgnakke 分布関数より棄却法を用いて求める。
  N2, O2 の自由度： $2$
  緩和値 $Zr = 5$ とする。
- 振動
  量子化された Larsen-Bolgnakke 分布関数より棄却法を用いて求める。
  緩和値 $Z_v$ は下記式より計算される。 \begin{align*} Z_v = ( C_1 / T^w ) e^{C_2 T^{-1/3}} \end{align*} ここで定数 $C_1$ と $C_2$ は分子種依存の定数である。

化学反応モデルChemical Reaction Model

全衝突エネルギーとアレニウス方程式の定数より衝突確率を計算する Total collision energy (TCE) model

解離 - 再結合
- 反応式
  N2 + M <---> N + N + M
  O2 + M <---> O + O + M
  NO + M <---> N + O + M
  M は第3体粒子である。
原子交換
- 反応式
  NO + O <---> N + O2
  N2 + O <---> N + NO