A physics-based-preconditioned quasi-Newton solver for energy- and momentum-conserving implicit hybrid particle-in-cell simulations
Abstract
Explicit hybrid particle-in-cell (PIC) schemes, with kinetic ions and fluid electrons, accumulate numerical heating over long simulations and are unstable to finite-grid instabilities. Energy-conserving implicit integration removes heating and instabilities, but may be inefficient due to the need of nonlinear iteration. We present a fast quasi-Newton solver for an implicit hybrid-PIC formulation. The nonlinear residual employs the electric field and electron pressure as nonlinear unknowns, and nonlinearly eliminates particles. The formulation is strictly energy conserving and, on Cartesian meshes, also momentum conserving. The nonlinear iteration is efficiently performed with a quasi-Newton solver that replaces the true Jacobian with an approximate one without spoiling conservation properties upon convergence. The approximate Jacobian replaces the particle ion linear response with an analytical mass matrix that is local and requires minimal storage. This approximate Jacobian is efficiently inverted with a physics-based multigrid preconditioner, adapted from Hall MHD. We verify the scheme linearly with whistler, kinetic Alfvén, and cylindrical normal modes, and nonlinearly with the GEM reconnection challenge. On GEM, explicit integration needs a forty-times-smaller time step, and nine times the wall-clock time, to reach an energy conservation comparable to that of the implicit solver. We apply the method to the simulation of a three-dimensional field-reversed configuration formation by external coils.
Related
- A fully implicit electromagnetic energy- and charge-conserving PIC model for simulations of high-β plasma
- A semi-implicit unified gas-kinetic wave-particle method with energy-conserving electromagnetic coupling for multiscale plasmas
- A θ-implicit hybrid particle-in-cell algorithm for kinetic simulations of magnetized plasmas
- An Implicit, Moment-Conserving Low-Rank Tensor Method for the Hybrid Kinetic-Ion Fluid-Electron Plasma Model
- A Semi-Implicit All-Mach Numerical Method for Implicit Large Eddy Simulations
BibTeX
@misc{kumar2026physicsbasedpreconditioned,
title = {{A physics-based-preconditioned quasi-Newton solver for energy- and momentum-conserving implicit hybrid particle-in-cell simulations}},
author = {Prabhat Kumar and Luis Chacón and William J. White and Lujain Almarhabi and Stephen E. Clark},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.08402},
eprint = {2611.08402},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.08402},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - A physics-based-preconditioned quasi-Newton solver for energy- and momentum-conserving implicit hybrid particle-in-cell simulations
AU - Prabhat Kumar
AU - Luis Chacón
AU - William J. White
AU - Lujain Almarhabi
AU - Stephen E. Clark
PY - 2026
DA - 2026/11/05
PB - arXiv
DO - 10.48550/arxiv.2611.08402
UR - https://arxiv.org/abs/2611.08402
AB - Explicit hybrid particle-in-cell (PIC) schemes, with kinetic ions and fluid electrons, accumulate numerical heating over long simulations and are unstable to finite-grid instabilities. Energy-conserving implicit integration removes heating and instabilities, but may be inefficient due to the need of nonlinear iteration. We present a fast quasi-Newton solver for an implicit hybrid-PIC formulation. The nonlinear residual employs the electric field and electron pressure as nonlinear unknowns, and nonlinearly eliminates particles. The formulation is strictly energy conserving and, on Cartesian meshes, also momentum conserving. The nonlinear iteration is efficiently performed with a quasi-Newton solver that replaces the true Jacobian with an approximate one without spoiling conservation properties upon convergence. The approximate Jacobian replaces the particle ion linear response with an analytical mass matrix that is local and requires minimal storage. This approximate Jacobian is efficiently inverted with a physics-based multigrid preconditioner, adapted from Hall MHD. We verify the scheme linearly with whistler, kinetic Alfvén, and cylindrical normal modes, and nonlinearly with the GEM reconnection challenge. On GEM, explicit integration needs a forty-times-smaller time step, and nine times the wall-clock time, to reach an energy conservation comparable to that of the implicit solver. We apply the method to the simulation of a three-dimensional field-reversed configuration formation by external coils.
ER -