An \mathcal{O}(N) Maxwell solver with improved numerical dispersion properties

Abstract

A Maxwell solver derived from finite element method with \mathcal{O}(N) computing cost is developed to improve the numerical dispersion properties in relativistic particle-in-cell (PIC) simulations. The correction of the dispersion relation of the electromagnetic wave is achieved using the neighboring cells via an iteration scheme without decomposing into Fourier modes. The local nature of the communication is ideally suited to massively parallel computer architectures. This Maxwell solver constrains the Numerical Cherenkov instability (NCI) for the ultra-relativistic drifting pair plasma in x direction to large wave vectors for two dimensional grid. The growth rate of NCI is suppressed by using the low pass filtering.

arXiv comment16 pages, 4 figures, submitted to Computer Physics Communications

References in corpus (8)

• Numerical instability due to relativistic plasma drift in EM-PIC simulations
• Numerical stability analysis of the Pseudo-Spectral Analytical Time-Domain PIC algorithm
• Controlling the Numerical Cerenkov Instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction
• Improved Numerical Cherenkov Instability Suppression in the Generalized PSTD PIC Algorithm
• A Systematic Approach to Numerical Dispersion in Maxwell Solvers
• Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm
• Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit
• On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability