[論文レビュー] A Dispersion Minimized Mimetic Method for Cold Plasma

Vrushali A. Bokil, Vitaliy Gyrya+1 more|arXiv (Cornell University)|Apr 5, 2016

Electromagnetic Simulation and Numerical Methods参考文献 1被引用数 3

ひとこと要約

本稿では、冷たいプラズマ中のマクスウェル方程式に対して、指数型時刻差分法（ETD）と新規のm-適応最適化を用いた分散最小化型ミメティック有限差分（MFD）法を提案する。空間的に一般化された質量アンローディングMFDスキームにm-適応を適用し、時間的にETDを用いることで、4次精度の数値分散誤差を達成し、標準的な2次精度スキームと比較して誤差を2桁小さくした。この結果は数値的および理論的に検証された。

ABSTRACT

In this paper we consider the lowest edge-based mimetic finite difference (MFD) discretization in space for Maxwell's equations in cold plasma on rectangular meshes. The method uses a generalized form of mass lumping that, on one hand, eliminates a need for linear solves at every iteration while, on the other hand, retains a set of free parameters of the MFD discretization. We perform an optimization procedure, called m-adaptation, that identified a set of free parameters that lead to the smallest numerical dispersion. The choice of the time stepping proved to be critical for successful optimization. Using exponential time differencing we were able to reduce the numerical dispersion error from second to fourth order of accuracy in mesh size. It was not possible to achieve this order of magnitude reduction in numerical dispersion error using the standard leapfrog time stepping. Numerical simulations independently verify our theoretical findings.

研究の動機と目的

• 長時間にわたる、電気的に大きな領域におけるシミュレーションにおいて、極めて重要な数値分散誤差を低減すること。
• 従来、真空中で成功したとされるm-適応技術を、追加の運動量を有する線形分散性媒体（冷たいプラズマ）へと拡張すること。これは、追加の力学的自由度が存在するため、非自明な問題である。
• ミメティック有限差分スキームの族における最適な自由パラメータを同定し、数値分散を最小化しながら明示的時刻積分を維持すること。
• 指数型時刻差分法（ETD）が、標準的なフロント・フロント時刻積分とは異なり、高次精度の分散低減を可能にすることを示すこと。
• 既存の4次精度スキームよりもスティンシルサイズを小さくすることで、計算効率を向上させる手法の開発

提案手法

• 冷たいプラズマ中のマクスウェル方程式の空間離散化に、最低次元のエッジベースミメティック有限差分（MFD）法を用いる。メッシュは長方形である。
• 一般化された質量アンローディングを適用し、各時刻ステップにおける連立一次方程式の解法を不要にし、最適化に向けた自由パラメータを保持する。
• m-適応手順により、フーリエ領域における分散関係解析を用いて、これらの自由パラメータを最適化し、数値分散を最小化する。
• 時間積分に指数型時刻差分法（ETD）を用いる。これにより、剛性のある線形項の正確な積分と高次精度が達成される。
• 得られた指数型時刻ミメティック有限差分（ETMFD）スキームは、標準的なETMFDスキームの2次精度誤差とは異なり、4次精度の数値分散誤差を達成する。
• この手法は、EおよびJ場のH場からの明示的スタッキングを維持し、半陰的時刻平均を回避する。

実験結果

リサーチクエスチョン

• RQ1 真空中から冷たいプラズマモデルへのm-適応技術の成功した拡張は可能か？（特に、極化電流を記述する補助微分方程式が含まれるため。）
• RQ2 なぜフロント・フロント時刻積分は、冷たいプラズマのMFD族において2次精度を超える分散低減を達成できないのか？一方でETDは成功するのか？
• RQ3 一般化された質量アンローディングMFDスキームにおける最適な自由パラメータの集合は何か？ これは冷たいプラズマにおける数値分散を最小化する。
• RQ4 ETMFD法は4次精度の数値分散誤差を達成するのか？ また、この事実は、シミュレーションによる独立的検証が可能か？
• RQ5 最適化されたETMFD法のスティンシルサイズは、文献に登場する他の4次精度法と比較してどうなるか？

主な発見

• 指数型時刻差分法（ETD）を用いるm-適応手順により、メッシュサイズの2次から4次精度への数値分散誤差低減が成功した。
• 理論的解析と数値シミュレーションの両方で、ETMFD法が4次精度の分散誤差収束を示した。
• ETMFDスキームにおける最適な自由パラメータおよび局所的質量行列は、真空中の場合と同一であり、モデル間で構造的一致性が示された。
• ETMFD法はコンパクトなスティンシルサイズを維持しており、高次精度の基本離散化に依存する他の4次精度法と比較して、計算上の利点がある。
• 標準的なフロント・フロント時刻積分では、冷たいプラズマモデルにおいてm-適応による高次精度分散低減が不可能であることが判明し、時刻積分法選択の重要性が浮き彫りになった。
• 特殊解の数値シミュレーションにより、電場および電流密度の両方のL²誤差において4次収束が確認され、理論的結果の妥当性が裏付けられた。