Lattice Boltzmann approaches to magnetohydrodynamics and electromagnetism

\newcommand{\Jv}{\mathbf{J}} \newcommand{\uv}{\mathbf{u}} \newcommand{\Bv}{\mathbf{B}} \newcommand{\Ev}{\mathbf{E}} \newcommand{\gv}{\mathbf{g}} We present a lattice Boltzmann approach for magnetohydrodynamics and electromagnetism that expresses the magnetic field using a discrete set of vector distribution functions $\gv_i$. The $\gv_i$ were first postulated to evolve according to a vector Boltzmann equation of the form $$ \partial_t \gv_i + \xi_i \cdot \nabla \gv_i = - \frac{1}{\tau} \left( \gv_i - \gv_i^{(0)} \right), $$ where the $\xi_i$ are a discrete set of velocities. The right hand side relaxes the $\gv_i$ towards some specified functions $\gv_i^{(0)}$ of the fluid velocity $\uv$, and of the macroscopic magnetic field given by $\Bv = \sum_i \gv_i$. Slowly varying solutions obey the equations of resistive magnetohydrodynamics. This lattice Boltzmann formulation has been used in large-scale (up to $1800^3$ resolution) simulations of magnetohydrodynamic turbulence. However, this is only the simplest form of Ohm's law. We may simulate more realistic extended forms of Ohm's law using more complex collision operators. A current-dependent relaxation time yields a current-dependent resistivity $\eta(|\nabla\times\Bv|)$, as used to model ``anomalous'' resistivity created by small-scale plasma processes. Using a \textit{hydrodynamic} matrix collision operator that depends upon the magnetic field $\Bv$, we may simulate Braginskii's magnetohydrodynamics, in which the viscosity for strains parallel to the magnetic field lines is much larger than the viscosity for strains in perpendicular directions. Changing the collision operator again, from the above vector Boltzmann equation we may derive the full set of Maxwell's equations, including the displacement current, and Ohm's law, $$ - \frac{1}{c^2} \partial_t \mathbf{E}+ \nabla \times \Bv = \mu_o \Jv, \quad \Jv = \sigma ( \mathbf{E} + \uv \times \Bv). $$ The original lattice Boltzmann scheme was designed to reproduce resistive magnetohydrodynamics in the non-relativistic limit. However, the kinetic formulation requires a system of first order partial differential equations with collision terms. This system coincides with the full set of Maxwell's equations and Ohm's law, so we capture a much wider range of electromagnetic phenomena, including electromagnetic waves.