[論文レビュー] Rubrics for Charge Conserving Current Mapping in Finite Element Particle in Cell Methods

Zane D. Crawford, Scott O’Connor+2 more|arXiv (Cornell University)|Jan 28, 2021

Plasma Diagnostics and Applications被引用数 4

ひとこと要約

本論文は、de Rham複体写像と離散Hodge作用素を活用することで、有限要素粒子法（FE-PIC）における電荷保存を保証する一般化された数学的基準を確立する。電流マッピングがエッジ基底と面基底の有限要素基底の双対性を尊重する場合、電荷保存が達成されることを示し、数値結果によって検証された、本質的に電荷保存を満たす新しいFE-PICスキームを提示する。

ABSTRACT

Modeling of kinetic plasmas using electromagnetic particle in cell methods (EM-PIC) is a problem that is well worn, in that methods developed have been used extensively both understanding physics and exploiting them for device design. EM-PIC tools have largely relied on finite difference methods coupled with particle representations of the distribution function. Refinements to ensure consistency and charge conservation have largely been an ad-hoc efforts specific to finite difference methods. Meanwhile, solution methods for field solver have grown by leaps and bounds with significant performance metrics compared to finite difference methods. Developing new EM-PIC computational schemes that leverage modern field solver technology means re-examining analysis framework necessary for self-consistent EM-PIC solution. In this paper, we prescribe general rubrics for charge conservation, demonstrate how these are satisfied in conventional finite difference PIC as well as finite element PIC, and prescribe a novel charge conserving finite element PIC. Our effort leverages proper mappings on to de-Rham sequences and lays a groundwork for understanding conditions that must be satisfied for consistency. Several numerical results demonstrate the applicability of these rubrics.

研究の動機と目的

• 電磁気的粒子法（EM-PIC）シミュレーションにおける電荷保存を保証する一般化された、ソルバーに依存しない基準を開発すること。
• 有限要素法における電流マッピングが電荷保存を保つために満たすべき数学的条件を同定すること。
• 有限差分PICから現代の有限要素場ソルバーへの電荷保存原理の拡張すること。
• FE-PICスキームにおける一貫性があり発散なしの電流定着を構築するための体系的枠組みを提供すること。
• 数値的検証を通じて、提案された基準の適用可能性と正確性を示すこと。

提案手法

• de Rham複体における離散微分形式の双対性を要請することで電荷保存を定式化し、電流マッピングがHodgeスター作用素および回転関係を尊重することを保証する。
• 電場と磁場にエッジ基底（1形式）および面基底（2形式）の有限要素基底関数を用いて、離散Maxwell方程式を導出する。
• 粒子軌道の線積分を用いて電流の自由度を定義し、Ampereの法則およびFaradayの法則と整合性を保つ。
• 離散Hodgeスター作用素を用いて原始メッシュと双対メッシュの間をマッピングし、正確な場と電流の補間を可能にする。
• de Rham系列を用いて、離散的な回転および発散作用素が連続的対応物と同一の代数的恒等式を満たすことを保証する。
• 時間経過にわたる電荷保存およびMaxwell方程式との整合性を確認する数値実験により、スキームを検証する。

実験結果

リサーチクエスチョン

• RQ1 有限要素PIC法における電流マッピングが電荷保存を保証するためには、どのような数学的条件を満たす必要があるか？
• RQ2 任意の有限要素場ソルバーに対して、電荷保存型電流定着を体系的に導出する方法は何か？
• RQ3 de Rham複体写像と離散Hodge作用素は、どのようにしてFE-PICにおける一貫性のある電荷保存を可能にするか？
• RQ4 提案されたFE-PICスキームは、従来のFDTD-PICと比較して、電荷保存性および整合性においてどのように異なるか？
• RQ5 提案された基準は、有限要素法を超えた他の場ソルバータイプに対しても一般化可能か？

主な発見

• 本論文は、de Rham複体を介してエッジ基底と面基底の有限要素基底関数の双対性を尊重する電流マッピングがなされれば、FE-PICにおける電荷保存が保証されることを確立した。
• 提案された手法により、各タイムステップで離散的連続の式が満たされ、不自然な電荷生成が排除された。
• 数値結果により、長時間にわたるシミュレーションにおいても、機械精度に近い電荷保存が維持されていることが確認された。
• 構造的グリッドに適用した場合、本手法は古典的なFDTD-PICの結果を回復し、既存のスキームとの整合性を検証した。
• de Rham系列と離散Hodge作用素の使用により、回転および発散作用素が連続的状態と同一の代数的恒等式を満たすことが保証された。
• この枠組みにより、後処理や発散クリーニングを必要としない、新たな本質的電荷保存型FE-PICスキームの構築が可能になった。