[論文レビュー] Continuum Kinetic Simulations of Plasma Sheaths and Instabilities

Petr Cagaš|arXiv (Cornell University)|Sep 17, 2018

Plasma Diagnostics and Applications参考文献 60被引用数 7

ひとこと要約

本稿では、プラズマシースと不安定性をシミュレートするための、新規の不連続ガレルキンに基づく連続的運動論的手法を提示する。この手法により、物理的電子放出に基づく自己自己整合的モデル化が可能となり、プラズマ-壁相互作用を再現できる。Landau減衰、2ストリーム不安定性およびWeibel不安定性の正確な再現が確認され、Weibel不安定性の非線形的飽和状態において、電場エネルギーが磁場エネルギーと同等の大きさに達することが判明した。一方、壁からの電子放出モデルの導入は、シース構造や電場分布に顕著な影響を及ぼすことが明らかになった。

ABSTRACT

A careful study of plasma-material interactions is essential to understand and improve the operation of devices where plasma contacts a wall. Key contributions of this work include (i) novel continuum kinetic algorithms with novel boundary conditions that directly discretize the Vlasov/Boltzmann equation using the discontinuous Galerkin method, (ii) fundamental studies of a plasma sheath physics with collisions, ionization, and physics-based wall emission, and (iii) theoretical and numerical studies of the linear growth and nonlinear saturation of the kinetic Weibel instability, including its role in plasma sheaths. The continuum kinetic algorithm has been shown to compare well with theoretical predictions of Landau damping of Langmuir waves and the two-stream instability. Benchmarks are also performed using the electromagnetic Weibel instability and excellent agreement is found between theory and simulation. The role of the electric field is significant during nonlinear saturation of the Weibel instability, which was not noted in previous studies. A focus is put on understanding plasma sheath physics which is essential for studying plasma-material interactions. Initial simulations are performed using a collisionless model to match classical sheath theory and the Bohm criterion. A collision operator and volumetric physics-based source terms are introduced and effects of heat flux are briefly discussed. Novel physics-based boundary conditions are developed and included in a general manner with the continuum kinetic algorithm for bounded plasma simulations. These are the first continuum kinetic simulations using self-consistent, wall emission boundary conditions with broad applicability across a variety of regimes.

研究の動機と目的

• 核融合、推進、宇宙システムに不可欠な高精度なプラズマ-壁相互作用の数値的シミュレーションの前進。
• 粒子-セル（PIC）法の制限を回避するため、不連続ガレルキン法を用いてVlasov-Boltzmann方程式を直接解く連続的運動論フレームワークの開発。
• 量子力学的モデルを用いた自己自己整合的で物理的根拠を持つ壁からの電子放出を組み込み、シースシミュレーションの現実性を向上。
• Weibel不安定性の非線形的飽和状態における電場の役割を解明し、従来の研究で無視されてきた点を特定。
• 正確なプラズマ-材料相互作用モデリングのための、シース形成、衝突効果、境界条件に関するベンチマーク品質のシミュレーションを確立。

提案手法

• 位相空間におけるVlasov-Boltzmann方程式の不連続ガレルキン（DG）離散化を用いて、連続的運動論的シミュレーションを実装。
• 高速かつ再現性のあるシミュレーションを実現するため、機械最適化されたモーダルDGカーネルを備えたGkeyll 2.0フレームワークを活用。
• 一般化された反射関数を導入し、鏡像反射や吸収に基づく反射を含む柔軟で自己自己整合的な壁モデリングを可能に。
• 体積内での衝突をモデル化し、シース領域における分布関数の整合性を保つためにBGK衝突演算子を適用。
• 自己自己整合的かつ物理的根拠を持つ電子放出モデル（Furman-PiviおよびBronold-Fehske）を実装し、壁からの電子放出をシミュレート。
• 線形運動論理論とのベンチマークを実施するため、シミュレーションから成長率および減衰率を抽出するためにスイーピングフィッティング手法を用いる。

実験結果

リサーチクエスチョン

• RQ1 連続的運動論的DG法は、Landau減衰や2ストリーム不安定性といった基本的運動論現象をどの程度正確に再現できるか？
• RQ2 Weibel不安定性の非線形的飽和段階において、電場の役割は何か？また、磁場エネルギーと比較してどうなるか？
• RQ3 自己自己整合的で物理的根拠を持つ電子放出モデルは、プラズマ-壁系におけるシース構造や電場プロファイルにどのような影響を及えるか？
• RQ4 衝突および熱フラックス効果は、現実的なプラズマシミュレーションにおけるシースダイナミクスをどの程度変化させるか？
• RQ5 一般化された反射関数は、計算ドメインを縮小しつつも物理的忠実性を保ちながら、正確で対称的なシースシミュレーションを可能にするか？

主な発見

• 連続的運動論的DG法は、2ストリーム不安定性の理論的成長率と0.3%の一致を達成し、粗い速度解像度でも数値的精度が保証された。
• Weibel不安定性の非線形段階において、電場エネルギーは磁場エネルギーと同等の大きさに達することがあり、これは以前に報告のなかった現象である。
• BronoldとFehske（2015）の電子吸収モデルを導入すると、理想の吸収壁と比較して壁近傍の電子密度が2倍に増加した。
• BronoldとFehskeに基づく境界条件を用いることで、シースにおける電場の大きさは約60%低減された。これは顕著な物理的影響を示している。
• 自己自己整合的電子放出モデリングにより、シースポテンシャルおよびイオン加速に測定可能な変化が生じ、現実のプラズマ-壁シミュレーションでは電子放出を無視できないことが示された。
• 本手法は、衝突なしのシミュレーションにおいてもBohm基準およびシース形成を正確に再現でき、イオンが壁から数個のDebye長さの距離でBohm速度に達した。