[論文レビュー] A spectral, quasi-cylindrical and dispersion-free Particle-In-Cell algorithm

Remi Lehé, Manuel Kirchen+3 more|Jul 16, 2015

Spectroscopy and Laser Applications被引用数 3

ひとこと要約

本稿では、円対称系における分散のない電磁波伝播を可能にする、ハンケル変換とフーリエ変換を組み合わせた新規なスペクトル的・準円筒座標系粒子法（PIC）アルゴリズムを提示する。径方向および軸方向におけるスペクトル法の活用により、偽の数値分散（特にゼロ次数値チチェレノフ効果）が排除され、レーザー・プラズマおよびビーム・プラズマシミュレーションにおいて計算効率を維持する。

ABSTRACT

We propose a spectral Particle-In-Cell (PIC) algorithm that is based on the combination of a Hankel transform and a Fourier transform. For physical problems that have close-to-cylindrical symmetry, this algorithm can be much faster than full 3D PIC algorithms. In addition, unlike standard finite-difference PIC codes, the proposed algorithm is free of numerical dispersion. This algorithm is benchmarked in several situations that are of interest for laser-plasma interactions. These benchmarks show that it avoids a number of numerical artifacts, that would otherwise affect the physics in a standard PIC algorithm - including the zero-order numerical Cherenkov effect.

研究の動機と目的

• 近似的に円対称な系に対して、3次元直交座標系PICシミュレーションの高い計算コストを低減すること。
• 標準の有限差分PICコードにおいて主要な誤差源である真空中の偽の数値分散を排除すること、特に相対論的粒子およびレーザーパルスに対して。
• 準円筒座標系アルゴリズムの効率性と、スペクトル（PSATD型）手法の分散フリーな精度を統合されたフレームワークで統合すること。
• 数値アーティファクト（例：数値チチェレノフ効果）を伴わず、レーザー駆動加速およびその他の相対論的プラズマ現象の高精度なシミュレーションを可能にすること。

提案手法

• 径方向依存性にハンケル変換、軸方向依存性にフーリエ変換を用いたスペクトルPICアルゴリズムを定式化する。
• 偏微分の有限差分近似を回避するため、スペクトル空間でマクスウェル方程式を解く擬似スペクトル的手法を適用する。
• 離散ハンケル変換（DHT）および逆離散ハンケル変換（IDHT）を用いて、径方向スペクトルモードと実空間グリッド点との間の変換を実施する。
• 特に非ゼロの軸対称モードに対して、直交性を保証し正確な再構成を可能にするために、変換行列に制約を課す。
• DHTおよびIDHT行列を事前に計算・逆行列化する行列形式を採用し、精度および安定性を維持する。
• 相対論的速さでさえも、偽の位相速度誤差が生じないよう、場の式を解く方法を適用する。

実験結果

リサーチクエスチョン

• RQ1 真空中に偽の数値分散が存在しない、準円筒座標系のスペクトルPICアルゴリズムを構築可能か？
• RQ2 スペクトル法（ハンケルおよびフーリエ変換）をPIC粒子プッシュャーとどのように組み合わせ、円対称系における精度を維持できるか？
• RQ3 提案されたアルゴリズムは、標準の有限差分PICコードと比較して、ゼロ次数値チチェレノフ効果をどの程度抑制するか？
• RQ4 非軸対称効果を完全に3次元的に再現しつつ、2次元シミュレーションと同等の計算効率を達成できるか？
• RQ5 非ゼロの方位モードに対して離散ハンケル変換を実装する際の数値的課題は何か？ それらをどのように緩和できるか？

主な発見

• 解析解とのベンチマークにより、真空中に偽の数値分散が完全に存在しないことが確認された。
• 標準のPICシミュレーションでよく知られたアーティファクトであるゼロ次数値チチェレノフ効果が、本手法で効果的に抑制された。
• レーザー・ウェーブフィールド加速のベンチマークにおいて、標準の有限差分PICコードが数値誤差を生じるのに対し、本アルゴリズムは正しい脱位長を正しく予測した。
• 適切な行列制約を課えた離散ハンケル変換の使用により、径方向モード全域で正確かつ安定した場の再構成が可能となった。
• 近似的に円対称な系に対して、完全な3次元直交座標系PICシミュレーションに比べて顕著な高速化が達成された。
• 本手法は $n=0$ および $n\neq0$ の両方の方位モードに対して頑健であり、変換行列に実験的制約を課すことで数値的安定性と精度が保証された。