约束损失与逃逸粒子分析

在磁约束聚变装置中，等离子体被强磁场约束。带电粒子绕磁力线做回旋运动，同时沿磁力线运动。理想情况下，粒子被限制在磁面上。但实际存在多种机制导致粒子离开约束区，称为约束损失。约束性能常用粒子约束时间 $\tau_p$ 和能量约束时间 $\tau_E$ 衡量。损失分析关注损失锥、损失分数、首次通过时间、损失位置等。

约束损失可分为两类：一类是粒子轨道本身偏离磁面，最终撞击壁面；另一类是碰撞或湍流导致跨场输运，使粒子逐渐扩散出去。对于高能粒子，如聚变α粒子、中性束注入离子和逃逸电子，其轨道偏移更大，损失风险更高。

主要损失机制

• 经典输运：库仑碰撞导致跨场扩散，扩散系数约为 $D \sim \rho_i^2 \nu$，其中 $\rho_i$ 为离子回旋半径，$\nu$ 为碰撞频率。
• 新经典输运：环向场不均匀导致捕获粒子，形成香蕉轨道。香蕉轨道宽度 $\Delta_b \sim q \rho_i / \sqrt{\epsilon}$，其中 $q$ 为安全因子，$\epsilon$ 为反纵横比。当碰撞频率与反弹频率相当时，损失显著增大。
• 湍流输运：漂移波、离子温度梯度模、捕获电子模等引起异常输运，通常远大于新经典输运。
• 磁流体不稳定性：撕裂模、磁岛、边缘局域模和破裂会破坏磁面，导致快速粒子损失。
• 磁纹波损失：环向场线圈数目有限，产生纹波磁场，使捕获粒子离开约束区。
• 高能粒子损失：聚变α粒子、中性束注入离子和逃逸电子因轨道偏移大，容易损失。

逃逸粒子产生与动力学

逃逸电子是约束损失分析中的重要对象。当电场加速超过碰撞阻力时，电子进入逃逸状态。临界速度 $v_c$ 满足：

$$eE = m_e \nu_c v_c$$

其中 $e$ 为电子电荷，$E$ 为电场，$m_e$ 为电子质量，$\nu_c$ 为碰撞频率。逃逸电子的产生主要有两种机制：

1. Dreicer机制：热电子尾部被电场加速，直接进入逃逸区。
2. 雪崩机制：已有逃逸电子通过大角度碰撞产生次级逃逸电子，导致逃逸电子数量指数增长。

相对论效应会降低临界速度，使雪崩阈值约为 $E/E_D \sim 0.01-0.1$，其中 $E_D$ 为Dreicer电场。逃逸电子能量可达数十MeV，方向性强，可能严重损坏第一壁。

分析方法与建模

约束损失与逃逸粒子分析通常采用以下方法：

• 导心近似：将快速回旋运动平均，追踪导心轨道，降低计算维度。
• 轨道追踪：直接积分粒子运动方程，判断粒子是否穿过最后闭合磁面。
• 蒙特卡洛碰撞：模拟库仑碰撞和辐射过程，统计损失概率。
• Fokker-Planck方程：描述速度空间分布函数演化：

$$\frac{\partial f}{\partial t} + \mathbf{v}\cdot\nabla f - \frac{e}{m}\mathbf{E}\cdot\frac{\partial f}{\partial \mathbf{v}} = C(f) + S$$

其中 $C(f)$ 为碰撞算子，$S$ 为源项。

• 损失面映射：计算粒子损失位置和损失时间，评估壁面热负荷。

简化导心轨道追踪伪代码如下：

for each particle:
    while t < t_max:
        v_parallel += (e * E / m) * dt
        v_perp = sqrt(v**2 - v_parallel**2)
        R = R0 + rho * cos(theta)
        Z = rho * sin(theta)
        if R > R_wall or Z > Z_wall:
            record_loss()
            break
        t += dt

诊断与缓解

实验上可通过硬X射线谱、电子回旋辐射、汤姆逊散射、损失探针和红外热像诊断逃逸电子与损失粒子。缓解措施包括：

• 增大电子密度，提高碰撞频率，抑制逃逸电子产生。
• 注入杂质，通过辐射能量耗散降低电场。
• 施加共振磁扰动，增强逃逸电子损失。
• 使用破碎弹丸或激光吹气，快速提高密度并耗散能量。

总结

约束损失与逃逸粒子分析是磁约束聚变研究的关键环节。理解经典输运、新经典输运、湍流输运和不稳定性等损失机制，建立导心轨道、蒙特卡洛和Fokker-Planck等模型，结合诊断与缓解手段，可有效提高等离子体约束性能并降低逃逸电子危害。对于入门者，建议从单粒子轨道和损失锥概念入手，逐步过渡到统计输运与自洽模拟。