Magnetic confinement fusion reactors feature numerous disparate time-scales, with the gyrofrequency induced by the strong magnetic field among the fastest. While asymptotic models of charged particle motion permit stepping over this scale, there is increasing concern that their ordering may break down in certain key regions of such reactors. This generates interest in asymptotic-preserving schemes that can recover the correct asymptotic limit when it is valid while still converging to full dynamics in the limit of small time-step. We present such a scheme which, in addition to the desired asymptotic properties, preserves the exact energy and charge conservation enjoyed by implicit particle-in-cell (PIC) methods developed over the last decade. We discuss the scheme’s theoretical development, an efficient and robust nonlinear solver, and numerical results for both single-particle test cases and self-consistent PIC simulations. This paper is co-authored with L. Chacón, O. Koshkarov, and G. Chen
Copyright. All Rights Reserved.