阿尔芬波与磁流体波

在等离子体物理中，磁流体力学（MHD）将等离子体视为导电流体，用密度、速度、压强和磁场描述宏观行为。理想 MHD 的线性波解给出三类基本模式：阿尔芬波、快磁声波和慢磁声波，它们统称磁流体波。理解这些波的色散关系、偏振和传播方向，是分析太阳风、磁层、托卡马克以及数值模拟稳定性的基础。
理想 MHD 方程组可写为：

$$
\frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{v})=0
$$

$$
\rho\left(\frac{\partial \mathbf{v}}{\partial t}+\mathbf{v}\cdot\nabla\mathbf{v}\right)
=-\nabla p+\mathbf{J}\times\mathbf{B}
$$

$$
\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B})
$$

$$
\nabla\times\mathbf{B}=\mu_0\mathbf{J},\qquad \nabla\cdot\mathbf{B}=0
$$

其中 $\rho$ 为质量密度，$\mathbf{v}$ 为速度，$p$ 为热压强，$\mathbf{B}$ 为磁场，$\mathbf{J}$ 为电流密度。若考虑绝热闭合，$p/\rho^\gamma=\text{const}$，并定义声速

$$
c_s=\sqrt{\frac{\gamma p}{\rho}}
$$

以及阿尔芬速度

$$
v_A=\frac{B_0}{\sqrt{\mu_0\rho_0}}
$$

将背景量取为均匀场 $\mathbf{B}_0$，并对小扰动线性化，可得到平面波解。设波矢 $\mathbf{k}$ 与 $\mathbf{B}_0$ 的夹角为 $\theta$，则磁流体波分为阿尔芬波和磁声波。

阿尔芬波：磁张力驱动的横向波

阿尔芬波是一种非压缩横波，其恢复力来自磁张力。理想 MHD 中，阿尔芬波的色散关系为：

$$
\omega^2=k^2v_A^2\cos^2\theta
$$

因此相速度为

$$
v_{\text{ph}}=\frac{\omega}{k}=v_A\cos\theta
$$

当波沿背景磁场传播时，$\theta=0$，相速度最大，等于 $v_A$；当波垂直背景磁场传播时，$\theta=90^\circ$，相速度为零。阿尔芬波的群速度沿 $\mathbf{B}_0$ 方向，扰动表现为速度场与磁场横向振荡，密度和压强几乎不变，满足

$$
\nabla\cdot\delta\mathbf{v}\approx 0,\qquad \delta\rho\approx 0,\qquad \delta p\approx 0
$$

一个典型估算：在托卡马克边缘，取 $B_0=2,\text{T}$，$n=10^{19},\text{m}^{-3}$，质子质量 $m_i=1.67\times10^{-27},\text{kg}$，则

$$
\rho_0=n m_i=1.67\times10^{-8},\text{kg/m}^3
$$

$$
v_A=\frac{2}{\sqrt{4\pi\times10^{-7}\times1.67\times10^{-8}}}
\approx 1.4\times10^7,\text{m/s}
$$

即约 $1.4\times10^4,\text{km/s}$。这说明强磁场、低密度等离子体中阿尔芬波速度很高，会显著限制显式 MHD 模拟的时间步长。

快磁声波与慢磁声波

磁声波是压缩波，恢复力同时包含磁压和热压。其色散关系为：

$$
\omega^2=\frac{k^2}{2}\left[v_A^2+c_s^2
\pm\sqrt{(v_A^2+c_s^2)^2-4v_A^2c_s^2\cos^2\theta}\right]
$$

其中取“$+$”号为快磁声波，取“$-$”号为慢磁声波。相速度平方为：

$$
v_{\text{fast,slow}}^2=\frac{1}{2}\left[v_A^2+c_s^2
\pm\sqrt{(v_A^2+c_s^2)^2-4v_A^2c_s^2\cos^2\theta}\right]
$$

主要特征如下：

• 快磁声波：通常为压缩波，速度大于 $v_A$ 和 $c_s$ 中的较大者；在垂直传播时速度达到 $\sqrt{v_A^2+c_s^2}$。
• 慢磁声波：速度小于 $v_A$ 和 $c_s$ 中的较小者；在垂直传播时速度趋于零。
• 平行传播时，快波速度为 $\max(v_A,c_s)$，慢波速度为 $\min(v_A,c_s)$，阿尔芬波速度为 $v_A$。
• 在低 $\beta$ 等离子体中，$c_s\ll v_A$，慢波常表现为沿磁场方向的声学型压缩扰动。

三种模式的对比

在建模中，这三类波具有不同意义。阿尔芬波常用于描述磁力线横向振荡、日冕加热和阿尔芬本征模；快磁声波决定显式 MHD 模拟的 CFL 条件；慢磁声波则与沿磁场的压缩扰动和结构形成有关。

在等离子体建模中的应用与注意事项

进行 MHD 数值模拟时，需要重点关注以下问题：

1. 时间步长：显式格式通常要求
   $$
   \Delta t\le \text{CFL}\frac{\Delta x}{v_{\max}}
   $$
   其中 $v_{\max}$ 常取快磁声波速度。强磁场区域会显著减小 $\Delta t$。
2. 边界条件：阿尔芬波沿磁力线传播，边界若不能正确吸收或反射，会引入非物理驻波。
3. 耗散与分辨率：理想 MHD 无耗散，但数值格式会引入数值耗散，可能改变慢波和阿尔芬波的振幅。
4. 动理学修正：当波频率接近离子回旋频率，或垂直波长接近离子回旋半径时，理想 MHD 失效，需要考虑动理学阿尔芬波等模型。

总之，阿尔芬波与磁声波构成了磁化等离子体宏观波动理论的骨架。掌握其色散关系和偏振性质，不仅有助于理解空间与聚变等离子体中的能量输运，也是构建可靠 MHD 数值模型的前提。