[論文レビュー] Application of the averaging method to the gyrokinetic plasma

Emmanuel Frénod|HAL (Le Centre pour la Communication Scientifique Directe)|Feb 26, 2007

Advanced Numerical Methods in Computational Mathematics参考文献 37被引用数 14

ひとこと要約

この論文は、強い磁場による高速振動を示す特異摂動ODE系に従うギャロアクイネティックプラズマ系に平均化法を適用する。粒子軌道の厳密な漸近展開をεの級数として確立し、各項が非振動的減少方程式の解に作用する周期的作用素であることを示し、核融合および分離応用における荷電粒子動力学の正確なモデル化を可能にする。

ABSTRACT

we show that the solution to an oscillatory-singularly perturbed ordinary differential equation may be asymptotically expanded into a sum of oscillating terms. Each of those terms writes as an oscillating operator acting on the solution to a non oscillating ordinary differential equation with an oscillating correction added to it. The expression of the non oscillating ordinary differential equations are defined by a recurrence relation. We then apply this result to problems where charged particles are submitted to large magnetic field.

研究の動機と目的

• 強い磁場に起因する高速振動を示す特異摂動ODE系の解の漸近的二尺度展開を正当化すること。
• プラズマ物理学におけるガイド中心および有限ラーミア半径領域の数学的基盤を提供すること。
• εの累乗の形で粒子軌道を体系的に展開すること。ここでεはサイクロトロン周波数の逆数を表す。
• 磁気的閉じ込め核融合および同位体分離分野におけるプラズマをシミュレートするための効率的な粒子格子法（PIC）の開発を可能にすること。
• 展開における遅いダイナミクスを支配する減少方程式の再帰関係を確立すること。

提案手法

• 時間および位相における周期的高速振動を示す特異摂動ODE系を扱うために平均化法を用いる。
• 解Xε(t)の形式的漸近展開をεの級数として構築し、各項Xⁱ(t, (t−s)/ε; x, s) が高速位相変数に関して2π周期的であることを保証する。
• 一次項X⁰は、高速ダイナミクスを平均化して得られる非振動的ODEを満たし、修正項は再帰的に生成される。
• 解は遅い（マクロな）成分と速い（サイクロトロン）成分に分解され、速い部分はベクトル場b(t, Z)に従う。
• 平均化フレームワークの有効性を保証するため、高速部分系∂Z/∂θ = b(t, Z)がすべてのtおよびzに対して2π周期的解Z(t, θ; z)を有することに依存する。
• 展開式を元のODEに代入し、εの累乗を一致させることで、高次項の再帰関係を導出する。

実験結果

リサーチクエスチョン

• RQ1 ε → 0 のとき、高速振動を示す特異摂動ODE系の解をεの累乗で体系的に展開する方法は何か？
• RQ2 強い磁場下における粒子軌道の漸近展開における一次項および高次項の構造はどのようなものか？
• RQ3 振動的かつ特異摂動を有するギャロアクイネティック系に平均化法を厳密に適用する方法は何か？
• RQ4 ガイド中心および有限ラーミア半径領域における遅いダイナミクスを支配する減少方程式は何か？
• RQ5 粒子運動の周期的補正項は、元々の周期的高速ダイナミクスからどのように生じるのか？

主な発見

• 解Xε(t; x, s) はεの累乗における厳密な漸近展開を有し、各項Xⁱ(t, (t−s)/ε; x, s) は高速位相変数に関して2π周期的である。
• 一次項X⁰は、高速ダイナミクスを平均化して得られる非振動的ODEを満たし、初期条件X⁰(s, 0; x, s) = x を満たす。
• i ≥ 1 についての高次項Xⁱは、平均化されたベクトル場の微分を含む再帰関係によって再帰的に生成される。
• 高速部分系∂Z/∂θ = b(t, Z) がすべてのtおよびzに対して2π周期的解を有することを仮定すれば、展開は有効である。
• この方法は、ギャロアクイネティックシミュレーションにおける平均化モデルの使用を正当化する。特に、磁気的閉じ込め核融合および同位体分離分野において有効である。
• このフレームワークは、高速のサイクロトロン運動と遅いガイド中心運動を分離することで、プラズマ物理学における正確かつ効率的なPIC法の構築に数学的基盤を提供する。