Improved n=1 Empirical Error Field Penetration Threshold Scaling with Ohmic and L-Mode Conventional Tokamak Plasma Discharges

本文基于扩展的传统欧姆和L模托卡马克放电数据库，提出了一种改进的n=1经验误差场穿透阈值标度律，该标度律具有更高的拟合精度和更低的误差，从而为未来传统托卡马克的工程公差与设计提供更可靠的依据。

要点

想象一下，托卡马克（一种旨在产生聚变能的装置）就像一个由强大磁场固定的、巨大的、不可见的超热气体漩涡。理想情况下，这个漩涡是完美对称的，就像一个光滑旋转的陀螺。然而，在现实世界中，固定它的磁铁并不完美。它们存在微小的倾斜、位移和缺陷。这些缺陷会产生“误差场”——微小的、不需要的磁场涟漪，可能扰乱气体的平滑旋转。

如果这些涟漪变得过强，它们可能导致漩涡形成一个“结”（即磁岛），并卡在原地。一旦卡住，整个系统就会崩溃，导致被称为“破裂”的突然停机。这对机器及其建造者来说都是坏消息。

问题：多强才算太强？
工程师需要知道确切的极限：这些磁缺陷能大到什么程度机器才会崩溃？ 如果设定的极限太低，他们就必须以不可能达到的精度建造机器，使其造价天文数字般高昂且建造缓慢。如果设定的极限太高，机器可能会崩溃。

多年来，科学家们一直试图建立一条“经验法则”（标度律），根据机器的大小和气体行为来预测这一极限。但旧规则有些不可靠，就像边缘模糊的地图。

解决方案：更清晰的地图
本文提出了一张新的、升级版的“地图”（经验标度律），它更加清晰可靠。以下是他们如何利用简单类比做到的：

• 清洗数据： 研究人员查阅了全球托卡马克（如 DIII-D、JET 和 KSTAR）过去实验的庞大数据库。他们决定只关注机器中特定的“天气”类型：“欧姆”模式和"L 模”条件。他们排除了"H 模”，因为那种状态就像一座坚固的堡垒——极难破坏，因此无助于我们理解机器的最薄弱点。通过聚焦于脆弱状态，他们找到了真正的危险区域。
• 添加新要素： 他们加入了来自两台特定机器的新数据：J-TEXT（体积较小，运行电流较低）以及更多来自 JET（体积巨大，与未来的 ITER 机器相似）的数据。这就像在汽车安全数据库中添加新的试驾记录。你需要在数据中包含小型汽车和巨型卡车，才能知道安全规则如何适用于未来可能建造的任何车辆。
• 更优的数学： 他们使用了一种更复杂的数学方法，来寻找机器尺寸、磁场强度、气体密度和电流之间的关系。他们发现，等离子体电流（流经气体的电量）是一个此前未充分考量的关键因素。

新发现
新的“经验法则”告诉我们：

1. 更高的密度是你的朋友： 将更多气体压缩进机器中，会使误差场更难引发崩溃。
2. 更大的机器具有惊人的韧性： 大型机器（如未来的 ITER）比我们之前认为的更能承受更大的磁缺陷。
3. 电流至关重要： 流经等离子体的电流大小会改变机器对这些误差的反应方式。

这对未来的意义
本文特别关注ITER项目，这是一个正在建设中的大型国际聚变实验。利用他们新的、更清晰的地图，研究人员运行了数百万次模拟（就像用略微不同的初始条件运行一百万次天气预报）。

结果： 他们发现 ITER 的状况比我们想象的要好得多。磁误差的“危险区域”比 ITER 预期存在的实际缺陷要远得多。

• 旧地图： 暗示 ITER 有相当大的几率会绊倒自己的鞋带（即发生锁定模态）。
• 新地图： 显示这种情况发生的几率微乎其微（在最可能的场景下，概率低于百万分之一）。

