跳到主要内容
Sandbox Physics
等离子体动理学 · 边学边做实验

L2 · 第 3 章

消失的波纹躲在哪里,网格又怎样把它放回来

密度波纹消退,是因为分布在速度方向折成了条纹,而且条纹以固定的速率变细。计算机把这些条纹存在速度网格上。条纹细过一个格子时,网格会把它误读成粗条纹,波纹便在一个只由网格决定的时刻重新出现。

推导

条纹以固定速率收紧

把波纹写成复数模式 eikxe^{ikx},记住真实的波纹是它的实部。自由流动下,eikxe^{ikx} 变成 eik(x−vt)=eikxe−iktve^{ik(x-vt)}=e^{ikx}e^{-iktv}。在固定位置看,分布在速度方向以波数 κ=kt\kappa = kt 振荡:条纹的密度和时间成正比。

密度要对速度积分。光滑的麦克斯韦分布乘上快速振荡再积分,结果几乎为零,振荡越快越接近零,这就是相混合。打开电场后,电子之间还会互相推,但条纹收紧的速率几乎不变。

eik(x−vt)=eikx e−iκv,κ(t)=kte^{ik(x-vt)}=e^{ikx}\,e^{-i\kappa v},\qquad \kappa(t)=kt

网格能存的条纹有极限

间距为 Δv\Delta v 的速度网格只在 vj=jΔvv_j=j\Delta v 处取样 e−iκve^{-i\kappa v}。当 κΔv\kappa\Delta v 达到 2π2\pi,每个取样值都回到 κ=0\kappa=0 时的数值。网格里存的又是原来那道光滑波纹,速度积分不再相消,电场就回来了。

回来的时刻叫复现时间,只取决于波数和速度间距,与等离子体的物理无关。它准时到来、不看物理,所以加密速度网格是最干净的检验:物理信号原地不动,复现会往后推。

k TR Δv=2π⟹TR=2πk Δvk\,T_R\,\Delta v=2\pi\qquad\Longrightarrow\qquad T_R=\frac{2\pi}{k\,\Delta v}

完整例题

例题:同一道波纹,两种网格

相空间 Lab 的速度范围是 −6-6 到 66。用 256256 格时,Δv=12/256≈0.047\Delta v=12/256\approx0.047,k=0.5k=0.5 的波纹要到 TR=2π/(0.5×0.047)≈268T_R=2\pi/(0.5\times0.047)\approx268 才复现,远在 70 个时间单位之外。t=40t=40 时一道条纹在速度方向宽 2π/(kt)≈0.312\pi/(kt)\approx0.31,仍占大约七格,网格分得清。

换成 3232 格,Δv=0.375\Delta v=0.375,TR≈33.5T_R\approx33.5。运行“粗网格”预设,电场在 t=33.5t=33.5 附近重新爬升,正好落在历史图的虚线复现标记上。复现之前测得的衰减率仍是 −0.1536-0.1536,和细网格相同,因为早期拟合窗口在网格失效之前就结束了。

三种会“回来”的信号

回来的电场可能是物理,可能是网格复现,也可能是速度边界的反射。Lab 的速度方向首尾相接,从上边出去的条纹会从下边进来;只要分布在边界附近可以忽略,这就无害。麦克斯韦分布的运行取 ±6\pm6,两束电子的运行取 ±9\pm9,都满足这一点。

每种嫌疑都有自己的检验。速度格数加倍,复现时间跟着加倍;间距不变而加宽速度范围,边界效应会变。物理信号经得起两种检验。项目的基准测试对第 6 章的回波做了这两项:用 128128、256256、512512 格,以及速度上限 66 和 1212,回波时刻都不动。

带着问题动手

先预测,再改变,最后解释

  1. 先算出 k=0.5k=0.5 的波纹在 64 格速度网格上的复现时间。
  2. 在相空间 Lab 里选 64 个速度格,看电场何时回来;再试 128 格。
  3. 用条纹频谱图解释这次回归:测得的条纹功率沿着金线走,直到条纹细到网格的极限。

三个相空间 Lab 用等离子体单位求解一维电子模型,每个 Lab 都能导出参数、模型版本和电场历史。

试着自己解释

1. k=0.5k=0.5、速度范围 [−6,6][-6,6] 分成 64 格,复现时间是多少?

展开答案与解释

Δv=12/64=0.1875\Delta v=12/64=0.1875,TR=2π/(0.5×0.1875)≈67T_R=2\pi/(0.5\times0.1875)\approx67,刚好落在 70 个时间单位以内。

2. 为什么复现集中在一个时刻,而不是慢慢出现?

展开答案与解释

对所有网格点都有 kTRvj=2πjkT_R v_j=2\pi j,每个取样同时回到起始相位。原先让波纹隐身的相消,在同一刻处处解除。

3. 同事在 64 格网格上看到 t=67t=67 时电场复活,说这是回波。用哪一个检验就能判定?

展开答案与解释

改用 128 格重跑。若是复现,会挪到 t≈134t\approx134;若是物理上的复活,仍停在 t=67t=67。

可以放在旁边读的来源

  1. R. Fitzpatrick, Plasma Physics lecture notes, University of Texas at Austin: Landau damping, resonant particles and the plasma dispersion function

    相连的三页推导麦克斯韦等离子体的线性响应、波速附近电子的作用,以及我们求介电根时用到的等离子体色散函数。

  2. C. Z. Cheng and G. Knorr, The integration of the Vlasov equation in configuration space, Journal of Computational Physics 22, 330 (1976)

    把一步 Vlasov 推进拆成位置平移和速度平移,通常引用这篇论文,本站的求解器也这样做。只核对了书目信息,求解器靠精确极限和守恒律来检验。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策