自由境界プラズマに対する非線形拡張MHDシミュレーション技法の構築

点击翻译按钮获取中文摘要

二流体プラズマ方程式（イオンと電子の流体方程式とマクスウェル方程式）に準中性条件を課した拡張MHD方程式は、通常のMHD（磁気流体力学）方程式 に比べるとホール効果と電子慣性効果を含んでおり、MHDでは扱えないプラズマの挙動も物理的に正しく扱えるモデルとして注目されている。特に、宇宙空間や実験室のプラズマは密度が低い真空領域に囲まれていることが多く、MHDではそこでアルヴェン波の伝搬速度が無限大に発散するので理論も数値計算も破綻してしまう。MHDでは電子慣性効果が働く電子の慣性長のスケールを小さいとして無視してしまうが、慣性長はプラズマ密度の二乗根に反比例するため、密度がほとんどゼロになっていくプラズマの周辺部ではむしろ電子慣性効果が支配的となるはずである。本研究では、拡張MHD方程式を解くことにより、このようなプラズマ領域と真空領域が共存する自由境界問題を物理的に正しく扱うことを目標とし、そのための理論構築と数値解析コードの開発を行っている。具体的に、円柱プラズマの軸方向に電圧をかけて放電プラズマを閉じ込めるZピンチ平衡を解析の対象とし、拡張MHD方程式を解くことで真空領域に囲まれたプラズマの平衡解が得られることを示した。また、開発中の数値シミュレーションコードでも、時間発展の結果、同じ平衡解に収束することが確認できた。プラズマの密度が厳密にゼロになる境界位置では、数学的に特異な解となるため数値的には正則化を行うことで安定な計算が可能になる。まだ単純な問題を扱ってはいるものの、拡張MHD方程式を解くことによって、プラズマ領域と真空領域を自己無撞着に解析できることを初めて示すことができ、より複雑な問題への応用が期待できる。