Related papers: Asymptotic-preserving deterministic particle metho…
In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker-Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source.…
In this paper, we introduce a novel variant of the CBO method that incorporates jumps according to an $\alpha$-stable stochastic process in a kinetic framework. This extension gives rise to nonlocal stochastic effects, which improve the…
We present an analytic solution of a simple set of equations that govern the expansion of boost-invariant plasmas of massless particles. These equations describe, approximately, the early time, collisionless regime, and the transition to…
Achieving large-scale kinetic modelling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications such as plasma etching, kinetic simulations are necessary to…
An added-mass partitioned (AMP) algorithm is described for solving fluid-structure interaction (FSI) problems coupling incompressible flows with thin elastic structures undergoing finite deformations. The new AMP scheme is fully…
We propose a novel non-compact, positivity-preserving scheme for linear non-divergence form parabolic equations. Based on the Feynman-Kac formula, the solution is expressed as a conditional expectation of an associated diffusion process.…
We present a novel flow-based kinetic approach, inspired by continuous normalizing flows, for plasma simulation that unifies the complementary strengths of direct Vlasov solvers and particle-based methods. By tracking the distribution…
This work explores the application of deep operator learning principles to a problem in statistical physics. Specifically, we consider the linear kinetic equation, consisting of a differential advection operator and an integral collision…
Differentiable programming, enabled by automatic differentiation (AD), provides a robust framework for gradient-based optimization in computational plasma physics. While optimization is often only used towards design, we demonstrate that it…
Introducing a reduced particle stiffness in discrete element method (DEM) allows for bigger time steps and therefore fewer total iterations in a simulation. Although this approach works well for dry non-adhesive particles, it has been shown…
In this paper, we propose and validate a two-species Multiscale model for a Poisson-Nernst-Planck (PNP) system, focusing on the correlated motion of positive and negative ions under the influence of a trap. Specifically, we aim to model…
We propose a low Mach number, Godunov-type finite volume scheme for the numerical solution of the compressible Euler equations of gas dynamics. The scheme combines Klein's non-stiff/stiff decomposition of the fluxes (J. Comput. Phys.…
The direct-forcing immersed boundary method (DF-IBM) algorithm previously developed by the authors is extended by coupling the Navier-Stokes equations with the Newton-Euler equations for rigid body dynamics within the DF-IBM framework. This…
Currently existing energy-stable parametric finite element methods for surface diffusion flow and other flows are usually limited to first-order accuracy in time. Designing a high-order algorithm for geometric flows that can also be…
The dynamics of electron-plasma waves are described at arbitrary collisionality by considering the full Coulomb collision operator. The description is based on a Hermite-Laguerre decomposition of the velocity dependence of the electron…
We develop deterministic particle schemes to solve non-local scalar conservation laws with congestion. We show that the discrete approximations converge to the unique entropy solution with an explicit rate of convergence under more general…
In this work, the asymptotic state of nonlinear Landau damping in one-dimensional plasma has been examined using a quasi-linear model and a second-order symplectic integrator. The dispersion relation of the plateau distribution function for…
Hamiltonian and action principle (HAP) formulations of plasma physics are reviewed for the purpose of explaining structure preserving numerical algorithms. Geometric structures associated with and emergent from HAP formulations are…
In this article, we design Asymptotic-Preserving Particle-In-Cell methods for the Vlasov-Maxwell system in the quasi-neutral limit, this limit being characterized by a Debye length negligible compared to the space scale of the problem.…
In this study, we present an optimal implicit algorithm specifically designed to accurately solve the multi-species nonlinear 0D-2V axisymmetric Fokker-Planck-Rosenbluth (FPR) collision equation while preserving mass, momentum, and energy.…