Skangaroo

Reduced Magnetohydrodynamics

Reduced Magnetyhdrodynamics (RMHD) is a set of four equations derived from the resistive MHD equations. It reduces the free paramameters of the equations down to four independent scalar variables and thus provides closure of MHD. It is important for investigations on linear and nonlinear plasma dynamics in Tokamak geometry and numerical simulations. The model was first developped by Strauß in 1976.

.

Outline of derivation

Multiscale expansion

Time $t$ and space ${\vec {x}}$ are expanded on multiple scales as

$t={\frac {1}{\varepsilon ^{2}}}T+\tau$ and

${\vec {x}}={\vec {X}}+\varepsilon {\vec {\xi }}$

within a small parameter $\varepsilon$ where $\tau$ and ${\vec {\xi }}$ denote a fast varying scale and $T$ and ${\vec {X}}$ a slowly varying scale. The magnetic field ${\vec {B}}$, electric field ${\vec {E}}$, current density ${\vec {J}}$, bulk fluid velocity ${\vec {V}}$, plamsa pressure $p$ and plasma density $\rho$ from the MHD model are exanded in $\varepsilon$ similar to:

${\vec {V}}={\vec {V}}_{0}({\vec {X}},T)+\sum \limits _{i=1}^{\infty }\varepsilon ^{i}{\vec {V}}_{i}({\vec {X}},T,{\vec {\xi }},\tau )$

This practically means that an equilibrium solution of the MHD model only varies on the slow scale whereas perturbations to this equilibrium can vary on both fast and slow scales. In particular are all fast derivatives of order $\varepsilon ^{0}$ quantaties equal to zero.

Flute Approximation

It is assumed that fast spatial variations are only possible in the plane perpendicular to the equilibrium magnetic field:

${\vec {b}}_{0}\cdot {\frac {\partial }{\partial {\vec {\xi }}}}=0$

This is the only additional assumption used and currently referred to as Flute Approximation.

Closure of the equations

The electromagnetic fields are expressed by an electrostatic potential $\phi$ and a vectorpotential ${\vec {A}}$

${\vec {B}}=\nabla \times {\vec {A}}$ and

${\vec {E}}=-\nabla \phi +{\frac {\partial }{\partial t}}{\vec {A}}$

Then all expanded quantaties are inserted into the resistive MHD equations. Since every order in epsilon is assumed to be an order of magnitude smaller than the previous ones, all terms of equal order in $\varepsilon$ have to balance each other. By going through the equations order by order in $\varepsilon$ some quantaties are found to be equal to zero or directly related to others. As independent variables remain the four following scalars:

• $P_{1}$, the lowest order pressure pertubation
• $\phi _{2}$, the second order electrostatic potential perturbation
• $\psi _{2}$, the component parallel to the magnetic field of the second order vectorpotential
• $V_{1\parallel }$, the lowest order velocity perturbation parallel to the magnetic Field

The evolution of these quantaties is governed by the following set of four coupled differential equations.

RMHD equations

