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

Chen, G.; Chacón, L.

doi:10.1016/j.cpc.2015.08.008

Title: 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

DOI:https://doi.org/10.1016/j.cpc.2015.08.008· OSTI ID:1329577

Chen, G. ^[1]; Chacón, L. ^[1]

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

For decades, the Vlasov–Darwin model has been recognized to be attractive for particle-in-cell (PIC) kinetic plasma simulations in non-radiative electromagnetic regimes, to avoid radiative noise issues and gain computational efficiency. However, the Darwin model results in an elliptic set of field equations that renders conventional explicit time integration unconditionally unstable. We explore a fully implicit PIC algorithm for the Vlasov–Darwin model in multiple dimensions, which overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. The finite-difference scheme for Darwin field equations and particle equations of motion is space–time-centered, employing particle sub-cycling and orbit-averaging. This algorithm conserves total energy, local charge, canonical-momentum in the ignorable direction, and preserves the Coulomb gauge exactly. An asymptotically well-posed fluid preconditioner allows efficient use of large cell sizes, which are determined by accuracy considerations, not stability, and can be orders of magnitude larger than required in a standard explicit electromagnetic PIC simulation. Finally, we demonstrate the accuracy and efficiency properties of the algorithm with various numerical experiments in 2D–3V.

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 National Nuclear Security Administration (NNSA)

Grant/Contract Number:: AC52-06NA25396

OSTI ID:: 1329577

Alternate ID(s):: OSTI ID: 1246662

Report Number(s):: LA-UR-15-21463

Journal Information:: Computer Physics Communications, Vol. 197, Issue C; ISSN 0010-4655

Publisher:: ElsevierCopyright Statement

Country of Publication:: United States

Language:: English

Citation Metrics:

Cited by: 56 works

Citation information provided by
Web of Science

References (44)

