Deepak Bhoriya,Dinshaw S. Balsara,Vladimir Florinski,Harish Kumar
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.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to more accurately simulate the behavior of plasma in the space plasma environment, especially in the case of plasma pressure anisotropy. The traditional magnetohydrodynamics (MHD) model assumes that the plasma pressure is isotropic, that is, the pressure is the same in all directions, which is reasonable in many cases. However, in the actual space plasma environment, such as the solar wind, the Earth's magnetosphere, and the magnetospheres of other planets, the plasma pressure is often anisotropic, that is, the pressure \( p_{\parallel} \) parallel to the magnetic field direction and the pressure \( p_{\perp} \) perpendicular to the magnetic field direction are not equal. The existence of this anisotropic pressure has an important impact on the dynamic behavior of the plasma, but the traditional MHD model cannot accurately describe these phenomena.
To overcome this limitation, Chew, Goldberger, and Low (CGL) proposed an extended MHD model, which allows the plasma to have an anisotropic pressure tensor. The CGL model can better describe the behavior of the plasma in the magnetic field by introducing two independent pressure variables \( p_{\parallel} \) and \( p_{\perp} \). However, the CGL model also brings some numerical challenges, including:
1. **Non - conservation form**: The CGL equations are not in conservation form, which makes it more difficult to use standard numerical methods (such as the finite - volume method).
2. **Stiff source terms**: The CGL equations contain stiff source terms. When the pressure anisotropy is large, these source terms will cause the numerical solution to be unstable.
3. **Loss of hyperbolicity**: If the difference between \( p_{\parallel} \) and \( p_{\perp} \) is too large, the CGL equations may lose hyperbolicity, making numerical solution very difficult.
To solve these problems, the paper proposes a new numerical model. Based on physical principles, this model maintains the hyperbolicity of the numerical system by introducing a physical pressure relaxation method based on scattering - angle scattering, and naturally tends to zero anisotropy in the high - β plasma limit. This method can not only handle magnetic fields of any strength, including cases where the magnetic field is close to zero, but also operate robustly in the weak - collision and strong - collision limits.
In summary, the main goal of this paper is to develop a numerical model that can effectively simulate the behavior of plasma with anisotropic pressure, thereby improving the modeling accuracy of the space plasma environment.