从线性增长读到 MRI 湍流的去留

打开白雪宁的 Hall-effect Controlled Gas Dynamics in Protoplanetary Disks: II. Full 3D Simulations toward the Outer Disk，arXiv:1409.2511v1，提交日期为 2014 年 9 月 8 日。固定使用这个版本；PDF 页眉另有较晚的内部草稿日期。阅读第 2 节，重点看印刷页 2—6 的 2.1—2.3 节。主图是印刷页 4 的 Fig.2，从一开始数也是 PDF 第 4 页。

再准备一份旁读材料：白雪宁与 James M. Stone 的 Local Study of Accretion Disks with a Strong Vertical Magnetic Field: Magnetorotational Instability and Disk Outflow，arXiv:1210.6661v2，日期为 2013 年 2 月 17 日。稍后用它的式 7、3.2.1 节和印刷页／PDF 第 10 页的表 2 比较净磁通。这次只核读指定段落和图表，没有重做两篇论文的模拟。

先找 Fig.2 改变了什么

从 Fig.1 下方的 2.3 节实验设置读起。这里用三维、无竖直分层、等温的剪切箱模拟，包含 Hall 项和双极扩散。“无分层”指没有加入竖直重力。箱体沿径向、轨道方向和竖直方向分别长 $4H$ 、 $4H$ 和 $2H$ ，其中 $H=c_s/\Omega$ ， $c_s$ 是声速， $\Omega$ 是局部角频率。

看曲线前先读图注。背景竖直磁场与转动方向相反，初始气体压强与竖直磁压之比为 $\beta_0=10^4$ ，Hall 长度为 $l_H=0.1H$ 。开始时使用理想 MHD，直到 $t=60\Omega^{-1}$ 才打开 Hall 项和双极扩散。此后的三条彩色曲线比较双极扩散耦合数 $\mathrm{Am}=100,10,1$ 。

开关打开时，气体已经处于湍流状态。这次实验问的是：已有的有限扰动能否持续？线性本征值问题问的则是：给指定背景一个任意小的扰动，它起初怎样变化？初态的区别会改变问题本身。

原创阅读示意图，区分论文中的Hall与双极扩散开关实验和另一个微小线性通道问题 — 原创阅读示意图，分开标出论文的开关实验与微小扰动问题；没有复制论文曲线。

读清坐标轴和参数

Fig.2 横轴的时间单位是 $\Omega^{-1}$ 。一整圈轨道周期为 $2\pi\Omega^{-1}$ ，所以开关约在 9.55 圈之后打开，不能把 60 直接读成 60 圈。纵轴是对数刻度的麦克斯韦应力。论文取磁导率、参考密度和声速都为一，曲线相当于把 $-\langle B_xB_y\rangle$ 除以参考气压 $P_0=\rho_0c_s^2$ 。纵向相等的距离表示相等的倍数变化；曲线上的起伏是时间变化，不是误差棒。

2.1 节定义了几种耦合量：

\begin{aligned} \mathrm{Am}&=\frac{\gamma_i\rho_i}{\Omega},\\ l_H&=\frac{v_A}{\omega_H},\qquad \chi=\frac{\omega_H}{\Omega}. \end{aligned}

其中 $\gamma_i\rho_i$ 是中性粒子与离子的碰撞频率， $\omega_H$ 是 Hall 频率， $v_A$ 是阿尔芬速度。在这套参数约定中， $\mathrm{Am}$ 越小，中性气体与离子的耦合相对于轨道时钟越弱。论文把 Hall 长度和 $\chi$ 定义为正数，再单独写明磁场极性。比较 Hall 实验时，光知道场强还不够，还得保留方向。

开关之后，红色的 $\mathrm{Am}=100$ 曲线仍保持明显应力，绿色与蓝色的 10、1 两组则大幅衰减。图注明确说明，这组背景没有线性不稳定的 MRI 模。这里展示的是既有非线性状态在指定条件下能否存活，“持续”指图中截至 $300\Omega^{-1}$ 的时段。它没有检验静止背景从任意微小种子出发，能否走到同一个状态。读这张图时，图注、开关时刻和初态必须放在一起。Bai 2014，Fig.2 与 2.3 节。

