[論文レビュー] New applications for the Boris Spectral Deferred Correction algorithm for plasma simulations

Kris Smedt, Daniel Ruprecht+3 more|arXiv (Cornell University)|Oct 15, 2021

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

ひとこと要約

本稿では、粒子-場（PIC）シミュレーションにおける高次粒子プッシュャとしてBorisスペクトル遅延補正（Boris-SDC）アルゴリズムの使用を提案し、標準的なBoris法と比較して精度が向上し、数値的ドリフトが低減することを示している。相対論的拡張版のBoris-SDCは、0.99cを超える速度領域でも高次収束を維持し、力がゼロの条件下でドリフトを示さない特性を示しており、高精度なプラズマシミュレーションにおいて優位性を発揮するが、現在の計算コストのトレードオフを伴う。

ABSTRACT

The paper investigates two new use cases for the Boris Spectral Deferred Corrections (Boris-SDC) time integrator for plasma simulations. First, we show that using Boris-SDC as a particle pusher in an electrostatic particle-in-cell (PIC) code can, at least in the linear regime, improve simulation accuracy compared with the standard second order Boris method. In some instances, the higher order of Boris-SDC even allows a much larger time step, leading to modest computational gains. Second, we propose a modification of Boris-SDC for the relativistic regime. Based on an implementation of Boris-SDC in the extsc{runko} PIC code, we demonstrate for a relativistic Penning trap that Boris-SDC retains its high order of convergence for velocities ranging from $0.5c$ to $>0.99c$. We also show that for the force-free case where acceleration from electric and magnetic field cancel, Boris-SDC produces less numerical drift than Boris.

研究の動機と目的

• 電磁場PICシミュレーションにおけるBoris-SDCを高次粒子プッシャとして評価すること。標準的なBoris法は2次精度である。
• PICコードにおける2次時間積分の制限を是正すること。これは、高次空間離散化が導入されていても、全体の精度を制限する要因である。
• Boris-SDCを相対論的領域に拡張し、相対論的ローレンツ方程式の高次積分を可能にすること。
• 力がゼロの状況（電場と磁場力が相殺する）における数値的ドリフトを評価すること。これは長期的シミュレーションの忠実性を検証するための重要なテストである。
• Boris-SDCの計算コストが、実際のプラズマ構造における精度向上に見合うかどうかを検討すること。

提案手法

• 元来単粒子動力学を目的として開発されたBoris-SDCアルゴリズムを、RUNKO PICコードフレームワーク内での粒子プッシャとして適応する。
• スペクトル遅延補正（SDC）を用いて、低次時間積分を繰り返し補正することで、配置ノード上で複数回のスイープを実行し、高次精度を達成する。
• 相対論的修正をBoris-SDCに実装するため、電磁場下での相対論的運動量とローレンツ力に対応するように速度更新ステップを再定式化する。
• 2つのテストケースに適用する：電荷-電場PICシミュレーションにおける2ストリーム不安定性とランダム減衰。精度と作業-精度効率を評価する。
• 0.5cから>0.99cの速度を有するペンリングトラップ設定で、相対論的Boris-SDCの収束次数と安定性を検証する。
• 力がゼロの状況（合計加速度がゼロ）における数値的ドリフトを評価する。ノード数と反復回数を変化させ、Boris、Vay、Boris-SDCを比較する。

実験結果

リサーチクエスチョン

• RQ1 Boris-SDCは、完全な静電場PICコードにおいて、標準的なBoris法よりも精度の高い高次粒子プッシャとして効果的であるか？
• RQ2 PICシミュレーションにおけるBoris-SDCの使用は、特に高次空間離散化と組み合わせた場合に、測定可能な計算効率の向上をもたらすか？
• RQ3 Boris-SDCは、高次収束性と安定性を保ったまま、相対論的領域に適切に一般化可能か？
• RQ4 力がゼロの条件下で、Boris-SDCの数値的ドリフト性能はいかがであり、長期的なエネルギー・運動量保存が重要な状況でどう機能するか？
• RQ5 特定のBoris-SDC設定がゼロまたはほぼゼロの数値的ドリフトを示す理由は何か？この現象は理論的に説明可能か？

主な発見

• Boris-SDCは非相対論的および相対論的領域の両方で高次収束を達成しており、観測された収束次数は四則演算の次数に一致し、理論的期待値に近い。
• 非相対論的状況では、Boris-SDCは同じ時間刻みで標準的なBoris法よりも位置誤差と速度誤差を低減するが、空間離散化誤差が支配的であるため、計算コストが効率の向上を相殺する。
• 相対論的ペンリングトラップシミュレーションでは、0.5cから>0.99cの速度範囲でBoris-SDCが高次収束を維持し、極端な条件下でも安定性を示している。
• 力がゼロの状況では、M=3ノードかつK=1反復のBoris-SDCは数値的ドリフトを示さず、M=5ノードかつK=1反復の設定では最小限のドリフト（BorisおよびVayの両方を下回る）を示した。
• 一部のBoris-SDC設定では完全にドリフトがゼロであったが、著者らはこの現象の理論的説明を欠いている。
• 誤差が10−3未満の範囲では、相対論的シミュレーションにおいて計算上の利点が観察された。これは、高精度領域における効率的利点の可能性を示唆している。