跳到正文
返回核聚变能源
至真探索 · 核聚变能源

核聚变湍流输运与陀螺动理学：从GENE多通道验证到E级GPU与W7-X带状流

2026-08-30

plain-核聚变湍流输运与陀螺动理学 2026-08-30 Sun 16:22 时间戳: 20260830T162228

---

核聚变研究现在最缺的不是实验数据，是预测能力。

JET托卡马克创了聚变能量产出纪录，W7-X仿星器拿了仿星器三乘积（密度×温度×约束时间）的世界第一，NIF首次点燃了聚变。但一个聚变电站到底该怎么设计，哪些参数组合能让等离子体待得最久、损失最少，这些问题至今没有可靠的"先算后做"的答案。Lawson判据说三乘积得过 3×10²¹ m⁻³·keV·s 这道坎，现在密度和温度都能摸到，唯一卡住的是约束时间 τE——而 τE 几乎完全由等离子体内部那些看不见的小漩涡决定。这些漩涡叫微湍流，尺度比离子回旋半径还小，却能把核心的热量像漏斗一样往边缘搬。

所以问题归结成一件具体的事：能不能用计算机把这些小漩涡算准？

2025年3月，一篇发在 Nature Communications 上的论文给出了目前为止最肯定的回答。Klara Höfler 和她在马克斯·普朗克等离子体物理研究所（IPP）加兴（Garching）分部的同事们，把一个叫 GENE 的模拟代码和 ASDEX Upgrade 托卡马克上的实验数据做了一次"多通道"对比——同时比对涨落幅度、波数谱、密度和温度之间的相位关系、径向关联长度——这些东西以前都是一个个单独验，这次一起上。结果：模拟定量复现了实验观测。这是头一回。

这篇论文的标题直白："Milestone in predicting core plasma turbulence: successful multi-channel validation of the gyrokinetic code GENE"（核心等离子体湍流预测的里程碑：陀螺动理学代码GENE的成功多通道验证）。Nature Communications 第16卷，文章号2558，DOI: 10.1038/s41467-025-56997-2。通讯作者 Klara Höfler 同时挂了 IPP 加兴和慕尼黑工业大学（TUM）自然科学院物理系的牌子。

要把这件事说清楚，得先回到一个更根本的问题：为什么算这些小漩涡这么难，为什么要专门发明一套叫"陀螺动理学"的理论来对付它们。

一、陀螺动理学：在MHD和粒子轨道之间找一条路

等离子体里的事，能用的描述工具大致分三档。

最粗的一档是磁流体力学（MHD），把等离子体当成一种导电流体。它能看到大尺度的平衡和不稳定性，但看不见那些回旋半径尺度上的精细结构——偏偏约束时间恰恰由这些小结构决定。MHD 像用广角镜看风景，能看到山脉走向，看不到石缝里的流水。

最细的一档是直接解 Vlasov-Maxwell 方程，跟踪每一个粒子在电磁场里的运动。理论上最准，实际上算不动。一团聚变等离子体里有多少粒子？10²⁰ 量级每立方米。任何超级计算机都追不动这个数。

中间这一档，就是陀螺动理学（gyrokinetics）。它的核心想法是：在一个强磁场里，带电粒子一边沿磁力线走，一边绕着磁力线做快速回旋运动。回旋频率很高，轨道半径（回旋半径 ρ）很小。如果我们要看的物理过程——比如那些驱动湍流的不稳定性——频率远低于回旋频率、尺度远大于回旋半径，那就可以把快速回旋那个自由度平均掉，只保留沿磁力线的运动和引导中心的漂移。粒子被等效成一个个绕磁力线转的"带电环"，回旋角度不再出现，五维相空间（三个空间坐标加平行速度和磁矩）代替了原来的六维。

这个降维操作让计算量降了好几个数量级，同时保留了有限回旋半径效应、朗道阻尼、磁捕获这些MHD里没有的动力学。这是妥协，但是个聪明的妥协。

"陀螺动理学"这个词最早在1960年代末出现，Rutherford、Frieman、Taylor、Hastie 一批人在搞线性漂移波的动理学描述。真正的突破是 1982 年，普林斯顿的 Edward A. Frieman 和 Liu Chen（陈骝）在 *Physics of Fluids* 25卷502页给出了第一套非线性陀螺动理学方程组。在此之前，陀螺动理学只能算线性问题；Frieman-Chen 证明非线性项也可以做回旋平均，保留任意 k⊥ρ 的有限回旋半径效应，而且不需要人为截断。这件事的意义在于：湍流本质上是非线性的，没有非线性方程组，就没有后面的所有模拟。

