Abstract:Coulomb collisions in particle simulations for weakly coupled plasmas are modeled by the Landau-Fokker-Planck equation, which is typically solved by Monte-Carlo (MC) methods. One of the main disadvantages of MC is the timestep accuracy constraint {\nu}{\Delta}t << 1 to resolve the collision frequency {\nu}. The constraint becomes extremely stringent for self-collisions in the presence of high charge state species and for inter-species collisions with large mass disparities (such as present in Inertial Confinement Fusion hohlraums), rendering long-time-scale simulations prohibitively expensive or impractical. To overcome these difficulties, we explore a hybrid Maxwellian-MC (HMMC) model for particle simulations. Specifically, we devise a collisional algorithm that describes weakly collisional species with particles, and highly collisional species and fluid components with Maxwellians. We employ the Lemons method for particle-Maxwellian collisions, enhanced with a more careful treatment of low-relative-speed particles, and a five-moment model for Maxwellian-Maxwellian collisions. Particle-particle binary collisions are dealt with classic Takizuka-Abe MC, which we extend to accommodate arbitrary particle weights to deal with large density disparities without compromising conservation properties. HMMC is strictly conservative and significantly outperforms standard MC methods in situations with large mass disparities among species or large charge states, demonstrating orders of magnitude improvement in computational efficiency. We will substantiate the accuracy and performance of the proposed method with several examples of varying complexity, including both zero-dimensional relaxation and one-dimensional transport problems, the latter using a hybrid kinetic-ion/fluid-electron model.

Reducing the computational complexity of implicit schemes in the modeling of kinetic inelastic collisions in a partially ionized plasma

An Implicit Lattice Boltzmann Method for Simulations of Compressible Plasma Kinetics

Direct Implicit and Explicit Energy-Conserving Particle-in-Cell Methods for Modeling of Capacitively-Coupled Plasma Devices

Importance of accurate consideration of the electron inertia in hybrid-kinetic simulations of collisionless plasma turbulence: 1. The 2D limit

A fast implicit solver for semiconductor models in one space dimension

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker-Planck-Boltzmann Model

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

Krylov-based Adaptive-Rank Implicit Time Integrators for Stiff Problems with Application to Nonlinear Fokker-Planck Kinetic Models

Sparse Data Structures for Efficient State-to-State Kinetic Simulations

High-order Implicit Particle-in-cell Method for Plasma Simulations at Solid Densities.

High-order implicit solver in conservative formulation for tokamak plasma transport equations

Explicit energy-conserving modification of relativistic PIC method

A Fully Explicit Integrator for Modeling Astrophysical Reactive Flows

An efficient implicit boundary integral solver for the Vlasov-Maxwell system

High order finite volume schemes with IMEX time stepping for the Boltzmann model on unstructured meshes

Effects of Artificial Collisions, Filtering, and Nonlocal Closure Approaches on Hermite-based Vlasov-Poisson Simulations

Numerical Investigation of Non-equilibrium Electron Effects on the Collisional Ionization Rate in the Collisional-Radiative Model

Relativistic Vlasov-Maxwell modelling using finite volumes and adaptive mesh refinement

A fully implicit, scalable, conservative nonlinear relativistic Fokker-Planck 0D-2P solver for runaway electron

A multiscale hybrid Maxwellian-Monte-Carlo Coulomb collision algorithm for particle simulations