[论文解读] 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

一句话总结

本文提出了一种全局守恒的有限元代码 DELiTE，用于轴对称磁流体动力学（MHD）模拟，采用线性三角形单元和基于矩阵的离散微分算子，确保质量、能量、环向通量和角动量的全局守恒。该代码成功模拟了紧凑等离子体的形成、磁悬浮和压缩过程，与实验诊断结果高度一致，特别是在通量守恒和压缩过程中环向磁场的演化方面表现优异。

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模拟中的数值不稳定性和非物理解问题。
• 在通用聚变公司的一项新实验中，对紧凑等离子体（CTs）的形成、悬浮和磁压缩过程进行建模与仿真。
• 通过引入绝缘壁和真空场耦合，实现对等离子体-壁相互作用和线圈感应场效应的精确模拟。
• 将数值模型与实验诊断进行验证，特别是偏滤器和环向磁场测量结果。

提出的方法

• DELiTE框架在非结构化三角形网格上使用线性有限元，以矩阵形式推导离散微分算子。
• 构造的离散算子满足拟度量性质，包括离散乘积法则和散度定理，确保在离散层面实现守恒定律。
• 通过用矩阵算子替换连续导数，对单流体双温MHD方程进行离散化，从而在全局范围内保持守恒性。
• 在绝缘区域中单独求解极向真空场，并与离子体区域耦合，以维持环向通量守恒。
• 采用显式时间推进，添加人工扩散以保证数值稳定性，但受限于时间步长和扩散系数的取值范围。
• 将等离子体MHD求解与外部线圈电流及真空场求解耦合，以模拟磁悬浮和压缩过程。

实验结果

研究问题

• RQ1 能否设计一种有限元MHD代码，使其在轴对称几何下全局守恒质量、能量、环向通量和角动量？
• RQ2 DELiTE代码在紧凑等离子体形成、悬浮和磁压缩过程中，对实验诊断的再现能力如何？
• RQ3 引入绝缘壁和真空场耦合对等离子体-壁相互作用和通量守恒有何影响？
• RQ4 为何模拟显示在压缩过程中极向磁通量守恒优于实验结果，其差异原因是什么？
• RQ5 该方案中的人工扩散对密度和磁场诊断的准确性有何影响？

主要发现

• DELiTE代码通过拟度量离散算子，成功实现了总能量、粒子数、环向通量和角动量的全局守恒。
• 模拟结果再现了实验观察：改进的线圈构型可减少等离子体-壁相互作用，表现为压缩过程中壁接触减少。
• 内探针处模拟的极向磁场在压缩过程中持续上升，与通量守恒一致，而实验中约145 μs时出现崩溃。
• 模拟的环向磁场（$B_\phi$）与实验平均$B_\phi$测量结果在定性上一致，特别是在电流转向低电感路径时。
• 该模型捕捉到了$B_\phi$随短路轴电流增加而上升的一般趋势，但轴对称建模无法再现环向不对称不稳定性。
• 代码的一阶精度和显式时间推进限制了模拟速度和稳定性，尤其在高扩散系数下，提示未来工作需采用隐式时间积分方法。