反应器高密度条件下小边缘局域模的台基起源与平行能量通量外推：ITER与SPARC的可行性验证

为应对未来聚变装置中边缘局域模(ELM)带来的强热流挑战，研究人员聚焦于高密度、高边缘碰撞频率下产生的小型ELM(简称小ELM)开展研究。通过结合DIII-D和ASDEX Upgrade的实验测量与BOUT++代码模拟，该工作首次揭示了小ELM起源于台基底部靠近分界面区域，并仅微弱扰动台基结构。基于此发现，对Eich的ELM峰值平行能量通量模型进行了ad hoc修正，成功复现了实验数据。首次将该模型外推至ITER和SPARC，预测的小ELM能量通量值低于钨单块表面的熔化阈值，为未来聚变装置的安全运行提供了重要依据。

想象一下，在追求“人造太阳”——可控核聚变的道路上，科学家们不仅要驯服上亿度的等离子体，还要解决一个看似微小却威力巨大的难题：边缘局域模(ELM)。在托卡马克（一种磁约束聚变装置）中，等离子体边缘的磁面会周期性地失稳崩塌，瞬间将大量能量和粒子抛向装置的第一壁，尤其是偏滤器靶板。对于未来的聚变电站（FPPs）如ITER和SPARC，未经控制的大型I型ELM所产生的热负荷足以损坏昂贵的钨制偏滤器组件。因此，寻找既能维持良好等离子体性能，又能产生“温和”ELM的等离子体运行区间，是实现聚变能源的关键挑战之一。近年来，一种在高边缘密度、高碰撞频率条件下产生的小型ELM（或称高碰撞频率小ELM）等离子体运行模式引起了广泛关注。与常见的大振幅I型ELM相比，这种小ELM具有频率高、幅度小、等离子体能量损失低的特点，且与部分脱靶偏滤器（一种通过辐射等方式分散热负荷的技术）的功率排废气解决方案兼容，展现出巨大潜力。然而，这种小ELM的物理起源、其对等离子体台基（边缘压力陡峭的区域）结构的具体影响，以及将其产生的峰值热负荷外推到ITER等下一代大型装置是否安全，都还存在不确定性。本文发表在《Journal of Plasma Physics》上的研究，正是为了系统回答这些问题，通过综合实验分析与数值模拟，首次为评估小ELM作为未来聚变装置可行运行模式的潜力，提供了关键的定量依据。

为开展这项研究，作者团队采用了多管齐下的技术方法。首先，基于DIII-D和ASDEX Upgrade托卡马克装置，建立了高密度小ELM实验案例数据库。关键诊断技术包括：1) 高速红外热成像(IRTV)，用于测量偏滤器靶板的峰值热流，时间分辨率足以分辨单个ELM事件；2) 汤姆逊散射(TS)系统和条件平均技术，用于获取ELM周期内的台基电子密度和温度剖面；3) 快速往复探针(RCP)，用于在刮削层(SOL)中测量ELM期间的电子温度和密度；4) 宽场束发射光谱(BES)，用于分析边缘区域的密度涨落。其次，利用BOUT++六场流体模型对代表性小ELM案例进行了线性和非线性模拟，以探究其稳定性、模式结构及其对台基的扰动。最后，基于实验和模拟揭示的物理图像，对Eich等人提出的ELM峰值平行能量通量经验模型进行了修正，并使用修正后的模型对DIII-D和ASDEX Upgrade的实验数据进行了首次系统性基准测试，进而外推到ITER和SPARC的参考情景。

研究结果主要包含以下几个方面：

3.1. 小ELM能量含量与功率平衡

通过分析等离子体储能信号，发现研究中的小ELM平均仅携带了等离子体储能的0.7%，远低于典型的I型ELM。通过比较储能下降与基于IRTV测量计算的偏滤器沉积能量，发现大约80%的ELM能量被沉积到偏滤器，两者在测量误差范围内基本一致。

3.2. 小ELM导致的台基结构动力学

