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

Chacon, L.; Chen, G.

doi:10.1016/j.jcp.2016.03.070

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

DOI:https://doi.org/10.1016/j.jcp.2016.03.070· OSTI ID:1255248

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

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

Here, we extend a recently proposed fully implicit PIC algorithm for the Vlasov–Darwin model in multiple dimensions (Chen and Chacón (2015) [1]) to curvilinear geometry. As in the Cartesian case, the approach is based on a potential formulation (Φ, A), and overcomes many difficulties of traditional semi-implicit Darwin PIC algorithms. Conservation theorems for local charge and global energy are derived in curvilinear representation, and then enforced discretely by a careful choice of the discretization of field and particle equations. Additionally, the algorithm conserves canonical-momentum in any ignorable direction, and preserves the Coulomb gauge ∇ • A = 0 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. We demonstrate the accuracy and efficiency properties of the algorithm with numerical experiments in mapped meshes in 1D-3V and 2D-3V.

View Accepted Manuscript (DOE)

View Accepted Manuscript (Publisher)

Cite

Export

Save

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

Sponsoring Organization:: USDOE

Grant/Contract Number:: AC52-06NA25396

OSTI ID:: 1255248

Alternate ID(s):: OSTI ID: 1348010

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

Journal Information:: Journal of Computational Physics, Vol. 316, Issue C; ISSN 0021-9991

Publisher:: ElsevierCopyright Statement

Country of Publication:: United States

Language:: English

Citation Metrics:

Cited by: 27 works

Citation information provided by
Web of Science

