Related papers: A multi-dimensional, energy- and charge-conserving…
We describe a charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This…
The advent of high-power Hall thrusters and the increasing interest towards their use as a primary propulsion system for various missions have given a new boost to the efforts aiming at self-consistent predictive modeling of this thruster…
Plasma dynamics is a multi-scale problem that involves many spatial and temporal scales. Turbulence connects the disparate scales in this system through a cascade that is established by nonlinear interactions. Most astrophysical plasma…
Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These…
Fusion energy offers the potential for the generation of clean, safe, and nearly inexhaustible energy. While notable progress has been made in recent years, significant challenges persist in achieving net energy gain. Improving plasma…
A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code,…
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…
We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in [Ricketson & Chac\'on, JCP, 2020] to include finite Larmor-radius (FLR) effects in the presence of electric-field…
The effects of noise in particle-in-cell (PIC) and Vlasov simulations of the Buneman instability in unmagnetized plasmas are studied. It is found that, in the regime of low drift velocity, the linear stage of the instability in PIC…
There are many interesting physical processes which involve the generation of high density plasmas in large volumes. However, when modeling these systems numerically, the large densities and volumes present a significant computational…
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…
We introduce a novel particle-in-Fourier (PIF) scheme that extends its applicability to non-periodic boundary conditions. Our method handles free space boundary conditions by replacing the Fourier Laplacian operator in PIF with a mollified…
In this article, we present an in-depth verification of the generalized electrostatic reduced-order particle-in-cell (PIC) scheme in a cross electric and magnetic field configuration representative of a radial-azimuthal section of a Hall…
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…
Room-temperature metals and semi-metals which consist of a gas of bound electrons in a near-continuum band structure can be classified as cold quantum plasmas. This insight suggests that Particle-in-Cell (PIC) simulations, traditionally…
The accurate prediction of occurrence and strength of kinetic instabilities in plasmas remains a significant challenge in nuclear fusion research. To accurately capture the plasmas dynamics one is required to solve the Vlasov equation for…
Modeling multi-scale collisionless magnetized processes constitutes an important numerical challenge. By treating electrons as a fluid and ions kinetically, the so-called hybrid Particle-In-Cell (PIC) codes represent a promising…
The Vlasov-Maxwell system of equations, which describes classical plasma physics, is extremely challenging to solve, even by numerical simulation on powerful computers. By linearizing and assuming a Maxwellian background distribution…
A Particle-in-Cell (PIC) numerical simulation of the electron Weibel instability is applied in a frame of Darwin (radiationless) approximation of the self-consistent fields of sparse plasma. As a result, we were able to supplement the…
This paper formulates a new particle-in-cell method for the Vlasov-Maxwell system. Under the Lorenz gauge condition, Maxwell's equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The…