Frieman-Chen 那套方程用了传统的多重尺度展开，管用但推导很繁，而且没有显式保证能量守恒和相空间体积守恒。1980年代中后期，Dubin 等人（1983，静电）和 Hahm 等人（1988，电磁）改用哈密顿/拉格朗日体系的 Lie 变换微扰方法重新推导，自然保住了守恒律，也更容易推到高阶。2007年，Alain Brizard（当时在 JAERI 东海，后到 Saint Michael's College, Vermont）和 T.S. Hahm（普林斯顿等离子体物理实验室 PPPL）在 *Reviews of Modern Physics* 79卷421–468页写了一篇综述："Foundations of nonlinear gyrokinetic theory"，把整个体系的来龙去脉和三大支柱——回旋中心哈密顿 Vlasov 方程、回旋动理学 Maxwell（Poisson-Ampère）方程、精确能量守恒律——讲透了。这篇综述至今是陀螺动理学领域的标准参考。

有了方程组，下一步就是写代码去解它。

二、GENE的谱系：从一个人的ETG代码到一个国际项目

GENE 的全名是 "Gyrokinetic Electromagnetic Numerical Experiment"，名字里就藏着野心——不是辅助工具，是"数值实验"。

1999年，Frank Jenko 在 IPP 加兴做博士后期间写下了 GENE 的第一版。当时的目标很具体：算电子温度梯度（ETG）驱动的湍流。2000年，Jenko 和 Bill Dorland（马里兰大学）、Manfred Kotschenreuther、B.N. Rogers 在 *Physics of Plasmas* 7卷1904页发表了 ETG 湍流的论文，在同期 Dorland 等人在 *Physical Review Letters* 85卷5579页发了姊妹篇。这两篇把 ETG 模式从一个理论猜想变成了可以用模拟验证的物理对象。GENE 就是从这批工作里长出来的。

Jenko 后来成了 IPP 加兴"托卡马克理论"部门的负责人，GENE 也从一个 ETG 专用脚本演变成一个覆盖多种不稳定性（ITG、TEM、ETG、KBM、微撕裂模）的通用平台。二十多年来，一支以 IPP 为核心的国际团队持续在开发，主要节点包括加兴的 IPP（Tokamak Theory 部门）、加兴的马克斯·普朗克计算与数据设施（MPCDF，前身是 RZG）、慕尼黑工业大学（TUM）、美国得克萨斯大学奥斯汀分校（UT Austin）、瑞士洛桑联邦理工（EPFL）的瑞士等离子体中心（SPC）。

GENE 的版本演化有几条主线。

最早是"局部"（flux-tube）版：算等离子体里一小条径向窄带，边界周期性处理，计算量相对小，适合做核心湍流的基础物理。2011年，Tobias Görler 带头把它扩展成"全局"版——径向上不再假设周期性，用有限差分代替谱方法，可以处理整个等离子体半径的剖面演化。论文发在 *Journal of Computational Physics* 230卷7053–7071页，作者里 Görler、Jenko、Told 来自 IPP 加兴的 Tokamak Theory 部门，Lapillonne、Brunner 来自 EPFL CRPP，Dannert 来自 MPCDF（当时叫 RZG），Merz 来自明斯特大学的研究组。全局版让 GENE 第一次能算和实验直接对应的剖面输运。

然后是仿星器版。托卡马克是轴对称的，仿星器不是——W7-X 的等离子体截面沿环向每个位置都在变形。把 GENE 搬到仿星器上，几何系数全得重做，这花了将近五年。2020年，Maurice Maurer（TUM 博士，导师 Jenko）作为第一作者在 *Journal of Computational Physics* 420卷109694页发表了 GENE-3D："A global gyrokinetic turbulence code for stellarators"。作者阵容庞大：Maurer、Bañón Navarro、Görler、Told、Jarema、Jenko 来自 IPP 加兴 Tokamak Theory，Dannert 来自 MPCDF，Restelli、Hindenlang 来自 IPP 加兴的"等离子体物理数值方法"（NMPP）部门，Merlo 当时挂在外部机构。GENE-3D 是第一个能在仿星器全三维几何下同时算离子和电子动力学（含磁场涨落）的陀螺动理学代码。

再往后是边缘版。等离子体从核心到壁有个过渡区，叫边缘和刮削层（edge and SOL），这里有 X 点、有开放磁面、有密度温度低到涨落幅度跟背景差不多大的区域，传统的沿磁面坐标在这里会出奇点。2021年，Dominik Michels（TUM 计算机学院博士，导师 Jenko 和 Sonnendrücker）在 *Computer Physics Communications* 264卷107986页发表了 GENE-X："A full-f gyrokinetic turbulence code based on the flux-coordinate independent approach"。合作者 Stegmeir、Ulbl、Jarema 来自 IPP 加兴 Tokamak Theory，Jenko 同时挂 IPP 和外部机构。GENE-X 用"磁通量坐标无关"（FCI）方法，在非结构化局部笛卡尔网格上算，能从磁轴一直算过分离面进到刮削层。它还有一个关键差别：用 full-f 而不是 δ-f。δ-f 假设涨落远小于背景，full-f 不做这个假设——边缘区涨落大到和背景同量级，必须用 full-f。

到这里，GENE 家族构成了一个从核心到边缘、从托卡马克到仿星器的完整工具链。但工具有了，问题没完：它能信吗？

三、多通道验证：把模拟和实验放到同一张表上

