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

JOURNAL OF COMPUTATIONAL PHYSICS

卷:258

Article

Chen, G.¹  Chacon, L.¹  Leibs, C. A.²  Knoll, D. A.¹  Taitano, W.³ 

[1] Los Alamos Natl Lab, Los Alamos, NM 87545 USA

[2] Univ Colorado, Boulder, CO 80309 USA

[3] Univ New Mexico, Albuquerque, NM 87131 USA

关键词: Electrostatic particle-in-cell;    Implicit methods;    Direct implicit;    Implicit moment;    Energy conservation;    Charge conservation;    Physics based preconditioner;    JFNK solver;   

DOI  :  10.1016/j.jcp.2013.10.052

来源: Elsevier

【 摘 要 】

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. (C) 2013 Elsevier Inc. All rights reserved.

【 授权许可】

Free