The Equations of Ideal Magnetohydrodynamics (MHD)


The ideal MHD equations are in terms of the eight variables in addition to time t which are conserved in a volume: 
density $\rho$, momentum, $\rho \mathbf u$, magnetic flux density $\mathbf B$ and energy density, $E$.
These equations expressed in conservation form are as follows:


Conservation of Mass:
\begin{equation}
\frac{\partial \rho}{\partial t}+\boldsymbol{\nabla}\cdot(\rho \mathbf u)=0
\end{equation}
Conservation of Momentum:
\begin{equation}
\frac{\partial}{\partial t}\left(
\rho \mathbf u
\right)
+\boldsymbol{\nabla}
\cdot
\left[
\rho \mathbf u \otimes \mathbf u
+\left(
p_{B}
\right)
\mathbf{{\overline {\overline I}}}
-\mathbf{B} \otimes \mathbf{B}
\right]
=0
\end{equation}
Conservation of Flux:
\begin{equation}
\frac{\partial \mathbf{B} }{\partial t}+
\boldsymbol{\nabla}
\cdot
(
\mathbf{u} \otimes \mathbf{B}-
\mathbf{B} \otimes \mathbf{u})
= 0
\end{equation}
Conservation of Energy:
\begin{equation}
\frac{\partial E}{\partial t}
+\boldsymbol{\nabla}
\cdot
\left[
\left( E + p_{B}\right) \mathbf{u}
-(\mathbf{u}\cdot\mathbf{B})\mathbf{B}
\right]  = 0
\end{equation}
where $E$ is the total energy, kinetic, internal and magnetic 
\begin{equation}
E =
\frac{1}{2}\rho |\mathbf{u}|^2 + 
\frac{p}{\gamma-1}+
\frac{1}{2} |\mathbf{B}|^2
\end{equation}
and $p_{B}$ is the sum of thermal and magnetic pressure:
\begin{equation}
p_{B} = p + \frac{1}{2}|\mathbf{B}|^2
\end{equation}
The units are chosen so that $\mathbf B$ absorbs a factor of $1/\sqrt {4 \pi}$.
The adiabatic index is $\gamma = 5/3$ for a monatomic gas throughout the simulations.
$\mathbf{I}$ is the identity tensor.
There are eight equations and nine variables. In order to close the system one more equation is required, the equation of state.
The equation of state is chosen to be the ideal gas equation.


```python

```

80871b6a724a04b566d628f4eaf4653a987166b6
ipynb
Jupyter Notebook
MHD Equations.ipynb
mhdproject/mhdvanleer
63d3a3f94ae2265dd78bc05b23fa3a6cbc950309
MIT
2018-12-31T17:40:48.000Z
2021-06-17T22:10:25.000Z
MHD Equations.ipynb
mhdproject/mhdvanleer
63d3a3f94ae2265dd78bc05b23fa3a6cbc950309
MIT
2017-07-12T18:45:00.000Z
2017-07-23T10:05:03.000Z
MHD Equations.ipynb
garethcmurphy/mhdvanleer
63d3a3f94ae2265dd78bc05b23fa3a6cbc950309
MIT
2017-11-03T19:47:55.000Z
2017-11-03T19:47:55.000Z
Qwen/Qwen-72B
1. YES
2. YES