Numerical Cherenkov instabilities in electromagnetic particle codes Godfrey, Brendan B. Journal of Computational Physics, Vol. 15, Issue 4 https://doi.org/10.1016/0021-9991(74)90076-X	journal	August 1974
Some Electromagnetic Plasma Simulation Methods and Their Noise Properties Langdon, A. Bruce Physics of Fluids, Vol. 15, Issue 6 https://doi.org/10.1063/1.1694042	journal	January 1972
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
An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov–Darwin particle-in-cell algorithm Chen, G.; Chacón, L. Computer Physics Communications, Vol. 185, Issue 10 https://doi.org/10.1016/j.cpc.2014.05.010	journal	October 2014
A self-consistent magnetostatic particle code for numerical simulation of plasmas Busnardo-Neto, J.; Pritchett, P. L.; Lin, A. T. Journal of Computational Physics, Vol. 23, Issue 3 https://doi.org/10.1016/0021-9991(77)90096-1	journal	March 1977
Hybrid simulations of quasineutral phenomena in magnetized plasma Byers, J. A.; Cohen, B. I.; Condit, W. C. Journal of Computational Physics, Vol. 27, Issue 3 https://doi.org/10.1016/0021-9991(78)90016-5	journal	June 1978
Low-frequency electromagnetic (Darwin) applications in plasma simulation Hewett, D. W. Computer Physics Communications, Vol. 84, Issue 1-3 https://doi.org/10.1016/0010-4655(94)90214-3	journal	November 1994
The Darwin Direct Implicit Particle-in-Cell (DADIPIC) Method for Simulation of Low Frequency Plasma Phenomena Gibbons, Matthew R.; Hewett, Dennis W. Journal of Computational Physics, Vol. 120, Issue 2 https://doi.org/10.1006/jcph.1995.1160	journal	September 1995
Nonlinear perturbative electromagnetic (Darwin) particle simulation of high intensity beams Lee, W. Wei-li; Qin, Hong; Davidson, Ronald C. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Vol. 464, Issue 1-3 https://doi.org/10.1016/S0168-9002(01)00109-7	journal	May 2001
Study of hot electron beam transport in high density plasma using 3D hybrid-Darwin code Taguchi, T.; Antonsen, T. M.; Mima, K. Computer Physics Communications, Vol. 164, Issue 1-3 https://doi.org/10.1016/j.cpc.2004.06.038	journal	December 2004
Simulations of electromagnetic effects in high-frequency capacitively coupled discharges using the Darwin approximation Eremin, Denis; Hemke, Torben; Brinkmann, Ralf Peter Journal of Physics D: Applied Physics, Vol. 46, Issue 8 https://doi.org/10.1088/0022-3727/46/8/084017	journal	February 2013
Darwin–Vlasov simulations of magnetised plasmas Schmitz, H.; Grauer, R. Journal of Computational Physics, Vol. 214, Issue 2 https://doi.org/10.1016/j.jcp.2005.10.013	journal	May 2006
Boundary conditions for the Darwin model Weitzner, H.; Lawson, W. S. Physics of Fluids B: Plasma Physics, Vol. 1, Issue 10 https://doi.org/10.1063/1.859058	journal	October 1989
An analysis of the Darwin model of approximation to Maxwell’s equations Degond, Pierre; Raviart, Pierre-Arnaud Forum Mathematicum, Vol. 4, Issue 4 https://doi.org/10.1515/form.1992.4.13	journal	January 1992
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
Development of a Consistent and Stable Fully Implicit Moment Method for Vlasov--Ampère Particle in Cell (PIC) System Taitano, William T.; Knoll, Dana A.; Chacón, Luis SIAM Journal on Scientific Computing, Vol. 35, Issue 5 https://doi.org/10.1137/120881385	journal	January 2013
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
Fluid preconditioning for Newton–Krylov-based, fully implicit, electrostatic particle-in-cell simulations Chen, G.; Chacón, L.; Leibs, C. A. Journal of Computational Physics, Vol. 258 https://doi.org/10.1016/j.jcp.2013.10.052	journal	February 2014
An efficient mixed-precision, hybrid CPU–GPU implementation of a nonlinearly implicit one-dimensional particle-in-cell algorithm Chen, G.; Chacón, L.; Barnes, D. C. Journal of Computational Physics, Vol. 231, Issue 16 https://doi.org/10.1016/j.jcp.2012.04.040	journal	June 2012
Implicit moment particle simulation of plasmas Mason, Rodney J. Journal of Computational Physics, Vol. 41, Issue 2 https://doi.org/10.1016/0021-9991(81)90094-2	journal	June 1981
Time-filtering particle simulations with ωpe Δt ⪢ 1 Denavit, J. Journal of Computational Physics, Vol. 42, Issue 2 https://doi.org/10.1016/0021-9991(81)90249-7	journal	August 1981
An implicit method for electromagnetic plasma simulation in two dimensions Brackbill, J. U.; Forslund, D. W. Journal of Computational Physics, Vol. 46, Issue 2 https://doi.org/10.1016/0021-9991(82)90016-X	journal	May 1982
CELEST1D: an implicit, fully kinetic model for low-frequency, electromagnetic plasma simulation Vu, H. X.; Brackbill, J. U. Computer Physics Communications, Vol. 69, Issue 2-3 https://doi.org/10.1016/0010-4655(92)90165-U	journal	March 1992
Direct implicit large time-step particle simulation of plasmas Langdon, A. Bruce; Cohen, Bruce I.; Friedman, Alex Journal of Computational Physics, Vol. 51, Issue 1 https://doi.org/10.1016/0021-9991(83)90083-9	journal	July 1983
Electromagnetic direct implicit plasma simulation Hewett, Dennis W.; Bruce Langdon, A. Journal of Computational Physics, Vol. 72, Issue 1 https://doi.org/10.1016/0021-9991(87)90075-1	journal	September 1987
Implicit particle simulation of electromagnetic plasma phenomena Kamimura, T.; Montalvo, E.; Barnes, D. C. Journal of Computational Physics, Vol. 100, Issue 1 https://doi.org/10.1016/0021-9991(92)90311-L	journal	May 1992
Elimination of electromagnetic radiation in plasma simulation: The Darwin or magneto inductive approximation Hewett, Dennis W. Space Science Reviews, Vol. 42, Issue 1-2 https://doi.org/10.1007/BF00218221	journal	October 1985
A hierarchy of approximate models for the Maxwell equations Raviart, Pierre-Arnaud; Sonnendrücker, Eric Numerische Mathematik, Vol. 73, Issue 3 https://doi.org/10.1007/s002110050196	journal	May 1996
A unified approach to the Darwin approximation Krause, Todd B.; Apte, A.; Morrison, P. J. Physics of Plasmas, Vol. 14, Issue 10 https://doi.org/10.1063/1.2799346	journal	October 2007
Hyperbolic Conservation Laws in Continuum Physics Dafermos, Constantine M. Grundlehren der mathematischen Wissenschaften https://doi.org/10.1007/978-3-662-22019-1	book	January 2000
The Origin of Spurious Solutions in Computational Electromagnetics Jiang, Bo-nan; Wu, Jie; Povinelli, L. A. Journal of Computational Physics, Vol. 125, Issue 1 https://doi.org/10.1006/jcph.1996.0082	journal	April 1996
One-Dimensional Plasma Model in the Presence of a Magnetic Field Hasegawa, Akira Physics of Fluids, Vol. 11, Issue 9 https://doi.org/10.1063/1.1692232	journal	January 1968
Divergence Correction Techniques for Maxwell Solvers Based on a Hyperbolic Model Munz, C. -D.; Omnes, P.; Schneider, R. Journal of Computational Physics, Vol. 161, Issue 2 https://doi.org/10.1006/jcph.2000.6507	journal	July 2000
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
An analytical particle mover for the charge- and energy-conserving, nonlinearly implicit, electrostatic particle-in-cell algorithm Chen, G.; Chacón, L. Journal of Computational Physics, Vol. 247 https://doi.org/10.1016/j.jcp.2013.04.002	journal	August 2013
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
The Darwin Model as a Tool for Electromagnetic Plasma Simulation Kaufman, Allan N. Physics of Fluids, Vol. 14, Issue 2 https://doi.org/10.1063/1.1693451	journal	January 1971
Nested Iteration and First-Order Systems Least Squares for a Two-Fluid Electromagnetic Darwin Model Leibs, C. A.; Manteuffel, T. A. SIAM Journal on Scientific Computing, Vol. 37, Issue 5 https://doi.org/10.1137/14097793X	journal	January 2015
Kinetic Alfvén waves and electron physics. I. Generation from ion-ion streaming Yin, L.; Winske, D.; Daughton, W. Physics of Plasmas, Vol. 14, Issue 6 https://doi.org/10.1063/1.2734950	journal	June 2007
An implicit “drift-Lorentz” mover for plasma and beam simulations Cohen, R. H.; Friedman, A.; Grote, D. P. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Vol. 606, Issue 1-2 https://doi.org/10.1016/j.nima.2009.03.083	journal	July 2009
A charge- and energy-conserving implicit, electrostatic particle-in-cell algorithm on mapped computational meshes Chacón, L.; Chen, G.; Barnes, D. C. Journal of Computational Physics, Vol. 233 https://doi.org/10.1016/j.jcp.2012.07.042	journal	January 2013
Numerical Simulation of the Weibel Instability in One and Two Dimensions Morse, R. L. Physics of Fluids, Vol. 14, Issue 4 https://doi.org/10.1063/1.1693518	journal	January 1971
A finite element formulation of the Darwin PIC model for use on unstructured grids Sonnendrücker, Eric; Ambrosiano, John J.; Brandon, Scott T. Journal of Computational Physics, Vol. 121, Issue 2 https://doi.org/10.1016/s0021-9991(95)90119-1	journal	October 1995
Assessment and Intervention Strategies for the Bilingual Language-Disordered Student Langdon, Henriette W. Exceptional Children, Vol. 50, Issue 1 https://doi.org/10.1177/001440298305000105	journal	September 1983

