[論文レビュー] Gyrokinetic Continuum Simulation of Turbulence in Open-Field-Line Plasmas

E. L. Shi|arXiv (Cornell University)|Aug 24, 2017

Magnetic confinement fusion research参考文献 184被引用数 12

ひとこと要約

本稿では、トカマクのスクラップ・オフ層におけるオープン磁場線での乱流をモデル化できる、エネルギー保存型不連続ガラーキン法と新規の導電性シース境界条件を用いた、最初のギャロアクティブラジエーション連続体シミュレーションコードを提示する。主な貢献は、安定でモーメント保存型のアルゴリズムであり、これにより境界領域におけるプラズマ乱流（電磁効果およびシース物理を含む）の正確なシミュレーションが可能になった。直線的およびヘリカルな磁場線幾何構造の両方での検証が行われた。

ABSTRACT

The properties of the boundary plasma in a tokamak are now recognized to play a key role in determining the achievable fusion power and the lifetimes of plasma-facing components. Accurate quantitative modeling and improved qualitative understanding of the boundary plasma ultimately require five-dimensional gyrokinetic turbulence simulations, which have been successful in predicting turbulence and transport in the core. The additional challenges of boundary-plasma simulation necessitate the development of new gyrokinetic codes or major modifications to existing core gyrokinetic codes. In this thesis, we develop the first gyrokinetic continuum code capable of simulating plasma turbulence on open magnetic field lines, which is a key feature of a tokamak scrape-off layer. In contrast to prior attempts at this problem, we use an energy-conserving discontinuous Galerkin discretization in space. To model the interaction between the plasma and the wall, we design conducting-sheath boundary conditions that permit local currents into and out of the wall. We start by designing spatially one-dimensional kinetic models of parallel SOL dynamics and solve these systems using novel continuum algorithms. By generalizing these algorithms to higher dimensions and adding a model for collisions, we present results from the first gyrokinetic continuum simulations of turbulence on two types of open-field-line systems. The first simulation features uniform and straight field lines, such as found in some linear plasma devices. The second simulation is of a hypothetical model we developed of the NSTX scrape-off layer featuring helical field lines. These developments comprise a major step towards a gyrokinetic continuum code for quantitative predictions of turbulence and transport in the boundary plasma of magnetic fusion devices.

研究の動機と目的

• オープン磁場線プラズマ、特にトカマクのスクラップ・オフ層における乱流をシミュレートできるギャロアクティブラジエーション連続体コードの開発を目的とする。
• 境界プラズマシミュレーションにおける大振幅揺動、電磁効果、複雑な磁場幾何構造の課題に対処することを目的とする。
• 正確な長時間シミュレーションを実現するため、エネルギー保存型およびモーメント保存型の数値スキームの実装を目的とする。
• 壁への局所的電流交換を許容する導電性シース境界条件を用いて、プラズマ-壁相互作用をモデル化することを目的とする。
• 線形装置およびNSTXに類似した幾何構造を代表する直線的およびヘリカルな磁場線配置の両方におけるコードの検証を目的とする。

提案手法

• 空間離散化に不連続ガラーキン有限要素法を採用し、位相空間におけるエネルギー保存および高次精度を確保した。
• 壁への局所的電流注入および抽出を許容する導電性シース境界条件を開発し、プラズマ-壁相互作用をモデル化した。
• 標準の多項式基底と比較して、分解能効率を向上させるために、速度空間に指数関数的重み付き多項式を導入した。
• 分布関数のすべてのモーメントを保存するように修正した不連続ガラーキンアルゴリズムを策定し、標準実装における欠陥を是正した。
• 一次元の並進方向SOL動的モデルを高次元に一般化し、現実的な乱流シミュレーションに適合する衝突モデルを統合した。
• 線形化された方程式を用いたギャロアクティブラジエーションフレームワークを採用し、イオンおよび電子の動的挙動を記述した。また、修正された並進電流モデルにより、シース電流効果を組み込んだ。

実験結果

リサーチクエスチョン

• RQ1 ギャロアクティブラジエーション連続体コードは、トカマクのスクラップ・オフ層におけるオープン磁場線での乱流を正確にシミュレートできるか？
• RQ2 導電性シース境界条件は、オープン磁場線系における乱流の安定性およびダイナミクスにどのように影響するか？
• RQ3 指数関数的重み付き多項式は、標準多項式と比較して、速度空間離散化の効率をどの程度向上させるか？
• RQ4 修正された不連続ガラーキンアルゴリズムは分布関数のモーメントを保存するか？また、なぜこれが長時間シミュレーションにおいて重要なのか？
• RQ5 シース電流の影響は、オープン磁場線構造における乱流モードの増幅率にどのような影響を及えるか？

主な発見

• エネルギー保存特性を持つ不連続ガラーキン法を用いた、オープン磁場線におけるギャロアクティブラジエーション連続体シミュレーションが、初めて成功裏に実現された。
• 導電性シース境界条件により、スクラップ・オフ層物理に不可欠なプラズマ-壁間の電流交換が現実的にモデル化された。
• 修正された不連続ガラーキンアルゴリズムは、標準実装とは異なり、分布関数のすべてのモーメントを正確に保存した。
• 速度空間における指数関数的重み付き多項式は、特に急峻な勾配を示す領域において、分布関数の分解能効率を顕著に向上させた。
• 強いシース電流極限（εₜ ≫ 1）において、乱流モードの増幅率は1/√εₜの要因で低下し、顕著な安定化効果が確認された。
• 直線的およびヘリカルな磁場線幾何構造におけるシミュレーション結果から、コードの頑健性および現実的なトカマク境界構造への適用可能性が確認された。