[论文解读] Nonlinear excitation of low-n harmonics in reduced MHD simulations of edge-localized modes

I. Krebs, M. Hoelzl+2 more|arXiv (Cornell University)|May 16, 2013

Magnetic confinement fusion research参考文献 2被引用 3

一句话总结

本研究利用JOREK代码中的简化MHD模拟，研究了在边缘局域模（ELMs）期间低n模态的非线性激发。研究提出了一种二次耦合模型，表明能量通过二阶非线性相互作用从主导模态传递至次主导的低n模态，特别是n = 1模，从而解释了其在ELMs早期非线性阶段于离子等离子体边界处显著增长及空间局域化的原因。

ABSTRACT

Nonlinear simulations of the early ELM phase based on a typical type-I ELMy ASDEX Upgrade discharge have been carried out using the reduced MHD code JOREK. The analysis is focused on the evolution of the toroidal Fourier spectrum. It is found that during the nonlinear evolution, linearly subdominant low-n Fourier components, in particular the n = 1, grow to energies comparable with linearly dominant harmonics. A simple model is developed, based on the idea that energy is transferred among the toroidal harmonics via second order nonlinear interaction. The simple model reproduces and explains very well the early nonlinear evolution of the toroidal spectrum in the JOREK simulations. Furthermore, it is shown for the n = 1 harmonic, that its spatial structure changes significantly during the transition from linear to nonlinearly driven growth. The rigidly growing structure of the linearly barely unstable n = 1 reaches far into the plasma core. In contrast, the nonlinearly driven n = 1 has a rigidly growing structure localized at the plasma edge, where the dominant toroidal harmonics driving the n = 1 are maximal and in phase. The presented quadratic coupling model might explain the recent experimental observation of strong low-n components in magnetic measurements [Wenninger et al., Non-linear magnetic perturbations during edge localized modes in TCV dominated by low n mode components, submitted to Nuclear Fusion].

研究动机与目标

• 理解托卡马克等离子体中类型-I ELMs早期阶段低n环向模态的非线性激发机制。
• 解释实验观测中磁测量结果中显著的低n成分（特别是n = 1），这些成分无法由线性理论预测。
• 研究从线性上勉强不稳定到非线性驱动增长的n = 1模态的转变过程。
• 构建一个能够捕捉ELMs期间环向模态之间能量转移基本物理机制的最小化模型。

提出的方法

• 使用JOREK代码对典型类型-I ELMy ASDEX Upgrade放电进行非线性简化MHD模拟。
• 分析环向傅里叶谱的时间演化，以追踪不同模态成分之间的能量分布。
• 基于环向模态之间的二阶非线性相互作用，建立二次耦合模型。
• 将模型预测结果与模拟结果进行对比，以验证能量转移机制。
• 研究n = 1模态在从线性到非线性增长阶段的空间结构演化。
• 识别非线性激发n = 1模态最有效的位置，以及导致该局域化的主导驱动模态及其同相特性。

实验结果

研究问题

• RQ1 尽管在理论上属于次主导模态，为何次主导的低n模态（特别是n = 1）在ELMs的早期非线性阶段仍能显著增长？
• RQ2 何种物理机制使得能量能从主导模态有效传递至n = 1等低n成分？
• RQ3 n = 1模态的空间结构如何从线性阶段过渡到非线性阶段？
• RQ4 简单的二次耦合模型能否再现观测到的环向傅里叶谱的非线性演化？
• RQ5 在等离子体中何处非线性激发n = 1模态最为有效？其局域化由什么因素决定？

主要发现

• 尽管在理论上属于线性次主导模态，n = 1模态在非线性阶段仍增长至与线性主导模态相当的能量水平。
• 二次耦合模型成功再现了环向谱的早期非线性演化，证实了通过二阶非线性相互作用实现的能量转移。
• n = 1模态的空间结构从线性增长阶段中延伸至等离子体核心的刚性增长形态，转变为非线性增长阶段中在等离子体边界处的局域化结构。
• 非线性驱动的n = 1模态局域化于等离子体边界，该区域主导驱动模态达到最大值且同相，从而实现了高效的能量传递。
• 该模型为近期实验观测中ELMs期间磁测量结果出现强低n成分的现象提供了物理解释。
• 空间结构的转变凸显了非线性耦合动力学在塑造ELMs最终模态结构中的关键作用。