[论文解读] A model for the fast evaluation of prompt losses of energetic ions in stellarators

J. L. Velasco, I. Calvo+5 more|arXiv (Cornell University)|Jun 10, 2021

Magnetic confinement fusion research参考文献 36被引用 22

一句话总结

本文提出了一种基于平均回旋漂移-漂移方程的快速、径向局域化模型，用于分类托卡马克中高能离子轨道为约束或瞬时损失。通过计算临界漂移参数 Γα，该模型相较于全蒙特卡罗模拟，能以高精度预测与构型相关的瞬时损失，从而实现高效集成至托卡马克优化工作流程中，以最小化高能离子损失。

ABSTRACT

A good understanding of the confinement of energetic ions in non-axisymmetric magnetic fields is key for the design of reactors based on the stellarator concept. In this work, we develop a model that, based on the radially-local bounce-averaged drift-kinetic equation, classifies orbits and succeeds in predicting configuration-dependent aspects of the prompt losses of energetic ions in stellarators. Such a model could in turn be employed in the optimization stage of the design of new devices.

研究动机与目标

• 开发一种计算高效的模型，用于预测非轴对称托卡马克磁构型中高能离子的瞬时损失。
• 基于回旋平均径向漂移，将粒子轨道分类为约束或瞬时损失，以实现对磁构型性能的快速评估。
• 在一系列 W7-X 平衡态（不同 β、镜像项和旋转变换）下，通过与高保真全轨道蒙特卡罗模拟（ASCOT）对比，验证该模型的准确性。
• 通过确保速度与精度，使该模型能够集成至托卡马克设计优化流程中，实现对构型的实时筛选。
• 探索该模型在 W7-X 之外的普适性，包括在螺旋形和准对称托卡马克中的应用，并评估其在指导全局模拟设置方面的潜力。

提出的方法

• 该模型使用回旋平均漂移-漂移方程，计算磁漂移的径向分量，重点聚焦于回旋平均径向漂移速度。
• 引入一个临界参数 Γα = (2/π) arctan(|∂αJ| / |∂sJ|)，用于量化切向与径向漂移分量的比值，从而决定轨道的约束性。
• 若 Γα > 0.5，则将轨道分类为瞬时损失，表明由于非零的回旋平均径向磁漂移，存在显著的径向漂移。
• 该方法在 KNOSOS 代码中实现，该代码严格处理与磁面相切的磁漂移，并已集成至 STELLOPT 优化套件中。
• 模型预测结果在 71 个自由边界 VMEC 平衡态（KJM 托卡马克构型）上与 ASCOT 蒙特卡罗模拟进行了对比，其 β 值（0–7%）和抛物线形压强分布各不相同。
• 该模型进一步扩展以通过 vE 引入径向电场，从而在漂移计算中考虑未来 E×B 效应和速度空间中的碰撞扩散。

实验结果

研究问题

• RQ1 是否可以使用一种径向局域化、回旋平均化的模型，在无需全轨道追踪的情况下，准确预测托卡马克中高能离子的瞬时损失？
• RQ2 在类似 W7-X 的托卡马克中，临界参数 Γα 与不同磁构型和 β 值下的瞬时损失分数之间存在何种相关性？
• RQ3 该模型在多大程度上能捕捉瞬时损失对径向位置（s）和磁构型参数（镜像项、旋转变换、β）的依赖关系？
• RQ4 该模型是否能识别出瞬时损失对 β 的非单调或阈值型行为，例如在 2–3% 附近存在最小损失区域？
• RQ5 在真实构型中，引入径向电场和碰撞效应后，该模型对瞬时损失的预测将如何变化？

主要发现

• 该模型在一系列 W7-X 平衡态中对瞬时损失分数的预测具有高精度，与 ASCOT 蒙特卡罗模拟在 s = 0.06、0.25 和 0.50 处产生的高能离子结果高度一致。
• 瞬时损失随 β 增加而近似线性减少，在最内层磁面（s = 0.06）处观察到约 2–3% β 时的最小斜率，表明可能存在一个优化的“甜点”区域。
• 该模型正确捕捉了在 s = 0.06 处低 β 对瞬时损失的不利影响，其中反磁效应可能加剧约束损失，而这一特征在 s = 0.25 处未被观察到。
• 该模型在不同磁构型下（包括镜像项和旋转变换的变化）的预测保持稳健，显示出其对构型的敏感性与强大的预测能力。
• 该模型揭示了超香蕉轨道（其特征为大的回旋平均径向漂移）是瞬时损失的主要来源，且当 Γα > 0.5 时，其出现时间被准确预测。
• 该模型的成功表明，回旋平均漂移-漂移方程可定量描述高能离子的类经典输运，通过聚焦于初始束缚轨道分类，实现更快的全局模拟。