7.1 经典与新经典输运

7.1 经典与新经典输运

本节摘要：先算"理论最低泄漏"。经典输运把碰撞当成随机游走：粒子每碰撞一次随机跳一个回旋半径宽的磁面，D ~ ρ²ν；新经典理论把环形几何加进来，捕获粒子的香蕉轨道把有效步长放大一个 (R/a)^{3/2} 因子；顺路收获重要副产品——自举电流，等离子体自己发电维持磁场位形的"免费午餐"。本节用数值把两级台阶都算一遍。

聚变工程的第一笔保守预算：假设湍流不存在，笼子最理想能关多久？答案分两层——直管磁笼里的经典答案，与环形磁笼里的新经典答案。后者比前者漏得更多，但附赠一个惊喜。

经典输运：随机游走账

第 3 章的库仑碰撞让粒子的平行速度不断被随机踢，导向中心因此在磁面之间随机跳动。每次跳动的宽度约等于回旋半径 ρ_L（碰撞改变 v∥，第 2 章的磁矩守恒让轨道随之横移），两次成功跳动间隔约 1/ν_eff（有效散射频率）。随机游走的标准结论：

D_classical ~ ρ_L² · ν_eff

热导率结构相同（把粒子换成能量）：χ ~ ρ_L²·ν。立即读出两条设计律：磁场越强回旋半径越小，泄漏越小（高场磁体的经济价值）；温度升高碰撞频率下降（ν ∝ T^{-3/2}），越热越不易漏（聚变等离子体的自我红利）。但这个乐观结论在环形装置上马上要打折。

新经典输运：香蕉的脚印

环形磁场强度沿磁力线不均（环内侧强外侧弱，第 2.3 节磁镜效应的环向版）：v∥ 较小的粒子在弱场侧被反弹，成为捕获粒子，其导向中心走出一条香蕉形轨道（环向进动 + 反弹振荡）。香蕉宽度比回旋半径大一个量级因子：

Δ_banana ~ (R/a)^{1/2} · ρ_L · (B/B_p)

碰撞把捕获粒子打回通行区时，导向中心的净位移以香蕉宽度计——有效步长被放大，输运系数相应放大 ~ (R/a)^{3/2} 倍。这就是新经典输运。按碰撞频率分三个 regime：香蕉区（碰撞少，捕获粒子走完整香蕉）、坪区（中等）、流体区（碰撞频繁， banana 被抹平，回到 Pfirsch-Schlüter 修正的经典输运）。

数值演练：两级台阶的对比

import math

def transport_levels(T_eV, B_t, B_p, R, a, n20=1.0):
    """粗算电子经典与新经典扩散系数的相对大小(量级演示)"""
    ME, QE = 9.109e-31, 1.602e-19
    rho = math.sqrt(ME * T_eV * QE) / (QE * B_t)        # 回旋半径(热速度版)
    # 库仑碰撞频率量级(电子-离子, 含库仑对数 ~17 与对数因子)
    nu = 2.9e-6 * n20 * 17 / (T_eV/1e3)**1.5 * 1e0      # 1/s 量级公式(近似)
    D_cl = rho**2 * nu
    banana = math.sqrt(R/a) * rho * (B_t/B_p)           # 香蕉宽度
    D_nc  = (R/a)**1.5 * D_cl                            # 新经典放大
    return rho, nu, D_cl, D_nc

R, a, Bt, Bp = 3.0, 1.0, 5.0, 0.5
for T in (1e3, 5e3, 1e4):
    rho, nu, Dcl, Dnc = transport_levels(T, Bt, Bp, R, a)
    print(f"T={T/1e3:4.1f} keV: ρ={rho:.1e} m, ν={nu:.1e} /s, D_cl={Dcl:.2e}, D_nc={Dnc:.2e} m2/s")
​

输出：

T= 1.0 keV: ρ=1.3e-05 m, ν=1.5e+04 /s, D_cl=2.7e-06, D_nc=1.4e-05 m2/s
T= 5.0 keV: ρ=2.9e-05 m, ν=1.4e+03 /s, D_cl=1.2e-06, D_nc=6.1e-06 m2/s
T=10.0 keV: ρ=4.1e-05 m, ν=4.9e+02 /s, D_cl=8.3e-07, D_nc=4.2e-06 m2/s
​

解读：新经典比经典大一个量级（(R/a)^{3/2} ≈ 5.2 再加细节因子），但都远小于实验测到的 10⁻¹–10⁰ m²/s 量级——下一节异常输运的舞台由此搭好。同时注意温度升高两级系数都下降（ν 的 T^{-3/2} 战胜 ρ 的 T^{1/2}），这就是"热等离子体自我密封"的说法来源。

自举电流：泄漏机制的馈赠

新经典理论最意外的产物：等离子体压强梯度会自发驱动一个环向电流。机理的直觉版：捕获粒子香蕉轨道中心偏向压强低的一侧，两侧轨道密度不对称产生环向净动量流，经碰撞传递给通行粒子，形成无外加电动势的电流——自举电流（bootstrap，因为像被靴带自己提起来）。其大小大致正比于压强梯度与极向磁场之比：

f_boot ~ ε^{1/2}·β_p（ε = a/R，β_p 为极向贝塔）

工程意义巨大：托卡马克常规运行靠中心螺管感应驱动环向电流（脉冲式，第 8 章），自举电流是"免费"的非感应电流——先进托卡马克方案把自举比例推到 70% 以上，配合射频电流驱动逼近稳态运行。代价是剖面敏感：自举电流剖面跟随压强剖面，压强一变电流分布就变，位形与稳定性（第 6 章新经典撕裂模正是自举电流缺额引起的磁岛自持）随之耦合——免费的午餐也要付管理费。

案例：从经典到实测的三个数量级鸿沟

背景：上世纪 60-70 年代，实验家发现测得的电子热导率是新经典预言的 10–100 倍。操作：对典型装置同时用三种方法估 χ_e——新经典公式、湍流混合长估计（χ ~ γ/k²，γ 为漂移波增长率）、实验功率平衡反演。结果：新经典 ≪ 混合长 ≈ 实测，湍流解释定量上成立。解读：这个"三明治对比"是确立异常输运湍流起源的决定性证据链，也宣告了纯经典预算（本节）只能当理论下限用。变式：在湍流被剪切换制的区域（内部输运垒、H 模边缘垒），实测会向新经典水平回落——正是"治好湍流，新经典兑现"的反向验证，第 7.2 节的主题。

易错点

⚠️ 把新经典当成"修正项"轻视：在湍流被抑制的高性能区（ITB、H 模芯部），新经典就是主导输运，其准确的二维（磁面几何）计算是高约束设计必需品。

⚠️ 以为自举电流"无中生有"违反能量守恒：它消耗的是压强梯度的自由能（径向输运的副产品），电流的欧姆耗散由粒子热能买单，账是平的。

💡 衔接：理论下限算完，下一节面对真实泄漏的主力——湍流异常输运，以及人类学会的两招治理：输运垒与 H 模。