[Paper Review] A Unified Gas Kinetic Scheme for Multi-scale Plasma Transport

Chang Liu, Kun Xu|arXiv (Cornell University)|Sep 17, 2016

Ionosphere and magnetosphere dynamics3 citations

TL;DR

This paper proposes a unified gas kinetic scheme (UGKS) for multi-scale plasma transport by combining Vlasov-BGK equations for electrons and ions with Maxwell's equations, enabling a seamless transition from kinetic to hydrodynamic plasma regimes. The scheme achieves a unified treatment of collisionless to highly collisional plasma dynamics through a time-evolving, scale-dependent flux reconstruction that couples particle transport, collisions, and electromagnetic fields without splitting.

ABSTRACT

A unified gas kinetic scheme (UGKS) for multi-scale and multi-component plasma transport is constructed. The current scheme is a direct modeling method, where the time evolution solutions from the Vlasov-BGK equations for both electron and ion, and the Maxwell equations are used to construct the scale-dependent plasma simulation. As a result, based on the modeling scales of mesh size and time step, the discretized governing equations for the whole plasma regimes are obtained. The UGKS takes into account the electron inertia, full electromagnetic field equations, and separate electron and ion evolution. The physics recovered in UGKS ranges from the kinetic Vlasov equation to the hydrodynamic magnetohydrodynamic (MHD) equations, with a unified treatment in all scales from the collisionless particle transport to the hydrodynamic wave interactions. The UGKS presents a plasma description which is more general than the Vlasov equation in the kinetic scale and all kinds of MHD equations in the hydrodynamic scale, such as Hall, Resistive, and Ideal Magnetohydrodynamics (MHD). All above single scale equations become the subsets of the UGKS. The key dynamics in UGKS is the non-splitting treatment of particle collision, acceleration, and transport in the construction of the numerical flux at cell interface, and this flux is evaluated from the scale and time evolving integral solution of the Vlasov-BGK model. The evolution of electromagnetic field and the plasma flow are coupled in a semi-implicit manner. The UGKS is able to give a physically reliable solution for plasma evolution from the collisionless limit to the highly collisional one.

Motivation & Objective

• To develop a single computational framework that captures plasma behavior across all scales, from kinetic to hydrodynamic regimes.
• To address the limitations of existing methods that require separate models for different plasma regimes, such as Vlasov for kinetic and MHD for fluid-like behavior.
• To unify the treatment of electron and ion dynamics, including inertia and electromagnetic field evolution, in a consistent numerical scheme.
• To eliminate the need for separate solvers by constructing a unified scheme that reduces to known equations (e.g., ideal MHD, Hall MHD, resistive MHD) as special cases.
• To ensure physical consistency across scales by coupling particle transport, collisions, and electromagnetic field evolution in a non-splitting flux formulation.

Proposed method

• The scheme is based on direct modeling of the Vlasov-BGK equations for electrons and ions, with time evolution solutions used to construct numerical fluxes at cell interfaces.
• The flux is computed from the scale- and time-dependent integral solution of the Vlasov-BGK model, ensuring accurate treatment of particle transport, collisions, and acceleration in a unified manner.
• Maxwell’s equations are solved simultaneously with the plasma equations in a semi-implicit coupling strategy to maintain field-particle consistency.
• The method inherently captures the transition from collisionless to collisional regimes by adjusting the mesh size and time step, which determine the effective scale of the solution.
• The scheme avoids operator splitting by treating particle collisions, acceleration, and transport together in the flux evaluation, preserving physical accuracy across scales.
• The formulation naturally reduces to known equations—such as ideal, resistive, and Hall MHD—when the system is in the hydrodynamic limit, ensuring consistency with established models.

Experimental results

Research questions

• RQ1 Can a single numerical scheme accurately simulate plasma transport across all scales, from kinetic to hydrodynamic regimes, without requiring model switching?
• RQ2 How can particle collisions, electromagnetic fields, and transport be consistently coupled in a non-splitting manner across multiple scales?
• RQ3 To what extent does the UGKS recover known plasma models (e.g., ideal MHD, Hall MHD, resistive MHD) as limiting cases?
• RQ4 Can the scheme maintain physical accuracy in both collisionless and highly collisional plasma conditions using the same algorithmic framework?
• RQ5 What is the role of the time-evolving integral solution of the Vlasov-BGK model in enabling scale-adaptive plasma simulations?

Key findings

• The UGKS provides a unified description of plasma dynamics that encompasses the Vlasov equation in the kinetic regime and all standard MHD models (ideal, resistive, Hall) in the hydrodynamic regime.
• The scheme successfully captures the transition from collisionless to collisional plasma behavior by adapting to the local mesh and time step scales.
• The non-splitting treatment of particle transport, collisions, and electromagnetic field evolution ensures physical consistency across all scales.
• The method maintains accuracy and stability in both the kinetic and hydrodynamic limits, with no need for separate solvers or model switching.
• The UGKS naturally reduces to known equations such as ideal MHD and resistive MHD in appropriate asymptotic limits, validating its consistency with established plasma models.
• The semi-implicit coupling of electromagnetic fields and plasma flow enables stable and physically reliable simulations over a wide range of plasma conditions.