References (45)

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 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
Bounded Multi-scale Plasma Simulation: Application to Sheath Problems Parker, S. E.; Friedman, A.; Ray, S. L. Journal of Computational Physics, Vol. 107, Issue 2 https://doi.org/10.1006/jcph.1993.1153	journal	August 1993
OSIRIS: A Three-Dimensional, Fully Relativistic Particle in Cell Code for Modeling Plasma Based Accelerators Fonseca, R. A.; Silva, L. O.; Tsung, F. S. Computational Science — ICCS 2002, p. 342-351 https://doi.org/10.1007/3-540-47789-6_36	book	January 2002
Three-dimensional deformable-grid electromagnetic particle-in-cell for parallel computers Wang, J.; Kondrashov, D.; Liewer, P. C. Journal of Plasma Physics, Vol. 61, Issue 3 https://doi.org/10.1017/S0022377899007552	journal	April 1999
A finite-volume particle-in-cell method for the numerical treatment of Maxwell-Lorentz equations on boundary-fitted meshes Munz, C. -D.; Schneider, R.; Sonnendr�cker, E. International Journal for Numerical Methods in Engineering, Vol. 44, Issue 4 https://doi.org/10.1002/(SICI)1097-0207(19990210)44:4<461::AID-NME510>3.0.CO;2-#	journal	February 1999
Body-fitted electromagnetic PIC software for use on parallel computers Eastwood, J. W.; Arter, W.; Brealey, N. J. Computer Physics Communications, Vol. 87, Issue 1-2 https://doi.org/10.1016/0010-4655(94)00165-X	journal	May 1995
High-order nodal discontinuous Galerkin particle-in-cell method on unstructured grids Jacobs, G. B.; Hesthaven, J. S. Journal of Computational Physics, Vol. 214, Issue 1 https://doi.org/10.1016/j.jcp.2005.09.008	journal	May 2006
Three-Dimensional Numerical Simulation of a 30-GHz Gyrotron Resonator With an Explicit High-Order Discontinuous-Galerkin-Based Parallel Particle-In-Cell Method Stock, Andreas; Neudorfer, Jonathan; Riedlinger, Marc IEEE Transactions on Plasma Science, Vol. 40, Issue 7 https://doi.org/10.1109/TPS.2012.2195509	journal	July 2012
Electron–Ion Coupling in Mesothermal Plasma Beam Emission: Full Particle PIC Simulations Wang, Joseph; Chang, Ouliang; Cao, Yong IEEE Transactions on Plasma Science, Vol. 40, Issue 2 https://doi.org/10.1109/TPS.2011.2179066	journal	February 2012
3D Hybrid Simulation Code Using Curvilinear Coordinates Bagdonat, T.; Motschmann, U. Journal of Computational Physics, Vol. 183, Issue 2 https://doi.org/10.1006/jcph.2002.7203	journal	December 2002
Electromagnetic particle-in-cell simulations on magnetic reconnection with adaptive mesh refinement Fujimoto, Keizo; Sydora, Richard D. Computer Physics Communications, Vol. 178, Issue 12 https://doi.org/10.1016/j.cpc.2008.02.010	journal	June 2008
CPIC: A Curvilinear Particle-in-Cell Code for Plasma–Material Interaction Studies Delzanno, Gian Luca; Camporeale, Enrico; Moulton, J. David IEEE Transactions on Plasma Science, Vol. 41, Issue 12 https://doi.org/10.1109/TPS.2013.2290060	journal	December 2013
Turbulent Transport Reduction by Zonal Flows: Massively Parallel Simulations Lin, Z. Science, Vol. 281, Issue 5384 https://doi.org/10.1126/science.281.5384.1835	journal	September 1998
Application of adaptive mesh refinement to PIC simulations in heavy ion fusion Vay, J. -L.; Friedman, A.; Grote, D. P. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Vol. 544, Issue 1-2 https://doi.org/10.1016/j.nima.2005.01.232	journal	May 2005
Density-Conserving Shape Factors for Particle Simulations in Cylindrical and Spherical Coordinates Ruyten, Wilhelmus M. Journal of Computational Physics, Vol. 105, Issue 2 https://doi.org/10.1006/jcph.1993.1070	journal	April 1993
Symmetric Spline Weighting for Charge and Current Density in Particle Simulation Verboncoeur, J. P. Journal of Computational Physics, Vol. 174, Issue 1 https://doi.org/10.1006/jcph.2001.6923	journal	November 2001
Correction factors for PIC accumulation on radial grids Larson, David J.; Hewett, Dennis W.; Langdon, A. Bruce Computer Physics Communications, Vol. 90, Issue 2-3 https://doi.org/10.1016/0010-4655(95)00104-N	journal	October 1995
A unified Monte Carlo interpretation of particle simulations and applications to non‐neutral plasmas Aydemir, A. Y. Physics of Plasmas, Vol. 1, Issue 4 https://doi.org/10.1063/1.870740	journal	April 1994
Electromagnetic gyrokinetic δf particle-in-cell turbulence simulation with realistic equilibrium profiles and geometry Chen, Yang; Parker, Scott E. Journal of Computational Physics, Vol. 220, Issue 2 https://doi.org/10.1016/j.jcp.2006.05.028	journal	January 2007
Localization schemes in 2D boundary-fitted grids Westermann, Thomas Journal of Computational Physics, Vol. 101, Issue 2 https://doi.org/10.1016/0021-9991(92)90008-M	journal	August 1992
Particle-in-cell simulations with moving boundaries—adaptive mesh generation Westermann, T. Journal of Computational Physics, Vol. 114, Issue 2 https://doi.org/10.1006/jcph.1994.1157	journal	October 1994
Algorithms for interpolation and localization in irregular 2D meshes Seldner, David; Westermann, Thomas Journal of Computational Physics, Vol. 79, Issue 1 https://doi.org/10.1016/0021-9991(88)90001-0	journal	November 1988
An arbitrary curvilinear-coordinate method for particle-in-cell modeling Fichtl, C. A.; Finn, J. M.; Cartwright, K. L. Computational Science & Discovery, Vol. 5, Issue 1 https://doi.org/10.1088/1749-4699/5/1/014011	journal	January 2012
Use of a Hybrid Code for Global-Scale Plasma Simulation Swift, Daniel W. Journal of Computational Physics, Vol. 126, Issue 1 https://doi.org/10.1006/jcph.1996.0124	journal	June 1996
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
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
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
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
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
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
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
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
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
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
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
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
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
Role of electron physics in slow mode shocks Daughton, William; Winske, Dan; Yin, Lin Journal of Geophysical Research: Space Physics, Vol. 106, Issue A11 https://doi.org/10.1029/2001JA000005	journal	November 2001
Dissipation in oblique slow shocks: DISSIPATION IN OBLIQUE SLOW SHOCKS Yin, L.; Winske, D.; Daughton, W. Journal of Geophysical Research: Space Physics, Vol. 110, Issue A9 https://doi.org/10.1029/2005JA011031	journal	September 2005
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
Laminar wave-train structure of collisionless magnetic slow shocks Coroniti, F. V. Nuclear Fusion, Vol. 11, Issue 3 https://doi.org/10.1088/0029-5515/11/3/008	journal	May 1971
The structure and evolution of slow mode shocks Winske, D.; Stover, Elmer K.; Gary, S. Peter Geophysical Research Letters, Vol. 12, Issue 5 https://doi.org/10.1029/GL012i005p00295	journal	May 1985
Electron kinetic effects in switch-off slow shocks Brackbill, J. U.; Vu, H. X. Geophysical Research Letters, Vol. 20, Issue 19 https://doi.org/10.1029/93GL02243	journal	October 1993

Cited By (2)

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
Metriplectic integrators for the Landau collision operator Kraus, Michael; Hirvijoki, Eero Physics of Plasmas, Vol. 24, Issue 10 https://doi.org/10.1063/1.4998610	journal	October 2017

Similar Records

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

Chen, G.; Chacón, L.

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

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

Chacon, Luis

An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system

Journal Article · Mon Jun 29 00:00:00 EDT 2020 · Journal of Computational Physics · OSTI ID:1255248

Anderson, S. E.; Taitano, W. T.; Chacón, L.; +1 more

Related Subjects

97 MATHEMATICS AND COMPUTING
70 PLASMA PHYSICS AND FUSION TECHNOLOGY
Mathematics
Magnetic Fusion Energy

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

Citation Formats

References (45)

Cited By (2)

Similar Records

Related Subjects