Asymptotic-preserving gyrokinetic implicit particle-orbit integrator for arbitrary electromagnetic fields (2312.00730v1)
Abstract: We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in [Ricketson & Chac\'on, JCP, 2020] to include finite Larmor-radius (FLR) effects in the presence of electric-field length-scales comparable to the particle gyro-radius (the gyro-kinetic limit). We introduce two modifications to the earlier scheme. The first is the explicit gyro-averaging of the electric field at the half time-step, along with an analogous modification to the current deposition, which we show preserves total energy conservation in implicit PIC schemes. The number of gyrophase samples is chosen adaptively, ensuring proper averaging for large timesteps, and the recovery of full-orbit dynamics in the small time-step limit. The second modification is an alternating large and small time-step strategy that ensures the particle trajectory samples gyrophases evenly. We show that this strategy relaxes the time-step restrictions on the scheme, allowing even larger speed-ups than previously achievable. We demonstrate the new method with several single-particle motion tests in a variety of electromagnetic field configurations featuring gyro-scale variation in the electric field. The results demonstrate the advertised ability to capture FLR effects accurately even when significantly stepping over the gyration time-scale.
- Multiscale gyrokinetics for rotating tokamak plasmas: fluctuations, transport and energy flows. Reports on Progress in Physics, 76(11):116201, 2013.
- Plasma physics via computer simulation. CRC press, 2004.
- Simulation of low-frequency, electromagnetic phenomena in plasmas. In Multiple time scales, pages 271–310. Elsevier, 1985.
- Foundations of nonlinear gyrokinetic theory. Reviews of modern physics, 79(2):421, 2007.
- “one size fits all” analytic solutions to the Grad–Shafranov equation. Physics of Plasmas, 17(3):032502, 2010.
- A charge-and energy-conserving implicit, electrostatic particle-in-cell algorithm on mapped computational meshes. Journal of Computational Physics, 233:1–9, 2013.
- Uniformly accurate methods for vlasov equations with non-homogeneous strong magnetic field. Mathematics of Computation, 88(320):2697–2736, 2019.
- Uniformly accurate methods for three dimensional vlasov equations under strong magnetic field with varying direction. SIAM Journal on Scientific Computing, 42(2):B520–B547, 2020.
- G. Chen and L. Chacón. An energy-and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov–Darwin particle-in-cell algorithm. Computer Physics Communications, 185(10):2391–2402, 2014.
- G. Chen and L. Chacón. A multi-dimensional, energy-and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm. Computer Physics Communications, 197:73–87, 2015.
- An energy-and charge-conserving, implicit, electrostatic particle-in-cell algorithm. Journal of Computational Physics, 230(18):7018–7036, 2011.
- An efficient mixed-precision, hybrid CPU–GPU implementation of a nonlinearly implicit one-dimensional particle-in-cell algorithm. Journal of Computational Physics, 231(16):5374–5388, 2012.
- Fluid preconditioning for Newton–Krylov-based, fully implicit, electrostatic particle-in-cell simulations. Journal of computational physics, 258:555–567, 2014.
- A semi-implicit, energy-and charge-conserving particle-in-cell algorithm for the relativistic vlasov-maxwell equations. arXiv preprint arXiv:1903.01565, 2019.
- Y. Chen and S. E. Parker. Particle-in-cell simulation with Vlasov ions and drift kinetic electrons. Physics of Plasmas, 16(5):052305, May 2009.
- D. Curreli. private communication.
- Kinetic particle simulations in a global toroidal geometry. Physics of Plasmas, 26(8):082507, Aug 2019.
- Nonlinear gyrokinetic equations. The Physics of fluids, 26(12):3524–3535, 1983.
- T Zh Esirkepov. Exact charge conservation scheme for particle-in-cell simulation with an arbitrary form-factor. Computer Physics Communications, 135(2):144–153, 2001.
- F. Filbet and L.M. Rodrigues. Asymptotically stable particle-in-cell methods for the vlasov–poisson system with a strong external magnetic field. SIAM Journal on Numerical Analysis, 54(2):1120–1146, 2016.
- F. Filbet and L.M. Rodrigues. Asymptotically preserving particle-in-cell methods for inhomogeneous strongly magnetized plasmas. SIAM Journal on Numerical Analysis, 55(5):2416–2443, 2017.
- Asymptotically preserving particle methods for strongly magnetized plasmas in a torus. Journal of Computational Physics, 480:112015, 2023.
- A fast implicit algorithm for highly magnetized charged particle motion. Open Plasma Physics Journal, 3:36–41, 2010.
- Nonlinear gyrokinetic theory for finite-beta plasmas. The Physics of fluids, 31(7):1940–1948, 1988.
- Large-stepsize integrators for charged-particle dynamics over multiple time scales. Numerische Mathematik, 151(3):659–691, 2022.
- Fully kinetic simulation of ion-temperature-gradient instabilities in tokamaks. Plasma, 1(1):105–118, May 2018.
- Ilon Joseph. Guiding center and gyrokinetic orbit theory for large electric field gradients and strong shear flows. Physics of Plasmas, 28(4), 2021.
- Fast nonlinear iterative solver for an implicit, energy-conserving, asymptotic-preserving charged-particle orbit integrator. Journal of Computational Physics, 459:111146, 2022.
- WW Lee. Gyrokinetic approach in particle simulation. The Physics of Fluids, 26(2):556–562, 1983.
- Conversion from gyrokinetic to full-orbit particles for plasma sheath boundary conditions. In APS Meeting Abstracts, 2018.
- Numerical error in electron orbits with large ωceδtsubscript𝜔𝑐𝑒𝛿𝑡\omega_{ce}\delta titalic_ω start_POSTSUBSCRIPT italic_c italic_e end_POSTSUBSCRIPT italic_δ italic_t. Journal of Computational Physics, 97(1):91–102, 1991.
- A fast, high-order solver for the Grad–Shafranov equation. Journal of Computational Physics, 243:28–45, 2013.
- L.F. Ricketson and L. Chacón. An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields. arXiv preprint arXiv:1904.09478, 2019.
- L.F. Ricketson and G. Chen. A pseudospectral implicit particle-in-cell method with exact energy and charge conservation. Computer Physics Communications, 291:108811, 2023.
- Exact energy-momentum conserving algorithms and symplectic schemes for nonlinear dynamics. Computer methods in applied mechanics and engineering, 100(1):63–116, 1992.
- B. J. Sturdevant. Fully kinetic ion models for magnetized plasma simulations. Doctoral Dissertation at CU Boulder, 2016.
- Low frequency fully kinetic simulation of the toroidal ion temperature gradient instability. Physics of Plasmas, 24(8):081207, Aug 2017.
- An implicit δf𝛿𝑓\delta fitalic_δ italic_f particle-in-cell method with sub-cycling and orbit averaging for lorentz ions. Journal of Computational Physics, 316:519 – 533, 2016.
- Development of a consistent and stable fully implicit moment method for Vlasov–Ampère particle in cell (PIC) system. SIAM Journal on Scientific Computing, 35(5):S126–S149, 2013.
- A new charge conservation method in electromagnetic particle-in-cell simulations. Computer Physics Communications, 156(1):73–85, 2003.
- J. Villasenor and O. Buneman. Rigorous charge conservation for local electromagnetic field solvers. Computer Physics Communications, 69(2-3):306–316, 1992.
- Accurate numerical solution of charged particle motion in a magnetic field. Journal of Computational Physics, 116(2):384–387, 1995.
- Geometric electrostatic particle-in-cell algorithm on unstructured meshes. 12 2020.