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.

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

A Mesh-Scale Monte Carlo Method for Modeling Triboelectric Charging of Chemically Identical Particles

An efficient numerical method for charged particle transport based on hybrid collision model and machine learning

Modeling of the non-Maxwellian response of DT plasmas to alpha particle transport in inertial confinement fusion (ICF) hotspot

Compensated Monte Carlo Collision Model for Particle-In-Cell Simulation in High-Pressure Plasmas

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

Application of the view factor model on the particle-in-cell and Monte Carlo collision code.

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

A Machine-Learning Augmented Mixed Method for Simulating Α Particle Transport in Inertial Confinement Fusion

A COLLISION MODEL IN PLASMA PARTICLE SIMULATIONS

Monte Carlo method for investigating population kinetics in non-local thermodynamic equilibrium plasmas

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

High-efficient Particle-in-cell/monte Carlo Model for Complex Solution Domain Andsimulation of Anode Layer Ion Source

A general-purpose Monte Carlo particle transport code based on inverse transform sampling for radiotherapy dose calculation

Multiscale simulation of neutral particle flows in the plasma edge

A simplified positron elastic scattering model for Monte Carlo calculations

Monte Carlo Method for Calculating Particle Radiation Detector Correction Factors

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

A Monte Carlo Simulation of Ion Transport in ECR Microwave Plasma

Numerical Simulation of Plasma and Electromagnetic Characteristics for an Ion Thruster Discharge Chamber