Particle pinch in global tokamak edge simulations

技术摘要：全球托卡马克边缘模拟中的粒子 pinch（粒子向内输运）现象

问题陈述
在磁约束聚变等离子体中，经常观察到“粒子 pinch”（即粒子相对于压力梯度的向内输运），这导致了中心密度剖面的峰值化。虽然这对于理解密度台阶（density pedestal）形成、密度极限以及 L-H 转变至关重要，但其潜在机制尚未得到完全理解，尤其是在边缘区域。现有的理论通常侧重于核心区，或在边缘模型中依赖于 ad hoc（特设）假设（例如固定核心密度或预设 pinch 速度）。此外，大多数边缘湍流模拟运行在短时标（亚毫秒级）的湍流时间尺度上，忽略了中性粒子动力学或长时标的剖面演化，因此无法捕捉在输运时间尺度上自洽发展的密度剖面。

方法论
作者采用了全局漂移气球（Global Drift Ballooning, GDB）代码，这是一种基于降低阶数的 Braginskii 方程的、通量驱动的、全局性的 3D 电磁湍流模型。模拟设置模仿了一个类似于内壁限制的 Alcator C-Mod 配置（单电荷氘等离子体，$B_0 = 5.3$ T）。

• 初始条件： 初始化为一个平坦的密度剖面（$n_0 = 0.8 \times 10^{20} \text{ m}^{-3}$）。
• 源项： 粒子注入源（$S_n$）通过高斯分布局部化在最后闭合磁面（LCFS）附近；热源（$S_{Te,i}$）则局部化在核心侧，以维持约 180 eV 的温度。
• 时间尺度： 模拟运行了约 10 ms，覆盖了输运时间尺度，允许密度剖面从平坦状态演化为准稳态。
• 分析： 研究利用了简化平板模型的线性稳定性分析、非线性控制模拟（关闭特定物理项）以及对径向粒子通量进行平衡态与涨落分量的详细分解。

关键结果与机制
模拟成功地在无需 ad hoc 假设的情况下重现了中心峰值的密度剖面，揭示了驱动向内粒子通量的两种截然不同的机制，它们在剖面演化的不同阶段占据主导地位：

1. 早期阶段（负 $\eta_\alpha$）：

   • 条件： 最初，密度梯度为正（凹陷剖面）而温度梯度为负，导致 $\eta_\alpha = L_n/L_{T\alpha} < 0$。
   • 机制： 漂移波湍流通过电子热扩散驱动向内的粒子通量。线性分析表明，当 $\eta_\alpha < 0$ 时，密度与势能涨落之间的相位差（$\delta_t$）变为正值，从而驱动向内的 $E \times B$ 通量。
   • 验证： 通过控制模拟关闭电子热动力学项（$\partial T_e/\partial t = 0$）后，系统趋于稳定且向内通量微乎其微，证实了电子热扩散是该阶段的主要驱动力。

2. 后期阶段（正 $\eta_\alpha$）：

   • 条件： 一旦密度剖面趋于平坦并开始向中心峰值化，$\eta_\alpha$ 变为正值。此时热扩散驱动减弱，且湍流通量通常转为向外。
   • 机制： 一种持续存在的非理想平衡态 $E \times B$ 通量驱动了中心峰值化。这种通量源于平衡态密度（$\bar{n}_1$）和静电势（$\bar{\phi}_1$）的极向不对称分量。
   • 物理本质： 虽然经典新古典理论（Pfirsch-Schlüter 机制）预测这些不对称性会在高阶项中相互抵消（由于玻尔兹曼电子响应），但模拟显示存在一个次级通量。广义欧姆定律中的非理想项——特别是平行电阻率、电子惯性和电磁感应——打破了玻尔兹曼关系。这产生了一个有限的极向相位差（$\delta_s \simeq -0.08\pi$），介于 $\bar{n}_1$ 与 $\bar{\phi}_1$ 之间，从而导致净向内通量。
   • 验证： 关闭驱动不对称性的横向热通量项后，峰值剖面发生弛豫，证实了 Pfirsch-Schlüter 不对称性与后期阶段 pinch 之间的联系。

主要贡献

• 自洽的剖面演化： 证明了通过通量驱动的全局边缘模拟可以自发地生成中心峰值的密度剖面（从平坦剖面开始），而无需固定边界条件或预设 pinch 速度。
• 双机制识别： 揭示了粒子 pinch 并非单一现象，而是从由湍流驱动的热扩散过程（早期阶段）向非理想平衡态 $E \times B$ 过程（后期阶段）的转变。
• 非理想平衡态通量： 确定了一种特定的机制，即电子力平衡中的非理想修正（电阻率、惯性、感应）打破了经典的玻尔兹曼响应，使得 Pfirsch-Schlüter 不对称性能驱动净径向通量。这补充而非矛盾于经典新古典理论。
• 通量分解： 提供了一种严谨的径向粒子通量分解方法，区分了湍流跨场输运、经典新古典项以及次级平衡态 $E \times B$ 通量。

意义与主张
本文声称这些结果解决了关于边缘密度剖面的长期疑问，特别是关于密度台阶的形成。通过展示向内粒子输运可以从基础等离子体物理（漂移波湍流和非理想平衡态效应）中自洽产生，而非依赖外部假设，本研究为边缘输运提供了更完整的图景。作者指出，尽管所使用的特定碰撞闭合（Braginskii）适用于高密度、低温度的边缘（如 C-Mod），但所识别的机制——极向对称性破缺和非玻尔兹曼响应——可以推广到包括仿星器和带有共振磁扰动在内的其他托卡马克配置中。这项工作强调了进行通量驱动、长时标模拟对于捕捉边缘粒子输运全复杂性的必要性。