Going Beyond the MHD Approximation: Physics-Based Numerical Solution of the CGL Equations

Abstract

We present a new numerical model for solving the Chew-Goldberger-Low system of equations describing a bi-Maxwellian plasma in a magnetic field. Heliospheric and geospace environments are often observed to be in an anisotropic state with distinctly different parallel and perpendicular pressure components. The CGL system represents the simplest leading order correction to the common isotropic MHD model that still allows to incorporate the latter's most desirable features. However, the CGL system presents several numerical challenges: the system is not in conservation form, the source terms are stiff, and unlike MHD it is prone to a loss of hyperbolicity if the parallel and perpendicular pressures become too different. The usual cure is to bring the parallel and perpendicular pressures closer to one another; but that has usually been done in an ad hoc manner. We present a physics-informed method of pressure relaxation based on the idea of pitch-angle scattering that keeps the numerical system hyperbolic and naturally leads to zero anisotropy in the limit of very large plasma beta. Numerical codes based on the CGL equations can, therefore, be made to function robustly for any magnetic field strength, including the limit where the magnetic field approaches zero. The capabilities of our new algorithm are demonstrated using several stringent test problems that provide a comparison of the CGL equations in the weakly and strongly collisional limits. This includes a test problem that mimics interaction of a shock with a magnetospheric environment in 2D.

Accepted in "The Astrophysical Journal (ApJ)"

References in corpus (14)

• An Unsplit Godunov Method for Ideal MHD via Constrained Transport
• MHD Turbulence Revisited
• Second Order Accurate Schemes for Magnetohydrodynamics With Divergence-Free Reconstruction
• Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
• Multidimensional HLLE Riemann solver; Application to Euler and Magnetohydrodynamic Flows
• The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas
• Shearing Box Simulations of the MRI in a Collisionless Plasma
• A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows
• Divergence-Free Reconstruction of Magnetic Fields and WENO Schemes for Magnetohydrodynamics
• The Distribution of Pressures in a Supernova-Driven Interstellar Medium. I. Magnetized Medium
• Magnetic field amplification and evolution in turbulent collisionless MHD: an application to the ICM
• A new framework for magnetohydrodynamic simulations with anisotropic pressure
• Limits on the ions temperature anisotropy in turbulent intracluster medium
• Gas Kinetic Schemes for Solving the Magnetohydrodynamic Equations with Pressure Anisotropy

Cited by in corpus (2)

• Two-Dimensional Kelvin-Helmholtz Instability with Anisotropic Pressure
• An Alternative Finite Difference WENO-like Scheme with Physical Constraint Preservation for Divergence-Preserving Hyperbolic Systems