| 研究課題/領域番号 |
21K03498
|
|---|
| 研究種目 |
基盤研究(C)
|
| 配分区分 | 基金 |
|---|
| 応募区分 | 一般 |
|---|
| 審査区分 |
小区分14010:プラズマ科学関連
|
|---|
| 研究機関 | 東北大学 |
|---|
研究代表者 |
廣田 真 東北大学, 流体科学研究所, 准教授 (40432900)
|
|---|
| 研究期間 (年度) |
2021-04-01 – 2024-03-31
|
| 研究課題ステータス |
交付 (2022年度)
|
| 配分額 *注記 |
1,300千円 (直接経費: 1,000千円、間接経費: 300千円)
2023年度: 390千円 (直接経費: 300千円、間接経費: 90千円)
2022年度: 520千円 (直接経費: 400千円、間接経費: 120千円)
2021年度: 390千円 (直接経費: 300千円、間接経費: 90千円)
|
|---|
| キーワード | 磁気流体力学 / 自由境界問題 / 二流体プラズマ / 数値シミュレーション |
|---|
| 研究開始時の研究の概要 |
磁気流体力学(MHD)は磁場閉じ込め核融合プラズマや宇宙プラズマの巨視的なスケールの現象を理解するために広く用いられている。しかし、プラズマが真空領域と接しているような自由境界問題、あるいはプラズマ密度がほぼゼロである領域を含む問題は現実に数多く存在するが、数値計算でMHD方程式を解くのが技術的に困難とされる。本研究ではMHD方程式において無視されている微視的スケールの物理(二流体効果)を取り入れた拡張MHD方程式を解くことにより、このような真空中のプラズマの運動を適切にシミュレーションをする技法を構築する。
|
| 研究実績の概要 |
二流体プラズマ方程式（イオンと電子の流体方程式とマクスウェル方程式）に準中性条件を課した拡張MHD方程式は、通常のMHD（磁気流体力学）方程式 に比べるとホール効果と電子慣性効果を含んでおり、MHDでは扱えないプラズマの挙動も物理的に正しく扱えるモデルとして注目されている。特に、宇宙空間や実験室のプラズマは密度が低い真空領域に囲まれていることが多く、MHDではそこでアルヴェン波の伝搬速度が無限大に発散するので理論も数値計算も破綻してしまう。MHDでは電子慣性効果が働く電子の慣性長のスケールを小さいとして無視してしまうが、慣性長はプラズマ密度の二乗根に反比例するため、密度がほとんどゼロになっていくプラズマの周辺部ではむしろ電子慣性効果が支配的となるはずである。
本研究では、拡張MHD方程式を解くことにより、このようなプラズマ領域と真空領域が共存する自由境界問題を物理的に正しく扱うことを目標とし、そのための理論構築と数値解析コードの開発を行っている。具体的に、円柱プラズマの軸方向に電圧をかけて放電プラズマを閉じ込めるZピンチ平衡を解析の対象とし、拡張MHD方程式を解くことで真空領域に囲まれたプラズマの平衡解が得られることを示した。また、開発中の数値シミュレーションコードでも、時間発展の結果、同じ平衡解に収束することが確認できた。プラズマの密度が厳密にゼロになる境界位置では、数学的に特異な解となるため数値的には正則化を行うことで安定な計算が可能になる。まだ単純な問題を扱ってはいるものの、拡張MHD方程式を解くことによって、プラズマ領域と真空領域を自己無撞着に解析できることを初めて示すことができ、より複雑な問題への応用が期待できる。
|
| 現在までの達成度 (区分) |
現在までの達成度 (区分)
2: おおむね順調に進展している
理由
拡張MHDにおいて、自由境界（プラズマと真空の界面）をもつ厳密な平衡解が存在することを示すことができた。ただし、界面において密度分布が滑らかではないことが数学的に予測されたので、プラズマの密度が十分に下がったところでは比熱比を調節することで解を滑らかにすることができることを発見した。これによって、界面の形成や移動の過程を数値的に安定に解くことができるようになった。
|
| 今後の研究の推進方策 |
自由境界をもつ平衡解を明らかにし、数値的にも再現することができたため、外部キンクモードのような不安定モードを加えて三次元シミュレーションを行う。まずは計算コードが安定に動作することを確認しなければならないが、計算コストが飛躍的に増加するため、反復法を用いたアルゴリズムで計算時間がどれくらい必要とされるのかが懸念事項と思われる。
|
