托卡马克破裂-ITER的达摩克利斯之剑
托卡马克破裂:ITER的达摩克利斯之剑
聚变堆里最吓人的东西,不是1.5亿度的等离子体,而是它"突然消失"的那一瞬间。一个正在运行的托卡马克,如果等离子体突然失去约束,几兆安培的电流会在几毫秒内猝灭,把相当于一百吨机车以五百公里时速撞墙的能量砸向第一壁——这就是"破裂"(disruption)。更糟的是,这种剧烈变化会把一小撮电子加速到接近光速,变成一束能熔化金属的"逃逸电子"。
ITER组织自己的口径是:2018年,理事会下设的科学与技术咨询委员会把等离子体破裂定性为"对ITER任务的严重威胁";到了2024年,官网直接用了"存在性威胁"这个词。需要说清楚的是,破裂不是"安全级"问题——ITER设计上能承受大约三千次破裂,真空室不会被炸穿——但每一次破裂都会损伤第一壁,而逃逸电流只要约300千安集中打在一个点上就能把壁熔化(近年的缓解研究对"局部损失"的耐受红线更严,有估算到约150千安;ITER官方还设了一条红线:逃逸电流超过2兆安即不可容忍),换一次壁就是装置可用率的大幅下降。所以准确的说法是:破裂是ITER的"头号运行/可用性风险",物理上不致命,运行上致命。
2025年,国际托卡马克物理活动(ITPA)的MHD、破裂与控制专题组发了一篇大综述,发表在Nuclear Fusion第65卷103001号(DOI 10.1088/1741-4326/ade7a0),通讯作者是ITER组织的I. Bandyopadhyay,作者名单约79人、横跨十几个国家——ITER组织、德国IPP、瑞士EPFL、美国哥伦比亚/PPPL/GA/MIT/ORNL、法国CEA、意大利那不勒斯、英国UKAEA、韩国NFRI、俄罗斯库尔恰托夫都有人署名。这篇综述是自2007年Hender那篇经典综述以来15年的总结,核心判断一句话:反应堆级装置对破裂和逃逸电流的容忍度远低于现有装置,必须完全避免或最小化。
那逃逸电子到底是什么?它怎么来的?
三条物理路径。第一条叫Dreicer产生,得名于1959年发表论文的物理学家H. Dreicer——注意,他当时在洛斯阿拉莫斯科学实验室,不是常被误传的贝尔实验室。机制是:当电场够强,超过电子所受库仑摩擦的峰值,麦克斯韦尾巴里那些高能电子就会被持续加速、一路狂奔。第二条叫雪崩,1997年Rosenbluth和Putvinski给出理论:一个高能逃逸电子通过大角碰撞把次级电子踢过临界动量,次级又踢出更多,形成指数爆炸。大装置里雪崩是主导——ITER一次大破裂里雪崩放大可以到e的50次方。第三条叫热尾:破裂时等离子体冷却太快,电子来不及碰撞热化,留下一截超热的"尾巴",随后感应电场把这截尾巴变成逃逸种子。
逃逸电子有多危险?JET的铍第一壁,建模显示在逃逸电子流密度2到12千安每平方米就会熔化。历史上T-4装置(1971-73)的钨限制器被熔穿过;Alcator C-Mod上逃逸束打在诊断电缆上喷出了大约2立方厘米的不锈钢和铜。ITER建模给出:集中的逃逸电流约300千安就能熔化壁面,扩散情形约1.9兆安。
那么怎么对付?这里就要讲缓解技术的演进。
早期方案叫MGI——大量气体注入。把高压氖或氘气往等离子体边缘一喷,靠辐射把能量耗散掉。但MGI有硬伤:气体在边缘就电离了,被磁场挡住到不了芯部,吸收率低,很难达到抑制逃逸电子所需的临界密度。不过要说明,ITPA联合实验(DIII-D、C-Mod、TEXTOR、FTU、KSTAR)实测产生逃逸所需的阈值电场比碰撞理论高约5到10倍——这意味着抑制所需的临界密度比理论(Rosenbluth密度)预测低,抑制逃逸可能比理论担心的容易一些,这个"临界密度"不是铁律。所以ITER选了下一代方案:SPI——碎裂弹丸注入。把低温固态弹丸加速到超音速(ITER设计超过每小时1800公里,JET实际150到600米每秒),撞上斜面碎裂成碎片,碎片比气体穿透深得多、能直达芯部(SPI冷锋能到q=1面,MGI只能到q=2)。SPI是橡树岭国家实验室(ORNL)的Baylor、Shiraki、Commaux等人开发的,先在DIII-D上验证,ITER最终选定它而放弃MGI。
ITER的破裂缓解系统(DMS)设计已经走完最终设计评审——2024年3月26到29日,主持人是通用原子能的George Sips,超三万页文档、一百多位专家评审通过。方案是24个赤道端口注入器加3个上端口注入器,弹丸大约"酒瓶塞"大小(直径28.5毫米、长57毫米),整个注入过程约12毫秒。注入策略用"交错方案"(这是建模场景里的参数,不是ITER官方硬件规范):先8发纯氢弹丸升高密度、压低电子温度以抑制逃逸种子生成,再3发掺氖弹丸靠辐射耗散能量。建模估计这套方案理论上能把逃逸电流限制在亚兆安量级(不同模型的估计在约0.24兆安到"仍不可接受"之间)。但要注意,瑞典Chalmers的Vallhagen等人2024年用DREAM模拟活化氘氚等离子体后给出的结论更悲观:即使交错注入,"也找不到可接受的逃逸电流"。首次等离子体注入目标2034年。
但别以为这就解决了。这里有个刺眼的内部矛盾:ITER的交错方案依赖约30帕斯卡的中性气压窗口,而JET实验验证过的"低Z良性终止"——用氘气二次注入让兆安级逃逸束安全消失(Paz-Soldan等人2021年那篇Nuclear Fusion)——发生在0.5到1帕斯卡。两个数量级的窗口冲突,说明抑制雪崩所需的路径和实验验证过的良性终止路径是两条不同的物理通道,还没有统一方案。DIII-D和TCV的实验明确显示中性气压窗口只有0.5到1帕,而ITER交错方案会产生30帕。逃逸电子缓解,至今没有"已解决"三个字。
模拟这边,三把大刀。M3D-C1(普林斯顿PPPL的Ferraro和Jardin开发)用高阶连续有限元,是SPARC破裂工程设计的核心工具;NIMROD(威斯康星大学主导)做全MHD;JOREK(德国IPP主导的国际社区代码)做约化MHD,欧洲和ITER深度使用。三者在垂直位移事件(VDE)上做过三方基准,壁力和晕电流实现了可量化的吻合——不过三种代码用了不同的物理模型(JOREK是约化MHD、另两个是全MHD)且对初值敏感,所以"吻合"不等于"一致"。
SPARC的非缓解破裂模拟是个标志性工作。2025年Physics of Plasmas第32卷082505号,作者是R. Datta、C. Clauser、N.M. Ferraro、C. Liu、R. Sweeney、R.A. Tinguely——机构是MIT PSFC加PPPL加CFS。这里要更正一个常见误解:V.A. Izzo不在这个作者名单里,她是另一篇用NIMROD做SPARC缓解模拟的论文的第一作者。这篇M3D-C1工作首次系统比较了不同的逃逸种子源(Dreicer、氚β衰变、康普顿散射),结论是Dreicer源主导,SPARC非缓解破裂会产生约7兆安的逃逸平台电流——8.7兆安的主等离子体电流,八成变成了一束相对论电子。要注意两点:一是M3D-C1是2D模拟、压住了3D随机场的逃逸损失,所以7兆安是保守上界;二是对"哪种种子源主导",不同代码结论不同——M3D-C1说Dreicer主导,而NIMROD和DREAM的SPARC研究认为热尾种子更关键。这束电子轰到第一壁上就是熔融和蒸发。
模拟与实验之间还有裂缝。RWTM(电阻壁撕裂模)那派预测ITER破裂的热猝灭时间长达70到125毫秒,比现役装置长得多,所以"ITER破裂可能比担心的更良性"——这跟传统评估相反,还没经过实验验证。而JOREK因为用日冕平衡杂质模型,给出的热猝灭时间比M3D-C1和NIMROD长两倍。2D轴对称模拟还漏掉了3D的MHD-逃逸电子相互作用——SPARC的3D模拟显示MHD活动反而增强逃逸电子产生。
AI能不能在破裂发生前拦住它?
这方面进展很快但门槛很高。2019年哈佛/普林斯顿的Kates-Harbeck、Svyatkovskiy和Tang在Nature上发了FRNN(融合循环神经网络)——注意卷页是Nature 568卷526到531页,不是常被记错的570卷528页。它用DIII-D和JET的数据训练,展示了跨装置迁移能力,能在破裂前30毫秒给出预警——这个量级的预警时间,正好是ITER缓解系统在考虑预警需求时关注的范围。2024年Nature上Seo、Kolemen等人更进一步:用深度强化学习训练一个控制器,在DIII-D上能提前最多300毫秒预测并主动规避撕裂模,这是从"预测"到"避免"的里程碑。
但工程门槛摆在那儿。MIT的Rea团队在DIII-D做的实时机器学习预测器(随机森林)虽然运行了900多炮,误报率却高到达不到ITER要求。EAST那边2025年报道了可解释的XGBoost预测:2094次破裂炮训练,30毫秒前真阳率96.7%、假阳率6.58%——数字很漂亮,但ITER的缓解系统需要"近零误报",因为误触发缓解本身就会损伤装置。更根本的问题:ITER没有历史破裂数据可以训练,跨装置外推的可靠性存疑。聚变AI最危险的失败模式,不是不工作,而是在有数据的区域表现够好、从而对没数据的区域产生虚假信心。
私营聚变的到来改变了整个时间表。SPARC(CFS/MIT,8.7兆安)用6阀MGI加被动逃逸缓解线圈;ARC(CFS商业堆)干脆把目标定为"每天容忍一次缓解后的破裂"——几十秒内重启、不断电——这是私营聚变的运行经济性逻辑,跟ITER"尽量别破"的哲学完全不同。英国STEP(球形托卡马克原型堆)那边,DREAM模拟显示满功率破裂会不可避免地产生接近全电流转化的逃逸电子束,两阶段SPI之后仍残余约13兆安——不可接受。前JET任务组长Guy Matthews给英国议会的证词说,STEP等离子体能量约为MAST-U的5000倍、JET的50倍,物理不确定性带来的失效风险极高。
中国在这张图里的位置,散落在这五个维度里,从没被整合过。
EAST(中科院等离子体所)做了启动期逃逸电子研究(判据E_loop/E_c大于10)、MGI/SPI缓解实验、双阈值场观测——2020年那篇Nuclear Fusion用实验证实了Dreicer和雪崩两个阈值电场,雪崩阈值约是Dreicer探测阈值的1.2倍,这是中国团队独立提出的机制级贡献。EAST还有可解释的AI破裂预测,30毫秒前真阳率96.7%。J-TEXT(华中科技大学)陈忠勇团队在小装置上实现了100%逃逸抑制、28兆安每秒的耗散速率——这是当时全球最高的耗散速率——还研发了电磁弹丸注入(EMPI)。中科大、华中科大、中科院等离子体所合作做了杂质辐射触发撕裂模破裂的非线性模拟(2023年Nuclear Fusion)。中国的装置类型覆盖相当全面——大装置实验(EAST)、小装置机制研究(J-TEXT)、非线性模拟(中科大/华科)、可解释AI预测(EAST)四件套都齐。但要诚实:欧洲的链其实同样完整甚至更完整——AUG、TCV、MAST-U还在运行,JOREK本身就是欧洲主导的国际代码,欧洲也有多组AI预测研究;韩国KSTAR也有大装置加AI。所以准确的说法是,中国在装置类型覆盖上进入了第一梯队,而不是"唯一"。
但中国也有差距。CFETR——中国自己的下一代聚变工程堆——虽然已有电流猝灭期逃逸电流模拟、氖MGI缓解模拟等初步工作,但尚未形成像ITER DMS那样集成的破裂缓解策略。中科大/华科的破裂模拟还没有进入M3D-C1/NIMROD/JOREK的国际三方基准体系。EAST那96.7%真阳率也没对接上"近零误报"的工程需求。
现在把镜头拉回那个达摩克利斯之剑。破裂这个从2007年Hender综述里的"物理问题",到2025年ITPA专刊定位的"存在性威胁",转折点是ITER运行需求和私营聚变几乎同时到来。三条线还在较劲:模拟线复现了Greenwald密度极限、预测了SPARC的7兆安逃逸平台,但RWTM"更良性"与JOREK"两倍长TQ"互相矛盾;理论线给出三机制,1D模型却系统性高估逃逸电流;AI线预测准确但没到工程门槛。缓解"是否已解决"仍未收敛——ITER选定的交错SPI与实验验证的低Z良性终止在物理通道上冲突。
开放问题很直白:2034年ITER首次等离子体注入之前,逃逸电子缓解能不能从"建模上限"变成"工程保证"?这个问题的答案,将决定聚变装置从"演示"走向"可用"的最后一公里。而JET这台唯一能产生并安全终止兆安级逃逸束的装置已经关停,实验验证面在收缩。剑还在头顶悬着。
---
参考文献
- I. Bandyopadhyay et al.(ITPA MHD, Disruption and Control Topical Group,约79位作者), *MHD, disruptions and control physics: Chapter 4 of the special issue: on the path to tokamak burning plasma operation*, Nuclear Fusion 65 (2025) 103001; DOI 10.1088/1741-4326/ade7a0
- T.C. Hender et al.(ITPA), *MHD stability, operational limits and disruptions*, Nuclear Fusion 47 (2007) S128-S202
- M. Lehnen et al.(ITER组织等), *Disruptions in ITER and strategies for their control and mitigation*, J. Nucl. Mater. 463 (2015) 39-48
- ITER技术报告 ITR-24-005(Disruption Management)
- R. Datta, C. Clauser, N.M. Ferraro, C. Liu, R. Sweeney, R.A. Tinguely(MIT PSFC/PPPL/CFS), *Coupled 2D MHD and runaway electron fluid simulations of SPARC disruptions*, Physics of Plasmas 32 (2025) 082505; DOI 10.1063/5.0272430
- V.A. Izzo, B. Stein-Lubrano, A. Battey, R. Sweeney, C. Hansen, R.A. Tinguely, *Disruption mitigation modeling for the SPARC tokamak*, Physics of Plasmas 32 (2025) 042507
- H. Bergström, S.-J. Liu, V. Bandaru, M. Hoelzl(IPP Garching/Chalmers), *Introduction of a 3D global non-linear full-f particle-in-cell model for runaway electrons in JOREK*, Plasma Phys. Control. Fusion 67 (2025) 035004; DOI 10.1088/1361-6587/adaee7
- C. Sommariva, E. Nardon, P. Beyer, M. Hoelzl, G.T.A. Huijsmans, D. van Vugt(Aix-Marseille/CEA/IPP/TU/e), *Test particles dynamics in the JOREK 3D non-linear MHD code*, Nuclear Fusion 58 (2018) 016043
- H. Dreicer(Los Alamos Scientific Lab), *Electron and Ion Runaway in a Fully Ionized Gas*, Phys. Rev. 115 (1959) 238
- M.N. Rosenbluth, S.V. Putvinski(ITER Joint Work Site), *Theory for avalanche of runaway electrons in tokamaks*, Nuclear Fusion 37 (1997) 1355
- P. Helander, H. Smith, T. Fülöp(CCFE/Chalmers), *Electron kinetics in a cooling plasma*, Physics of Plasmas 11 (2004) 5704
- M. Hoppe, O. Embreus, T. Fülöp(Chalmers), *DREAM: A fluid-kinetic framework for tokamak disruption runaway electron simulations*, Comput. Phys. Commun. 268 (2021) 108098
- C. Paz-Soldan, C. Reux et al.(Columbia/CEA/IPP/ITER/GA/ORNL/PPPL/UKAEA), *A novel path to runaway electron mitigation via deuterium injection and current-driven MHD instability*, Nuclear Fusion 61 (2021) 116058
- C. Reux et al.(JET), *Demonstration of Safe Termination of Megaampere Relativistic Electron Beams in Tokamaks*, Phys. Rev. Lett. 126 (2021) 175001
- D. Shiraki, N. Commaux, L.R. Baylor et al.(ORNL), *Dissipation of post-disruption runaway electron plateaus by shattered pellet injection in DIII-D*, Nuclear Fusion 58 (2018) 056006
- S. Jachmich et al.(ITER组织/JET), *Shattered pellet injection experiments at JET in support of the ITER disruption mitigation system design*, Nuclear Fusion 62 (2022) 026012
- L.R. Baylor et al.(ORNL), *Design and performance of shattered pellet injection systems for JET and KSTAR*, Nuclear Fusion 61 (2021) 106001
- N.M. Ferraro, S.C. Jardin et al.(PPPL), *The M3D-C1 approach to simulating 3D 2-fluid magnetohydrodynamics*, J. Phys.: Conf. Ser. 125 (2008) 012044
- J. Kates-Harbeck, A. Svyatkovskiy, W. Tang(Harvard/Princeton/PPPL), *Predicting disruptive instabilities in controlled fusion plasmas through deep learning*, Nature 568 (2019) 526-531; DOI 10.1038/s41586-019-1116-4
- J. Seo et al. + E. Kolemen(Princeton/PPPL/中央大学), *Avoiding fusion plasma tearing instability with deep reinforcement learning*, Nature 626 (2024) 746-751; DOI 10.1038/s41586-024-07024-9
- C. Rea, K.J. Montes, R.S. Granetz, R.A. Tinguely(MIT PSFC), *A real-time machine learning-based disruption predictor in DIII-D*, Nuclear Fusion 59 (2019) 096016
- N.W. Eidietis et al.(ITPA), *The ITPA disruption database*, Nuclear Fusion 55 (2015) 063030
- S. Ratynskaia et al.(KTH/IPP/CEA等), *Runaway electron-induced plasma facing component damage in tokamaks*, Plasma Phys. Control. Fusion 68 (2026) 053501
- 陈忠勇 等(华中科技大学), *Overview of runaway current suppression and dissipation on J-TEXT tokamak*, Plasma Science and Technology 24 (2022) 124009
- X. Zhu, L. Zeng(曾龙), Y. Liang et al.(ASIPP), *Observation of two threshold fields for runaway-electron generation in tokamaks*, Nuclear Fusion 60 (2020) 084002
- S.B. Zhao et al.(ASIPP), *Characteristics of plasma disruption mitigation achieved by MGI and SPI on EAST*, Nuclear Fusion 65 (2025) 016048
- S. Zeng, P. Zhu, R. Zhou, D.F. Escande(中科大/HUST/ASIPP), *Triggering of tearing instability by impurity radiation*, Nuclear Fusion 63 (2023) 016026
- Y.L. Yu, Z.Y. Chen(华中科技大学J-TEXT), *Optimal design and experimental testing of EMPI system*, Fusion Eng. Des. 208 (2024) 114701
- H.R. Strauss, B.C. Lyons, M. Knölker(HRS Fusion/GA), *Resistive wall tearing mode disruptions in DIII-D and ITER*, Physics of Plasmas 29 (2022) 112508
- O. Vallhagen, T. Fülöp et al.(Chalmers), *Runaway electron dynamics in ITER disruptions with shattered pellet injection*, Nuclear Fusion 64 (2024) 086033; arXiv:2401.14167
- Liu et al.(ASIPP), *An interpretable disruption predictor on EAST using XGBoost and SHAP*, Nuclear Fusion 65 (2025) 086035