磁流体平衡方程推导
在等离子体物理中，磁流体力学（Magnetohydrodynamics, 简称 MHD）模型是研究高温等离子体宏观行为的基础框架。当等离子体处于相对稳定的状态，即随时间变化极其缓慢时，我们可以忽略动力学演化过程，将其视为处于平衡态。磁流体平衡的分析是理解托卡马克、仿星器等磁约束聚变装置中等离子体约束性能的前提。本文将系统地推导磁流体平衡方程，并深入探讨其物理意义与应用。
磁流体力学模型将等离子体视为导电流体，其行为由流体力学方程与麦克斯韦方程组耦合决定。在推导平衡方程之前，我们首先回顾理想 MHD 的基本方程组：
- 动量方程：描述等离子体受力运动
$$ \rho \frac{d\mathbf{v}}{dt} = -\nabla p + \mathbf{J} \times \mathbf{B} + \rho \mathbf{g} $$ - 麦克斯韦-安培定律（忽略位移电流）：
$$ \mu_0 \mathbf{J} = \nabla \times \mathbf{B} $$ - 高斯磁定律：
$$ \nabla \cdot \mathbf{B} = 0 $$
其中，$\rho$ 为流体质量密度，$\mathbf{v}$ 为流体速度，$p$ 为等离子体热压强，$\mathbf{J}$ 为电流密度，$\mathbf{B}$ 为磁感应强度，$\mu_0$ 为真空磁导率。在磁约束等离子体中，重力 $\rho \mathbf{g}$ 通常远小于电磁力，可以忽略不计。
平衡态假设
所谓“平衡态”，是指系统的宏观物理量不随时间发生变化，即 $\partial / \partial t = 0$。同时，在静态平衡假设下，等离子体的宏观流动速度为零或小到可以忽略，即 $\mathbf{v} = 0$。
将上述假设代入动量方程，等号左侧的惯性项 $\rho \frac{d\mathbf{v}}{dt}$ 变为零。此时，动量方程退化为一个力平衡方程，表明等离子体内部的压力梯度被洛伦兹力完全抵消。
平衡方程推导
基于上述平衡态假设，动量方程简化为：
$$ -\nabla p + \mathbf{J} \times \mathbf{B} = 0 $$
为了更清晰地分析磁场与压强的关系，我们需要将洛伦兹力 $\mathbf{J} \times \mathbf{B}$ 用磁场 $\mathbf{B}$ 显式地表示出来。利用安培定律将 $\mathbf{J}$ 替换，可得：
$$ \mathbf{J} \times \mathbf{B} = \frac{1}{\mu_0} (\nabla \times \mathbf{B}) \times \mathbf{B} $$
利用矢量恒等式 $(\nabla \times \mathbf{B}) \times \mathbf{B} = (\mathbf{B} \cdot \nabla)\mathbf{B} - \nabla \left( \frac{B^2}{2} \right)$，洛伦兹力可以展开为两项：
$$ \mathbf{J} \times \mathbf{B} = \frac{1}{\mu_0} (\mathbf{B} \cdot \nabla)\mathbf{B} - \nabla \left( \frac{B^2}{2\mu_0} \right) $$
将此展开式代回平衡方程中，整理得到：
$$ \nabla p + \nabla \left( \frac{B^2}{2\mu_0} \right) = \frac{1}{\mu_0} (\mathbf{B} \cdot \nabla)\mathbf{B} $$
这就是磁流体平衡方程的一般矢量形式。等号左侧代表等离子体热压强梯度与磁压强梯度的叠加，等号右侧代表磁张力。
物理意义解析
从推导出的平衡方程可以看出，磁场对等离子体施加了两种不同性质的力：
- 磁压强：由 $-\nabla \left( \frac{B^2}{2\mu_0} \right)$ 项体现。磁压强 $p_B = \frac{B^2}{2\mu_0}$ 具有与流体热压强相同的量纲，它表现为磁场对自身产生的向内挤压的趋势。磁压强抵抗外部磁场进入等离子体内部。
- 磁张力：由 $\frac{1}{\mu_0} (\mathbf{B} \cdot \nabla)\mathbf{B}$ 项体现。类似于绷紧的橡皮筋产生的恢复力，磁张力沿着磁力线方向作用，抵抗磁力线的弯曲。
在平衡态下，等离子体的热压强与磁压强、磁张力达到力学平衡。为了衡量磁场约束等离子体热压强的效率，通常引入**比压（Beta, $\beta$）**参数：
$$ \beta = \frac{p}{B^2/(2\mu_0)} $$
$\beta$ 值越高，说明等离子体的热压强在总压强中占比越大，磁场的利用率越高；但 $\beta$ 值过高往往会导致磁流体不稳定性的发生。
一维圆柱等离子体平衡示例
为了更直观地理解平衡方程，我们考察最简单的一维圆柱等离子体平衡（如直线箍缩放电）。假设系统具有柱对称性，物理量仅依赖于半径 $r$，即 $\partial / \partial \theta = \partial / \partial z = 0$。
在柱坐标系 $(r, \theta, z)$ 下，设磁场仅有轴向分量 $B_z(r)$ 和角向分量 $B_\theta(r)$，即 $\mathbf{B} = B_\theta \hat{\theta} + B_z \hat{z}$。由于对称性，磁力线无弯曲，磁张力项 $(\mathbf{B} \cdot \nabla)\mathbf{B}$ 在径向的投影为零。
此时，平衡方程退化为压强梯度与磁压强梯度在径向的平衡：
$$ \frac{dp}{dr} + \frac{d}{dr}\left( \frac{B_\theta^2 + B_z^2}{2\mu_0} \right) = 0 $$
对上式从 $0$ 到等离子体边界 $a$ 进行积分，并假设边界处热压强 $p(a) = 0$，边界外侧角向磁场 $B_\theta(a) = \mu_0 I / (2\pi a)$（$I$ 为等离子体总电流），轴向磁场 $B_z(a) = B_0$，可得：
$$ \langle p \rangle = \frac{B_\theta^2(a)}{2\mu_0} + \frac{B_0^2 - B_z^2(0)}{2\mu_0} $$
其中 $\langle p \rangle$ 为平均压强。该结果表明，等离子体内部的压强由两部分磁压强差来支撑：一是表面角向磁压强（箍缩效应），二是内部与外部轴向磁压强之差（磁塞效应）。
总结
磁流体平衡方程推导揭示了磁约束等离子体实现宏观静止的力学本质：等离子体热压强梯度必须被磁压强梯度与磁张力完全平衡。通过将洛伦兹力分解为磁压强与磁张力，我们能够更深刻地理解磁场约束等离子体的物理机制。无论是简单的圆柱构型，还是复杂的环构型（如托卡马克），磁流体平衡方程都是进行装置设计、参数优化以及稳定性分析不可或缺的理论基石。