В магнитогидродинамике (МГД) плазма рассматривается как проводящая жидкость, взаимодействующая с магнитным полем. Основные уравнения МГД включают:
$$ \rho \left(\frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v}\cdot\nabla)\mathbf{v}\right) = -\nabla p + \mathbf{J}\times \mathbf{B} + \rho \mathbf{g} + \mathbf{F}_{\text{дисс}}, $$
где ρ — плотность плазмы, v — скорость, p — давление, J — ток плотности, B — магнитное поле, g — ускорение силы тяжести, Fдисс — диссипативные силы (вязкость, резистивность).
$$ \frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{v} \times \mathbf{B}) - \nabla \times (\eta \mathbf{J}), $$
где η — магнитная диффузия ($\eta = \frac{1}{\sigma \mu_0}$), σ — электропроводность, μ0 — магнитная проницаемость вакуума.
$$ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0, $$
$$ \frac{\partial}{\partial t}\left(\frac{\rho v^2}{2} + \frac{B^2}{2\mu_0} + \epsilon\right) + \nabla \cdot \mathbf{S} = 0, $$
где ϵ — внутренняя энергия, S — поток энергии (включая магнитный и кинетический вклад).
МГД волны описывают малые возмущения в плазме под действием магнитного поля. Их свойства зависят от ориентации волнового вектора относительно магнитного поля и параметров среды.
$$ v_A = \frac{B_0}{\sqrt{\mu_0 \rho}}, $$
где B0 — величина основного магнитного поля.
Особенности:
$$ v_{f,s}^2 = \frac{1}{2} \left(v_A^2 + c_s^2 \pm \sqrt{(v_A^2 + c_s^2)^2 - 4 v_A^2 c_s^2 \cos^2\theta}\right), $$
где $c_s = \sqrt{\gamma p/\rho}$ — скорость звука в плазме, θ — угол между магнитным полем и направлением волны, vf — скорость fast-волны, vs — скорость slow-волны.
Неустойчивость возникает, когда возмущения в плазме растут с течением времени, а не затухают. Они играют ключевую роль в термоядерных реакторах и астрофизике.
Неустойчивость Релея–Тейлора: Возникает на границе двух слоев с разной плотностью при наличии ускорения, направленного от тяжелого слоя к легкому.
$$ \gamma = \sqrt{g k \frac{\rho_2 - \rho_1}{\rho_2 + \rho_1}}, $$
где k — волновое число, ρ2 > ρ1, g — ускорение.
Неустойчивость Кельвина–Гельмгольца: Возникает при сдвиге скоростей двух слоев плазмы.
$$ \gamma = k \frac{\Delta v (\rho_1 \rho_2)^{1/2}}{\rho_1 + \rho_2}. $$
Токовая неустойчивость (длинные токовые столбы, токи в токамаке): Связана с накоплением энергии магнитного поля, которая может высвободиться при нарушении равновесия.
Неустойчивость Тейлора–Кулона (sausage и kink): Волны деформируют токовую трубку: