Related papers: Finite spatial-grid effects in energy-conserving p…
We perform a scaling and performance portability study of the particle-in-cell scheme for plasma physics applications through a set of mini-apps we name "Alpine", which can make use of exascale computing capabilities. The mini-apps are…
A fundamental task in particle-in-cell (PIC) simulations of plasma physics is solving for charged particle motion in electromagnetic fields. This problem is especially challenging when the plasma is strongly magnetized due to numerical…
A conservative scheme has been formulated and verified for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from drift kinetic equation is used to…
Depending on the physical conditions involved the beam plasma systems may reveal new unstable regimes triggered by the wave instabilities of different nature. We show through linear theory and numerical simulations the existence of an…
In this paper, we investigate the kinetic stability of classical, collisional plasma - that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared to the length scale $L$ over which fields and bulk…
Here we present the mathematical and numerical details of a general hybrid model for plasmas. All grid quantities are stored at cell centers on the grid. The most common discretization of the fields in PIC solvers is to have the electric…
Reduced-order particle-in-cell (PIC) scheme is a novel modeling approach that enables computationally efficient electrostatic kinetic simulations of plasma. In our previous publications, we demonstrated the potentials of a preliminary…
Split form schemes for Euler and Navier-Stokes equations are useful for computation of turbulent flows due to their better robustness. This is because they satisfy additional conservation properties of the governing equations like kinetic…
The efficient numerical solution of many kinetic models in plasma physics is impeded by the stiffness of these systems. Exponential integrators are attractive in this context as they remove the CFL condition induced by the linear part of…
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…
Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC)…
We introduce $\texttt{A}$strophysical $\texttt{H}$ybrid-$\texttt{K}$inetic simulations with the $\texttt{flASH}$ code ($\texttt{AHKASH}$) -- a new Hybrid particle-in-cell (PIC) code developed within the framework of the multi-physics code…
Particle in Cell (PIC) simulations are a widely used tool for the investigation of both laser- and beam-driven plasma acceleration. It is a known issue that the beam quality can be artificially degraded by numerical Cherenkov radiation…
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…
When a charged particle moves through a plasma at a speed much higher than the thermal velocity of the plasma, it is subjected to the force of the electrostatic field induced in the plasma by itself and loses its energy. This process is…
The effect of a uniform magnetic field on particle transport in low-pressure radio frequency (RF) capacitively coupled plasma (CCP) has been studied using a particle-in-cell (PIC) model. Three distinct regimes of plasma behavior can be…
In this work, we propose and develop efficient and accurate numerical methods for solving the Kirchhoff-Love plate model in domains with complex geometries. The algorithms proposed here employ curvilinear finite-difference methods for…
The Energy Conserving Semi-Implicit Method (ECSIM) introduced by Lapenta (2017) has many advantageous properties compared to the classical semi-implicit and explicit PIC methods. Most importantly, energy conservation eliminates the growth…
The equations of electrostatic drift kinetics are observed to possess a symmetry associated with their intrinsic scale invariance. Under the assumptions of spatial periodicity, stationarity, and locality, this symmetry implies a particular…
The growth and saturation of Buneman-type instabilities is examined with a particle-in-cell (PIC) simulation for parameters that are representative for the foreshock region of fast supernova remnant (SNR) shocks. A dense ion beam and the…