对TS数据的条件平均分析表明，小ELM对台基结构的扰动非常微弱，其电子压力剖面在ELM期间和ELM间隙几乎重叠。相对变化(Δp_e/p_e)的最大值（约0.2）出现在靠近分界面的台基底部区域。相比之下，I型ELM则导致台基大面积坍塌，Δp_e/p_e在台基顶部区域可达0.6。BES测量也显示，在典型的小ELM放电中，归一化密度涨落的峰值位于分界面附近，支持了增强的边缘湍流活动图像。

3.3. 小ELM期间台基稳定性的BOUT++模拟

使用BOUT++代码对DIII-D的一个小ELM案例进行模拟。线性稳定性分析表明，台基底部区域对高环向模数(n)的电阻性气球模不稳定。线性模态结构显示，压力的扰动从归一化磁通Ψ_n~ 0.98延伸到分界面外(Ψ_n~ 1.01)。非线性模拟进一步证实，模态在台基底部发展并饱和，只压平了该区域的压力梯度，并未深入台基顶部，这与实验观测到的小ELM仅微弱扰动台基的结果一致。

3.4. 通过SOL局域测量推断小ELM的台基起源

利用OMP附近快速RCP测量的SOL电子压力数据，并通过指数衰减函数外推至分界面，推断出小ELM起源于Ψ_n~ 0.99到Ψ_n= 1.0之间的狭窄区域，即台基底部/分界面附近。这与之前对I型ELM的研究（通常起源于台基中部/顶部）形成鲜明对比。

3.5. 用实验数据基准测试修正的ELM能量通量模型

研究人员测试了Eich的ELM峰值平行能量通量模型。如果像处理I型ELM一样，采用台基顶部的电子压力(p_e)作为输入，模型会严重高估（超过8倍）DIII-D小ELM的实验测量值。然而，基于前述实验和模拟结果，将模型中的p_e修改为取自台基底部(Ψ_n= 0.995)的值后，模型成功复现了DIII-D和ASDEX Upgrade的小ELM实验数据，平均偏差在40%以内，大部分数据点落在两倍偏差范围内。这为模型的可靠性提供了支持。

3.6. 将峰值ELM能量通量外推至SPARC和ITER

利用上述验证过的修正模型，首次对ITER (15 MA D-T情景) 和SPARC装置的高密度小ELM峰值平行能量通量进行了外推估计。对于ITER，若假设p_e取自分界面，预测值约为0.43 MJ m^-2；若p_e取自稍靠内的Ψ_n= 0.98，预测值约为2.3 MJ m^-2。对于SPARC，使用分界面参数预测值约为0.7 MJ m^-2。这些预测值均低于ITER钨单块偏滤器靶板表面因ELM导致熔化的名义阈值（对应平行能量通量约为12 MJ m^-2），初步表明该运行模式在材料耐受性方面可能是可接受的。

结论与讨论部分对全文发现进行了总结。本研究的核心在于通过实验与模拟相结合，明确了高密度小ELM起源于台基底部/分界面附近，仅微弱扰动台基结构。基于此物理图像修正的ELM能量通量模型，能够较好地复现当前装置（DIII-D, ASDEX Upgrade）的小ELM实验数据。利用该模型首次对ITER和SPARC进行的预测显示，小ELM产生的峰值能量通量可能低于钨靶板的熔化阈值。这标志着在评估小ELM作为下一代聚变装置可行运行模式的道路上迈出了重要一步。作者在讨论中也指出了外推中存在的不确定性和注意事项，例如磁场线入射角、靶板前缘耐受性、以及ELM能量需与ELM间隙热流叠加考虑等。未来的研究需要结合更多装置和诊断的数据来进一步验证和完善模型，并重点考察在该小ELM模式下，同时实现等离子体良好约束和偏滤器完全脱靶的可行性。这项研究为未来聚变装置应对功率排废气挑战提供了一个有希望且经过初步定量评估的潜在解决方案。