Direct simulation Monte Carlo法

直接シミュレーションモンテカルロ法（DSMC法）

直接シミュレーションモンテカルロ法（DSMC法）は、希薄気体の流れを扱う数値的手法で、航空工学の名誉教授であるGraeme A. Birdによって開発されました。この方法は、ボルツマン方程式を解決するために確率論的アプローチを用いるものです。DSMC法は、特に有限クヌーセン数を持つ流れのシミュレーションに効果的で、様々な産業や研究で広く利用されています。たとえば、スペースシャトルの再突入時の空気力学から、ターボ分子ポンプや真空放電に関する解析、さらにはMEMS（微小電子機械システム）の設計に至るまで、適用範囲は多岐にわたります。

希薄気体とDSMC法

DSMC法の基本的な前提として、気体分子の平均自由行程が代表的な長さスケールと同程度、あるいはそれ以上であることが挙げられます。これは、クヌーセン数が1以上である場合を想定しています。超音速や極超音速の流れでは、この流れの希薄度がクヌーセン数とマッハ数の積、もしくはTsienのパラメーターによって特徴づけられます。これらの希薄流れにおいて、ナビエ・ストークス方程式は不正確になることが多く、ボルツマン方程式がより適切に流れを表現します。DSMC法では、シミュレーション粒子を使ってボルツマン方程式に従う流れをモデル化します。また、これを連続流領域まで拡張するモデルもあり、ナビエ・ストークス方程式との比較も示されています。

サンプル粒子

希薄気体ではあっても、現実の気体分子の数は膨大です。すべてをシミュレーションすることは不可能であり、したがってサンプル粒子と呼ばれる仮想の粒子が利用されます。これらのサンプル粒子は、実際の分子の特性を代表する存在で、その数は「重み」として表されます。DSMC法では、重みの小さいサンプル粒子を多数扱うことで、シミュレーションの精度を高めますが、その分処理時間も増加するというトレードオフがあります。非定常流の特性を正確に再現するためには、物理的な時間と空間の変化に基づいてシミュレーションを進めなければなりません。

衝突

DSMC法では、分子の移動と衝突が分離して考えられ、衝突フェーズは平均衝突時間より短い時間間隔で処理されます。分子間の衝突と表面との衝突は、確率的なモデルによって計算されます。一般的には、剛体球モデルや可変剛体球モデルなど、さまざまなモデルが衝突の解析に使用されます。これらのモデルは、衝突確率の計算や衝突する粒子の選定に役立ちます。

DSMCソフトウェアとツール

DSMC法には、いくつかのシミュレーションツールが存在します。Bird教授によるオリジナルのプログラムであるDS1V、DS2V、DS3Vは、GUIを備えており、シミュレーションの構築と後処理に役立ちます。商用パッケージとしては、Bernd Rose博士が開発したPI-DSMCや、サンディア国立研究所のSPARTA、NASAが制作したMAPなどがあります。これらのプログラムは、それぞれ異なる特徴や機能を持ち、研究者やエンジニアのニーズに応えています。

様々な国や企業によって開発された製品もあり、特に国内では、日本のウェーブフロント社が提供するDSMC-Neutralsや、韓国のMetariver Technology社によるSAMADII/SCiVなどがあります。これらのツールは、各分野におけるDSMC法の利用を促進しています。

参考資料

最後に、DSMC法についてさらに深く知りたい方には、Graeme A. Birdの著作『The DSMC Method』をおすすめします。これは、この方法について詳細に解説した教科書です。DSMC法の理解を更に深めるために、関連する教材やチュートリアルも活用すると良いでしょう。

もう一度検索