Related papers: On numerical errors to the fields surrounding a re…
A recent proof-of-principle study proposes a nonlinear electrostatic implicit particle-in-cell (PIC) algorithm in one dimension (Chen, Chacon, Barnes, J. Comput. Phys. 230 (2011) 7018). The algorithm employs a kinetically enslaved…
We introduce a new electrostatic particle-in-cell algorithm capable of using large timesteps compared to particle gyro-period under a uniform external magnetic field. The algorithm extends earlier electrostatic fully implicit PIC…
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…
The recently developed energy conserving semi-implicit method (ECsim) for PIC simulation is applied to multiple scale problems where the electron-scale physics needs to be only partially retained and the interest is on the macroscopic or…
Recently proposed modifications of the standard particle-in-cell (PIC) method resolve long-standing limitations such as exact preservation of physically conserved quantities and unbiased ensemble down-sampling. Such advances pave the way…
In this paper, we develop an asymptotic-preserving and energy-conserving (APEC) Particle-In-Cell (PIC) algorithm for the Vlasov-Maxwell system. This algorithm not only guarantees that the asymptotic limiting of the discrete scheme is a…
We propose a spectral Particle-In-Cell (PIC) algorithm that is based on the combination of a Hankel transform and a Fourier transform. For physical problems that have close-to-cylindrical symmetry, this algorithm can be much faster than…
In a previous paper, we developed a new particle-in-cell method for the Vlasov-Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a…
When supersonic plasma flows collide, many physical processes contribute to the morphology of the resulting shock. One of these processes is the acceleration of non-thermal ions, which will, eventually, reach relativistic speeds and become…
We present a new technique for transferring momentum and velocity between particles and grid with Particle-In-Cell (PIC) calculations which we call Affine-Particle-In-Cell (APIC). APIC represents particle velocities as locally affine,…
With the increases in computational power and advances in machine learning, data-driven learning-based methods have gained significant attention in solving PDEs. Physics-informed neural networks (PINNs) have recently emerged and succeeded…
This paper provides a rigorous analysis of boundary element methods for the magnetic field integral equation on Lipschitz polyhedra. The magnetic field integral equation is widely used in practical applications to model electromagnetic…
It is demonstrated that the inherent structure of Boris pusher may result in the accumulation of errors in numerical integration when estimating the works of electric fields on particles. This, in turn, leads to an incorrect estimations of…
Partial differential equations with distributional sources---in particular, involving (derivatives of) delta distributions---have become increasingly ubiquitous in numerous areas of physics and applied mathematics. It is often of…
Over the past decades, multiple gyrokinetic codes have shown to be able to simulate turbulence and associated transport in the core of Tokamak devices. However, their application to the edge and scrape-off layer (SOL) region presents…
In order to perform physically faithful particle-in-cell (PIC) simulations, the Gauss's law stands as a critical requirement, since its violation often leads to catastrophic errors in long-term plasma simulations. This work proposes a novel…
The origin of the Finite Grid Instability (FGI) is studied by resolving the dynamics in the 1D electrostatic Particle-In-Cell (PIC) model in spectral domain at the single particle level and at the collective motion level. The spectral…
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…
A simulation package employing a Particle in Cell (PIC) method is developed to study the high current beam transport and the dynamics of plasmas. This package includes subroutines those are suited for various planned projects at University…
A method of solving Maxwell equations in a vicinity of a multipole particle (moving along an arbitrary trajectory) is proposed. The method is based on a geometric construction of a trajectory-adapted coordinate system, which simplifies…