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

Stanier, A.; Chacón, L.; Chen, G.

doi:10.1016/j.jcp.2018.09.038

Title: A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm

Journal Article · Thu Sep 20 00:00:00 EDT 2018 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2018.09.038· OSTI ID:1477697

^[1];

^[1]; Chen, G. ^[1]

Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

The quasi-neutral hybrid model with kinetic ions and fluid electrons is a promising approach for bridging the inherent multi-scale nature of many problems in space and laboratory plasmas. A novel, implicit, particle-in-cell based scheme for the hybrid model is derived for fully 3D electromagnetic problems with multiple ion species, which features global mass, momentum and energy conservation. The scheme includes sub-cycling and orbit-averaging for the ions, with cell-centered finite differences and implicit midpoint time advance. To reduce discrete particle noise, the scheme allows arbitrary-order shape functions for the particle-mesh interpolations and the application of conservative binomial smoothing. The algorithm is verified for a number of test problems to demonstrate the correctness of the implementation, the unique conservation properties, and the favorable stability properties of the new scheme. In particular, there is no indication of unstable growth of the finite-grid instability for a population of cold ions drifting through a uniform spatial mesh, in a set-up where several commonly used non-conservative schemes are highly unstable.

View Accepted Manuscript (DOE)

View Accepted Manuscript (Publisher)

Cite

Export

Save

Research Organization:: Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)

Sponsoring Organization:: USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); USDOE National Nuclear Security Administration (NNSA)

Grant/Contract Number:: AC52-06NA25396

OSTI ID:: 1477697

Alternate ID(s):: OSTI ID: 1756197

Report Number(s):: LA-UR-18-21708

Journal Information:: Journal of Computational Physics, Vol. 376; ISSN 0021-9991

Publisher:: ElsevierCopyright Statement

Country of Publication:: United States

Language:: English

Citation Metrics:

Cited by: 11 works

Citation information provided by
Web of Science

References (56)

The link between shocks, turbulence, and magnetic reconnection in collisionless plasmas Karimabadi, H.; Roytershteyn, V.; Vu, H. X. Physics of Plasmas, Vol. 21, Issue 6 https://doi.org/10.1063/1.4882875	journal	June 2014
Role of Ion Kinetic Physics in the Interaction of Magnetic Flux Ropes Stanier, A.; Daughton, W.; Chacón, L. Physical Review Letters, Vol. 115, Issue 17 https://doi.org/10.1103/PhysRevLett.115.175004	journal	October 2015
Low frequency fully kinetic simulation of the toroidal ion temperature gradient instability Sturdevant, B. J.; Chen, Y.; Parker, S. E. Physics of Plasmas, Vol. 24, Issue 8 https://doi.org/10.1063/1.4999945	journal	August 2017
Quasineutral hybrid simulation of macroscopic plasma phenomena Harned, Douglas S. Journal of Computational Physics, Vol. 47, Issue 3 https://doi.org/10.1016/0021-9991(82)90094-8	journal	September 1982
Electromagnetic ion beam instabilities: Comparison of one- and two-dimensional simulations Winske, D.; Quest, K. B. Journal of Geophysical Research, Vol. 91, Issue A8 https://doi.org/10.1029/JA091iA08p08789	journal	January 1986
Pegasus: A new hybrid-kinetic particle-in-cell code for astrophysical plasma dynamics Kunz, Matthew W.; Stone, James M.; Bai, Xue-Ning Journal of Computational Physics, Vol. 259 https://doi.org/10.1016/j.jcp.2013.11.035	journal	February 2014
Current Advance Method and Cyclic Leapfrog for 2D Multispecies Hybrid Plasma Simulations Matthews, Alan P. Journal of Computational Physics, Vol. 112, Issue 1 https://doi.org/10.1006/jcph.1994.1084	journal	May 1994
Re-forming supercritical quasi-parallel shocks: 1. One- and two-dimensional simulations Thomas, V. A.; Winske, D.; Omidi, N. Journal of Geophysical Research, Vol. 95, Issue A11 https://doi.org/10.1029/JA095iA11p18809	journal	January 1990
Magnetic self-organization in Hall-dominated magnetorotational turbulence Kunz, Matthew W.; Lesur, Geoffroy Monthly Notices of the Royal Astronomical Society, Vol. 434, Issue 3 https://doi.org/10.1093/mnras/stt1171	journal	July 2013
Use of a hybrid code to model the Earth's magnetosphere Swift, Daniel W. Geophysical Research Letters, Vol. 22, Issue 3 https://doi.org/10.1029/94GL03082	journal	February 1995
A global method of solving the electron-field equations in a zero-inertia-electron-hybrid plasma simulation code Hewett, D. W. Journal of Computational Physics, Vol. 38, Issue 3 https://doi.org/10.1016/0021-9991(80)90155-2	journal	December 1980
A robust method for handling low density regions in hybrid simulations for collisionless plasmas Amano, Takanobu; Higashimori, Katsuaki; Shirakawa, Keisuke Journal of Computational Physics, Vol. 275 https://doi.org/10.1016/j.jcp.2014.06.048	journal	October 2014
A second-order semi-implicit method for hybrid simulation Cheng, Jianhua; Parker, Scott E.; Chen, Yang Journal of Computational Physics, Vol. 245 https://doi.org/10.1016/j.jcp.2013.03.017	journal	July 2013
An implicit δf particle-in-cell method with sub-cycling and orbit averaging for Lorentz ions Sturdevant, Benjamin J.; Parker, Scott E.; Chen, Yang Journal of Computational Physics, Vol. 316 https://doi.org/10.1016/j.jcp.2016.04.036	journal	July 2016
An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm Chen, G.; Chacón, L.; Barnes, D. C. Journal of Computational Physics, Vol. 230, Issue 18 https://doi.org/10.1016/j.jcp.2011.05.031	journal	August 2011
The energy conserving particle-in-cell method Markidis, Stefano; Lapenta, Giovanni Journal of Computational Physics, Vol. 230, Issue 18 https://doi.org/10.1016/j.jcp.2011.05.033	journal	August 2011
Energy-conserving numerical approximations for Vlasov plasmas Lewis, H. Ralph Journal of Computational Physics, Vol. 6, Issue 1 https://doi.org/10.1016/0021-9991(70)90012-4	journal	August 1970
A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm Chen, G.; Chacón, L. Computer Physics Communications, Vol. 197 https://doi.org/10.1016/j.cpc.2015.08.008	journal	December 2015
A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions Chacón, L.; Chen, G. Journal of Computational Physics, Vol. 316 https://doi.org/10.1016/j.jcp.2016.03.070	journal	July 2016
Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme Squire, J.; Qin, H.; Tang, W. M. Physics of Plasmas, Vol. 19, Issue 8 https://doi.org/10.1063/1.4742985	journal	August 2012
Variational formulation of particle algorithms for kinetic plasma simulations Evstatiev, E. G.; Shadwick, B. A. Journal of Computational Physics, Vol. 245 https://doi.org/10.1016/j.jcp.2013.03.006	journal	July 2013
Hamiltonian splitting for the Vlasov–Maxwell equations Crouseilles, Nicolas; Einkemmer, Lukas; Faou, Erwan Journal of Computational Physics, Vol. 283 https://doi.org/10.1016/j.jcp.2014.11.029	journal	February 2015
Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov–Maxwell equations Qin, Hong; Liu, Jian; Xiao, Jianyuan Nuclear Fusion, Vol. 56, Issue 1 https://doi.org/10.1088/0029-5515/56/1/014001	journal	December 2015
GEMPIC: geometric electromagnetic particle-in-cell methods Kraus, Michael; Kormann, Katharina; Morrison, Philip J. Journal of Plasma Physics, Vol. 83, Issue 4 https://doi.org/10.1017/S002237781700040X	journal	July 2017
Nonphysical noises and instabilities in plasma simulation due to a spatial grid Okuda, Hideo Journal of Computational Physics, Vol. 10, Issue 3 https://doi.org/10.1016/0021-9991(72)90048-4	journal	December 1972
On energy and momentum conservation in particle-in-cell plasma simulation Brackbill, J. U. Journal of Computational Physics, Vol. 317 https://doi.org/10.1016/j.jcp.2016.04.050	journal	July 2016
Gridless finite-size-particle plasma simulation Vlad, G.; Briguglio, S.; Fogaccia, G. Computer Physics Communications, Vol. 134, Issue 1 https://doi.org/10.1016/S0010-4655(00)00191-0	journal	February 2001
Finite grid instability and spectral fidelity of the electrostatic Particle-In-Cell algorithm Huang, C. -K.; Zeng, Y.; Wang, Y. Computer Physics Communications, Vol. 207 https://doi.org/10.1016/j.cpc.2016.05.021	journal	October 2016
Finite-Grid Instability in Quasineutral Hybrid Simulations Rambo, P. W. Journal of Computational Physics, Vol. 118, Issue 1 https://doi.org/10.1006/jcph.1995.1086	journal	April 1995
Hybrid simulations of magnetic reconnection with kinetic ions and fluid electron pressure anisotropy Le, A.; Daughton, W.; Karimabadi, H. Physics of Plasmas, Vol. 23, Issue 3 https://doi.org/10.1063/1.4943893	journal	March 2016
Demonstration of Anisotropic Fluid Closure Capturing the Kinetic Structure of Magnetic Reconnection Ohia, O.; Egedal, J.; Lukin, V. S. Physical Review Letters, Vol. 109, Issue 11 https://doi.org/10.1103/PhysRevLett.109.115004	journal	September 2012
The island coalescence problem: Scaling of reconnection in extended fluid models including higher-order moments Ng, Jonathan; Huang, Yi-Min; Hakim, Ammar Physics of Plasmas, Vol. 22, Issue 11 https://doi.org/10.1063/1.4935302	journal	November 2015
A non-staggered, conservative, , finite-volume scheme for 3D implicit extended magnetohydrodynamics in curvilinear geometries Chacón, L. Computer Physics Communications, Vol. 163, Issue 3 https://doi.org/10.1016/j.cpc.2004.08.005	journal	November 2004
Rigorous charge conservation for local electromagnetic field solvers Villasenor, John; Buneman, Oscar Computer Physics Communications, Vol. 69, Issue 2-3 https://doi.org/10.1016/0010-4655(92)90169-Y	journal	March 1992
Charge conservation in electromagnetic PIC codes; spectral comparison of Boris/DADI and Langdon-Marder methods Mardahl, P. J.; Verboncoeur, J. P. Computer Physics Communications, Vol. 106, Issue 3 https://doi.org/10.1016/S0010-4655(97)00094-5	journal	November 1997
Exact charge conservation scheme for Particle-in-Cell simulation with an arbitrary form-factor Esirkepov, T. Zh. Computer Physics Communications, Vol. 135, Issue 2 https://doi.org/10.1016/S0010-4655(00)00228-9	journal	April 2001
Jacobian-free Newton–Krylov methods: a survey of approaches and applications Knoll, D. A.; Keyes, D. E. Journal of Computational Physics, Vol. 193, Issue 2 https://doi.org/10.1016/j.jcp.2003.08.010	journal	January 2004
An optimal, parallel, fully implicit Newton–Krylov solver for three-dimensional viscoresistive magnetohydrodynamics Chacón, L. Physics of Plasmas, Vol. 15, Issue 5 https://doi.org/10.1063/1.2838244	journal	May 2008
A Flexible Inner-Outer Preconditioned GMRES Algorithm Saad, Youcef SIAM Journal on Scientific Computing, Vol. 14, Issue 2 https://doi.org/10.1137/0914028	journal	March 1993
Orbit-averaged implicit particle codes Cohen, Bruce I.; Freis, Robert P.; Thomas, Vincent Journal of Computational Physics, Vol. 45, Issue 3 https://doi.org/10.1016/0021-9991(82)90108-5	journal	March 1982
Accurate Numerical Solution of Charged Particle Motion in a Magnetic Field Vu, H. X.; Brackbill, J. U. Journal of Computational Physics, Vol. 116, Issue 2 https://doi.org/10.1006/jcph.1995.1037	journal	February 1995
Heuristic stability theory for finite-difference equations Hirt, C. W. Journal of Computational Physics, Vol. 2, Issue 4 https://doi.org/10.1016/0021-9991(68)90041-7	journal	June 1968
Curvature-compensated convective transport: SMART, A new boundedness- preserving transport algorithm Gaskell, P. H.; Lau, A. K. C. International Journal for Numerical Methods in Fluids, Vol. 8, Issue 6 https://doi.org/10.1002/fld.1650080602	journal	June 1988
A linear dispersion relation for the hybrid kinetic-ion/fluid-electron model of plasma physics Told, D.; Cookmeyer, J.; Astfalk, P. New Journal of Physics, Vol. 18, Issue 7 https://doi.org/10.1088/1367-2630/18/7/075001	journal	July 2016
Numerical simulation of nonoscillatory mirror waves at the Earth's magnetosheath Price, Channon P.; Swift, Daniel W.; Lee, Lou-Chuang Journal of Geophysical Research, Vol. 91, Issue A1 https://doi.org/10.1029/JA091iA01p00101	journal	January 1986
Ion anisotropy instabilities in the magnetosheath Gary, S. Peter; Fuselier, Stephen A.; Anderson, Brian J. Journal of Geophysical Research: Space Physics, Vol. 98, Issue A2 https://doi.org/10.1029/92JA01844	journal	February 1993
Observation of mirror waves downstream of a quasi-perpendicular shock Hubert, D.; Perche, C.; Harvey, C. C. Geophysical Research Letters, Vol. 16, Issue 2 https://doi.org/10.1029/GL016i002p00159	journal	February 1989
Geospace Environmental Modeling (GEM) Magnetic Reconnection Challenge Birn, J.; Drake, J. F.; Shay, M. A. Journal of Geophysical Research: Space Physics, Vol. 106, Issue A3 https://doi.org/10.1029/1999JA900449	journal	March 2001
A fully nonlinear characteristic method for gyrokinetic simulation Parker, S. E.; Lee, W. W. Physics of Fluids B: Plasma Physics, Vol. 5, Issue 1 https://doi.org/10.1063/1.860870	journal	January 1993
A Particle-in-cell Method with Adaptive Phase-space Remapping for Kinetic Plasmas Wang, B.; Miller, G. H.; Colella, P. SIAM Journal on Scientific Computing, Vol. 33, Issue 6 https://doi.org/10.1137/100811805	journal	January 2011
Current sheet bifurcation and collapse in electron magnetohydrodynamics Zocco, A.; Chacón, L.; Simakov, Andrei N. Physics of Plasmas, Vol. 16, Issue 11 https://doi.org/10.1063/1.3264102	journal	November 2009
Fluid vs. kinetic magnetic reconnection with strong guide fields Stanier, A.; Simakov, Andrei N.; Chacón, L. Physics of Plasmas, Vol. 22, Issue 10 https://doi.org/10.1063/1.4932330	journal	October 2015
A 2D high-β Hall MHD implicit nonlinear solver Chacón, L.; Knoll, D. A. Journal of Computational Physics, Vol. 188, Issue 2 https://doi.org/10.1016/S0021-9991(03)00193-1	journal	July 2003
On magnetic reconnection regimes and associated three-dimensional asymmetries: Hybrid, Hall-less hybrid, and Hall-MHD simulations Karimabadi, H. Journal of Geophysical Research, Vol. 109, Issue A9 https://doi.org/10.1029/2004ja010478	journal	January 2004
Magnetic field and density fluctuations at perpendicular supercritical collisionless shocks Winske, D.; Quest, K. B. Journal of Geophysical Research, Vol. 93, Issue A9 https://doi.org/10.1029/ja093ia09p09681	journal	January 1988
A Robust Method for Handling Low Density Regions in Hybrid Simulations for Collisionless Plasmas Amano, Takanobu; Higashimori, Katsuaki; Shirakawa, Keisuke arXiv https://doi.org/10.48550/arxiv.1406.6613	text	January 2014

Cited By (1)

PolyPIC: the Polymorphic-Particle-in-Cell Method for Fluid-Kinetic Coupling Markidis, Stefano; Olshevsky, Vyacheslav; Sishtla, Chaitanya Prasad arXiv https://doi.org/10.48550/arxiv.1807.05183	text	January 2018

Figures / Tables (5)

Similar Records

A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes

Journal Article · Tue Mar 15 00:00:00 EDT 2022 · Journal of Computational Physics · OSTI ID:1477697

Stanier, Adam John; Chacon, Luis

A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm

Journal Article · Tue Aug 11 00:00:00 EDT 2015 · Computer Physics Communications · OSTI ID:1477697

Chen, G.; Chacón, L.

A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions

Journal Article · Tue Apr 19 00:00:00 EDT 2016 · Journal of Computational Physics · OSTI ID:1477697

Chacon, L.; Chen, G.

Related Subjects

97 MATHEMATICS AND COMPUTING
70 PLASMA PHYSICS AND FUSION TECHNOLOGY
hybrid
plasma
particle-in-cell
implicit
conservative
finite grid instability