Related papers: An Arbitrary Curvilinear Coordinate Method for Par…
This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma…
Applying gyrokinetic simulations for theoretical turbulence and transport studies to the plasma edge and scrape-off layer (SOL) presents significant challenges. To in particular account for steep density and temperature gradients in the…
In this article, we introduce a novel dimensionality reduction formulation for the Poisson's equation in the Vlasov-Poisson system that yields a reduced-order particle-in-cell scheme. This scheme allows a remarkable reduction in the…
High-resolution simulations of particle-based kinetic plasma models typically require a high number of particles and thus often become computationally intractable. This is exacerbated in multi-query simulations, where the problem depends on…
In this study, a novel method for simulating plasma dynamics using parallel programming has been developed. The equations based on Particle-in-Cell (PIC) method were utilized and adapted for this purpose. We utilized 35 processors from…
We develop an algorithm for bosonic path integral molecular dynamics (PIMD) simulations with periodic boundary conditions (PBC) that scales quadratically with the number of particles. Path integral methods are a powerful tool to simulate…
We develop a method for simulating colloidal suspensions using multiparticle collision dynamics (MPCD) with a discrete particle model represented as a rigid body. The key steps for incorporating the rigid-body constraints are to thermalize…
The PIC model relies on two building blocks. The first stems from the capability of computers to handle only up to $\sim10^{10}$ particles, while real plasmas contain from $10^4$ to $10^{20}$ particles per Debye sphere: a coarse-graining…
Three-dimensional (3D) particle-in-cell (PIC) simulations are used to investigate the interaction of ultrahigh intensity lasers ($> 10^{20}$ W/cm$^{-2}$) with matter at overcritical densities. Intense laser pulses are shown to penetrate up…
The conductor-like polarization model (C-PCM) with switching/Gaussian smooth discretization is a widely used implicit solvation model in chemical simulations. However, its application in quantum mechanical calculations of large-scale…
This paper presents a novel approach for simulating plasma instabilities in tokamak plasmas using the piecewise field-aligned finite element method in combination with the particle-in-cell method. Our method traditionally aligns the…
With an implicit Particle-in-cell/Monte Carlo model, capacitively coupled plasmas are studied in two-dimensional and axisymmetric geometry. Self-bias dc voltage effects are self-consistently considered. Due to finite length effects, the…
The acceleration of protons and electrons in a reconnecting current sheet (RCS) is simulated with a particle-in-cell (PIC) 2D3V code for the proton-to-electron mass ratio of 100. The electro-magnetic configuration forming the RCS…
For many plasma physics problems, three-dimensional and kinetic effects are very important. However, such simulations are very computationally intensive. Fortunately, there is a class of problems for which there is nearly azimuthal symmetry…
Particle-mesh methods, such as the particle-in-cell (PIC) method, cannot retain exact pairwise interaction at sub-cell scales. For dense nonneutral relativistic electron bunches, this makes it difficult to accurately capture the…
A two-dimensional particle-in-cell code for simulation of low-frequency electromagnetic processes in laboratory plasmas has been developed. The code uses the Darwin method omitting the electromagnetic wave propagation. The Darwin method…
Large linear systems with sparse, non-symmetric matrices arise in the modeling of Markov chains or in the discretization of convection-diffusion problems. Due to their potential to solve sparse linear systems with an effort that is linear…
We present a robust and accurate numerical method for the anisotropic diffusion equation in curvilinear coordinates. This study extends the recent work [Muir et al., Computer Physics Communications, 2025] for solving the anisotropic…
The numerical solution of high dimensional Vlasov equation is usually performed by particle-in-cell (PIC) methods. However, due to the well-known numerical noise, it is challenging to use PIC methods to get a precise description of the…
Hamiltonian particle-based simulations of plasma dynamics are inherently computationally intensive, primarily due to the large number of particles required to obtain accurate solutions. This challenge becomes even more acute in many-query…