Related papers: An Explicitly Solvable Energy-Conserving Algorithm…
Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions,…
The traditional explicit electrostatic momentum-conserving Particle-in-cell algorithm requires strict resolution of the electron Debye length to deliver numerical accuracy. The explicit electrostatic energy-conserving Particle-in-Cell…
In this paper, we introduce and discuss an exactly energy-conserving Particle-in-Cell method for arbitrary curvilinear coordinates. The flexibility provided by curvilinear coordinates enables the study of plasmas in complex-shaped domains…
We report a new particle in cell (PIC) method based on the semi-implicit approach. The novelty of the new method is that unlike any of its semi-implicit predecessors at the same time retains the explicit computational cycle and conserves…
We propose an energy-stable parametric finite element method (ES-PFEM) to discretize the motion of a closed curve under surface diffusion with an anisotropic surface energy $\gamma(\theta)$ -- anisotropic surface diffusion -- in two…
Based on the previously developed Energy Conserving Semi Implicit Method (ECsim) code, we present its cylindrical implementation, called ECsim-CYL, to be used for axially symmetric problems. The main motivation for the development of the…
We propose a novel Skew Gradient Embedding (SGE) framework for systematically reformulating thermodynamically consistent partial differential equation (PDE) models-capturing both reversible and irreversible processes-as generalized gradient…
The Energy Conserving semi-implicit method (ECsim), presented by Lapenta in 2017, is a Particle in Cell (PIC) algorithm for the simulation of plasmas. Energy conservation is achieved within a semi-implicit formulation that does not require…
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…
Coulomb collision is a fundamental diffusion process in plasmas that can be described by the Landau-Fokker-Planck (LFP) equation or the stochastic differential equation (SDE). While energy and momentum are conserved exactly in the LFP…
Electromagnetic wave scattering from planar dielectric films deposited on one-dimensional, randomly rough, perfectly conducting substrates is studied by numerical simulations for both p- and s-polarization. The reduced Rayleigh equation,…
We benchmark and analyze the error of energy conservation (EC) scheme for Particle in cell/Monte-Carlo Couple (PIC/MCC) algorithms by a radio frequency discharging simulation. The plasma heating behaviors and electron distributing functions…
In this paper, exponential energy-preserving methods are formulated and analysed for solving charged-particle dynamics in a strong and constant magnetic field. The resulting method can exactly preserve the energy of the dynamics. Moreover,…
The Energy Conserving Descent (ECD) algorithm was recently proposed (De Luca & Silverstein, 2022) as a global non-convex optimization method. Unlike gradient descent, appropriately configured ECD dynamics escape strict local minima and…
The computational bottleneck of molecular dynamics is the pairwise additive long-range interactions between particles. The random batch Ewald (RBE) method provides a highly efficient and superscalable solver for long-range interactions, but…
The Lorenz--Mie formulation of electromagnetic scattering by a homogeneous, isotropic, dielectric-magnetic sphere was extended to incorporate topologically insulating surface states characterized by a surface admittance $\gamma$.…
A new approach using a hyperbolic-equation system (HES) is proposed to solve for the electron fluids in quasi-neutral plasmas. The HES approach avoids treatments of cross-diffusion terms which cause numerical instabilities in conventional…
A new computational algorithm, the discrete singular convolution (DSC), is introduced for computational electromagnetics. The basic philosophy behind the DSC algorithm for the approximation of functions and their derivatives is studied.…
The $T$-matrix, often obtained with Waterman's extended boundary condition method (EBCM), is a widely-used tool for fast calculations of electromagnetic scattering by particles. Here we investigate the quasistatic or long-wavelength limit…
Studying single-particle dynamics over many periods of oscillations is a well-understood problem solved using symplectic integration. Such integration schemes derive their update sequence from an approximate Hamiltonian, guaranteeing that…