[論文レビュー] A globally conservative finite element MHD code and its application to the study of compact torus formation, levitation and magnetic compression

Carl Dunlea, Ivan Khalzov|arXiv (Cornell University)|Jul 31, 2019

Magnetic confinement fusion research参考文献 16被引用数 4

ひとこと要約

本稿では、線形三角形要素と行列形式の離散微分作用素を用いて、質量、エネルギー、トロイダル磁束、角運動量の全領域的保存を保証する、軸対称磁気流体力学（MHD）のグローバルに保存的な有限要素コード、DELiTEを提示する。このコードは、コンパクトトロイダスの形成、磁気浮上、圧縮を効果的にモデル化できており、特に磁束保存と圧縮中のトロイダル磁場の時間発展に関して、実験的診断と強い一致を示している。

ABSTRACT

The DELiTE (Differential Equations on Linear Triangular Elements) framework was developed for spatial discretisation of partial differential equations on an unstructured triangular grid in axisymmetric geometry. The framework is based on discrete differential operators in matrix form, which are derived using linear finite elements and mimic some of the properties of their continuous counterparts. A single-fluid two-temperature MHD code is implemented in this framework. The inherent properties of the operators are used in the code to ensure global conservation of energy, particle count, toroidal flux, and angular momentum. The code was applied to study a novel experiment in which a compact torus (CT), produced with a magnetized Marshall gun, is magnetically levitated off an insulating wall and then magnetically compressed through the action of currents in the levitation/compression coils located outside the wall. We present numerical models for CT formation, levitation, and magnetic compression, and comparisons between simulated and experimental diagnostics.

研究の動機と目的

• 軸対称系において、質量、エネルギー、トロイダル磁束、角運動量の内在的全領域的保存を実現する有限要素MHDコードの開発。
• 離散的作用素レベルに保存則を埋め込むことで、複雑なMHDシミュレーションにおける数値不安定性や物理的に不自然な解の発生を是正。
• General Fusionでの新規実験を想定し、コンパクトトロイダス（CT）の形成、浮上、磁気的圧縮をモデル化・シミュレート。
• 絶縁壁と真空磁場結合を組み込むことで、プラズマ-壁相互作用およびコイル由来磁場効果の正確なシミュレーションを可能に。
• 主にポロイダルおよびトロイダル磁場測定に基づく実験的診断と比較して、数値モデルの妥当性を検証。

提案手法

• DELiTEフレームワークは、非構造的三角形メッシュ上に線形有限要素を用い、行列形式の離散微分作用素を導出。
• 離散的作用素は、離散積の法則および発散定理を満たすように構築され、離散レベルでの保存則の保証を実現。
• 単一流体二温度MHD方程式は、連続的微分を行列作用素に置き換えることで離散化され、保存則がグローバルに保持される。
• 絶縁領域におけるポロイダル真空磁場は別個に解き、プラズマ領域に結合することで、トロイダル磁束保存を維持。
• 明示的時間積分が用いられ、数値的安定性を確保するための人工的拡散が追加され、時間ステップサイズおよび拡散係数の範囲が制限される。
• プラズマMHD解と外部コイル電流および真空磁場解を結合することで、磁気浮上および圧縮をモデル化。

実験結果

リサーチクエスチョン

• RQ1 軸対称幾何における有限要素MHDコードとして、質量、エネルギー、トロイダル磁束、角運動量の全領域的保存を実現できるか？
• RQ2 コンパクトトロイダスの形成、浮上、磁気的圧縮の過程で、DELiTEコードは実験的診断をどれほど正確に再現できるか？
• RQ3 絶縁壁と真空磁場結合を組み込むことで、プラズマ-壁相互作用および磁束保存にどのような影響が生じるか？
• RQ4 圧縮過程において、シミュレーションではポロイダル磁束保存が実験よりも良好に得られるが、その原因は何か？
• RQ5 スキームに組み込まれた人工的拡散は、密度および磁場診断の精度にどのように影響を与えるか？

主な発見

• DELiTEコードは、ミメティックな離散的作用素を用いることで、全エネルギー、粒子数、トロイダル磁束、角運動量のグローバルな保存を成功裏に実現。
• シミュレーションは、コイル構成の改善がプラズマ-壁相互作用を低減することを実験的観察と一致させ、圧縮過程における壁接触の減少を示した。
• 内側のプローブにおけるシミュレートされたポロイダル磁場は、圧縮過程でも上昇を続けるが、これは磁束保存と整合的であり、実験では約145 μsで観察された崩壊とは対照的である。
• シミュレートされたトロイダル磁場（$B_\phi$）は、実験的平均$B_\phi$測定値と定性的に一致しており、特に自己誘導インダクタンスの低いパスに電流が分散する際の挙動で顕著である。
• モデルは、クローバードシャフト電流の増加に伴う$B_\phi$上昇の一般的なトレンドを捉えているが、軸対称モデルではトロイダル非対称不安定性を再現できない。
• コードの1次精度および明示的時間積分により、シミュレーションの速度と安定性が制限され、特に高拡散係数では顕著である。今後は、暗黙的時間積分の導入が望ましいと考えられる。