磁化プラズマの巨視的運動を記述する磁気流体力学的(MHD)モデルについて,これまで計算不可能だった定常状態の計算を可能にする新しい計算法(疑似アニーリング)を開発した.本研究課題では,低ベータ簡約化MHDに疑似アニーリングを適用し,円柱プラズマの内部に磁気島があったりヘリカル変形した定常状態が,系のエネルギー極小状態として求められることをシミュレーションで実証した.さらに,高ベータ簡約化MHDに拡張し,トロイダルプラズマの定常状態計算にも成功した.また本研究を通じ,ハミルトン系の先進的な数値計算法が必要となり,荷電粒子の運動も含め支配方程式のハミルトン構造を保持した数値計算法の開発も行った.
|