A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm
Abstract
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.
- Authors:
-
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States). Theoretical Division
- Publication Date:
- Research Org.:
- Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
- Sponsoring Org.:
- USDOE National Nuclear Security Administration (NNSA)
- OSTI Identifier:
- 1329577
- Alternate Identifier(s):
- OSTI ID: 1246662
- Report Number(s):
- LA-UR-15-21463
Journal ID: ISSN 0010-4655
- Grant/Contract Number:
- AC52-06NA25396
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Computer Physics Communications
- Additional Journal Information:
- Journal Volume: 197; Journal Issue: C; Journal ID: ISSN 0010-4655
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 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
Citation Formats
Chen, G., and Chacón, L. A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm. United States: N. p., 2015.
Web. doi:10.1016/j.cpc.2015.08.008.
Chen, G., & Chacón, L. A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm. United States. https://doi.org/10.1016/j.cpc.2015.08.008
Chen, G., and Chacón, L. Tue .
"A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm". United States. https://doi.org/10.1016/j.cpc.2015.08.008. https://www.osti.gov/servlets/purl/1329577.
@article{osti_1329577,
title = {A multi-dimensional, energy- and charge-conserving, nonlinearly implicit, electromagnetic Vlasov–Darwin particle-in-cell algorithm},
author = {Chen, G. and Chacón, L.},
abstractNote = {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.},
doi = {10.1016/j.cpc.2015.08.008},
journal = {Computer Physics Communications},
number = C,
volume = 197,
place = {United States},
year = {Tue Aug 11 00:00:00 EDT 2015},
month = {Tue Aug 11 00:00:00 EDT 2015}
}
Web of Science
Works referenced in this record:
Numerical Cherenkov instabilities in electromagnetic particle codes
journal, August 1974
- Godfrey, Brendan B.
- Journal of Computational Physics, Vol. 15, Issue 4
Some Electromagnetic Plasma Simulation Methods and Their Noise Properties
journal, January 1972
- Langdon, A. Bruce
- Physics of Fluids, Vol. 15, Issue 6
The energy conserving particle-in-cell method
journal, August 2011
- Markidis, Stefano; Lapenta, Giovanni
- Journal of Computational Physics, Vol. 230, Issue 18
An energy- and charge-conserving, nonlinearly implicit, electromagnetic 1D-3V Vlasov–Darwin particle-in-cell algorithm
journal, October 2014
- Chen, G.; Chacón, L.
- Computer Physics Communications, Vol. 185, Issue 10
A self-consistent magnetostatic particle code for numerical simulation of plasmas
journal, March 1977
- Busnardo-Neto, J.; Pritchett, P. L.; Lin, A. T.
- Journal of Computational Physics, Vol. 23, Issue 3
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
Low-frequency electromagnetic (Darwin) applications in plasma simulation
journal, November 1994
- Hewett, D. W.
- Computer Physics Communications, Vol. 84, Issue 1-3
The Darwin Direct Implicit Particle-in-Cell (DADIPIC) Method for Simulation of Low Frequency Plasma Phenomena
journal, September 1995
- Gibbons, Matthew R.; Hewett, Dennis W.
- Journal of Computational Physics, Vol. 120, Issue 2
Nonlinear perturbative electromagnetic (Darwin) particle simulation of high intensity beams
journal, May 2001
- 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
Study of hot electron beam transport in high density plasma using 3D hybrid-Darwin code
journal, December 2004
- Taguchi, T.; Antonsen, T. M.; Mima, K.
- Computer Physics Communications, Vol. 164, Issue 1-3
Simulations of electromagnetic effects in high-frequency capacitively coupled discharges using the Darwin approximation
journal, February 2013
- Eremin, Denis; Hemke, Torben; Brinkmann, Ralf Peter
- Journal of Physics D: Applied Physics, Vol. 46, Issue 8
Darwin–Vlasov simulations of magnetised plasmas
journal, May 2006
- Schmitz, H.; Grauer, R.
- Journal of Computational Physics, Vol. 214, Issue 2
Boundary conditions for the Darwin model
journal, October 1989
- Weitzner, H.; Lawson, W. S.
- Physics of Fluids B: Plasma Physics, Vol. 1, Issue 10
An analysis of the Darwin model of approximation to Maxwell’s equations
journal, January 1992
- Degond, Pierre; Raviart, Pierre-Arnaud
- Forum Mathematicum, Vol. 4, Issue 4
An energy- and charge-conserving, implicit, electrostatic particle-in-cell algorithm
journal, August 2011
- Chen, G.; Chacón, L.; Barnes, D. C.
- Journal of Computational Physics, Vol. 230, Issue 18
Development of a Consistent and Stable Fully Implicit Moment Method for Vlasov--Ampère Particle in Cell (PIC) System
journal, January 2013
- Taitano, William T.; Knoll, Dana A.; Chacón, Luis
- SIAM Journal on Scientific Computing, Vol. 35, Issue 5
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
Fluid preconditioning for Newton–Krylov-based, fully implicit, electrostatic particle-in-cell simulations
journal, February 2014
- Chen, G.; Chacón, L.; Leibs, C. A.
- Journal of Computational Physics, Vol. 258
An efficient mixed-precision, hybrid CPU–GPU implementation of a nonlinearly implicit one-dimensional particle-in-cell algorithm
journal, June 2012
- Chen, G.; Chacón, L.; Barnes, D. C.
- Journal of Computational Physics, Vol. 231, Issue 16
Implicit moment particle simulation of plasmas
journal, June 1981
- Mason, Rodney J.
- Journal of Computational Physics, Vol. 41, Issue 2
Time-filtering particle simulations with ωpe Δt ⪢ 1
journal, August 1981
- Denavit, J.
- Journal of Computational Physics, Vol. 42, Issue 2
An implicit method for electromagnetic plasma simulation in two dimensions
journal, May 1982
- Brackbill, J. U.; Forslund, D. W.
- Journal of Computational Physics, Vol. 46, Issue 2
CELEST1D: an implicit, fully kinetic model for low-frequency, electromagnetic plasma simulation
journal, March 1992
- Vu, H. X.; Brackbill, J. U.
- Computer Physics Communications, Vol. 69, Issue 2-3
Direct implicit large time-step particle simulation of plasmas
journal, July 1983
- Langdon, A. Bruce; Cohen, Bruce I.; Friedman, Alex
- Journal of Computational Physics, Vol. 51, Issue 1
Electromagnetic direct implicit plasma simulation
journal, September 1987
- Hewett, Dennis W.; Bruce Langdon, A.
- Journal of Computational Physics, Vol. 72, Issue 1
Implicit particle simulation of electromagnetic plasma phenomena
journal, May 1992
- Kamimura, T.; Montalvo, E.; Barnes, D. C.
- Journal of Computational Physics, Vol. 100, Issue 1
Elimination of electromagnetic radiation in plasma simulation: The Darwin or magneto inductive approximation
journal, October 1985
- Hewett, Dennis W.
- Space Science Reviews, Vol. 42, Issue 1-2
A hierarchy of approximate models for the Maxwell equations
journal, May 1996
- Raviart, Pierre-Arnaud; Sonnendrücker, Eric
- Numerische Mathematik, Vol. 73, Issue 3
A unified approach to the Darwin approximation
journal, October 2007
- Krause, Todd B.; Apte, A.; Morrison, P. J.
- Physics of Plasmas, Vol. 14, Issue 10
Hyperbolic Conservation Laws in Continuum Physics
book, January 2000
- Dafermos, Constantine M.
- Grundlehren der mathematischen Wissenschaften
The Origin of Spurious Solutions in Computational Electromagnetics
journal, April 1996
- Jiang, Bo-nan; Wu, Jie; Povinelli, L. A.
- Journal of Computational Physics, Vol. 125, Issue 1
One-Dimensional Plasma Model in the Presence of a Magnetic Field
journal, January 1968
- Hasegawa, Akira
- Physics of Fluids, Vol. 11, Issue 9
Divergence Correction Techniques for Maxwell Solvers Based on a Hyperbolic Model
journal, July 2000
- Munz, C. -D.; Omnes, P.; Schneider, R.
- Journal of Computational Physics, Vol. 161, Issue 2
Accurate Numerical Solution of Charged Particle Motion in a Magnetic Field
journal, February 1995
- Vu, H. X.; Brackbill, J. U.
- Journal of Computational Physics, Vol. 116, Issue 2
An analytical particle mover for the charge- and energy-conserving, nonlinearly implicit, electrostatic particle-in-cell algorithm
journal, August 2013
- Chen, G.; Chacón, L.
- Journal of Computational Physics, Vol. 247
Rigorous charge conservation for local electromagnetic field solvers
journal, March 1992
- Villasenor, John; Buneman, Oscar
- Computer Physics Communications, Vol. 69, Issue 2-3
The Darwin Model as a Tool for Electromagnetic Plasma Simulation
journal, January 1971
- Kaufman, Allan N.
- Physics of Fluids, Vol. 14, Issue 2
Nested Iteration and First-Order Systems Least Squares for a Two-Fluid Electromagnetic Darwin Model
journal, January 2015
- Leibs, C. A.; Manteuffel, T. A.
- SIAM Journal on Scientific Computing, Vol. 37, Issue 5
Kinetic Alfvén waves and electron physics. I. Generation from ion-ion streaming
journal, June 2007
- Yin, L.; Winske, D.; Daughton, W.
- Physics of Plasmas, Vol. 14, Issue 6
An implicit “drift-Lorentz” mover for plasma and beam simulations
journal, July 2009
- 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
A charge- and energy-conserving implicit, electrostatic particle-in-cell algorithm on mapped computational meshes
journal, January 2013
- Chacón, L.; Chen, G.; Barnes, D. C.
- Journal of Computational Physics, Vol. 233
Numerical Simulation of the Weibel Instability in One and Two Dimensions
journal, January 1971
- Morse, R. L.
- Physics of Fluids, Vol. 14, Issue 4
The energy conserving particle-in-cell method
journal, August 2011
- Markidis, Stefano; Lapenta, Giovanni
- Journal of Computational Physics, Vol. 230, Issue 18
A finite element formulation of the Darwin PIC model for use on unstructured grids
journal, October 1995
- Sonnendrücker, Eric; Ambrosiano, John J.; Brandon, Scott T.
- Journal of Computational Physics, Vol. 121, Issue 2
The Darwin Model as a Tool for Electromagnetic Plasma Simulation
journal, January 1971
- Kaufman, Allan N.
- Physics of Fluids, Vol. 14, Issue 2
Assessment and Intervention Strategies for the Bilingual Language-Disordered Student
journal, September 1983
- Langdon, Henriette W.
- Exceptional Children, Vol. 50, Issue 1
Works referencing / citing this record:
Vlasov-Maxwell system with physical reasonableness constraints
journal, May 2019
- Lin, H.; Liu, C. P.
- Plasma Research Express, Vol. 1, Issue 2
ViDA: a Vlasov–DArwin solver for plasma physics at electron scales
journal, October 2019
- Pezzi, Oreste; Cozzani, Giulia; Califano, Francesco
- Journal of Plasma Physics, Vol. 85, Issue 5
Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code
journal, May 2017
- Hager, Robert; Lang, Jianying; Chang, C. S.
- Physics of Plasmas, Vol. 24, Issue 5
Interpolation-free particle simulation
journal, January 2020
- Lin, H.; Liu, C. P.
- Laser and Particle Beams, Vol. 38, Issue 1
A Post-Newtonian Expansion Including Radiation Damping for a Collisionless Plasma
journal, September 2019
- Bauer, Sebastian
- Journal of Nonlinear Science, Vol. 30, Issue 1
Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration
journal, May 2018
- Ghosh, Debojyoti; Dorf, Mikhail A.; Dorr, Milo R.
- Journal of Scientific Computing, Vol. 77, Issue 2
A Semi-implicit Particle-in-cell Expanding Box Model Code for Fully Kinetic Simulations of the Expanding Solar Wind Plasma
journal, January 2019
- Elena Innocenti, Maria; Tenerani, Anna; Velli, Marco
- The Astrophysical Journal, Vol. 870, Issue 2
A Post-Newtonian Expansion Including Radiation Damping for a Collisionless Plasma
text, January 2019
- Bauer, Sebastian
- Unpublished
Kinetic Simulation of Collisional Magnetized Plasmas with Semi-Implicit Time Integration
preprint, January 2017
- Ghosh, Debojyoti; Dorf, Mikhail A.; Dorr, Milo R.
- arXiv
Electromagnetic full-$f$ gyrokinetics in the tokamak edge with discontinuous Galerkin methods
text, January 2019
- Mandell, N. R.; Hakim, A.; Hammett, G. W.
- arXiv