skip to main content

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

A recent proof-of-principle study proposes an energy- and charge-conserving, nonlinearly implicit electrostatic particle-in-cell (PIC) algorithm in one dimension [9]. 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 as 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.
Authors:
 [1] ;  [1] ;  [2] ;  [1] ;  [3]
  1. Los Alamos National Laboratory, Los Alamos, NM 87545 (United States)
  2. University of Colorado Boulder, Boulder, CO 80309 (United States)
  3. University of New Mexico, Albuquerque, NM 87131 (United States)
Publication Date:
OSTI Identifier:
22230859
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 258; Other Information: Copyright (c) 2013 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICAL METHODS AND COMPUTING; 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ACCELERATION; ALGORITHMS; APPROXIMATIONS; CHARGE CONSERVATION; CONVERGENCE; ENERGY CONSERVATION; EQUATIONS; FLUIDS; ION ACOUSTIC WAVES; MASS; NONLINEAR PROBLEMS; SIMULATION