[論文レビュー] Theoretical analysis of the saturation phase of the $1/1$ energetic-ion-driven resistive interchange mode

J. Varela, D. A. Spong+5 more|arXiv (Cornell University)|Jan 3, 2022

Magnetic confinement fusion research参考文献 40被引用数 15

ひとこと要約

本研究では、FAR3dコードを用いた非線形ギャロア流体シミュレーションを用い、大型ヘリカル装置（LHD）における1/1のエネルギー粒子駆動型抵抗性磁気インターフェースモード（EIC）の飽和期をモデル化した。実験で観測された主要な特徴を再現した。具体的には、エネルギー粒子βの増加に伴うモード振幅および幅の増大、励起閾値、9.4 kHzのバースト発生、高調波成分、二次モード結合に起因する内向き伝播、および圧力プロファイル回復に連動した安定化が再現された。

ABSTRACT

The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. A set of non linear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The non linear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP $\beta$ increases, the EP beta threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the non linear destabilization of the 3/4 and 2/3 EPMs, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.

研究の動機と目的

• 大型ヘリカル装置（LHD）における1/1のエネルギー粒子駆動型抵抗性磁気インターフェースモード（EIC）の非線形的飽和期をシミュレートすること。
• 実験的に観測されたEICの多段階的進化、特にバースト発生およびモード伝播を再現すること。
• ヘリカルトラップエネルギー粒子（EPs）がEICの不安定化および飽和に果たす役割を調査すること。
• モード誘発圧力フラットネスとその後のRIC再活性化との間のフィードバック機構を分析すること。
• シミュレーションモデルが実験周波数、モード構造、およびβ閾値行動をどれだけ正確に再現できるかを検証すること。

提案手法

• 熱プラズマに対しては簡略化MHDモデル、エネルギー粒子（EPs）に対してはギャロア流体モデルを用い、解析的TAE成長率と整合する閉じ込め関係を導入したFAR3dギャロア流体コードを採用。
• 熱プラズマには簡略化MHDモデル、EPsにはギャロア流体モデルを用い、解析的TAE成長率と一致する閉じ込め関係を採用。
• ヘリカルトラップEPsの進行運動を含めるために、平均化されたドリフト速度作用素を修正。
• 線形波-粒子共鳴効果（ランダム減衰/増幅）およびジオデシックアコースティック波結合のための並進運動量応答を組み込み。
• LHDショット116190の固定境界平衡をVMECコードから得て、実験プロファイルに一致するパラメータを設定。
• 径方向には有限差分法、トロイダルおよびポリオダル角度にはフーリエ展開を用い、(Δt)³の精度を達成するための半陽解法・陽解法2段階スキームを採用。

実験結果

リサーチクエスチョン

• RQ1 ヘリカルトラップエネルギー粒子が存在する中で、1/1 EICはどのように飽和期に至るか？
• RQ2 LHD実験で観測されたバースト発生（周波数上昇および複雑な径方向モード構造を伴う）を引き起こすメカニズムは何か？
• RQ3 2/3および3/4のEPMsのような二次モードが、バースト期における1/1 EICの内向き伝播にどのように寄与するか？
• RQ4 圧力プロファイルのフラットネスと回復が、RICの安定化およびEICバースト後のRIC再活性化に果たす役割は何か？
• RQ5 FAR3dコードは、実験的に測定された1/1 EICの周波数、振幅、モード構造をどの程度正確に再現できるか？

主な発見

• シミュレーションは、実験測定と整合する9.4 kHzの1/1 EICバースト発生を再現した。
• バースト期に2/2および3/3の高調波が不安定化し、r/a = [0.6, 0.85]範囲で観測された複雑な径方向モード構造と一致した。
• 1/1 EICの内向き伝播は、3/4および2/3のEPMsとの非線形結合に起因し、バースト期に共鳴重ね合わせが生じた。
• 安定化段階は、EICに起因する圧力勾配のフラットネスが消失した後に圧力勾配回復が生じ、RICの再活性化を可能にする。
• EP βの1/1 EIC励起閾値が正確に再現され、βが上昇するに従いモード振幅および幅が増大した。
• 本モデルは、4段階の全過程の進化を捉えた：RIC発生（第I段階）、EIC増幅（第II段階）、周波数ジャンプを伴うバースト（第III段階）、RIC再活性化を伴う島形成（第IV段階）