Fluid preconditioning for Newton-Krylov-based, fully implicit, electrostatic particle-in-cell simulations

Abstract

A recent proof-of-principle study proposes an energy- and charge-conserving, nonlinearly implicit electrostatic particle-in-cell (PIC) algorithm in one dimension [Chen et al, J. Comput. Phys., 230 (2011) 7018]. The algorithm in the reference employs an unpreconditioned Jacobian-free Newton-Krylov method, which ensures nonlinear convergence at every timestep (resolving the dynamical timescale of interest). Kinetic enslavement, which is one key component of the algorithm, not only enables fully implicit PIC a practical approach, but also allows preconditioning the kinetic solver with a fluid approximation. This study proposes such a preconditioner, in which the linearized moment equations are closed with moments computed from particles. Effective acceleration of the linear GMRES solve is demonstrated, on both uniform and non-uniform meshes. The algorithm performance is largely insensitive to the electron-ion mass ratio. Numerical experiments are performed on a 1D multi-scale ion acoustic wave test problem.

21 pages, 8 figures

References in corpus (1)

• An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov-Darwin particle-in-cell algorithm

Cited by in corpus (14)

• A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov-Darwin particle-in-cell algorithm
• An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov-Darwin particle-in-cell algorithm
• An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields
• Arbitrarily high order Convected Scheme solution of the Vlasov-Poisson system
• ECsim-CYL: Energy Conserving Semi-Implicit particle in cell simulation in axially symmetric cylindrical coordinates
• An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system
• Verification of a Fully Implicit Particle-in-Cell Method for the $v_{\parallel }$ Formalism of Electromagnetic Gyrokinetics in the XGC Code
• Accelerating the convergence of Newton's method for nonlinear elliptic PDEs using Fourier neural operators
• An implicit, conservative electrostatic particle-in-cell algorithm for paraxial magnetic nozzles
• Implicit temporal discretization and exact energy conservation for particle methods applied to the Poisson-Boltzmann equation
• The shape of convection in 2D and 3D global simulations of stellar interiors
• Asymptotic-preserving Particle-In-Cell methods for the Vlasov-Maxwell system near quasi-neutrality
• Machine learning accelerated particle-in-cell plasma simulations
• Bridging kinetic plasma descriptions and single fluid models