核心结论
这篇文章不仅仅是在说“聚变很难”。它为工程师提供了一把更自信、更精确的尺子，用来测量机器的公差。因为新规则显示机器对磁误差更具鲁棒性，工程师可能不需要以如此极端、昂贵的精度来制造磁铁。这可以在保持机器安全的同时节省时间和金钱。

简而言之：他们利用更好的数据和更聪明的数学，清理并修正了一张关于磁安全极限的模糊、令人困惑的地图，并发现聚变发电厂的未来比我们之前认为的更安全、更可实现。

技术摘要

以下是论文《改进的 n=1 经验误差场穿透阈值标度律：基于欧姆和 L 模常规托卡马克等离子体放电》的详细技术总结。

1. 问题陈述

托卡马克需要精确的轴对称性才能稳定运行。然而，制造缺陷（线圈倾斜、位移和未对准）会产生非轴对称的误差场（EFs）。如果这些场超过临界穿透阈值，它们可能触发磁重联，形成磁岛，并导致磁岛“锁定”到器壁上。这将引发锁定模，随后导致等离子体破裂。

预测该阈值对于以下方面至关重要：

• 工程公差：确定建造期间允许的线圈未对准程度。
• 成本与进度：过于严格的公差会增加建造成本和时间；过于宽松的公差则可能导致装置失效。
• 未来装置：像ITER和SPARC这样的项目具有极高的存储能量，使得破裂特别危险。

先前的经验标度律存在以下问题：

• 拟合质量低：决定系数（$R^2$）低至 0.38。
• 不确定性高：ITER 的预测不确定性高达 48%。
• 数据不一致：混合包含了 H 模、L 模和欧姆等离子体、球形托卡马克（NSTX）以及不同装置间变化的本征误差场建模水平。
• 参数空间有限：缺乏大型大半径装置（如 ITER）和特定密度区域的数据。

2. 方法论

作者通过优化数据库和拟合方法，开发了一种更新的经验标度律。

A. 数据库优化

• 区域选择：研究将数据集限制在常规（非球形）托卡马克的欧姆和 L 模运行区域。排除了 H 模，因为它对误差场穿透更具抵抗力，而作者旨在模拟“最危险”的区域。排除了饱和欧姆约束（SOC）放电，以避免混淆密度依赖性。
• 装置包含：
  • 新增：来自J-TEXT的新数据（提供低电流/密度点和大半径尺寸）以及来自JET的扩展数据集（21 次放电，对比之前的 4 次，涵盖大半径和高等离子体电流）。
  • 排除：NSTX（球形托卡马克）和H 模放电。
• 标准化：
  • DIII-D：从分组平衡近似转变为对每次放电使用 SURFMN 代码进行单独的平衡重构，以更准确地模拟本征误差场。
  • J-TEXT：使用TokaMaker代码重构平衡（因为缺乏常规磁重构），确保与用于所有其他装置的GPEC（广义扰动平衡代码）兼容。
  • 指标：使用了主导模重叠（$\delta$），这是一种通过 GPEC 计算、考虑等离子体响应且与坐标系无关的指标。

B. 标度律方法论

• 回归技术：采用前向选择最小二乘回归来确定最佳对数线性幂律。
• 参数选择：扩展了参数空间，超越了以往研究，明确包含对等离子体电流（$I_p$）、环向场（$B_T$）、大半径（$R_0$）、**密度（$n_e$）和归一化压力（$\beta_n/\ell_i$）**的依赖关系。
• 验证：使用条件数检测多重共线性（在保留 5 个变量时停止，以保持条件数 < 20），并使用核密度估计（KDE）比较普通最小二乘法（OLS）与加权最小二乘法（WLS），以避免对数据丰富区域的过采样。
• 不确定性量化：假设锁定模识别时间存在 10% 的基准不确定性，通过蒙特卡洛模拟传播，生成预测的概率密度函数（PDFs）和累积分布函数（CDFs）。

3. 主要贡献

1. 新标度律：推导了两个新的标度方程（方程 9 用于 OLS，方程 10 用于 WLS），统计指标显著改善。
   • 方程 9（OLS）：$\delta \propto (\beta_n/\ell_i)^{0.25} I_p^{-0.97} R_0^{1.88} n_e^{0.77} B_T^{0.20}$
   • 方程 10（WLS）：$\delta \propto (\beta_n/\ell_i)^{0.13} I_p^{-1.01} R_0^{1.57} n_e^{0.56} B_T^{0.30}$
2. 包含显式$I_p$依赖：与依赖隐式依赖关系的先前标度律不同，新模型明确包含等离子体电流，揭示了较低的等离子体电流会增加对误差场的抵抗力。
3. 环向场依赖关系的转变：$B_T$的标度关系从先前模型中的强负依赖转变为小的正依赖，这是由包含$I_p$和移除球形托卡马克数据所驱动的。
4. 鲁棒性测试：证明了该标度律在移除单个装置（如 DIII-D、J-TEXT、C-MOD）时具有鲁棒性，而“最小”装置子集（如 JET+KSTAR+C-MOD）则导致严重的过拟合。

4. 结果

• 拟合质量改善：新标度律的$R^2$达到0.63（OLS）和0.66（WLS），相比之前的数值（0.38–0.58）有显著改善。
• 不确定性降低：
  • 约化卡方：从 2020 年模型的 17.0 改善至12.99。
  • 平均对数绝对误差：从 1.336 改善至1.306（相当于将预测不确定性从$\pm33.6\%$降低至$\pm30.6\%$）。
• 参数依赖关系：
  • 密度（$n_e$）：确认了正依赖关系（$n_e^{0.56-0.77}$），与许多单个装置一致，但不同于某些理论模型。
  • 大半径（$R_0$）：强正依赖关系（$R_0^{1.57-1.88}$），表明大型装置更具抵抗力。
  • 等离子体电流（$I_p$）：强负依赖关系（$I_p^{-0.97}$ 至 $-1.01$），表明较高的电流使等离子体更容易发生锁定。
• 装置特异性异常值：
  • NSTX：新标度律显著低估了 NSTX 的误差场耐受性（预测$\delta \sim 10^{-5}$，而实际约为$\sim 10^{-4}$），这可能是由于其球形性质和高$q_{95}$剖面所致，证实它应从常规托卡马克标度律中排除。
  • C-MOD：来自 C-MOD 的高密度数据对于捕捉非单调密度标度行为至关重要。

5. 对 ITER 的意义和影响

• ITER 预测：
  • 与先前模型相比，新标度律预测 ITER 的误差场穿透阈值要高得多。
  • ITER 最可能的本征$\delta$：$0.38 \times 10^{-4}$。
  • 预测阈值（方程 10）：$4.01 \times 10^{-4}$（带 1 个标准差的不确定性范围）。
  • 锁定概率：计算得出 ITER 在基线场景下发生锁定的概率极低（对于最可能的本征场为0.0001%，即使对于最大可能的本征场也仅为22.3%）。
• 工程影响：
  • 结果表明，ITER 目前的线圈对准工程公差是保守且充足的。
  • 该标度律为设计未来装置（如 SPARC）和优化现有机器运行提供了更可靠的工具，通过基于本征误差场的蒙特卡洛模拟来预测模锁定的概率。
• 未来工作：作者建议增加来自JT-60SA（大半径）和ASDEX-U（高$\beta_n$）的数据将进一步完善标度律。此外，开发基于具有明确旋转测量的放电的标度律可能会进一步提高预测能力，因为旋转是误差场物理中的主导因素。

总之，本文建立了一种更稳健、物理一致且统计可靠的 n=1 误差场穿透经验标度律，显著增强了人们对下一代常规托卡马克设计和运行的信心。