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

doi:10.1016/j.jcp.2016.03.070

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

Abstract

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.

Authors:

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

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

Publication Date:: 2016-04-19

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

Sponsoring Org.:: USDOE

OSTI Identifier:: 1255248

Alternate Identifier(s):: OSTI ID: 1348010

Report Number(s):: LA-UR-15-27639
Journal ID: ISSN 0021-9991

Grant/Contract Number:: AC52-06NA25396

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 316; Journal Issue: C; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

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

Citation Formats


                    Chacon, L., and Chen, G. A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions.  United States: N. p., 2016. 
        Web.  doi:10.1016/j.jcp.2016.03.070.

Copy to clipboard


                    Chacon, L., & Chen, G. A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions.  United States.  doi:10.1016/j.jcp.2016.03.070.

Copy to clipboard


                    Chacon, L., and Chen, G. Tue .  
        "A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions".  United States.  doi:10.1016/j.jcp.2016.03.070.  https://www.osti.gov/servlets/purl/1255248.

Copy to clipboard


                    
@article{osti_1255248,

  title        = {A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions},

  author       = {Chacon, L. and Chen, G.},

  abstractNote = {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.},

  doi          = {10.1016/j.jcp.2016.03.070},

  journal      = {Journal of Computational Physics},

  number       = C,

  volume       = 316,

  place        = {United States},

  year         = {2016},

  month        = {4}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

DOI: 10.1016/j.jcp.2016.03.070

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 5 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Works referencing / citing this record:

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

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

Similar Records in DOE PAGES and OSTI.GOV collections:

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

Journal Article Chen, G. ; Chacón, L. - Computer Physics Communications

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,more »« less
Cited by 8
DOI: 10.1016/j.cpc.2015.08.008

Full Text Available
Fully implicit Particle-in-cell algorithms for multiscale plasma simulation

Conference Chacon, Luis

The outline of the paper is as follows: Particle-in-cell (PIC) methods for fully ionized collisionless plasmas, explicit vs. implicit PIC, 1D ES implicit PIC (charge and energy conservation, moment-based acceleration), and generalization to Multi-D EM PIC: Vlasov-Darwin model (review and motivation for Darwin model, conservation properties (energy, charge, and canonical momenta), and numerical benchmarks). The author demonstrates a fully implicit, fully nonlinear, multidimensional PIC formulation that features exact local charge conservation (via a novel particle mover strategy), exact global energy conservation (no particle self-heating or self-cooling), adaptive particle orbit integrator to control errors in momentum conservation, and canonical momenta (EM-PICmore »« less
Full Text Available
Exact charge and energy conservation in implicit PIC with mapped computational meshes

Journal Article Chen, Guangye ; Barnes, D. C. - Journal of Computational Physics

This paper discusses a novel fully implicit formulation for a one-dimensional electrostatic particle-in-cell (PIC) plasma simulation approach. Unlike earlier implicit electrostatic PIC approaches (which are based on a linearized Vlasov Poisson formulation), ours is based on a nonlinearly converged Vlasov Amp re (VA) model. By iterating particles and fields to a tight nonlinear convergence tolerance, the approach features superior stability and accuracy properties, avoiding most of the accuracy pitfalls in earlier implicit PIC implementations. In particular, the formulation is stable against temporal (Courant Friedrichs Lewy) and spatial (aliasing) instabilities. It is charge- and energy-conserving to numerical round-off for arbitrary implicitmore »« less
Energy-conserving perfect-conductor boundary conditions for an implicit, curvilinear Darwin particle-in-cell algorithm

Journal Article Chacón, L. ; Chen, G. - Journal of Computational Physics

We proposemore »« less
DOI: 10.1016/j.jcp.2019.04.032

Full Text Available
Energy-Charge-conserving, Implicit, Electrostatic Particle-in-Cell Algorithm

Journal Article Chen, Guangye ; Chacon, Luis ; Barnes, Daniel C - Journal of Computational Physics

This paper discusses a novel fully implicit formulation for a one-dimensional electrostatic particle-in-cell (PIC) plasma simulation approach. Unlike earlier implicit electrostatic PIC approaches (which are based on a linearized Vlasov-Poisson formulation), ours is based on a nonlinearly converged Vlasov-Ampere (VA) model. By iterating particles and fields to a tight nonlinear convergence tolerance, the approach features superior stability and accuracy properties, avoiding most of the accuracy pitfalls in earlier implicit PIC implementations. In particular, the formulation is stable against temporal (Courant-Friedrichs-Lewy) and spatial (aliasing) instabilities. It is charge- and energy-conserving to numerical round-off for arbitrary implicit time steps (unlike the earliermore »« less
DOI: 10.1016/j.jcp.2011.05.031

Similar Records

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

Abstract

Citation Formats

Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code journal, May 2017

Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code journal, May 2017

Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code
journal, May 2017

Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code
journal, May 2017