DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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:
ORCiD logo [1]; ORCiD logo [1]
  1. 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
  • DOI: 10.1063/1.4999945

Finite-Grid Instability in Quasineutral Hybrid Simulations
journal, April 1995


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
  • DOI: 10.1006/jcph.1997.5738

High-Resolution Hybrid Simulations of Kinetic Plasma Turbulence at Proton Scales
journal, October 2015


A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm
journal, January 2019


Role of Ion Kinetic Physics in the Interaction of Magnetic Flux Ropes
journal, October 2015


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
  • DOI: 10.1029/JA095iA11p18809

Controlling self-force errors at refinement boundaries for AMR-PIC
journal, February 2010


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
  • DOI: 10.1063/1.4932330

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
  • DOI: 10.1063/1.5062853

Current Advance Method and Cyclic Leapfrog for 2D Multispecies Hybrid Plasma Simulations
journal, May 1994


Jacobian-free Newton–Krylov methods: a survey of approaches and applications
journal, January 2004


Temperature gradient driven heat flux closure in fluid simulations of collisionless reconnection
journal, June 2018


Two-way coupled MHD-PIC simulations of magnetic reconnection in magnetic island coalescence
journal, May 2018


Magnetic field and density fluctuations at perpendicular supercritical collisionless shocks
journal, January 1988


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
  • DOI: 10.1063/1.5058140

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
  • DOI: 10.1029/2021JA029125

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
  • DOI: 10.1002/jgra.50497

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
  • DOI: 10.1063/1.4882875

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
  • DOI: 10.1016/j.cpc.2010.12.033

Subproton-Scale Cascades in Solar wind Turbulence: Driven Hybrid-Kinetic Simulations
journal, April 2016


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
  • DOI: 10.1063/1.4935302

Self-focusing of local plasma currents
journal, September 1965


Low-noise electromagnetic and relativistic particle-in-cell plasma simulation models
journal, September 1999


The coalescence instability in collisionless plasmas
journal, October 1992

  • Pritchett, P. L.
  • Physics of Fluids B: Plasma Physics, Vol. 4, Issue 10
  • DOI: 10.1063/1.860392

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
  • DOI: 10.1002/2017JA024180

Hybrid simulation codes with application to shocks and upstream waves
journal, October 1985


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
  • DOI: 10.1063/1.4971231

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
  • DOI: 10.1016/j.jcp.2014.06.048

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
  • DOI: 10.1063/1.5125014

Origin of Mercury’s double magnetopause: 3D hybrid simulation study with A.I.K.E.F.
journal, March 2012


A global method of solving the electron-field equations in a zero-inertia-electron-hybrid plasma simulation code
journal, December 1980


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
  • DOI: 10.1063/1.3026714

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
  • DOI: 10.1063/5.0054394

Electromagnetic ion beam instabilities: Comparison of one- and two-dimensional simulations
journal, January 1986


Hybrid simulations of quasineutral phenomena in magnetized plasma
journal, June 1978


Coalescence instability of magnetic islands
journal, January 1977

  • Finn, J. M.; Kaw, P. K.
  • Physics of Fluids, Vol. 20, Issue 1
  • DOI: 10.1063/1.861709

Magnetically Self-Focussing Streams
journal, June 1934


Global gyrokinetic particle-in-cell simulations of internal kink instabilities
journal, December 2012

  • Mishchenko, Alexey; Zocco, Alessandro
  • Physics of Plasmas, Vol. 19, Issue 12
  • DOI: 10.1063/1.4769379

Orbit-averaged implicit particle codes
journal, March 1982


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
  • DOI: 10.1063/1.5110512

Heuristic stability theory for finite-difference equations
journal, June 1968


A cancellation problem in hybrid particle-in-cell schemes due to finite particle size
journal, November 2020


A multidimensional quasineutral plasma simulation model
journal, November 1978


Use of a hybrid code to model the Earth's magnetosphere
journal, February 1995

  • Swift, Daniel W.
  • Geophysical Research Letters, Vol. 22, Issue 3
  • DOI: 10.1029/94GL03082

A Flexible Inner-Outer Preconditioned GMRES Algorithm
journal, March 1993

  • Saad, Youcef
  • SIAM Journal on Scientific Computing, Vol. 14, Issue 2
  • DOI: 10.1137/0914028

An unsupervised machine-learning checkpoint-restart algorithm using Gaussian mixtures for particle-in-cell simulations
journal, July 2021


MHD-kinetic hybrid code based on structure-preserving finite elements with particles-in-cell
journal, May 2021


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
  • DOI: 10.1063/1.4976712

Quasineutral hybrid simulation of macroscopic plasma phenomena
journal, September 1982


Flux Pileup in Collisionless Magnetic Reconnection: Bursty Interaction of Large Flux Ropes
journal, July 2011


Pegasus: A new hybrid-kinetic particle-in-cell code for astrophysical plasma dynamics
journal, February 2014