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

doi:10.1016/j.jcp.2018.09.038

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

Abstract

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.

Authors:

^[1];

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

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

Publication Date:: 2018-09-20

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

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

OSTI Identifier:: 1477697

Alternate Identifier(s):: OSTI ID: 1756197

Report Number(s):: LA-UR-18-21708
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: 376; 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; hybrid; plasma; particle-in-cell; implicit; conservative; finite grid instability

Citation Formats


                    Stanier, A., Chacón, L., and Chen, G. A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm.  United States: N. p., 2018. 
        Web.  doi:10.1016/j.jcp.2018.09.038.

Copy to clipboard


                    Stanier, A., Chacón, L., & Chen, G. A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm.  United States.  doi:10.1016/j.jcp.2018.09.038.

Copy to clipboard


                    Stanier, A., Chacón, L., and Chen, G. Thu .  
        "A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm".  United States.  doi:10.1016/j.jcp.2018.09.038.  https://www.osti.gov/servlets/purl/1477697.

Copy to clipboard


                    
@article{osti_1477697,

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

  author       = {Stanier, A. and Chacón, L. and Chen, G.},

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

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

  journal      = {Journal of Computational Physics},

  number       = ,

  volume       = 376,

  place        = {United States},

  year         = {2018},

  month        = {9}

}

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.2018.09.038

Other availability

Search WorldCat to find libraries that may hold this journal

Figures / Tables:

Figure 1: Comparison study between an explicit leapfrog-based momentum-conserving scheme (red, see Appendix D.1), a predictor-corrector method (green, see Appendix D.2), and the new implicit conservative algorithm (blue). Top panel: Particle phase space distribution (x_p, v_xp) at t = 8.8Δx/C_s (400 steps). Second panel: Change in total momentum vs timemore »

All figures and tables (5 total)

Save / Share:

Export Metadata

Save to My Library

Figures / Tables found in this record:

Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.

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
A curvilinear, fully implicit, conservative electromagnetic PIC algorithm in multiple dimensions

Journal Article Chacon, L. ; Chen, G. - Journal of Computational Physics

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. Anmore »« less
Cited by 5
DOI: 10.1016/j.jcp.2016.03.070

Full Text Available
Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Journal Article Barnes, D. C. ; Chacón, L. - Computer Physics Communications

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not freemore »« less
DOI: 10.1016/j.cpc.2020.107560
An efficient mixed-precision, hybrid CPU-GPU implementation of a nonlinearly implicit one-dimensional particle-in-cell algorithm

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

Recently, a fully implicit, energy- and charge-conserving particle-in-cell method has been developed for multi-scale, full-f kinetic simulations [G. Chen, et al., J. Comput. Phys. 230, 18 (2011)]. The method employs a Jacobian-free Newton-Krylov (JFNK) solver and is capable of using very large timesteps without loss of numerical stability or accuracy. A fundamental feature of the method is the segregation of particle orbit integrations from the field solver, while remaining fully self-consistent. This provides great flexibility, and dramatically improves the solver efficiency by reducing the degrees of freedom of the associated nonlinear system. However, it requires a particle push per nonlinearmore »« less
DOI: 10.1016/j.jcp.2012.04.040
Fluid preconditioning for Newton–Krylov-based, fully implicit, electrostatic particle-in-cell simulations

Journal Article Chen, G., E-mail: gchen@lanl.gov ; Chacón, L. ; Leibs, C. A. ; ... - Journal of Computational Physics

A recent proof-of-principle study proposes an energy- and charge-conserving, nonlinearly implicit electrostatic particle-in-cell (PIC) algorithm in one dimension [9]. The algorithm in the reference employs an unpreconditioned Jacobian-free Newton–Krylov method, which ensures nonlinear convergence at every timestep (resolving the dynamical timescale of interest). Kinetic enslavement, which is one key component of the algorithm, not only enables fully implicit PIC as a practical approach, but also allows preconditioning the kinetic solver with a fluid approximation. This study proposes such a preconditioner, in which the linearized moment equations are closed with moments computed from particles. Effective acceleration of the linear GMRES solvemore »« less
DOI: 10.1016/J.JCP.2013.10.052

Similar Records