Related papers: A structure-preserving particle discretisation for…
This paper contributes new insights into discretizing Coulomb collisions in kinetic plasma models. Building on the previous works [Carrillo et al. J. Comp. Phys. X 7:100066 (2020), Hirvijoki and Burby Phys. Plasmas 27(8):082307 (2020)], I…
This paper proposes a novel numerical integrator for modeling multispecies Coulomb collisions in kinetic plasmas. The proposed scheme provides an energy-, momentum-, and positivity-preserving particle discretization of the nonlinear Landau…
Collisions are crucial in governing particle and energy transport in plasmas confined in a magnetic mirror trap. Modern gyrokinetic codes model transport in magnetic mirrors, but some utilize approximate model collision operators. This…
Spectral approximation based on Hermite-Fourier expansion of the Vlasov-Poisson model for a collisionless plasma in the electro-static limit is provided, by including high-order artificial collision operators of Lenard-Bernstein type. These…
Generalising the work of Lenard and Bernstein, we introduce a new, fully relativistic model to describe collisional plasmas. Like the Fokker-Planck operator, this equation represents velocity diffusion and conserves particle number.…
Energy conserving particle-in-cell schemes are constructed for a class of reduced relativistic Vlasov--Maxwell equations of laser-plasma interaction. Discrete Poisson equation is also satisfied by the numerical solution. Specifically,…
The relaxation of a weakly collisional plasma, which is of fundamental interest to laboratory and astrophysical plasmas, can be described by the Boltzmann-Poisson equations with the Lenard-Bernstein collision operator. We perform a…
We present a novel structure-preserving framework for solving the Vlasov-Poisson-Landau system of equations using a particle in cell (PIC) discretization combined with discrete gradient time integrators. The Vlasov-Poisson-Landau system is…
We present a novel family of particle discretisation methods for the nonlinear Landau collision operator. We exploit the metriplectic structure underlying the Vlasov-Maxwell-Landau system in order to obtain disretisation schemes that…
The micro-macro (mM) decomposition approach is considered for the numerical solution of the Vlasov--Poisson--Lenard--Bernstein (VPLB) system, which is relevant for plasma physics applications. In the mM approach, the kinetic distribution…
A variational method is used to derive a self-consistent macro-particle model for relativistic electromagnetic kinetic plasma simulations. Extending earlier work [E. G. Evstatiev and B. A. Shadwick, J. Comput. Phys., vol. 245, pp. 376-398,…
We propose an explicit particle method for the Vlasov-Fokker-Planck equation that conserves energy at the fully discrete level. The method features two key components: a deterministic and conservative particle discretization for the…
The Landau collision integral is often considered the gold standard in the context of kinetic plasma simulation due to its conservative properties, despite challenges involved in its discretization. The primary challenge when implementing…
Ionization is critical in the formation and evolution of plasma dynamics; collisional ionization, in particular, is an often overlooked source of electrons when dealing with laser-plasma interactions. Ionization plays a crucial role in…
We present a generalized, data-driven collisional operator for one-component plasmas, learned from molecular dynamics simulations, to extend the collisional kinetic model beyond the weakly coupled regime. The proposed operator features an…
The Balescu-Lenard equation from plasma physics is widely considered to include a highly accurate correction to Landau's fundamental collision operator. Yet so far it has seen very little mathematical study. We perform an extensive…
This research note documents new developments regarding finite-element discretizations of the relativistic Beliaev-Budker Coulomb collision operator and the nonrelativistic Landau operator. Where energy conservation in a finite-element…
Ideally, binary-collision algorithms conserve kinetic momentum and energy. In practice, the finite size of collision cells and the finite difference in the particle locations affect the conservation properties. In the present work, we…
We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by…
In this paper, we propose and implement a structure-preserving stochastic particle method for the Landau equation. The method is based on a particle system for the Landau equation, where pairwise grazing collisions are modeled as diffusion…