模拟代码最怕一件事：自洽。自己跑出来漂亮的结果，换个参数就塌了；或者跟某个实验量对上了，跟另一个量又对不上。做验证（validation）的核心，就是看模拟在多少个独立的实验可观测量上能同时站住脚。

Greenwald（2010, Phys. Plasmas 17, 058101）和 Terry 等人（2008）给聚变验证立过规矩：不能只看全局量（比如时空平均热流），还得看涨落本身的特征——幅度、波数谱、交叉相位。Höfler 等人这次做的，正是把这些规矩一次性兑现。

实验在 ASDEX Upgrade 上做。ASDEX Upgrade 是 IPP 加兴从1991年开始运行的托卡马克，德国最大的聚变装置，R=1.65m，a=0.5m，全钨第一壁——这是它给 ITER 和未来电站探路的一个关键贡献，证明了钨壁在反应堆条件下可行。这次实验用的是氘低约束模式（L-mode）放电，磁场 −2.537 T，等离子体电流 0.87 MA，辅助加热包括 0.87 MW 中性束和 1.07 MW 电子回旋共振加热。选 L-mode 是刻意的：在 AUG 的 H-mode 下，归一化回旋半径 ρ* 的值会大到违反局部陀螺动理学的假设。用了上单零磁构型，避免跳到 H-mode。

设计了两个场景，区别在电子温度梯度——一个陡（steep），一个平（flat）——这样能改变湍流驱动条件。两套场景在两个径向位置都测了，结果一致，论文只展示一个。

关键在"多通道"到底多在哪里。他们在同一个放电上，用两种独立诊断同时测了电子温度涨落和密度涨落的全套特征：

CECE（关联电子回旋发射）测温度涨落谱。DBS（多普勒背散射）测密度涨落的波数谱。然后还要测密度和温度之间的交叉相位 α_nT——涨落不是各自随便动，它们之间的相位关系藏着湍流的物理。再测径向关联长度，看涨落结构在径向延伸多远。

模拟这边，GENE 拿实验测到的平衡和动理学剖面做输入，跑非线性湍流模拟，输出涨落场。但涨落场不能直接跟诊断信号比——诊断看到的是微波和等离子体相互作用后的信号，中间隔着一层物理。所以要用"合成诊断"（synthetic diagnostics）：对模拟输出的涨落场跑一遍和真实诊断相同的测量物理（射线追踪、二维全波分析），得到"如果诊断去测这个模拟，它会看到什么"。这样模拟和实验才在同一个观察面上比。

结果是什么？电子和离子的面积分热流跟输运分析（带马尔可夫链蒙特卡洛不确定度）对上。温度涨落谱对上。密度涨落波数谱对上。交叉相位对上。径向关联长度对上。论文里逐个画图，实验是实线带误差棒，模拟合成诊断是虚线或叉号，叠在一起。

Höfler 和合作者特别强调，这不是"某一个量对上了"的验证，是"迄今最大数量湍流量同时对比"。他们把过去的若干验证工作列了张表（论文 Fig.1a），每一项用了哪些可观测量，这次的工作把这些通道全占了。结论很直接：陀螺动理学代码已经成熟到可以作为预测核心湍流的可靠工具。

作者名单本身就是个国际网络。IPP 加兴出主力（Höfler 兼通讯作者、Görler、Happel、Molina、Bergmann、Conway、David、Fischer、Jenko、McDermott、Stroth），TUM 自然科学院物理系挂了 Höfler 和 Stroth，斯图加特大学界面过程工程与等离子体技术研究所出了 Lechte（做合成诊断的全波分析），MIT 等离子体科学与聚变中心出了 Bielajew、Denk、White，EPFL 瑞士等离子体中心挂了 Molina，通用原子（San Diego）挂了 Denk，巴黎综合理工的等离子体物理实验室出了 Hennequin。这是一个典型的"实验+理论+模拟+合成诊断"闭环团队，每一环都得有人。

这个里程碑不是凭空蹦出来的。Görler、White、Jenko 等人之前就在 DIII-D 和 AUG 上做过大量对比，发现过"离子热流低估近一个量级"的所谓 shortfall 问题，后来靠加全物理（碰撞、电磁、真实几何）逐步缩小。Höfler 这篇是这个长期努力的一次集中兑现。

四、W7-X上的带状流故事：当饱和机制换了人

托卡马克上验证通了，仿星器呢？

Wendelstein 7-X（W7-X）是 IPP 在格赖夫斯瓦尔德（Greifswald）的超导仿星器，2015年底首次产生等离子体。它的卖点是用三维磁场的"数值优化"把新经典输运压下去——Beidler 等人 2021 年在 *Nature* 596卷221页证明了这个优化确实生效：在创纪录的放电里，新经典损失占加热功率的 30%，而如果换一个没优化的磁场，新经典损失会超过加热功率，物理上不可能。换句话说，W7-X 观测到的那些剖面只有在低新经典损失的磁场里才可能出现。

新经典压下去了，剩下的就是湍流。W7-X 头几次实验就发现，总输运比新经典预测高，湍流是主犯。GENE-3D 正好派上用场。