Cited By (10)

Vlasov-Maxwell system with physical reasonableness constraints Lin, H.; Liu, C. P. Plasma Research Express, Vol. 1, Issue 2 https://doi.org/10.1088/2516-1067/ab19f1	journal	May 2019
ViDA: a Vlasov–DArwin solver for plasma physics at electron scales Pezzi, Oreste; Cozzani, Giulia; Califano, Francesco Journal of Plasma Physics, Vol. 85, Issue 5 https://doi.org/10.1017/s0022377819000631	journal	October 2019
Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code Hager, Robert; Lang, Jianying; Chang, C. S. Physics of Plasmas, Vol. 24, Issue 5 https://doi.org/10.1063/1.4983320	journal	May 2017
Interpolation-free particle simulation Lin, H.; Liu, C. P. Laser and Particle Beams, Vol. 38, Issue 1 https://doi.org/10.1017/s0263034619000806	journal	January 2020
A Post-Newtonian Expansion Including Radiation Damping for a Collisionless Plasma Bauer, Sebastian Journal of Nonlinear Science, Vol. 30, Issue 1 https://doi.org/10.1007/s00332-019-09580-1	journal	September 2019
Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration Ghosh, Debojyoti; Dorf, Mikhail A.; Dorr, Milo R. Journal of Scientific Computing, Vol. 77, Issue 2 https://doi.org/10.1007/s10915-018-0726-6	journal	May 2018
A Semi-implicit Particle-in-cell Expanding Box Model Code for Fully Kinetic Simulations of the Expanding Solar Wind Plasma Elena Innocenti, Maria; Tenerani, Anna; Velli, Marco The Astrophysical Journal, Vol. 870, Issue 2 https://doi.org/10.3847/1538-4357/aaf1be	journal	January 2019
A Post-Newtonian Expansion Including Radiation Damping for a Collisionless Plasma Bauer, Sebastian Unpublished https://doi.org/10.13140/rg.2.2.23041.10087	text	January 2019
Kinetic Simulation of Collisional Magnetized Plasmas with Semi-Implicit Time Integration Ghosh, Debojyoti; Dorf, Mikhail A.; Dorr, Milo R. arXiv https://doi.org/10.48550/arxiv.1707.08247	preprint	January 2017
Electromagnetic full-$f$ gyrokinetics in the tokamak edge with discontinuous Galerkin methods Mandell, N. R.; Hakim, A.; Hammett, G. W. arXiv https://doi.org/10.48550/arxiv.1908.05653	text	January 2019

Similar Records

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:1329577

Chacon, L.; Chen, G.

Low frequency, electrodynamic simulation of kinetic plasmas with the DArwin Direct Implicit Particle-In-Cell (DADIPIC) method

Thesis/Dissertation · Thu Jun 01 00:00:00 EDT 1995 · OSTI ID:1329577

Gibbons, Matthew Richard

Fully implicit Particle-in-cell algorithms for multiscale plasma simulation

Conference · Thu Jul 16 00:00:00 EDT 2015 · OSTI ID:1329577

Chacon, Luis

Related Subjects

97 MATHEMATICS AND COMPUTING
70 PLASMA PHYSICS AND FUSION TECHNOLOGY
Implicit particle-in-cell
Energy conservation
Charge conservation
Canonical momentum conservation
Vlasov–Darwin
Multi-scale
JFNK
Physics-based preconditioner

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

Citation Formats

References (44)

Cited By (10)

Similar Records

Related Subjects