[論文レビュー] Electromagnetic gyrokinetic instabilities in STEP

Daniel Kennedy, M. Giacomin+5 more|arXiv (Cornell University)|Jul 4, 2023

Magnetic confinement fusion researchPhysics and Astronomy被引用数 3

ひとこと要約

本研究では、STEP球形トカマク設計における電磁的マイクロ不安定性について線形ギャロアクイネティクス解析を行い、イオン温度勾配および捕獲電子によって駆動されるハイブリッドキネティックバルーンイングモード（KBM）がイオンバイノーマルスケールで支配的であることが特定された。衝突項を含むマイクロテアリングモード（MTM）は準主徴的だが有意であり、明確な径方向スケールを有する。結果は複数コードによるベンチマークと解像度収束性の検証により妥当性が確認され、次に続く非線形研究の信頼性が確保された。

ABSTRACT

We present herein the results of a linear gyrokinetic analysis of electromagnetic microinstabilites in the conceptual high-$β$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production, https://step.ukaea.uk). We examine a range of flux surfaces between the deep core and the pedestal top for two candidate flat-top operating points of the prototype device. Local linear gyrokinetic analysis is performed to determine the type of microinstabilities that arise under these reactor-relevant conditions. We find that the equilibria are dominated at ion binormal scales by a hybrid version of the Kinetic Ballooning Mode (KBM) instability that has significant linear drive contributions from the ion temperature gradient and from trapped electrons, while collisional Microtearing Modes (MTMs) are sub-dominantly also unstable at similar binormal scales. The hybrid-KBM and MTM exhibit very different radial scales. We study the sensitivity of these instabilities to physics parameters, and discuss potential mechanisms for mitigating them. The results of this investigation are compared to a small set of similar conceptual reactor designs in the literature. A detailed benchmark of the linear results is performed using three gyrokinetic codes; alongside extensive resolution testing and sensitivity to numerical parameters providing confidence in the results of our calculations, and paving the way for detailed nonlinear studies in a companion article.

研究の動機と目的

• 高β、炉規模、低アスペクト比の球形トカマク（STEP）平衡状態における電磁的マイクロ不安定性の特性を評価すること。
• 2つの候補となるSTEP運転点における、コアからペデスタル領域までのフラックスサーフェスにわたる、支配的および準主徴的不安定性を同定すること。
• 将来の非線形シミュレーションに向け、複数のギャロアクイネティクスコードを用いた数値収束性および解像度要件を評価すること。
• 3つのギャロアクイネティクスコード（GS2、CGYRO、GENE）間での結果のベンチマークを実施し、線形安定性予測の妥当性を検証すること。
• 次に続く非線形乱流シミュレーションの基礎を築くために、解像度およびパrameter感度を確立すること（論文IIの補足研究として）

提案手法

• 2つの候補となるSTEP平衡状態（STEP-EC-HDおよびSTEP-EB-HD）を用いて、GS2コードを用いた局所的線形ギャロアクイネティクス解析を実施する。
• コアからペデスタル領域をカバーする複数のフラックスサーフェス（Ψₙ = 0.49 から 0.95）においてシミュレーションを実施する。
• 線形成長率およびモード周波数のクロスバリデーションとベンチマークのため、3つのギャロアクイネティクスコード（GS2、CGYRO、GENE）を用いる。
• 数値的安定性を確保するため、径方向、並進方向（バルーンイング空間）、および速度空間グリッドにおける解像度収束テストを厳密に実施する。
• 衝突モデル（ローレンツ、スガマ、完全線形化フォッカー・プランク）および粒子種数（2〜5）に対する感度を調査し、特に成長率安定性との関係を検討する。
• 非線形シミュレーションに必要な最小解像度要件を特定するために、詳細なパラメータスキャン（nθ、nperiod、npass、nε）を実施する。

実験結果

リサーチクエスチョン

• RQ1 高β、炉に相当するSTEP平衡状態において、イオンバイノーマルスケールでの支配的電磁的マイクロ不安定性は何か？
• RQ2 フラックスサーフェスの階層にわたる、ハイブリッドキネティックバルーンイングモード（KBM）と衝突項を含むマイクロテアリングモード（MTM）の成長率および径方向スケールは、どのように比較されるか？
• RQ3 非線形シミュレーションの準備として、正確な線形安定性予測に必要な数値的解像度はどの程度か？
• RQ4 線形成長率およびモード周波数は、衝突演算子の選択およびプラズマ種の数にどの程度感度を示すか？
• RQ5 GS2、CGYRO、GENEの複数コードベンチマークにより、不安定性予測の妥当性はどの程度確認されるか？

主な発見

• 平衡状態は、イオンバイノーマルスケール（k_yρ_i ≲ 1）において、イオン温度勾配および捕獲電子の両方による顕著な線形駆動を持つハイブリッドキネティックバルーンイングモード（KBM）によって支配されている。
• 衝突項を含むマイクロテアリングモード（MTM）は準主徴的だが、同様のバイノーマルスケールで不安定であり、ハイブリッド-KBMとは異なり、はるかに広い径方向スケールを有する。
• 低k_yモードではnθ ≥ 32で解像度収束が達成されるが、k_yρ_s > 0.4のモードではnθ ≥ 64を要するため、高波数モードではより高い解像度が必要であることが示された。
• 速度空間解像度は、k_yスペクトル全体にわたり成長率および周波数に弱い依存性を示し、npassおよびnεのテスト範囲で収束が確認された。
• 線形スペクトルは衝突モデルの選択に対して僅かな感度を示すが、CGYROでスガマ演算子を使用した場合、ローレンツまたは完全線形化フォッカー・プランクモデルと比較して成長率が約20％低くなる傾向がある。
• GS2、CGYRO、GENE間の複数コードベンチマークにより、成長率およびモード周波数に強く一致が確認され、非線形シミュレーションへの応用に向けた結果の妥当性が検証された。