[論文レビュー] Direct simulation Monte Carlo schemes for Coulomb interactions in plasmas

Giacomo Dimarco, Russell Caflisch+1 more|arXiv (Cornell University)|Oct 1, 2010

High-Energy Particle Collisions Research被引用数 12

ひとこと要約

この論文は、ボビレフ＝ナンブ近似を用いて、ボルツマン散乱作用素の近似を適用することで、希薄ガスの直接シミュレーション・モンテカルロ（DSMC）スキームを、長距離相互作用を示すプラズマに一般化する。この手法により、時間ステップを大きくとることができ、拒否ベースの衝突サンプリングを回避し、元のBNスキームと同等の精度を維持しながら、計算効率を向上させ、収束速度を向上させ、統計的ゆらぎを低減する。これは、バード型の定式化において顕著である。

ABSTRACT

We consider the development of Monte Carlo schemes for molecules with Coulomb interactions. We generalize the classic algorithms of Bird and Nanbu-Babovsky for rarefied gas dynamics to the Coulomb case thanks to the approximation introduced by Bobylev and Nanbu (Theory of collision algorithms for gases and plasmas based on the Boltzmann equation and the Landau-Fokker-Planck equation, Physical Review E, Vol. 61, 2000). Thus, instead of considering the original Boltzmann collision operator, the schemes are constructed through the use of an approximated Boltzmann operator. With the above choice larger time steps are possible in simulations; moreover the expensive acceptance-rejection procedure for collisions is avoided and every particle collides. Error analysis and comparisons with the original Bobylev-Nanbu (BN) scheme are performed. The numerical results show agreement with the theoretical convergence rate of the approximated Boltzmann operator and the better performance of Bird-type schemes with respect to the original scheme.

研究の動機と目的

• 長距離クーロン相互作用を有するプラズマにおける、従来のボルツマンに基づくスキームが計算的に非現実的であるため、効率的なモンテカルロ手法の不足に応える。
• 古典的ボビレフ＝ナンブ（BN）スキームの限界を克服する。このスキームは非常に小さな時間ステップを必要とし、時間離散化誤差に苦しむ。
• 元々短距離力の希薄ガスに用いられるバードおよびナンブのDSMCアルゴリズムを、近似された散乱作用素を介してクーロン相互作用に一般化する。
• 提案されたスキームが理論的収束速度を維持しながら、数値的安定性を向上させ、統計的分散を低減することを示す。
• 時間離散化と衝突近似を分離することで、時間ステップを大きくとることができ、コストの高い受容・拒否手順を回避する。

提案手法

• 小角度散乱の極限において有効な弱い二体衝突をモデル化するため、ボルツマン散乱作用素のボビレフ＝ナンブ近似を、クーロン系に適応する。
• 近似された散乱作用素に基づく直接シミュレーション・モンテカルロ（DSMC）スキームを構築し、各時間ステップで全粒子が衝突するが、拒否サンプリングを伴わない。
• 時間積分と衝突近似を分離し、近似誤差を制御するパラメータτを用いることで、精度を保ちつつ、より大きな時間ステップを可能にする。
• 比較のため、バード型とボビレフ＝ナンブ型の両方のスキームを実装し、バード型はより安定した時間離散化手法を用いる。
• 近似された散乱作用素の数値的収束性を検証するため、τ → 0 のときの理論的収束率O(τ)をテストする。
• 統計的分散と収束行動を評価するため、高密度粒子数（N = 2×10⁶）と多数の実現（M = 10–1500）を用いた基準解を用いる。

実験結果

リサーチクエスチョン

• RQ1 希薄ガスの古典的DSMCアルゴリズムは、ボビレフ＝ナンブ近似を用いることで、クーロン相互作用を示すプラズマに一般化可能か？
• RQ2 提案されたDSMCスキームは、近似された散乱作用素の理論的収束率を達成するか？また、元のBNスキームと比較してどのように異なるか？
• RQ3 粒子数と時間ステップサイズを変化させた場合、バード型とBN型スキームの統計的ゆらぎと数値的安定性にはどのような差が生じるか？
• RQ4 時間離散化誤差はBNスキームの性能にどのような影響を及えるか？また、バード型スキームはこの問題をどのように軽減するか？
• RQ5 近似された散乱作用素のおかげで、精度を損なわずにDSMCシミュレーションにおける時間ステップを大きくとることができるか？

主な発見

• 近似された散乱作用素は、τ → 0 のとき、期待されるO(τ)の収束速度を示し、理論的基盤の妥当性を確認する。
• バード型スキームは、特に大きなτのとき、時間離散化誤差が小さいため、元のBNスキームよりも統計的ゆらぎが小さい。
• τが小さい場合（例：τ = 0.0625）、バード型とBN型スキームは同じ解に収束する。これは、近似誤差が小さい極限において両者が等価であることを示している。
• 両スキームの収束率r(τ)は、τを小さくするにつれて漸近的に1に近づき、理論的収束挙動が正当化される。
• 提案された手法により、従来のスキームよりも大きな時間ステップが可能となり、元のボルツマン散乱積分に必要な小さな時間ステップを回避できる。
• このスキームにより、各時間ステップで全粒子が衝突するようになり、計算効率が向上し、高価な受容・拒否ステップが不要になる。