---
title: "从线性增长读到 MRI 湍流的去留"
description: "按原论文图号读懂Hall与双极扩散实验，区分线性增长、湍流持续和净磁通对输运的影响。"
canonical_url: "https://sandboxphysics.com/zh/stories/magnetorotational-instability/read-beyond-linear-mri/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/magnetorotational-instability/read-beyond-linear-mri/"
last_modified: "2026-11-12"
---

L4 · 方法指南

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

Sandbox Physics · 2026-11-12

线性 MRI 增长率描述微小扰动最初怎样变化。要问盘能否持续输运角动量，还要知道气体从什么状态出发、有多少磁通穿过盘，以及气体与磁场耦合得多紧。一个很有帮助的论文实验，是先让湍流形成，再改变磁流体方程中的物理过程。

打开白雪宁的 [Hall-effect Controlled Gas Dynamics in Protoplanetary Disks: II. Full 3D Simulations toward the Outer Disk，arXiv:1409.2511v1](https://arxiv.org/pdf/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](https://arxiv.org/pdf/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与双极扩散开关实验和另一个微小线性通道问题](https://sandboxphysics.com/_astro/07-reading-map.zh.desktop.CxcfhTPS.svg)

原创阅读示意图，分开标出论文的开关实验与微小扰动问题；没有复制论文曲线。

## 读清坐标轴和参数

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 节](https://arxiv.org/pdf/1409.2511v1)。

再抬头看 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](https://arxiv.org/pdf/1210.6661v2)。

“饱和”指非线性演化限制了原来的增长，并形成后期的统计状态。线性方程会继续给出指数增长。要得到表中的长期 $\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](https://arxiv.org/pdf/1409.2511v1)）；Bai与Stone的arXiv1210.6661v2，式7、3.1.3／3.2.1节及表2（[固定PDF](https://arxiv.org/pdf/1210.6661v2)）。

## 动手检验预测

### [分开比较初期增长与持续性](https://sandboxphysics.com/zh/labs/magnetorotational-instability/)

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

## 继续阅读

磁张力怎样放大吸积盘中的扰动 · 阅读路径

### 磁张力怎样放大吸积盘中的扰动

1. [文章 **吸积盘的角动量由谁向外送？**](https://sandboxphysics.com/zh/stories/magnetorotational-instability/who-carries-angular-momentum/)
2. [文章 **旋转小区域里的四个实场**](https://sandboxphysics.com/zh/stories/magnetorotational-instability/four-fields-in-a-rotating-patch/)
3. [文章 **先算增长，再跟着能量检查**](https://sandboxphysics.com/zh/stories/magnetorotational-instability/growth-stress-and-energy/)
4. [文章 **竖直网格为什么会改变线性通道**](https://sandboxphysics.com/zh/stories/magnetorotational-instability/why-a-vertical-mesh-changes-a-channel/)
5. [文章 **怎样检验增长、传输与剪切功**](https://sandboxphysics.com/zh/stories/magnetorotational-instability/checking-growth-and-shear-work/)
6. [文章 **从保存的场数据检查 MRI 计算**](https://sandboxphysics.com/zh/stories/magnetorotational-instability/check-saved-mri-fields/)
7. [当前页面 **从线性增长读到 MRI 湍流的去留**](https://sandboxphysics.com/zh/stories/magnetorotational-instability/read-beyond-linear-mri/)