[論文レビュー] A model for the fast evaluation of prompt losses of energetic ions in stellarators

J. L. Velasco, I. Calvo+5 more|arXiv (Cornell University)|Jun 10, 2021

Magnetic confinement fusion research参考文献 36被引用数 22

ひとこと要約

本稿では、ストラルレーターにおけるエネルギーの高いイオンの軌道が閉じ込められているか即時損失であるかを分類する、高速で径方向に局所的なモデルを提示する。バウンス平均化された拡散力学方程式に基づき、臨界拡散パラメータ Γα を計算することで、完全なモンテカルロシミュレーションと比較して高い精度で構造依存の即時損失を予測可能となり、ストラルレーター最適化ワークフローへの効率的な統合を可能にする。

ABSTRACT

A good understanding of the confinement of energetic ions in non-axisymmetric magnetic fields is key for the design of reactors based on the stellarator concept. In this work, we develop a model that, based on the radially-local bounce-averaged drift-kinetic equation, classifies orbits and succeeds in predicting configuration-dependent aspects of the prompt losses of energetic ions in stellarators. Such a model could in turn be employed in the optimization stage of the design of new devices.

研究の動機と目的

• 非軸対称なストラルレーター磁場配置におけるエネルギーの高いイオンの即時損失を予測する計算的に効率的なモデルの開発。
• バウンス平均化された径方向拡散に基づき、粒子軌道を閉じ込め・即時損失に分類することで、磁場配置の性能を迅速に評価可能にする。
• W7-Xのさまざまな平衡状態（β、ミラー項、回転変換が変化する）において、高精度な全軌道モンテカルロシミュレーション（ASCOT）と照合してモデルを検証する。
• リアルタイムの構成スクリーニングに適した速度と正確性を確保することで、モデルをストラルレーター設計最適化パイプラインに統合可能にする。
• W7-Xを越えた一般性の検証を含め、ヘリオトロンおよび準対称ストラルレーターへの応用を検討し、グローバルシミュレーションの設定を支援する可能性を評価する。

提案手法

• モデルはバウンス平均化された拡散力学方程式を用い、磁場拡散の径方向成分を計算し、バウンス平均化された径方向拡散速度に注目する。
• モデルはパラメータ Γα = (2/π) arctan(|∂αJ| / |∂sJ|) を導入し、接線方向と径方向の拡散成分の比を定量化し、軌道の閉じ込めを決定する。
• Γα > 0.5 である場合、バウンス平均化された径方向磁場拡散が非ゼロであるため、軌道を即時損失と分類する。
• 本手法はKNOSOSコードに実装されており、フラックス面に接する磁場拡散を厳密に取り扱い、STELLOPT最適化スイートに統合されている。
• モデルの予測は、β（0–7%）と放物線的圧力プロファイルが変化するKJMストラルレーター構成の71個の自由境界VMEC平衡状態において、ASCOTモンテカルロシミュレーションと照合された。
• 拡散計算に vE を用いた径方向電場の導入により、将来の E×B 効果および速度空間における衝突拡散の組み込みが可能になる。

実験結果

リサーチクエスチョン

• RQ1 全軌道追跡を伴わずに、径方向に局所的でバウンス平均化されたモデルが、ストラルレーターにおけるエネルギーの高いイオンの即時損失を正確に予測できるか？
• RQ2 W7-Xに類似したストラルレーターにおいて、臨界パラメータ Γα が異なる磁場配置および β 値における即時損失率とどの程度相関するか？
• RQ3 β、径方向位置（s）、磁場配置パラメータ（ミラー項、回転変換、β）の変化に伴う即時損失の依存性をどの程度モデルが捉えられるか？
• RQ4 βに伴う即時損失の依存性に非単調的またはしきい値的挙動（例：2–3% の周辺に最小損失領域）が現れるかどうかをモデルが特定できるか？
• RQ5 径方向電場および衝突効果の組み込みが、現実的な構成における即時損失予測にモデルの予測をどのように変更するか？

主な発見

• モデルはW7-Xのさまざまな平衡状態において、エネルギーの高いイオンが s = 0.06、0.25、0.50 で生成された場合に、ASCOTモンテカルロシミュレーションと強く一致する高い精度で即時損失率を予測した。
• 即時損失はβの増加に伴い概ね線形に減少し、内側のフラックス面（s = 0.06）では2–3% β付近に最小勾配が観察された。これは最適化の最適領域である可能性を示唆する。
• モデルはs = 0.06における低βが即時損失に悪影響を及ぼすことを正しく捉えており、これはダイアモグネティック効果が閉じ込めを悪化させる可能性を示しており、s = 0.25 では観察されない特徴である。
• ミラー項や回転変換の変化を含むさまざまな磁場配置においても、モデルの予測は安定しており、構成依存性と予測能力の両方を示している。
• モデルは、バウンス平均化された径方向拡散が大きなスーパーバナナ軌道が即時損失の主因であることを明らかにした。また、Γα > 0.5 によってその発生が良好に予測された。
• モデルの成功は、バウンス平均化された拡散力学方程式がエネルギーの高いイオンのニューコンダクティブ輸送を定量的に記述できることを示唆し、初期の閉じ込め軌道分類に注目することで、より高速なグローバルシミュレーションが可能になる。