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.

A COLLISION MODEL IN PLASMA PARTICLE SIMULATIONS

Particle-resolved Direct Numerical Simulation of Collisions of Bidisperse Inertial Particles in a Homogeneous Isotropic Turbulence

A Coulomb Collision Model Based on Collision Cross Section and Its Effects on Kinetic Characteristics of Atmospheric Thermal Microplasmas

A Kinetic Model for Non-Local Thermodynamic Equilibrium Plasmas

Monte Carlo simulations of charged particle transport in plasmas: A fast single-scattering model

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

Particle simulation of collision dynamics for ion beam injection into a rarefied gas

Simulating Particle Collision Process Based on Monte Carlo Method

A tutorial overview on the angular scattering models of Electron-Neutral, Ion-Neutral, Neutral-Neutral, and Coulomb Collisions in Monte Carlo collision modeling on low temperature plasma

A Numerical Model For Ion Charge Distribution Of Plasmas In Collisional Radiative Steady State

Finite Element Modeling of the Transient Heat Conduction Between Colliding Particles

A Nonlinear, Conservative, Entropic Fokker-Planck Model for Multi-Species Collisions

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

Enhanced Collisional Losses from a Magnetic Mirror Using the Lenard-Bernstein Collision Operator

Collision System Size Scan of Collective Flows in Relativistic Heavy-Ion Collisions

A Fluid-Kinetic Hybrid Electron Model for Electromagnetic Simulations

Verification of a Monte Carlo binary collision model for simulating elastic and inelastic collisions in particle-in-cell simulations

Study of discrete-particle effects in a one-dimensional plasma simulation with the Krook type collision model

An Extended Model for Ion Charge State Distribution of Plasmas in Collisional Radiative Steady State

Improved numerical simulation model for nuclear reaction rate calculations in high-speed plasma collisions