Abstract:In this paper we present a customized finite-difference-time-domain (FDTD) Maxwell solver for the particle-in-cell (PIC) algorithm. The solver is customized to effectively eliminate the numerical Cerenkov instability (NCI) which arises when a plasma (neutral or non-neutral) relativistically drifts on a grid when using the PIC algorithm. We control the EM dispersion curve in the direction of the plasma drift of a FDTD Maxwell solver by using a customized higher order finite difference operator for the spatial derivative along the direction of the drift ($\hat 1$ direction). We show that this eliminates the main NCI modes with moderate $\vert k_1 \vert$, while keeps additional main NCI modes well outside the range of physical interest with higher $\vert k_1 \vert$. These main NCI modes can be easily filtered out along with first spatial aliasing NCI modes which are also at the edge of the fundamental Brillouin zone. The customized solver has the possible advantage of improved parallel scalability because it can be easily partitioned along $\hat 1$ which typically has many more cells than other directions for the problems of interest. We show that FFTs can be performed locally to current on each partition to filter out the main and first spatial aliasing NCI modes, and to correct the current so that it satisfies the continuity equation for the customized spatial derivative. This ensures that Gauss' Law is satisfied. We present simulation examples of one relativistically drifting plasmas, of two colliding relativistically drifting plasmas, and of nonlinear laser wakefield acceleration (LWFA) in a Lorentz boosted frame that show no evidence of the NCI can be observed when using this customized Maxwell solver together with its NCI elimination scheme.

Modified Piecewise Linear Current Density Recursive Convolution Finite-Difference Timedomain Method for Anisotropic Magnetised Plasmas

One-step Leapfrog ADI-FDTD Method for Anisotropic Magnetized Plasma

Improved FDTD simulation for predicting transmission characteristics of graphene and solid plasma composite with a biased magnetic field

One-step leapfrog ADI-FDTD method for simulating electromagnetic wave propagation in general dispersive media.

Development of a Novel FDTD (2, 4)-Compatible Conformal Scheme for Electromagnetic Computations of Complex Curved PEC Objects

The Penetrated and Scattered Field Distributions in the Presence of a Non-Uniform Cold Magnetised Plasma Column Surrounded By Multiple Strips: Tm Wave Incidence

Numerical Dispersion Optimized Leapfrog Adi-Fdtd Method And Its Application For Capturing Surface Current Distributions Of Complex Objects Illuminated By An Electromagnetic Pulse (Emp)

An Area-Averaged Conformal Technique for the Electromagnetic Simulations of Debye Medium

The Penetrated and Scattered Field Distributions in the Presence of A Nonuniform Cold Magnetized Plasma Column Surrounded by Multiple Apertures: Te-Wave Incidence

A modified high-order FDTD conformal technique for modeling perfectly conducting objects

Analysis of electromagnetic wave propagation and scattering characteristics of plasma shealth via high order ADE-ADI FDTD

Improved HIE-FDTD Method for Simulating Thin Layer Structure in Linear Dispersive Media

Newton-ADE FDTD Method for Time-Varying Plasma

Combined Piecewise Linear Recursive Convolution-Bilinear Transform Implementation of the CFS-PML for Unmagnetized Plasma

A 2-D Axisymmetric Magneto-Hydrodynamic Model of a DC Arc Plasma Torch and Its Solution Methodology

Tridiagonalized WLP-FDTD Implementation for Wave Propagation in Magnetized Plasma

Perfectly Matched Layer for Accurate FDTD for Anisotropic Magnetized Plasma

Efficient Implementation of the CPML in 3-D Hybrid Implicit-Explicit FDTD Method

Novel Scheme for GPU‐Accelerated Finite‐Difference Time‐Domain Simulation of Electromagnetic Wave Interaction With Magnetic Plasma

A Diagnostic Method of Radial Electron Density Distribution Based on Microwave Transmission of Cylindrical Plasma along Multipath

Controlling the Numerical Cerenkov Instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction