重联的流体图景在哪里不再够用

从指定电流片的小扰动开始

打开Loureiro、Schekochihin与Cowley的astro-ph/0703631v1。从方程3附近的平衡与匹配条件，读到方程18的内外解匹配。平衡系数由原文的磁场剖面决定。

把完整片长 $L_{\rm CS}$ 、半片厚 $\delta_{\rm CS}$ 与参考流动率 $\Gamma_0$ 写在纸边。无量纲波数 $\kappa$ 衡量沿薄片的波长；薄片展开还限制这些量的组合。

增长率描述指定薄片附近的小扰动。它还没有给出非线性重联速率、耀斑持续时间或磁岛离开薄层的速度。

\Gamma_0=\frac{2v_A}{L_{\rm CS}},\quad\kappa=\frac{k v_A}{\Gamma_0},\quad\epsilon=\frac{2\delta_{\rm CS}}{L_{\rm CS}}=\sqrt{\frac{2}{S_{\rm full}}}

这套完整片长定义，与前章按半片长定义的Lundquist数有固定换算。

方程18是一个含Euler伽马函数 $\Gamma$ 的标量匹配关系。这个数学函数与流动率 $\Gamma_0$ 是两回事。参考程序在 $0<\Lambda<1$ 的物理分支求根，用对数伽马函数的绝对值和已知符号避免数值溢出。

程序先由平衡匹配求出约 $\alpha=1.30693$ 的系数，再随波数解增长率、寻找峰值。改变薄度后，可复得增长率对完整片长Lundquist数的四分之一次方，以及最快波数的八分之三次方。

原图4的竖线标出展开条件开始失效的位置，自算图也区分超域的形式延拓。 $\epsilon=10^{-3}$ 时，峰值处 $\kappa\epsilon\simeq0.338$ ；取 $10^{-6}$ 时约为 $0.0601$ 。小于一，并不自动意味着渐近小量，所以求根的小数位不能照搬成物理精度。

-\frac\pi8(\kappa\alpha)^{1/3}\Lambda^{5/4}\frac{\Gamma[(\Lambda^{3/2}-1)/4]}{\Gamma[(\Lambda^{3/2}+5)/4]}=\frac{2\alpha^2}{\kappa\epsilon}

固定版本的方程18。所选分支的分子伽马函数为负，前面的负号再翻一次，左边就为正。

\frac\gamma{\Gamma_0}=\Lambda\alpha^{2/3}\kappa^{2/3},\qquad\gamma_{\max}/\Gamma_0\propto S_{\rm full}^{1/4},\quad\kappa_{\max}\propto S_{\rm full}^{3/8}

这些是带有指定平衡和展开条件的线性渐近预测。

python -W error scripts/astronomy/ac28/run_reference.py
python -W error scripts/astronomy/ac28/test_reference.py

再打开Zhou等人2019年发表于Physical Review Research的论文。方程1与2给出周期不可压缩约化MHD；数值研究从一列磁岛开始，与我们默认的双尺度受应力磁势不同。

先看图1：颜色表示电流，等值线表示磁通，随时间可见并合后的结构变大。图2画不同初始磁岛尺度Lundquist数下的总磁能。主图时间除以重联时间尺度，插图才用代码时间，这两种横轴不能混着比较。

论文实际采用 $2048^2$ 到 $8192^2$ 网格，磁Prandtl数为一，初始磁岛尺度Lundquist数从125到1786。文中明确说，这些运行不预期出现次级等离子体团，目视结果也与此相符。这是磁岛并合研究。

图3将拟合的时间幂指数与层级模型比较，后面的谱和缩放讨论再检验能量怎样移向较大尺度。读指数时，要同时读初态、分辨率与拟合区间。单个受应力结构的128网格画面，还没有复算这项统计结果。

最后打开Burch等人2016年的Science论文。观测发生在2015年10月16日的地球磁层顶。读图3至5时，先对齐电场、磁场、电子速度和分布各自的时间轴，再解释某一幅面板。

电子分布的采样间隔为三十毫秒。速度空间里的新月形，表示不同速度上有多少电子，并不是等离子体照片里一弯发亮的气体。要把航天器时间序列换成空间结构，还要知道边界运动与穿越几何。

比较离子、电子矩给出的电流与多航天器磁旋度给出的电流。电子参考系中的非理想电场为 $\mathbf E+\mathbf v_e\times\mathbf B$ ，其做功指标取电流与该电场的点积。指标出现峰值，并不单独确定局部熵产生率。

到这些尺度，压力张量、电子惯性和不同粒子种类的相对运动都进入问题。这个二维流体例子用电阻率联系非理想电场与电流；研究电子新月形还需要各速度上的粒子分布。

\mathbf J=en(\mathbf v_i-\mathbf v_e),\qquad\mathbf E'_e=\mathbf E+\mathbf v_e\times\mathbf B

这里近似电中性，离子带单个正电荷。

e

为正元电荷，

n

为共同数密度；非相对论参考系跟随电子整体运动。

D_e=\mathbf J\cdot\mathbf E'_e

电子参考系场做功指标，与标量碰撞加热或实测X射线强度都不同。

真正的光球足点计算，需要规定边界运动，并记录由此进入的磁能。要研究耀斑何时起始，还须有平衡、储能阶段及得到支持的触发机制。从一个已经受应力的周期场开始，还没有完成这条路线。

要预言X射线亮度，就补上发光气体和仪器响应；要预言加速粒子，就求分布或轨道及其能量交换。三维几何还会改变磁连接与允许的模态。每一项延伸都要找到新输入、新观测量与新检验。

研究日冕加热或太阳风时，也可以这样提出问题：先写清实际测了什么、来自哪里、在哪个尺度，再选择能求出它的模型。这三篇论文提供了一条可行的阅读次序，从受控线性薄片，走到充分分辨的流体动力学，再走到电子尺度观测。