[论文解读] Linear MHD stability studies with the STARWALL code

P. Merkel, E. Strumberger|arXiv (Cornell University)|Aug 20, 2015

Magnetic confinement fusion research参考文献 19被引用 20

一句话总结

本论文介绍了STARWALL代码，这是一种用于具有复杂多连通导电结构的3D理想与阻性平衡态下线性MHD稳定性分析的数值工具。它通过与OPTIM耦合，实现了对阻性壁模（RWMs）的鲁棒反馈稳定化，并引入了Lüst-Martensen项以精确建模垂直模，通过ASDEX Upgrade和ITER构型的模拟展示了其在详细本征模和电流响应分析方面的能力。

ABSTRACT

The STARWALL/CAS3D/OPTIM code package is a powerful tool to study the linear MHD stability of 3D, ideal equilibria in the presence of multiply-connected ideal and/or resistive conducting structures, and their feedback stabilization by external currents. Robust feedback stabilization of resistive wall modes can be modelled with the OPTIM code. Resistive MHD studies are possible combining STARWALL with the linear, resistive 2D CASTOR code as well as nonlinear MHD simulations combining STARWALL with the JOREK code. In the present paper, a detailed description of the STARWALL code is given and some of its applications are presented to demonstrate the methods used. Conducting structures are treated in the thin wall approximation and depending on their complexity they are discretized by a spectral method or by triangular finite elements. As an example, a configuration is considered consisting of an ideal plasma surrounded by a vacuum domain containing a resistive wall and bounded by an external wall. Ideal linear MHD modes and resistive wall modes in the presence of multiply-connected walls are studied. In order to treat the vertical mode self-consistently the STARWALL code has been completed by adding the so-called Luest-Martensen term generated by a constant normal displacement of the plasma. The appendix contains the computation of the 2D Fourier transform of singular inductance integrals, and the derivation of an asymptotic expansion for large Fourier harmonics.

研究动机与目标

• 开发一个灵活且高精度的数值框架，用于研究具有复杂导电结构的3D托卡马克构型中的线性MHD稳定性。
• 通过与STARWALL联合使用的OPTIM代码，实现对阻性壁模（RWMs）的鲁棒反馈稳定化。
• 通过引入Lüst-Martensen项，扩展对轴对称不稳定性中垂直模的处理。
• 支持与非线性和阻性MHD代码（如JOREK和CASTOR3D）的耦合，以实现全面的稳定性分析。
• 为下一代聚变装置（如ITER）中RWMs和垂直位移事件（VDEs）的研究提供计算基础。

提出的方法

• STARWALL采用薄壁近似，并通过谱方法或三角形有限元对导电结构进行离散化。
• 它在变分框架下求解理想MHD能量原理，通过边界积分方法引入真空和壁的贡献。
• 该代码计算3D平衡态下的本征模和增长率，通过与线性阻性2D CASTOR代码耦合，纳入阻性壁效应。
• 通过实现Lüst-Martensen项，对具有恒定法向离子位移的轴对称（n=0）模进行真空响应建模。
• 反馈稳定化通过OPTIM代码建模，该代码通过同时考虑不稳定和稳定模来优化线圈电压。
• STARWALL已集成到包含CAS3D（理想MHD）、NEMEC（3D平衡）和JOREK（非线性MHD）的更大代码套件中，支持多尺度模拟。

实验结果

研究问题

• RQ1 如何对具有多连通阻性和理想导电壁的3D平衡态实现线性MHD稳定性的精确计算？
• RQ2 壁几何复杂性对阻性壁模（RWMs）增长率和结构有何影响？
• RQ3 引入Lüst-Martensen项如何改善轴对称构型中垂直模稳定性的建模？
• RQ4 在实际托卡马克壁构型下，反馈控制在多大程度上能稳定RWMs？最优控制如何实现？
• RQ5 STARWALL与非线性和阻性MHD代码（如JOREK、CASTOR3D）的耦合如何提升对破裂和VDEs的预测能力？

主要发现

• STARWALL成功计算了具有阻性和理想壁的3D平衡态下的本征模和增长率，与解析基准结果高度一致。
• 引入Lüst-Martensen项后，能够精确建模轴对称垂直模，ASDEX Upgrade构型下n=0、n=1和n=2的模增长率分别为26 s⁻¹、67 s⁻¹和314 s⁻¹。
• 多模态（n=0,1,2）阻性壁表面电流因壁几何不对称性而表现出明显的耦合与分裂，本征值与单模解唯一关联。
• CAS3D中非等距径向网格在仅使用931个点的情况下，即可在共振面（如q=3/2）附近实现高分辨率，优于需约40,000个点的等距网格。
• 通过OPTIM对RWMs的反馈稳定化建模具有鲁棒性，该方法在优化线圈电压的同时避免了对稳定模的意外激发。
• 与JOREK耦合可实现对VDEs、破裂和QH模的非线性研究，初步结果表明阻性壁和环向电流效应的行为一致。