Exact and Locally Implicit Source Term Solvers for Multifluid-Maxwell Systems

Abstract

Recently, a family of models that couple multifluid systems to the full Maxwell equations draw a lot of attention in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they naturally retain non-ideal effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, etc. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid kinetic processes like plasma oscillation and electron cyclotron motion, even when they are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic forms solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with $S$ species, the locally implict algorithm involves inverting a local $3S+3$ matrix only, thus is very efficient. The performance can be further elevated by using the direct update formulas to skip null calculations. Benchmarks illustrated the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and large-scale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.

References in corpus (20)

• Athena: A New Code for Astrophysical MHD
• BOUT++: a framework for parallel plasma fluid simulations
• Comparison of multi-fluid moment models with Particle-in-Cell simulations of collisionless magnetic reconnection
• Global Ten-Moment Multifluid Simulations of the Solar Wind Interaction with Mercury: From the Planetary Conducting Core to the Dynamic Magnetosphere
• Continuum Kinetic and Multi-Fluid Simulations of Classical Sheaths
• Electron Physics in 3D Two-Fluid Ten-Moment Modeling of Ganymede's Magnetosphere
• A High-Order Relativistic Two-Fluid Electrodynamic Scheme with Consistent Reconstruction of Electromagnetic Fields and a Multidimensional Riemann Solver for Electromagnetism
• Entropy Stable Numerical Schemes for Two-Fluid Plasma Equations
• The island coalescence problem: scaling of reconnection in extended fluid models including higher-order moments
• New closures for more precise modeling of Landau damping in the fluid framework
• Relativistic Two-fluid Simulations of Guide Field Magnetic Reconnection
• Temperature gradient driven heat flux closure in fluid simulations of collisionless reconnection
• A Six-moment Multi-fluid Plasma Model
• Divergence-free Approximate Riemann Solver for the Quasi-neutral Two-fluid Plasma Model
• Magnetohydrodynamic motion of a two-fluid plasma
• Drift instabilities in thin current sheets using a two-fluid model with pressure tensor effects
• A new framework for magnetohydrodynamic simulations with anisotropic pressure
• Gaussian-Moment Relaxation Closures for Verifiable Numerical Simulation of Fast Magnetic Reconnection in Plasma
• Using the maximum entropy distribution to describe electrons in reconnecting current sheets
• An extended MHD study of the 16 October 2015 MMS diffusion region crossing

Cited by in corpus (11)

• Global Ten-Moment Multifluid Simulations of the Solar Wind Interaction with Mercury: From the Planetary Conducting Core to the Dynamic Magnetosphere
• Data-driven, multi-moment fluid modeling of Landau damping
• Data-driven modeling of Landau damping by physics-informed neural networks
• Machine-learning heat flux closure for multi-moment fluid modeling of nonlinear Landau damping
• Magnetic Field Reconstruction for a Realistic Multi-Point, Multi-Scale Spacecraft Observatory
• Numerical Study of Magnetic Island Coalescence Using Magnetohydrodynamics With Adaptively Embedded Particle-In-Cell Model
• Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling
• Fluid Simulations of Three-Dimensional Reconnection that Capture the Lower-Hybrid Drift Instability
• Experimental study of Alfvén wave reflection from an Alfvén-speed gradient relevant to the solar coronal holes
• A Parallel-Kinetic-Perpendicular-Moment Model for Magnetized Plasmas
• Surrogate Modeling of Landau Damping with Deep Operator Networks