拖曳以转 · 滚轮或两指开合以缩放
方载磁力线1,360条…
大幅之阿尔芬子(Alfvénon),即余折返磁力线探览(Switchback Field-Line Explorer)中之折返模型,今以三维示之:拖之可转。磁力线凡1,360条,各依其局部之角 θ = cos⁻¹(Bx/|B|) 着色;θ 者,与背景场 B₀ = x̂(箭头所示)之夹角也:场直处为钢蓝,60° 为品红,90° 为橙。过 90°,则场反折于己,是为折返,色续转为金:折者逾130条,最甚者 θ ≈ 145°。线之弯逾 25° 者,绘为实管,故每一折皆可于三维中循之。全本 Plotly 探览之页重59 MB;此图所载0.5 MB。

三维阿尔芬孤波

帕克太阳探测器察得,初生之太阳风中,满布大幅之阿尔芬式起伏,其状孤立。其为磁场之骤转,甚者全然反向(即所谓“折返”),处于本自平静之日冕磁场中,|B| 几近恒定,质子喷流独向一侧。教科书所载之非线性阿尔芬波,即球偏振阿尔芬波,以 B 之均值为其背景场 B₀。然于飞船所测之数据,此均值无定:虽同一太阳风流之内亦变,阿尔芬速度随之而变。

场强恒定之圆:背景矢量 b0 向上,波场 b1 与总场 b 皆止于圆上,速度 u1 沿 b1 而绘
|B| 恒定之约束:波场 b₁ 使总场 b = b₀ + b₁ 摆转于定径之球面上;流速 u₁ 与 b₁ 共线。出自 Huang et al. (2026a)。

吾等于 Huang et al. (2026a) 别立一法,书阿尔芬波为理想磁流体(MHD)之孤波解:局部之扰动,乘于未受扰之远场之上,|B| 恒定。欲于三维构此场,非易事也。须严守无散,且磁力线皆开放;|B| 须几近恒定;远场须不为所动。吾等以迭代之法构之:谱域之 Helmholtz–Hodge 投影以去其散度,|B| 之归一,二者交替而行,迭代乃收敛。名其所得曰阿尔芬子(Alfvénon)。

场偏转角 θ 之半透明等值面层层相套,于单位盒中成一局部之肾形阿尔芬子;灰色磁力线径直穿之
128³ 网格上之阿尔芬子。层层相套之面示 θ = cos⁻¹(Bx/|B|),即局部之场与未受扰之背景场 B₀ = x̂ 之夹角;其极值 ≈45°。灰色磁力线沿 x̂ 出入此结构:扰动局于一处,|B| 几近恒定,磁力线无不开放。出自 Huang et al. (2026a)。

以伪谱代码 LAPS 直接作磁流体模拟(β = 0.1,不设显式之黏性与电阻),见阿尔芬子以阿尔芬速度行进,持其形,仅有轻微之弛豫。扰动之三分量皆保阿尔芬性,总能量守恒至约 10⁻⁹。其徐徐变形,主因相混合:阿尔芬子之内,VA 之变约1%。密度之微小起伏,生于所构之场之小疵,如 |B| 未尽恒定。阿尔芬子于理想磁流体之下,其行止如非线性之阿尔芬式孤立解。

四行子图,时刻 0、5、10、100:一维之场与速度剖面,θ、阿尔芬速度与密度之切片
阿尔芬子之磁流体模拟,t = 0、5、10、100(各行)。(a) 过中心之 B 与 u 剖面,平移至同一位置;(b) x–y 中平面之 θ;(c) 阿尔芬速度 VA = |B|/√ρ;(d) 密度 ρ。注意 (c)、(d) 之色标,所跨仅 ±0.5%。出自 Huang et al. (2026a)。

阿尔芬子相撞 初果

二小阿尔芬子,屡屡相撞

二小阿尔芬子(θmax ≈ 45°)于 2×1×1 周期盒中迎面相撞,再三不已。模拟用伪谱磁流体代码 LAPS,网格 512×256×256。阿尔芬子独处时,|B| 几近恒定。二者相叠之处,|B| 偏离 1 数个百分点,其扰动化为压缩波之波前,向外扩展,辐射而去。左上:过中心沿 x 之 B(实线)与 u(虚线)。余图:三中平面之 |B|,t = 0 至 5.8。

湍流之级串,须反向传播之阿尔芬波相撞;阿尔芬子使纯净之相撞实验可行。孤立之阿尔芬子无压缩分量(δ|B| ≈ 0,δρ ≈ 0),传播之中始终如是。故二者相遇之时,凡所现之压缩能量,皆相撞所生。其理甚简:阿尔芬子各自 |B| 恒定,相叠则不然。相叠之处,量级为 b₊·b₋ 之磁压失衡激发磁声波,携去阿尔芬能量之一份,其量级为 (δb/B₀)²。帕克太阳探测器近日所见振幅甚大,于此振幅,此份之量级为一。

吾等方以 GPU 模拟迎面之相撞以验之:激波捕捉,双精度,网格至 1024×512×512,能量守恒优于 10⁸ 分之一。

<10⁻⁸每一双精度模拟之能量相对误差
5–6%二小阿尔芬子初遇既毕,内能之增幅
≈10%二者之一为大幅(θmax ≈ 164°)时之增幅
5迄今相撞模拟之数,自 512×256×256 至 1024×512×512

能量之转移随振幅而增;波包经相遇而出,受损而未毁。今方行者数事:

  • 相配之单波包对照与同向传播对照,以别相撞之加热与数值之加热;
  • 振幅之扫描,下至寻常波之极限;
  • 异形阿尔芬子之相撞;
  • 相撞所发快波、慢波之普查。

所期首用,在日冕之闭合磁环:足点反射,使反向传播之波为常。有灵台研究器灵,助行此番模拟,并录每一模拟之始末。

论文

  • Huang, Z., Velli, M., Shi, C., & Ding, Y. 2026, “Solitary Alfvén Waves,” ApJL 1006, L20. doi · arXiv:2512.02292
  • Huang, Z., Velli, M., & Ding, Y. 2025, “What are Switchbacks?” 穿越内嵌旋转间断之阿尔芬子,所得之磁场反转类乎折返。 arXiv:2512.12585
  • Ding, Y., Velli, M., Huang, Z., et al. 2026, “Spherically Polarized Alfvén Waves and the Gosling Boost.” arXiv:2608.05091

代码

  • helmholtz-skill:Helmholtz 交替投影求解器,以求无散、模恒定之三维场(亦见 Helmholtz_algo_public)
  • alfvenon-pluto-cuda:PLUTO 之 NVIDIA GPU 移植,已验证,开源,为阿尔芬子诸实验而作
  • LAPS:洛城加州大学(UCLA)之伪谱、MPI 并行三维霍尔磁流体(Hall-MHD)代码