这里有个有意思的物理。

在托卡马克里，离子温度梯度（ITG）湍流饱和的主要机制是"带状流"（zonal flow）。带状流是一种 m=0、n=0 的径向电场剪切流——湍流自己产生的副产品，反过来剪切湍流涡旋，把涡旋切碎，压低输运。这像湍流自己给自己踩刹车。带状流强，输运就被压低；带状流弱，湍流就肆虐。2000年 Dimits 等人在 *Physics of Plasmas* 7卷969页做的著名 Cyclone Base Case 比较，揭示了一个叫"Dimits shift"的现象：在临界梯度之上还有一段非线性临界梯度，在线性不稳定但非线性输运几乎为零——正是带状流撑出来的。这段低输运的"Dimits 态"由准静态的带状流楼梯（E×B staircase）维持，涡旋在强剪切区被剪，转折点处温度梯度被抹平。Ivanov（牛津，2022 EPS 会议）和 Li、Pueschel、Whelan、Terry（Nuclear Fusion 63, 2023）等人后来把 Dimits shift 的机制推到了模共振层面。

但在 W7-X 上，带状流弱了。

Proll、Xanthopoulos、Helander、Plunk、Faber、Görler、Smith、Pueschel 这群人（2020 arXiv:2002.10484；正式版 *Journal of Plasma Physics* 88卷905880112, 2022）用 GENE 算 W7-X 的 TEM（捕获电子模）湍流，发现一个反直觉的事：W7-X 的 TEM 热流比准对称仿星器 HSX 和托卡马克 DIII-D 都低——这是优化的功劳，因为 W7-X 近似满足"最大 J"性质，线性上 TEM 在很大参数空间稳定。但饱和机制不一样：托卡马克里 TEM 主要靠带状流饱和，W7-X 里带状流弱得多，得靠别的机制。

作者机构需要逐个核实。Proll 当时在埃因霍温理工大学核聚变科学技术组，Xanthopoulos、Helander、Plunk、Smith 在 IPP 格赖夫斯瓦尔德分部，Faber 在威斯康星大学麦迪逊分校的 HSX 等离子体实验室（同时挂物理系），Görler 在 IPP 加兴，Pueschel 在 UT Austin 聚变研究所。这是一张典型的"仿星器国际联合"表。

带状流为什么在 W7-X 弱？和三维几何有关。在托卡马克轴对称场里，带状流可以沿磁面均匀地形成；W7-X 的三维性让不同磁面上的曲率和捕获粒子分布不一样，带状流的"组织"能力被打散了。结果是：W7-X 不靠带状流踩刹车，而是靠几何本身（最大 J）在线性阶段就不让 TEM 起来。两条路，殊途都同归于低输运。

González-Jerez 等人（2022, *J. Plasma Physics* 88卷905880310）后来做了一个漂亮的基准测试：把新写的仿星器代码 stella 和 GENE 在 W7-X 几何里跑同一套 ITG/TEM 线性稳定性、线性带状流响应、非线性 ITG 热流，逐个对。作者是 CIEMAT 马德里国家聚变实验室的 González-Jerez、García-Regaña、Calvo，IPP 格赖夫斯瓦尔德的 Xanthopoulos、Alcusón、Geiger，IPP 加兴的 Bañón-Navarro，牛津大学 Rudolf Peierls 理论物理中心的 Barnes 和 Parra。仿星器陀螺动理学代码以前太少，没人敢信任何一家的结果，有了 stella 和 GENE 的互校，才有了可以往下做的基础。

更晚近的物理在往两个方向走。

一个是快粒子路线。Di Siena、Bañón Navarro、Jenko（2020, *Physical Review Letters*, arXiv:2003.00288）提出：在优化仿星器里，用离子回旋共振加热（ICRF）产生的超热离子，可以和 ITG 模发生共振波-粒相互作用，把湍流抑制掉最多 80%。Di Siena 当时同时挂 UT Austin 和 IPP 加兴，Bañón Navarro 和 Jenko 在 IPP 加兴。这是把托卡马克上发现的快离子抑制湍流机制"搬"到仿星器上的第一次理论证明，而且预言了 W7-X 上做这个实验的条件。后来 Bañón Navarro、Di Siena、Velasco（CIEMAT）、Wilms、Merlo（UT Austin）、Windisch（IPP 格赖夫斯瓦尔德）、LoDestro 和 Parker（LLNL）、Jenko 在 2023 年 *Nuclear Fusion* 上更进一步：用 GENE-3D 耦合新经典输运代码 KNOSOS 和一维输运求解器 TANGO，第一次从第一性原理出发自洽地预测 W7-X 的温度剖面，解释了 ECRH 放电里离子温度被"钳"在 ~1.5 keV 上不去的机制——ECRH 造成大的电子-离子温度比，放大 ITG 湍流，反过来限制离子温度。

