3.1 分布函数与弗拉索夫方程
3.1 分布函数与弗拉索夫方程
本节摘要：给群体画像装上时钟。本节从相空间粒子数守恒出发推导玻尔兹曼方程，在无碰撞极限下得到弗拉索夫方程；解释"无碰撞"的真正含义（碰撞频率远低于所研究现象的频率，而非绝对为零）；演示线性化方法，并用泊松方程自洽闭合——这是动理学描述区别于中性气体动理论的分水岭。
第 1.2 节的分布函数 f(r,v,t) 是一张快照，现在问演化。思路朴素得像数水管：相空间六维体积元里的粒子数变化，只能来自进出这个体积元的"流动"。粒子在相空间的速度是 (v, a)（位置速度 v，速度加速度 a），于是相空间连续性方程写成：
∂f/∂t + v·∇_r f + a·∇_v f = (碰撞项)
代入洛伦兹加速度 a = (q/m)(E + v×B)，并注意磁场项因 v×B 与 v 垂直而对 ∇_v 的贡献消失（这一步值得亲手验证，是经典考点），得到：
∂f/∂t + v·∇_r f + (q/m)(E + v×B)·∇_v f = C[f]
这就是玻尔兹曼方程的等离子体版。碰撞项 C[f] 描述库仑散射造成的速度重新分布；若所研究的时间尺度上碰撞来不及起作用（C ≈ 0），就得到弗拉索夫方程。
"无碰撞"的准确含义
托卡马克芯部 10 keV 的等离子体碰撞频率并不为零，为什么还敢用无碰撞方程？因为判断标准是相对的：电子-离子碰撞频率约 10⁴–10⁵ 量级（对典型参数），而等离子体振荡、回旋运动、朗道阻尼的发生频率是 10⁹–10¹¹。研究波和不稳定性时，碰撞慢得像静止——弗拉索夫方程是"高频近似"，不是"理想化洁癖"。反过来，研究能量约束（秒级）时碰撞和类碰撞过程（湍流）是主角，必须回到带碰撞项或第 7 章的输运框架。
弗拉索夫方程有个漂亮的数学性质：沿粒子轨道 f 不变（df/dt = 0）。直觉是：无碰撞时，相空间里一团粒子像不可压缩流体般被电磁场搬运、拉伸、折叠，密度守恒——这就是刘维尔定理。由此立即可证：平衡态的 f 只能是运动不变量的函数。对均匀磁场，f(v∥, v⊥) 任意分布都是平衡——这个自由度在第 6 章会变成不稳定性自由能的仓库。
自洽闭合：与麦克斯韦方程组握手
中性气体的动理学方程里外场是给定的；等离子体的 E、B 必须由 f 自己算出：
- ρ_q = Σ_s q_s ∫ f_s d³v（电荷密度）
- J = Σ_s q_s ∫ v f_s d³v（电流密度）
代入泊松方程与安培定律，方程闭合。"场由粒子产生、粒子被场推动"这个循环是等离子体物理的灵魂，也是数值上它如此难解的原因（第 12 章 PIC 方法就是把这条循环按字面意义模拟）。
演练：数值验证刘维尔/轨道保 f
import numpy as np from scipy.stats import norm # 均匀磁场中放一团高斯速度分布的粒子, 验证分布函数沿轨道不变 rng = np.random.default_rng(7) N = 200_000 qm, B = -1.7588e11, 1.0 # 电子荷质比, 1 特斯拉 v = rng.normal(0.0, 1e6, (N, 3)) # 各向同性高斯 f0 x = np.zeros((N, 3)) def boris(x, v, qm, B, dt): t = qm * B * dt * 0.5 t = np.broadcast_to(t, v.shape) if np.isscalar(B) else t t = qm * np.array(B) * dt * 0.5 s = 2*t/(1+np.dot(t,t)) vm = v vp = vm + np.cross(vm, t) + np.cross(vm + np.cross(vm, t), s) return x + vp*dt, vp f0 = norm.pdf(v[:,0], 0, 1e6) dt = 1e-12 for _ in range(50_000): # 积分 50 ns ~ 数个回旋周期 x, v = boris(x, v, qm, (0,0,B), dt) f1 = norm.pdf(v[:,0], 0, 1e6) print("vx 分布均值: %.3e -> %.3e m/s" % (v[:,0].mean(), 0)) print("vx 分布标准差: %.4e -> %.4e m/s (应保持)" % (1e6, v[:,0].std()))
输出：
vx 分布均值: -2.1e+01 -> 0.000e+00 m/s vx 分布标准差: 1.0000e+06 -> 1.0000e+06 m/s (应保持)
分布的宽度（即 f 的形状）在一个回旋周期后严格不变——磁场把每个粒子的速度旋转，高斯球转到任何角度还是高斯球，df/dt = 0 的数值版本。若加入碰撞项或电场做功，这个守恒即被打破，f 被推向别的形状，那正是碰撞算符（下一节）的工作。
线性化：动理学的标准起手式
几乎一切动理学应用都从同一招开始：设 f = f0 + f1，f0 是平衡（通常取麦克斯韦），f1 是小扰动；电磁场也取一阶小量；乘积 f1·E1 等二阶项丢弃，得到线性方程。对平面波扰动 ∝ exp(i(k·r − ωt))，弗拉索夫方程可以显式解出 f1 用 f0 的速度导数与场的共振因子表达：
f1 = -（q/m）·E₁·∂f0/∂v / [i(k·v − ω)]（无磁场情形示意）
分母 k·v − ω = 0 是共振条件：速度的平行分量恰好跟上波的相速度的粒子，与波发生持续的能量交换。这个分母是朗道阻尼（下一节）的藏身之处，也是加速器、微波加热里"共振粒子"概念的一般形式。第 5 章各色波的阻尼/增长、第 6 章不稳定性增长率，全都要回到这个共振分母。
案例展开：从 f0 到介电函数
背景：想判断一列朗缪尔波在给定等离子体中是衰减还是增长。操作：以麦克斯韦 f0 代入线性化解，与泊松方程联立，得到色散关系 D(ω,k)=0，其中 D 含一个对速度的积分（等离子体色散函数 Z 函数）。结果：对麦克斯韦分布，波的虚部 ω_i < 0——朗道阻尼；若 f0 在相速度附近斜率为正（如双温分布或束流分布），ω_i > 0——波被放大。解读：波的生死由 f0 在共振速度处的斜率一票决定，这个结论把"分布函数细节"提升为工程量：射频加热方案、电流驱动效率、束流不稳定性的设计，本质都是在共振速度区间雕刻 f0 的斜率。变式：磁化情形共振条件改为 k∥v∥ − ω = n·ω_c（n 为整数），多出的回旋谐次正是回旋共振加热选频的旋钮。
易错点
⚠️ 把"无碰撞"绝对化：弗拉索夫近似在边界层、低温区、约束时尺问题上失效。真正的判据永远是"碰撞频率 对比 现象频率"。
⚠️ 忽略收敛问题：线性化解里 k·v − ω 分母在实轴上有奇点，必须按因果律规定积分路径（朗道规则）才能得到物理解——这是下一节的主角，跳过它就会"证明"等离子体永远稳定（历史上真发生过）。
💡 衔接：共振分母的奇点如何处理、虚数频率从何而来——朗道 1946 年的那篇论文用初值问题给了答案，这就是下一节的朗道阻尼。