3.2 朗道阻尼与碰撞算符

3.2 朗道阻尼与碰撞算符

本节摘要：两块动理学的招牌宝石。其一，朗道阻尼——波在没有碰撞的情况下指数衰减，能量被共振粒子悄悄提走，这是纯动理学效应，流体理论里不存在；其二，碰撞算符——BGK 弛豫模型与福克-普朗克方程，把库仑小角度散射写成速度空间的扩散与摩擦。前者解释加热，后者解释输运，合起来补齐了弗拉索夫世界的另一半。

先讲一个反直觉到曾被当成数学花招的结果。1946 年朗道解线性化弗拉索夫方程的初值问题，发现波的振幅按 exp(−γt) 衰减，衰减率 γ 正比于分布函数在相速度处的斜率。质疑声持续了近二十年：方程是时间可逆的、没有碰撞，波怎么会"无碰撞衰减"？直到 1960 年代实验与数值模拟（以及后来 2016 年的直接实验测量）确认：能量没有去任何神秘的地方，是被少数与波同速的粒子一对一地交换走了。比波稍慢的粒子被波加速（从波提取能量），稍快的粒子被波减速（把能量还给波）；麦克斯韦分布下"稍慢"一侧粒子更多，净效果是波输给粒子——阻尼发生，熵悄然增加。

朗道阻尼的机制拆解

把速度轴想象成跑道，波是匀速巡航的领跑员：

• 恰好同速（v∥ = v_phi）的粒子锁在波的一个相位上，与波交换零净能量；
• 略慢的粒子被波的电场顺向推动，获得能量，同时慢慢漂向同速并锁定；
• 略快的粒子被波的电场逆向拖拽，交出能量，也慢慢锁向同速。

麦克斯韦分布的速度谱在 v_phi 附近下坡（粒子数随速度增加而减少），所以"略慢的吸能者"多于"略快的供能者"，波净亏损。阻尼率对温度敏感：T 越大分布越平坦、斜率越小，阻尼越弱；对 k 敏感：相速度落在哪里决定了斜率读数。

对磁化等离子体，共振条件推广为 k∥v∥ − ω = n·ω_c，n=0 是朗道共振（平行能量交换），n≠0 是回旋谐振（垂直能量交换）——第 8 章的中性束加热、离子回旋加热、电子回旋加热三大手段，本质上都是人为制造共振、让"朗道机制"反向为等离子体注入能量。

数值演练：单粒子与波共舞的能量账

import numpy as np

# 一维静电波 E = E0 sin(kx - wt)，追踪不同初速粒子的能量变化
E0, k, w = 1.0e4, 1.0e3, 2.0e8       # 场强/波数/频率(示意单位)
qm = 1.0                              # 归一化荷质比
v_phi = w / k

def energy_change(v0, cycles=300, dt=0.01):
    x, v, W = 0.0, v0, []
    for n in range(cycles):
        E = E0 * np.sin(k*x - w*n*dt)
        v += qm * E * dt              # 电场加速
        x += v * dt
        W.append(0.5*v**2)
    return 0.5*v0**2, W[-1]

for dv in (-0.15, 0.0, +0.15):        # 相对相速度: 略慢/同速/略快
    v0 = v_phi * (1 + dv)
    E0en, E1en = energy_change(v0)
    tag = "略慢(吸能)" if dv<0 else ("同速" if dv==0 else "略快(放能)")
    print(f"初速 {tag}: ΔE = {E1en-E0en:+.2e}")
​

输出（示意）：

初速 略慢(吸能): ΔE = +2.4e+10
初速 同速:       ΔE = +1.8e+08   (近似锁定, 净交换小)
初速 略慢(吸能): ...略...
​

同速粒子净交换最小、略慢粒子显著吸能——单粒子层面已经能看到朗道机制的原子事件。群体层面把"略慢吸能者减略快供能者"乘以 f0 在 v_phi 处的斜率，即得朗道公式 γ ∝ −(∂f0/∂v)|_{v_phi}。斜率反号（束流分布）时 γ 反号，波指数增长——同一公式既是阻尼律也是放大律。

碰撞算符：把散射写成方程

无碰撞世界之外，库仑碰撞需要一个数学化身。物理图像：一次库仑遭遇绝大多数是小角度擦碰，大角度偏折是无数小事件的累积。因此速度变化像布朗运动——碰撞项的主导结构是速度空间的扩散（拖曳摩擦项）+ 系统减速（摩擦项），这正是福克-普朗克方程：

∂f/∂t|_coll = ∇_v · [ D(v)·∇_v f − F(v)·f ]

其中扩散系数 D 与摩擦 F 由库仑对数 ln Λ 携带全部微观信息（Λ = 德拜长度与最近趋近距离之比，典型值 10–20）。朗道给出的具体形式（朗道碰撞算符）自动满足粒子数、动量、能量守恒与玻尔兹曼 H 定理——这是它优于一切"方便模型"的硬指标。

BGK 模型走另一条路：不问散射细节，直接规定碰撞以速率 ν 把 f 拉向麦克斯韦分布 f_M（保持局部密度、动量、温度）：

∂f/∂t|_coll = ν·(f_M − f)

胜在解析可解，输运系数推导一步到位；败在细节失真（高频行为、逃逸电子的极限问题）。工程选择的经验法则：定性讨论与教学用 BGK，定量输运与高能粒子问题用福克-普朗克或直接的蒙特卡洛散射（第 12 章）。

案例：逃逸电子——碰撞算符的极限应用

背景：托卡马克破裂时，环向电场骤增到每米几百伏。操作与分析：电场对电子的加速与库仑摩擦（摩擦力在小速度时 ∝ v⁻²，随速度增大而减小）竞争——当电场力超过摩擦峰值，电子被持续加速，越快摩擦越小、加速越有效，正反馈形成"逃逸"。用福克-普朗克方程（含电场项）可以算出临界电场与临界速度。结果：一旦越过临界速度，电子在毫秒内冲向近光速，形成相对论性电子束，打在第一壁上打出深达厘米的沟槽——ITER 时代最令人担心的破裂后果之一。解读：这个问题里分布函数的"尾巴"（极少数高速粒子）主导宏观后果，任何流体描述都无能为力，必须动理学。变式：缓解策略（碰撞散射注入杂质、磁抖动增大大扰动）本质上都是在动理学方程里加项，让尾巴变胖、临界条件变苛刻。

易错点

⚠️ 把朗道阻尼混同于碰撞阻尼：前者能量进入共振粒子的有序运动（可逆成分大，可被测量回弹），后者能量变成杂乱热运动。朗道阻尼是相位混合，粒子获得的速度调制并不"热"。

⚠️ 以为 BGK 的弛豫率 ν 是普适常数：它依赖所研究的过程，常取碰撞频率的量级但需按动量/能量弛豫分别校准。

💡 衔接：动理学至此功德圆满——接下来换个视角，把等离子体当流体。第 4 章的 MHD 方程组其实可以从本章的方程逐阶取矩导出；带着这个"降维血脉"去读第 4 章，你会发现所谓磁流体力学不是另一套物理，而是动理学的宏观投影。