另一个是微撕裂模（MTM）。Cu-Castillo、Bañón Navarro、Merlo、Reimold、……、Jenko 和 W7-X 团队（2026, *Physical Review Letters* 137卷035102，2026年7月14日发表）在 W7-X 一个湍流被压低、密度梯度陡、温度梯度中等的放电里，用 GENE 算出微撕裂模主导输运。W7-X 的准等向性和近最大 J 配置稳住了 ITG 和密度梯度驱动的 TEM，中等碰撞性和低磁剪切又给 MTM 让出了生长空间。他们在密度梯度扫描里看到，实验观测的阈值处输运大幅下降——正是从 ITG 主导到 MTM 主导的转变。作者里 Cu-Castillo、Bañón Navarro、Merlo、Jenko 在 IPP 加兴，Reimold 等在 IPP 格赖夫斯瓦尔德。

还有多尺度。Merlo、Bañón Navarro、Görler、Jenko（IPP 加兴）和 Wilms（Proxima Fusion，慕尼黑）在 2025 年用 GPU 版 GENE 做了 W7-X 的第一次多尺度湍流模拟——同时算离子尺度和电子尺度。在 ITG 主导的电静电区，电子尺度涨落会侵蚀带状流，放大离子尺度输运；反过来离子尺度湍流压低 ETG。而在电磁 MTM 条件下，ETG 的各向同性限制了它的抑制效应，MTM 得以存活。跨尺度耦合对 W7-X 的输运预测至关重要。

W7-X 上的故事整体是：几何优化把新经典压下去了，湍流接管了限制角色；带状流不再是主角，几何本身的线性稳定性和新的非线性机制（快粒子共振、MTM 转变、多尺度耦合）来决定输运水平。GENE 在这里的角色从"解释"在往"预测"走——2023 年的剖面预测论文已经做到了自洽推演。

五、E级计算与GPU：从解释到预测的那道坎

"解释"和"预测"听起来像一步之遥，计算上隔着一整个数量级。

解释实验，是把实验测到的剖面喂给代码，看代码算出来的涨落特征和实验对不对得上。预测，是不知道结果，先算一遍，告诉实验会看到什么。后者对计算量的要求高得多：要算长时间演化、要扫参数空间、要做整器建模（whole-device modeling），核心到壁一把抓。

Jenko 的话说得很直白：有了 E 级 GENE，才第一次从"解释实验"走到"预测实验"，目标是给 ITER 造一个"数字孪生"——一个虚拟的 ITER 等离子体，先算后做，省掉大量试错。

2022年底，欧盟委员会在"面向 E 级时代的卓越中心"招标下批了一个项目：Plasma-PEPSC（Pushing Flagship Plasma Simulation Codes to Tackle Exascale-Enabled Grand Challenges via Performance Optimisation and Codesign），由瑞典斯德哥尔摩 KTH 领导，IPP、MPCDF、TUM 和另外八个欧洲伙伴参与，总经费超过 700 万欧元，其中给 GENE 代码升级的是 214 万欧元，德国联邦教育研究部出一半。2023年1月1日启动，为期四年。

E 级计算机不是"把老代码拿到新机器上跑"。E 级机（10¹⁸ 次运算每秒）的性能跃升主要靠新硬件架构——大量 GPU 加 CPU——代码不重新适配，根本吃不到算力。Jenko 用了个比喻：十个文件夹十个人能比一个人快十倍，但突然给你一万人，只有完全重新分工才可能更快。这叫协同设计（co-design），硬件团队和代码团队一起设计。

GPU 移植这条线，GENE 和 GENE-X 各有进展。

GENE 的 GPU 版由 Kai Germaschewski（新罕布什尔大学，当时）、Tilman Dannert（MPCDF）、Gabriele Merlo（UT Austin Oden 研究所）、Frank Jenko（IPP 加兴）、Stephane Ethier（PPPL）、Ed D'Azevedo（ORNL）、Amitava Bhattacharjee（PPPL）等人在 2021 年完成。论文发在 *Physics of Plasmas* 28卷062501页，DOI: 10.1063/5.0046327。他们自己写了一个叫 gtensor 的库来支持性能可移植的 GPU 实现。在美国 Summit 超级计算机的单节点上，GPU 版相对纯 CPU 版快了将近 15 倍，典型 GPU 核达到峰值带宽的约 90%。瓶颈在多 GPU 通信——上千块 GPU 时通信成本开始吃掉收益。这个项目挂在美国 DOE 的 ECP（Exascale Computing Project）下面，叫 WDMapp（whole-device modeling application）。

GENE-X 的 GPU 版走的是另一条技术路线。Jordy Trilaksono、Philipp Ulbl（IPP 加兴）、Jeremy J. Williams（KTH 电气工程与计算机科学学院，计算科学与技术组）、Carl-Martin Pfeiler、Martin Finkbeiner、Tilman Dannert……Jenko 在 2025 年 PASC 会议（Platform for Advanced Scientific Computing Conference，ACM 出版）上发了论文，用 OpenACC 和 OpenMP offload 把 GENE-X 的 Fortran 算子封装到 C++17 层、再卸载到 GPU。做了单 GPU 屋顶线分析和多 GPU 扩展分析，主要算子性能显著提升，瓶颈同样转移到 GPU 间通信。

