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.

Energy-Preserving Coupling of Explicit Particle-In-Cell with Monte Carlo Collisions

Numerical heating in particle-in-cell simulations with Monte Carlo binary collisions

On the Energy Conservation Electrostatic Particle-In-Cell/Monte Carlo Simulation: Benchmark and Application to the Radio Frequency Discharges

Stochastic Weighted Particle Control for Electrostatic Particle-in-cell Monte Carlo Collision Simulations in an Axisymmetric Coordinate System

Parallelization And Optimization Of Electrostatic Particle-In-Cell/Monte-Carlo Coupled Codes As Applied To Rf Discharges

An Asymptotic-Preserving and Energy-Conserving Particle-In-Cell Method for Vlasov-Maxwell Equations

Explicit energy-conserving modification of relativistic PIC method

Particle-in-cell Simulation of Transport and Energy Deposition of Intense Proton Beams in Solid-State Materials

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

An explicit, energy-conserving particle-in-cell scheme

Numerical Approach to Predict Particle Breakage in Dense Flows by Coupling Multiphase Particle-In-Cell and Monte Carlo Methods

Moment-preserving Monte-Carlo Coulomb collision method for particle codes

Development of Two Asymptotic-Preserving Particle-in-Cell Methods in Quasi-Neutral Plasma Simulation

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

Advances in the Implementation of the Exactly Energy Conserving Semi-Implicit (ECsim) Particle in Cell Method

Numerical Thermalization in One- and Two-Dimensional Particle-In-Cell Simulations with Monte-Carlo Collisions

Implicit Temporal Discretization and Exact Energy Conservation for Particle Methods Applied to the Poisson–Boltzmann Equation

Computational Speed-Up Techniques for a PIC-MCC Computer Model for Use in Modeling the Plasma in an Ion Engine Discharge Chamber

Exact local conservation of energy in fully implicit PIC algorithms

Numerical study on wave attenuation via 1D fully kinetic electromagnetic particle-in-cell simulations