[論文レビュー] Adaptive Multi-Dimensional Particle In Cell

Giovanni Lapenta|ArXiv.org|Jun 4, 2008

Electrohydrodynamics and Fluid Dynamics被引用数 3

ひとこと要約

本稿では、空間グリッドの適応と粒子再配置を組み合わせることで、マルチスケールプラズマ現象を効率的にシミュレートできる、適応的多次元粒子法（PIC）を提案する。事後誤差指標を用いたグリッド細分化と動的粒子投入により、衝突なしシャocks やダスト帯電といった関心領域の精度を維持しながら計算コストを削減する。非適応的シミュレーションと比較して、数値的安定性と統計的精度の向上が実証された。

ABSTRACT

Kinetic Particle In Cell (PIC) methods can extend greatly their range of applicability if implicit time differencing and spatial adaption are used to address the wide range of time and length scales typical of plasmas. For implicit differencing, we refer the reader to our recent summary of the implicit moment PIC method implemented in our CELESTE3D code [G. Lapenta, Phys. Plasmas, 13, 055904 (2006)]. Instead, the present document deals with the issue of PIC spatial adaptation. Adapting a kinetic PIC code requires two tasks: adapting the grid description of the fields and moments and adapting the particle description of the distribution function. Below we address both issues. First, we describe how grid adaptation can be guided by appropriate measures of the local accuracy of the solution. Based on such information, grid adaptation can be obtained by moving grid points from regions of lesser interest to regions of higher interest or by adding and removing points. We discuss both strategies. Second, we describe how to adapt the local number of particles to reach the required statistical variance in the description of the particle population. Finally two typical applications of adaptive PIC are shown: collisionless shocks and charging of small bodies immersed in a plasma.

研究の動機と目的

• プラズマにおける広範な時間スケールおよび長さスケールを解像する上で、標準的な明示的PIC法の計算非効率性を解消すること。
• 局所的勾配や動的粒子分布を捉えることができない固定グリッド・固定粒子のPICシミュレーションの制限を克服すること。
• グリッド適応と粒子再配置を統合したフレームワークを構築し、暗黙的PIC法の精度を向上させ、ノイズを低減すること。
• 宇宙環境における衝突なしシャocks やダスト帯電といった複雑なプラズマ系の正確なシミュレーションを可能にすること。
• 適応的解像度と粒子制御を通じて、並列実装におけるエネルギー保存則と負荷分散を確保すること。

提案手法

• 局所的解の精度に基づいてグリッド適応を誘導するため、事後誤差指標を用いる。高解像度を要する領域を特定する。
• 2つのグリッド適応戦略を採用する：グリッド移動（高関心領域へポイントを移動）とグリッド細分化（必要に応じてポイント追加）。
• 統計的精度を維持しノイズを低減するために、細胞あたりの粒子数を動的に調整する粒子再配置を適用する。
• ダスト粒子の境界条件をモデル化するため、感受性および減衰項を含む埋め込み境界法を用いる。
• 精度と計算コストのバランスを取るために、移動メッシュ法（MMA）と適応メッシュ細分化（AMR）を実装する。
• 安定性制約を緩和し、適切な場合に大きな時間ステップと粗いグリッドを許容するために、暗黙的時間差分法を統合する。

実験結果

リサーチクエスチョン

• RQ1 多次元PICシミュレーションにおいて、どのように局所的精度指標に基づいてグリッド適応を誘導できるか？
• RQ2 適応的PIC法において、グリッド移動とグリッド細分化のどちらがより効果的か？
• RQ3 時間的に変化する粒子密度領域において、どのように粒子再配置が統計的精度を維持できるか？
• RQ4 適応的グリッドおよび粒子制御が、並列PICコードにおけるエネルギー保存則と負荷分散に与える影響は何か？
• RQ5 適応的PIC法は、衝突なしシャックス やダスト帯電といった複雑なプラズマ現象のシミュレーションをどのように改善するか？

主な発見

• 事後誤差指標に基づくグリッド適応により、シャックスの前線やシースなど勾配が急な領域の解像度が顕著に向上する。
• 粒子再配置により、細胞あたりの粒子数が一定に保たれ、低密度領域における統計的ノイズの増加が防止される。
• 適応的グリッドと粒子制御の組み合わせにより、ダスト帯電の安定的かつ正確なシミュレーションが可能となり、定常状態での電荷蓄積が達成された。
• 粒子再配置が行われない場合、ダスト周辺の細胞あたりの粒子数は時間経過とともに減少し、精度が低下し、平衡状態への収束が不可能になる。
• 固定グリッド・固定粒子のシミュレーションと比較して、特に空間的・時間的スケール分離が顕著な系において、適応的PIC法は計算効率と精度の両方を向上させる。