再抬头看 Fig.1 做一次交叉检查。它的横坐标把磁场极性与 Hall 耦合数的倒数组合起来，纵坐标是 $kv_{A0}/\Omega$ ，颜色给出以 $\Omega$ 为单位的线性增长率。左右两幅分别固定 $\mathrm{Am}=100$ 和 1。上图的颜色回答小扰动问题，下图的时间曲线回答持续性问题。不必从图片描出数据点，也能读清两者各自问了什么。

比较净磁通前，先看压强分母

翻到 2013 年论文的式 7。它把应力沿高度积分，再按压强积分归一化：

\alpha=\frac{\int \overline{T}_{xy}(z)\,dz} {\int \overline{\rho}(z)c_s^2\,dz}.

横线表示水平平均，表中还做了时间平均。这和某一个高度的应力不是同一个量。那篇论文使用有竖直分层、带净竖直磁通的剪切箱，因此密度结构与竖直区域也进入比较。

做一道短算题：表 2 中 B4 的麦克斯韦与雷诺两列相加，得到 $0.061+0.014=0.075$ ；B2 则为 $0.92+0.12=1.04$ 。这些是表中舍入后的数值，对应初始 $\beta_0$ 分别为 $10^4$ 和 $10^2$ 。若初始气压相同，后者的竖直场强是前者的十倍，因为磁压随场强的平方变化，不能误读成场强增加一百倍。表中例子说明，输运水平会随外加净磁通改变。3.2.1 节还指出，强磁场 B2 的积分值受箱体竖直范围和整体磁场条件影响，应按那组设置理解它。Bai 与 Stone 2013，式 7、3.2.1 节及表 2。

“饱和”指非线性演化限制了原来的增长，并形成后期的统计状态。线性方程会继续给出指数增长。要得到表中的长期 $\alpha$ ，需要计算后期的非线性状态。

用线性 Lab 回答线性问题

MRI Lab 跟随理想、不可压缩、水平匀质的通道扰动。控件改变剪切 $q$ 和磁张力频率 $K=kv_A$ ；显示的 $K/\Omega$ 使用 $\Omega=1$ 的旋转参考。速度随高度按正弦变化，水平磁扰动按余弦变化，流体竖直速度为零。小振幅场把增长、应力和背景剪切供能连在了一起。

这里的应力按 $\rho L^2\Omega_{\rm ref}^2$ 归一化， $L$ 为选定的高度尺度。不可压缩模型没有用于组成论文 $\alpha$ 的热压强，也没有 Hall 极性、双极扩散或饱和控件。把整个背景磁场反向，理想线性问题的增长率保持相同；Bai 的非理想模型则需要区分极性。

先做可操作的线性比较：固定开普勒剪切，依次比较参考张力、最快增长张力和导致振荡的强张力，解释各组中的能量来源。再回论文写下它多用了哪些输入：非线性初态、气体压强、磁场极性和非理想耦合。这些输入说明两种计算之间具体差在哪里。

带着问题继续读

如果把 Fig.2 的开关移到最初时刻，只凭原图注能断定后来的结果吗？不能，因为初值实验已经变了。接下来的计算需要先说明种子扰动，再检查不同初态会不会走向相同的后期行为。

若两篇论文报告不同应力，先比较平均体积、压强分母、外加净磁通和平均时段，再看竖直重力、边界和非理想项。研究分辨率收敛时，应比较相同物理设置。2013 年论文的 3.1.3 节讨论了饱和水平尚未解决的收敛问题；分辨出一个线性波长，与让后期统计应力收敛，是两项不同的检验。

下次读 MRI 图时，在旁边记下五件事：版本和页码，初态与固定参数，坐标轴和归一化，产生比较的那次改变，以及图中量能支持的一个结论。这样，增长率图、衰减时间序列和饱和输运表就各有清楚的位置。

来源

核读范围：Bai的arXiv1409.2511v1，第2节及Fig.1—2（固定PDF）；Bai与Stone的arXiv1210.6661v2，式7、3.1.3／3.2.1节及表2（固定PDF）。

动手检验预测

分开比较初期增长与持续性

先预测同一磁张力下，刚性转动与开普勒剪切中的微小通道会不会增长。
在旋转参考中固定张力频率为0.6，比较剪切参数0与1.5。
读取初期增长或振荡，以及两种应力。
指出剪切能的来源，再列出论文持续性实验额外使用的初态与非理想输入。