还有一个 GENE-X 的速度空间谱方法新进展。2025 年 PASC 一个小型研讨会上报告：GENE-X 里新实现了一种速度空间谱方法，和原来的欧拉（网格）版 GENE-X 结果高度一致（已用实验测量验证过），但速度空间分辨率可以降低 50 倍，带来显著加速。这是把边缘湍流模拟从"百万 CPU 小时、跑几周"往"能用于反应堆场景"推的一步。

把这几条线放一起看：GENE 有 GPU 版做核心湍流，GENE-3D 有仿星器三维几何，GENE-X 有边缘和刮削层的 full-f FCI 能力，加上 E 级项目和 GPU 移植，整个工具链在往"整器建模"——从磁轴到壁一次算完——靠拢。Höfler 那篇多通道验证，是这条链上核心段的一次硬验证；W7-X 上的一系列工作，是仿星器段的验证和发现；GENE-X 的边缘工作，是最后一块拼图。拼齐了，才谈得上给 ITER 或者以后的电站做数字孪生。

六、还在路上的事

这篇综述看到这里，有几件事值得点出来。

第一，验证不是一次完成的。Höfler 那篇是 ASDEX Upgrade 上的 L-mode 两个径向位置，特定磁场和电流。H-mode、更高参数、不同装置（JET、DIII-D、ITER）上的多通道验证还需要继续做。论文结尾特意说会把数据公开，供进一步的验证。科学的诚实就在这里：一次成功不是终点，是把门槛抬高，让下一次验证更严。

第二，仿星器上的物理和托卡马克系统性地不同。W7-X 上带状流弱、几何本身做线性稳定、靠快粒子共振和 MTM 转变来定输运——这些发现意味着不能把托卡马克上的经验直接搬过去。GENE-3D 在这里从解释往预测走，2023 年剖面预测是第一次从第一性原理自洽推演 W7-X 温度剖面。

第三，计算物理和实验物理之间的"翻译"——合成诊断——已经成了和代码本身同等重要的一环。没有 Lechte 做的全波分析、没有 CECE 和 DBS 的测量体积精确标定（用 Torbeam 和 ECRad），多通道对比无从谈起。诊断物理的精度是验证可信度的硬上限。

第四，E 级不是万能药。GPU 移植后瓶颈转移到通信、内存带宽，谱方法这类算法革新和代码重构同样重要。Plasma-PEPSC 项目四年期结束后能不能真的造出 ITER 的数字孪生，取决于硬件、算法、物理三头并进。

第五，作者机构的国际网络本身就是这个故事的一部分。IPP 加兴和格赖夫斯瓦尔德是轴心，但 MIT、EPFL、UT Austin、Oxford、CIEMAT、PPPL、ORNL、KTH、LLNL、UNH、Proxima Fusion……每一项关键进展都是一个跨大西洋、跨欧亚的联合体做出来的。GENE 代码本身是开源的，使用它的机构覆盖德、美、瑞士、法、英、荷、意、日、印、中、韩。聚变模拟不是任何一家的私活。

---

