A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes
Abstract
We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.
- Authors:
-
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Publication Date:
- Research Org.:
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC). Advanced Scientific Computing Research (ASCR); USDOE National Nuclear Security Administration (NNSA)
- OSTI Identifier:
- 1859901
- Report Number(s):
- LA-UR-21-29944
Journal ID: ISSN 0021-9991; TRN: US2305382
- Grant/Contract Number:
- 89233218CNA000001; AC52-06NA25396
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Computational Physics
- Additional Journal Information:
- Journal Volume: 459; Journal ID: ISSN 0021-9991
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; hybrid; curvilinear; plasma; particle-in-cell; space weather; fusion
Citation Formats
Stanier, Adam John, and Chacon, Luis. A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes. United States: N. p., 2022.
Web. doi:10.1016/j.jcp.2022.111144.
Stanier, Adam John, & Chacon, Luis. A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes. United States. https://doi.org/10.1016/j.jcp.2022.111144
Stanier, Adam John, and Chacon, Luis. Tue .
"A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes". United States. https://doi.org/10.1016/j.jcp.2022.111144. https://www.osti.gov/servlets/purl/1859901.
@article{osti_1859901,
title = {A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes},
author = {Stanier, Adam John and Chacon, Luis},
abstractNote = {We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.},
doi = {10.1016/j.jcp.2022.111144},
journal = {Journal of Computational Physics},
number = ,
volume = 459,
place = {United States},
year = {Tue Mar 15 00:00:00 EDT 2022},
month = {Tue Mar 15 00:00:00 EDT 2022}
}
Works referenced in this record:
Low frequency fully kinetic simulation of the toroidal ion temperature gradient instability
journal, August 2017
- Sturdevant, B. J.; Chen, Y.; Parker, S. E.
- Physics of Plasmas, Vol. 24, Issue 8
Finite-Grid Instability in Quasineutral Hybrid Simulations
journal, April 1995
- Rambo, P. W.
- Journal of Computational Physics, Vol. 118, Issue 1
Hybrid Simulations of the Effects of Energetic Particles on Low-Frequency MHD Waves
journal, September 1997
- Belova, E. V.; Denton, R. E.; Chan, A. A.
- Journal of Computational Physics, Vol. 136, Issue 2
High-Resolution Hybrid Simulations of Kinetic Plasma Turbulence at Proton Scales
journal, October 2015
- Franci, Luca; Landi, Simone; Matteini, Lorenzo
- The Astrophysical Journal, Vol. 812, Issue 1
A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm
journal, January 2019
- Stanier, A.; Chacón, L.; Chen, G.
- Journal of Computational Physics, Vol. 376
Role of Ion Kinetic Physics in the Interaction of Magnetic Flux Ropes
journal, October 2015
- Stanier, A.; Daughton, W.; Chacón, L.
- Physical Review Letters, Vol. 115, Issue 17
Re-forming supercritical quasi-parallel shocks: 1. One- and two-dimensional simulations
journal, January 1990
- Thomas, V. A.; Winske, D.; Omidi, N.
- Journal of Geophysical Research, Vol. 95, Issue A11
Controlling self-force errors at refinement boundaries for AMR-PIC
journal, February 2010
- Colella, Phillip; Norgaard, Peter C.
- Journal of Computational Physics, Vol. 229, Issue 4
Fluid vs. kinetic magnetic reconnection with strong guide fields
journal, October 2015
- Stanier, A.; Simakov, Andrei N.; Chacón, L.
- Physics of Plasmas, Vol. 22, Issue 10
Wavelet methods for studying the onset of strong plasma turbulence
journal, December 2018
- Le, A.; Roytershteyn, V.; Karimabadi, H.
- Physics of Plasmas, Vol. 25, Issue 12
Current Advance Method and Cyclic Leapfrog for 2D Multispecies Hybrid Plasma Simulations
journal, May 1994
- Matthews, Alan P.
- Journal of Computational Physics, Vol. 112, Issue 1
A non-staggered, conservative, , finite-volume scheme for 3D implicit extended magnetohydrodynamics in curvilinear geometries
journal, November 2004
- Chacón, L.
- Computer Physics Communications, Vol. 163, Issue 3
Jacobian-free Newton–Krylov methods: a survey of approaches and applications
journal, January 2004
- Knoll, D. A.; Keyes, D. E.
- Journal of Computational Physics, Vol. 193, Issue 2
Temperature gradient driven heat flux closure in fluid simulations of collisionless reconnection
journal, June 2018
- Allmann-Rahn, F.; Trost, T.; Grauer, R.
- Journal of Plasma Physics, Vol. 84, Issue 3
Two-way coupled MHD-PIC simulations of magnetic reconnection in magnetic island coalescence
journal, May 2018
- Makwana, Kirit; Keppens, Rony; Lapenta, Giovanni
- Journal of Physics: Conference Series, Vol. 1031
Magnetic field and density fluctuations at perpendicular supercritical collisionless shocks
journal, January 1988
- Winske, D.; Quest, K. B.
- Journal of Geophysical Research, Vol. 93, Issue A9
A computational model for exploring particle acceleration during reconnection in macroscale systems
journal, January 2019
- Drake, J. F.; Arnold, H.; Swisdak, M.
- Physics of Plasmas, Vol. 26, Issue 1
Astrophysical Explosions Revisited: Collisionless Coupling of Debris to Magnetized Plasma
journal, September 2021
- Le, Ari; Winske, Dan; Stanier, Adam
- Journal of Geophysical Research: Space Physics, Vol. 126, Issue 9
A new 3-D spherical hybrid model for solar wind interaction studies: 3-D SPHERICAL HYBRID MODEL
journal, August 2013
- Dyadechkin, S.; Kallio, E.; Jarvinen, R.
- Journal of Geophysical Research: Space Physics, Vol. 118, Issue 8
The link between shocks, turbulence, and magnetic reconnection in collisionless plasmas
journal, June 2014
- Karimabadi, H.; Roytershteyn, V.; Vu, H. X.
- Physics of Plasmas, Vol. 21, Issue 6
A.I.K.E.F.: Adaptive hybrid model for space plasma simulations
journal, April 2011
- Müller, Joachim; Simon, Sven; Motschmann, Uwe
- Computer Physics Communications, Vol. 182, Issue 4
Subproton-Scale Cascades in Solar wind Turbulence: Driven Hybrid-Kinetic Simulations
journal, April 2016
- Cerri, S. S.; Califano, F.; Jenko, F.
- The Astrophysical Journal, Vol. 822, Issue 1
The island coalescence problem: Scaling of reconnection in extended fluid models including higher-order moments
journal, November 2015
- Ng, Jonathan; Huang, Yi-Min; Hakim, Ammar
- Physics of Plasmas, Vol. 22, Issue 11
Self-focusing of local plasma currents
journal, September 1965
- Fadeev, V. M.; Kvabtskhava, I. F.; Komarov, N. N.
- Nuclear Fusion, Vol. 5, Issue 3
Low-noise electromagnetic and relativistic particle-in-cell plasma simulation models
journal, September 1999
- Sydora, R. D.
- Journal of Computational and Applied Mathematics, Vol. 109, Issue 1-2
The coalescence instability in collisionless plasmas
journal, October 1992
- Pritchett, P. L.
- Physics of Fluids B: Plasma Physics, Vol. 4, Issue 10
Structure and Properties of the Foreshock at Venus
journal, October 2017
- Omidi, N.; Collinson, G.; Sibeck, D.
- Journal of Geophysical Research: Space Physics, Vol. 122, Issue 10
Hybrid simulation codes with application to shocks and upstream waves
journal, October 1985
- Winske, D.
- Space Science Reviews, Vol. 42, Issue 1-2
Hybrid simulations of a parallel collisionless shock in the large plasma device
journal, December 2016
- Weidl, Martin S.; Winske, Dan; Jenko, Frank
- Physics of Plasmas, Vol. 23, Issue 12
A robust method for handling low density regions in hybrid simulations for collisionless plasmas
journal, October 2014
- Amano, Takanobu; Higashimori, Katsuaki; Shirakawa, Keisuke
- Journal of Computational Physics, Vol. 275
Three-dimensional stability of current sheets supported by electron pressure anisotropy
journal, October 2019
- Le, A.; Stanier, A.; Daughton, W.
- Physics of Plasmas, Vol. 26, Issue 10
Origin of Mercury’s double magnetopause: 3D hybrid simulation study with A.I.K.E.F.
journal, March 2012
- Müller, Joachim; Simon, Sven; Wang, Yung-Ching
- Icarus, Vol. 218, Issue 1
A global method of solving the electron-field equations in a zero-inertia-electron-hybrid plasma simulation code
journal, December 1980
- Hewett, D. W.
- Journal of Computational Physics, Vol. 38, Issue 3
Electrostatic mode associated with the pinch velocity in reversed field pinch simulations
journal, December 2008
- Delzanno, Gian Luca; Chacón, Luis; Finn, John M.
- Physics of Plasmas, Vol. 15, Issue 12
Bursty magnetic reconnection at the Earth's magnetopause triggered by high-speed jets
journal, September 2021
- Ng, J.; Chen, L. -J.; Omelchenko, Y. A.
- Physics of Plasmas, Vol. 28, Issue 9
Electromagnetic ion beam instabilities: Comparison of one- and two-dimensional simulations
journal, January 1986
- Winske, D.; Quest, K. B.
- Journal of Geophysical Research, Vol. 91, Issue A8
Hybrid simulations of quasineutral phenomena in magnetized plasma
journal, June 1978
- Byers, J. A.; Cohen, B. I.; Condit, W. C.
- Journal of Computational Physics, Vol. 27, Issue 3
A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm
journal, December 2015
- Chen, G.; Chacón, L.
- Computer Physics Communications, Vol. 197
Coalescence instability of magnetic islands
journal, January 1977
- Finn, J. M.; Kaw, P. K.
- Physics of Fluids, Vol. 20, Issue 1
Magnetically Self-Focussing Streams
journal, June 1934
- Bennett, Willard H.
- Physical Review, Vol. 45, Issue 12
Global gyrokinetic particle-in-cell simulations of internal kink instabilities
journal, December 2012
- Mishchenko, Alexey; Zocco, Alessandro
- Physics of Plasmas, Vol. 19, Issue 12
Orbit-averaged implicit particle codes
journal, March 1982
- Cohen, Bruce I.; Freis, Robert P.; Thomas, Vincent
- Journal of Computational Physics, Vol. 45, Issue 3
Hybrid particle-in-cell simulations of laser-driven plasma interpenetration, heating, and entrainment
journal, November 2019
- Higginson, D. P.; Amendt, P.; Meezan, N.
- Physics of Plasmas, Vol. 26, Issue 11
Heuristic stability theory for finite-difference equations
journal, June 1968
- Hirt, C. W.
- Journal of Computational Physics, Vol. 2, Issue 4
A cancellation problem in hybrid particle-in-cell schemes due to finite particle size
journal, November 2020
- Stanier, A.; Chacón, L.; Le, A.
- Journal of Computational Physics, Vol. 420
A multidimensional quasineutral plasma simulation model
journal, November 1978
- Hewett, D. W.; Nielson, C. W.
- Journal of Computational Physics, Vol. 29, Issue 2
Use of a hybrid code to model the Earth's magnetosphere
journal, February 1995
- Swift, Daniel W.
- Geophysical Research Letters, Vol. 22, Issue 3
A Flexible Inner-Outer Preconditioned GMRES Algorithm
journal, March 1993
- Saad, Youcef
- SIAM Journal on Scientific Computing, Vol. 14, Issue 2
An unsupervised machine-learning checkpoint-restart algorithm using Gaussian mixtures for particle-in-cell simulations
journal, July 2021
- Chen, G.; Chacón, L.; Nguyen, T. B.
- Journal of Computational Physics, Vol. 436
MHD-kinetic hybrid code based on structure-preserving finite elements with particles-in-cell
journal, May 2021
- Holderied, Florian; Possanner, Stefan; Wang, Xin
- Journal of Computational Physics, Vol. 433
The role of guide field in magnetic reconnection driven by island coalescence
journal, February 2017
- Stanier, A.; Daughton, W.; Simakov, Andrei N.
- Physics of Plasmas, Vol. 24, Issue 2
Quasineutral hybrid simulation of macroscopic plasma phenomena
journal, September 1982
- Harned, Douglas S.
- Journal of Computational Physics, Vol. 47, Issue 3
A generalized quasi-neutral fluid-particle hybrid plasma model and its application to energetic-particle-magnetohydrodynamics hybrid simulation
journal, August 2018
- Amano, Takanobu
- Journal of Computational Physics, Vol. 366
Flux Pileup in Collisionless Magnetic Reconnection: Bursty Interaction of Large Flux Ropes
journal, July 2011
- Karimabadi, H.; Dorelli, J.; Roytershteyn, V.
- Physical Review Letters, Vol. 107, Issue 2
Pegasus: A new hybrid-kinetic particle-in-cell code for astrophysical plasma dynamics
journal, February 2014
- Kunz, Matthew W.; Stone, James M.; Bai, Xue-Ning
- Journal of Computational Physics, Vol. 259