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

Stanier, A.; Chacón, L.; Chen, G.

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.  https://doi.org/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.  https://doi.org/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

https://doi.org/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 conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes

Journal Article Stanier, Adam John ; Chacon, Luis - Journal of Computational Physics

We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven tomore »« less
https://doi.org/10.1016/j.jcp.2022.111144
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
https://doi.org/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
https://doi.org/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
https://doi.org/10.1016/j.cpc.2020.107560

Full Text Available
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
https://doi.org/10.1016/j.jcp.2012.04.040

Similar Records