引用文献

  1. K. Höfler, T. Görler, T. Happel, C. Lechte, P. Molina, M. Bergmann, R. Bielajew, G.D. Conway, P. David, S.S. Denk, R. Fischer, P. Hennequin, F. Jenko, R.M. McDermott, A.E. White, U. Stroth, and the ASDEX Upgrade Team. "Milestone in predicting core plasma turbulence: successful multi-channel validation of the gyrokinetic code GENE." *Nature Communications* 16, 2558 (2025). DOI: 10.1038/s41467-025-56997-2
  2. 通讯作者 Klara Höfler：Max Planck Institute for Plasma Physics, Boltzmannstr. 2, Garching, Germany；Technical University of Munich, TUM School of Natural Sciences, Physics Department, James-Franck-Str. 1, Garching, Germany
  3. Carsten Lechte：University of Stuttgart, Institute of Interfacial Process Engineering and Plasma Technology, Pfaffenwaldring 31, Stuttgart, Germany
  4. Rachel Bielajew, Severin S. Denk, Anne E. White：Plasma Science and Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
  5. Denk 兼：General Atomics, San Diego, California, USA
  6. Pedro Molina 兼：Ecole Polytechnique Fédérale de Lausanne (EPFL), Swiss Plasma Center (SPC), Lausanne, Switzerland
  7. Pascale Hennequin：Laboratoire de Physique des Plasmas, Ecole Polytechnique, Palaiseau, France
  8. 其余 Görler, Happel, Molina, Bergmann, Conway, David, Fischer, Jenko, McDermott, Stroth：Max Planck Institute for Plasma Physics, Boltzmannstr. 2, Garching, Germany
  1. E.A. Frieman and L. Chen. "Nonlinear gyrokinetic equations for low-frequency electromagnetic waves in general plasma equilibria." *Physics of Fluids* 25, 502 (1982).
  1. A.J. Brizard and T.S. Hahm. "Foundations of nonlinear gyrokinetic theory." *Reviews of Modern Physics* 79, 421–468 (2007). DOI: 10.1103/RevModPhys.79.421
  2. Brizard：Department of Chemistry and Physics, Saint Michael's College, Colchester, Vermont 05439, USA（当时兼 JAERI, Tokai, Japan）
  3. Hahm：Plasma Physics Laboratory, Princeton University, Princeton, New Jersey 08543, USA
  1. F. Jenko, W. Dorland, M. Kotschenreuther, and B.N. Rogers. "Electron temperature gradient driven turbulence." *Physics of Plasmas* 7, 1904–1910 (2000). DOI: 10.1063/1.874014
  2. GENE 代码的首篇应用论文
  1. W. Dorland, F. Jenko, M. Kotschenreuther, and B.N. Rogers. "Electron temperature gradient turbulence." *Physical Review Letters* 85, 5579 (2000).
  1. T. Görler, X. Lapillonne, S. Brunner, T. Dannert, F. Jenko, F. Merz, and D. Told. "The global version of the gyrokinetic turbulence code GENE." *Journal of Computational Physics* 230(18), 7053–7071 (2011). DOI: 10.1016/j.jcp.2011.05.034
  2. Görler, Jenko, Told：Max Planck Institute for Plasma Physics, Tokamak Theory (TOK), Garching
  3. Lapillonne, Brunner：Ecole Polytechnique Fédérale de Lausanne, CRPP, Lausanne
  4. Dannert：Max Planck Computing and Data Facility (MPCDF，前身 RZG), Garching
  5. Merz：Helmholtz University Research Group 'Theory and Ab Initio Simulation of Plasma Turbulence'（明斯特青年研究组）, IPP
  1. M. Maurer, A. Bañón Navarro, T. Dannert, M. Restelli, F. Hindenlang, T. Görler, D. Told, D. Jarema, G. Merlo, and F. Jenko. "GENE-3D: A global gyrokinetic turbulence code for stellarators." *Journal of Computational Physics* 420, 109694 (2020). DOI: 10.1016/j.jcp.2020.109694
  2. Maurer 博士论文：Technische Universität München, 2020
  3. Maurer, Bañón Navarro, Görler, Told, Jarema, Jenko：IPP Tokamak Theory
  4. Dannert：MPCDF
  5. Restelli, Hindenlang：IPP Numerical Methods in Plasma Physics (NMPP)
  1. D. Michels, A. Stegmeir, P. Ulbl, D. Jarema, and F. Jenko. "GENE-X: A full-f gyrokinetic turbulence code based on the flux-coordinate independent approach." *Computer Physics Communications* 264, 107986 (2021). DOI: 10.1016/j.cpc.2021.107986
  2. Michels 博士论文：Technische Universität München, Fakultät für Informatik, 2021（导师 F. Jenko, E. Sonnendrücker）
  3. 全员：IPP Tokamak Theory (TOK), Garching
  1. K. Germaschewski, B. Allen, T. Dannert, M. Hrywniak, J. Donaghy, G. Merlo, F. Jenko, S. Ethier, A. Bhattacharjee, and E. D'Azevedo. "Toward exascale whole-device modeling of fusion devices: Porting the GENE gyrokinetic microturbulence code to GPU." *Physics of Plasmas* 28, 062501 (2021). DOI: 10.1063/5.0046327
  2. Germaschewski：University of New Hampshire, Durham
  3. Dannert：Max Planck Computing and Data Facility (MPCDF), Garching
  4. Merlo, Jenko：The University of Texas at Austin Oden Institute / IPP Garching
  5. Ethier, Bhattacharjee：Princeton Plasma Physics Laboratory (PPPL)
  6. D'Azevedo：Oak Ridge National Laboratory (ORNL)
  1. J. Trilaksono, P. Ulbl, J.J. Williams, C.-M. Pfeiler, M. Finkbeiner, T. Dannert, …, F. Jenko. "OpenACC and OpenMP-Accelerated Fortran/C++ Gyrokinetic Fusion Code GENE-X for Heterogeneous Architectures." Proceedings of the Platform for Advanced Scientific Computing Conference (PASC 2025), ACM, 2025.
  2. Trilaksono, Ulbl, Pfeiler, Finkbeiner, Dannert, Jenko：Max Planck Institute for Plasma Physics, Garching
  3. Williams：KTH, School of Electrical Engineering and Computer Science (EECS), Computational Science and Technology (CST)
  1. EUROfusion 项目 "Plasma-PEPSC"（Pushing Flagship Plasma Simulation Codes to Tackle Exascale-Enabled Grand Challenges via Performance Optimisation and Codesign）。KTH Stockholm 领导，IPP、MPCDF、TUM 及八个欧洲伙伴参与，2023年1月启动，四年期，总经费超 700 万欧元，GENE 部分 214 万欧元。
  1. A.M. Dimits et al. "Comparisons and physics basis of tokamak transport models and turbulence simulations." *Physics of Plasmas* 7, 969 (2000). DOI: 10.1063/1.874036
  2. Cyclone Base Case 比较，"Dimits shift" 的来源
  1. J.H.E. Proll, P. Xanthopoulos, P. Helander, G.G. Plunk, B.J. Faber, T. Görler, H.M. Smith, M.J. Pueschel. "Turbulence mitigation in maximum-J stellarators with electron-density gradient." *Journal of Plasma Physics* 88(1), 905880112 (2022). arXiv:2002.10484
  2. Proll：Science and Technology of Nuclear Fusion, Department of Applied Physics, Eindhoven University of Technology
  3. Xanthopoulos, Helander, Plunk, Smith：Max-Planck-Institut für Plasmaphysik, Greifswald
  4. Faber：HSX Plasma Lab, University of Wisconsin-Madison
  5. Görler：IPP Garching
  6. Pueschel：Institute for Fusion Studies, The University of Texas at Austin
  1. A. González-Jerez, P. Xanthopoulos, J.M. García-Regaña, I. Calvo, J. Alcusón, A. Bañón-Navarro, M. Barnes, F.I. Parra, J. Geiger. "Electrostatic gyrokinetic simulations in Wendelstein 7-X geometry: benchmark between the codes stella and GENE." *Journal of Plasma Physics* 88, 905880310 (2022). DOI: 10.1017/S0022377822000393
  2. González-Jerez, García-Regaña, Calvo：Laboratorio Nacional de Fusión, CIEMAT, Madrid, Spain
  3. Xanthopoulos, Alcusón, Geiger：Max-Planck-Institut für Plasmaphysik, Greifswald
  4. Bañón-Navarro：Max-Planck-Institut für Plasmaphysik, Garching
  5. Barnes, Parra：Rudolf Peierls Centre for Theoretical Physics, University of Oxford
  1. A. Di Siena, A. Bañón Navarro, F. Jenko. "Turbulence suppression by energetic particle effects in modern optimized stellarators." *Physical Review Letters* (2020). arXiv:2003.00288
  2. Di Siena：The University of Texas at Austin + Max Planck Institute for Plasma Physics, Garching
  3. Bañón Navarro, Jenko：IPP Garching
  1. A. Bañón Navarro, A. Di Siena, J.L. Velasco, F. Wilms, G. Merlo, T. Windisch, L.L. LoDestro, J.B. Parker, F. Jenko. "First-principles based plasma profile predictions for optimized stellarators." *Nuclear Fusion* (2023).
  2. Bañón Navarro, Di Siena, Wilms, Jenko：IPP Garching
  3. Velasco：Laboratorio Nacional de Fusión, CIEMAT, Madrid
  4. Merlo：The University of Texas at Austin
  5. Windisch：IPP Greifswald
  6. LoDestro, Parker：Lawrence Livermore National Laboratory (LLNL)
  1. H. Cu-Castillo, A. Bañón Navarro, G. Merlo, F. Reimold, T. Romba, O. Ford, …, J. Geiger, A. Goodman, A. Zocco, F. Jenko, and the W7-X Team. "Significance of Microtearing Turbulence in Turbulence-Reduced High-Density-Gradient Plasmas in Wendelstein 7-X." *Physical Review Letters* 137, 035102 (2026). DOI: 10.1103/nd3p-yb9t
  2. Cu-Castillo, Bañón Navarro, Merlo, Jenko：IPP Garching
  3. Reimold, Romba, Ford, Geiger, Goodman, Zocco 等：IPP Greifswald
  1. C.D. Beidler et al. (W7-X Team). "Wendelstein 7-X 的新经典输运优化验证." *Nature* 596, 221 (2021). DOI: 10.1038/s41586-021-03687-w
  2. Craig Beidler：IPP Greifswald, Stellarator Theory Division
  1. G. Merlo, A. Bañón Navarro, T. Görler, F. Jenko, F. Wilms. "Multiscale turbulence in stellarators." (2025)
  2. Merlo, Bañón Navarro, Görler, Jenko：Max Planck Institute for Plasma Physics, Garching
  3. Wilms：Proxima Fusion, Munich
  1. P.G. Ivanov. "Zonally dominated dynamics and the transition to strong turbulence in ion-scale plasma turbulence." EPS 2022 Conference, O5.104.
  2. Ivanov：Rudolf Peierls Centre for Theoretical Physics, University of Oxford + EURATOM/UKAEA Fusion Association, Culham Science Centre
  1. P.-Y. Li, M.J. Pueschel, P.W. Terry, G.G. Whelan. "On the role of mode resonances in regulating zonal-flow-moderated plasma microturbulence." *Nuclear Fusion* 63(2), 026039 (2023). DOI: 10.1088/1741-4326/acae37
  1. ASDEX Upgrade 装置：Max Planck Institute for Plasma Physics, Boltzmannstraße 2, Garching, Germany。1991年首次等离子体。R=1.65m, a=0.5m, 全钨第一壁。
  1. Wendelstein 7-X 装置：Max Planck Institute for Plasma Physics, Greifswald, Germany。2015年12月首次等离子体。50 块超导非平面磁线圈。