A Deep Dive into the Distribution Function: Understanding Phase Space Dynamics with Continuum Vlasov-Maxwell Simulations

05/27/2020
by   James Juno, et al.
0

In collisionless and weakly collisional plasmas, the particle distribution function is a rich tapestry of the underlying physics. However, actually leveraging the particle distribution function to understand the dynamics of a weakly collisional plasma is challenging. The equation system of relevance, the Vlasov-Maxwell-Fokker-Planck (VM-FP) system of equations, is difficult to numerically integrate, and traditional methods such as the particle-in-cell method introduce counting noise into the distribution function. In this thesis, we present a new algorithm for the discretization of VM-FP system of equations for the study of plasmas in the kinetic regime. Using the discontinuous Galerkin (DG) finite element method for the spatial discretization and a third order strong-stability preserving Runge-Kutta for the time discretization, we obtain an accurate solution for the plasma's distribution function in space and time. We both prove the numerical method retains key physical properties of the VM-FP system, such as the conservation of energy and the second law of thermodynamics, and demonstrate these properties numerically. These results are contextualized in the history of the DG method. We discuss the importance of the algorithm being alias-free, a necessary condition for deriving stable DG schemes of kinetic equations so as to retain the implicit conservation relations embedded in the particle distribution function, and the computational favorable implementation using a modal, orthonormal basis in comparison to traditional DG methods applied in computational fluid dynamics. Finally, we demonstrate how the high fidelity representation of the distribution function, combined with novel diagnostics, permits detailed analysis of the energization mechanisms in fundamental plasma processes such as collisionless shocks.

READ FULL TEXT
03/12/2022

Energy conserving particle-in-cell methods for relativistic Vlasov–Maxwell equations of laser-plasma interaction

Energy conserving particle-in-cell schemes are constructed for a class o...
10/05/2021

Inference and De-Noising of Non-Gaussian Particle Distribution Functions: A Generative Modeling Approach

The particle-in-cell numerical method of plasma physics balances a trade...
08/08/2022

A Stable Mimetic Finite-Difference Method for Convection-Dominated Diffusion Equations

Convection-diffusion equations arise in a variety of applications such a...
07/25/2021

A new family of second order convergent weakly-compressible SPH schemes

Despite the many advances in the use of weakly-compressible smoothed par...
07/22/2021

Conservative DG Method for the Micro-Macro Decomposition of the Vlasov-Poisson-Lenard-Bernstein Model

The micro-macro (mM) decomposition approach is considered for the numeri...
04/19/2022

Globally time-reversible fluid simulations with smoothed particle hydrodynamics

This paper describes an energy-preserving and globally time-reversible c...
05/31/2019

Improving the resolution of CryoEM single particle analysis

We present a new 3D refinement method for CryoEM single particle analysi...