${\begin{aligned}&{\frac {\partial \psi _{2}}{\partial \tau }}+\left({\vec {b}}_{0}\cdot {\frac {\partial }{\partial {\vec {X}}}}-{\frac {{\vec {b}}_{0}}{B_{0}}}\cdot {\frac {\partial \psi _{2}}{\partial {\vec {\xi }}}}\right)\phi _{2}={\frac {\eta }{\mu _{0}}}{\frac {\partial ^{2}\psi _{2}^{2}}{\partial {\vec {\xi }}^{2}}}&(1)\\&\rho _{0}\left({\frac {\partial }{\partial \tau }}+{\vec {V}}_{E}\cdot {\frac {\partial }{\partial {\vec {\xi }}}}\right)V_{1\parallel }=-\left({\vec {b}}_{0}\cdot {\frac {\partial }{\partial {\vec {X}}}}-{\frac {{\vec {b}}_{0}}{B_{0}}}\cdot {\frac {\partial \psi _{2}}{\partial {\vec {\xi }}}}\times {\frac {\partial }{\partial {\vec {\xi }}}}\right)P_{1}+{\frac {{\vec {b}}_{0}}{B_{0}}}\cdot {\frac {\partial \psi _{2}}{\partial {\vec {x}}}}\times {\frac {\partial P_{0}}{\partial {\vec {X}}}}+\mu _{\perp }{\frac {\partial ^{2}V_{1\parallel }}{\partial {\vec {\xi }}^{2}}}&(2)\\&\rho _{0}\left({\frac {\partial }{\partial \tau }}+{\vec {V}}_{E}\cdot {\frac {\partial }{\partial {\vec {\xi }}}}\right){\frac {\partial ^{2}\phi _{2}}{\partial {\vec {\xi }}^{2}}}=2B_{0}{\vec {b}}_{0}\times {\vec {\kappa }}_{0}\cdot {\frac {\partial P_{1}}{\partial {\vec {\xi }}_{\perp }}}-{\frac {B_{0}^{2}}{\mu _{0}}}\left(b_{0}\cdot {\frac {\partial }{\partial {\vec {X}}}}-{\frac {{\vec {b}}_{0}}{B_{0}}}\cdot {\frac {\partial \psi _{2}}{\partial {\vec {\xi }}}}\times {\frac {\partial }{\partial {\vec {\xi }}}}\right){\frac {\partial ^{2}\psi _{2}}{\partial {\vec {\xi }}^{2}}}+\mu _{\perp }{\frac {\partial ^{4}\phi _{2}}{\partial {\vec {\xi }}^{4}}}&(3)\\&\left(1+{\frac {5}{3}}{\frac {\mu _{0}P_{0}}{B_{0}^{2}}}\right)\left({\frac {\partial P_{1}}{\partial \tau }}+{\vec {V}}_{E}\cdot {\frac {\partial P_{1}}{\partial {\vec {\xi }}}}+{\vec {V}}_{E}\cdot {\frac {\partial P_{0}}{\partial {\vec {X}}}}\right)={\frac {5}{3}}P_{0}\left[{\frac {2}{B_{0}}}{\vec {b}}_{0}\cdot \kappa _{0}\times {\frac {\partial \phi _{2}}{\partial {\vec {X}}}}-\left({\vec {b}}_{0}\cdot {\frac {\partial }{\partial {\vec {X}}}}-{\frac {{\vec {b}}_{0}}{B_{0}}}\cdot {\frac {\partial \psi _{2}}{\partial \xi }}\times {\frac {\partial }{\partial {\vec {\xi }}}}\right)V_{1\parallel }+V_{1\parallel }{\vec {b}}_{0}\cdot {\frac {\partial \ln B_{0}}{\partial {\vec {X}}}}+{\frac {\eta }{B_{0}^{2}}}{\frac {\partial ^{2}P_{1}}{\partial {\vec {\xi }}^{2}}}\right]+\chi {\frac {\partial ^{2}P_{1}}{\partial {\vec {\xi }}^{2}}}&(4)\end{aligned}}$

Here it is assumed that an equilibrium solution of the MHD model is known and thus all quantaties denoted by a $0$ subscript. ${\vec {\kappa }}_{0}=({\vec {b}}_{0}\cdot \partial /\partial {\vec {X}}){\vec {b}}_{0}$ is the curvature vector of the 0th order magnetic field. $\mu _{0}$ is the magnetic constant and $\eta$ the plasma's resistivity which is assumed to be small. $\mu _{\perp }$ is the leading order viscosity and $\chi$ the heat diffusivity. ${\vec {V}}_{E}$ is the electric drift. The particle density $\rho _{0}$ in these equations is a passive scalar and thus can be assumed constant.

In the above set of equations only fast time derivatives occur. Hence they solely describe fast dynamics along the slowly varying equilbrium magnetic field. Thus Magnetosonic waves, which constraine the computational speed of the original MHD model are eliminated.

Low Beta RMHD

In the limit of low Beta the equations simplify further. By normalizing the prefactor $P_{0}$ on the right hand side of equation (4) with a characteristic magnetic pressure it gets the plasma beta as a prefactor. By assuming an ordering of $\beta \sim \varepsilon$ everything but the diffusion term on the rhs can be ignored. Then equation (4) simplifies to

$\left({\frac {\partial }{\partial \tau }}+{\vec {V}}_{E}\cdot {\frac {\partial }{\partial {\vec {\xi }}}}\right)P_{1}+{\vec {V}}_{E}\cdot {\frac {\partial P_{0}}{\partial {\vec {X}}}}=\chi {\frac {\partial ^{2}P_{1}}{\partial {\vec {\xi }}^{2}}}\quad (4b)$.

Now only equation (2) contains the velocity perturbation and it decouples from the rest of equations. Thus low beta RMHD contains only the equations (1), (3) and (4b) for the three variables P₁, Φ₂ and ψ₂.

References

https://web2.ph.utexas.edu/~morrison/84POF_morrison.pdf