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

Full Record
Other Related Research

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) (SC-21); USDOE National Nuclear Security Administration (NNSA)

OSTI Identifier:: 1477697

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 (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

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

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

Journal Article Chacón, L. ; Chen, G. - Journal of Computational Physics
Cited by 5
DOI: 10.1016/j.jcp.2016.03.070
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
3D electrostatic gyrokinetic electron and fully kinetic ion simulation of lower-hybrid drift instability of Harris current sheet

Journal Article Wang, Zhenyu ; Lin, Yu ; Wang, Xueyi ; ... - Physics of Plasmas

The eigenmode stability properties of three-dimensional lower-hybrid-drift-instabilities (LHDI) in a Harris current sheet with a small but finite guide magnetic field have been systematically studied by employing the gyrokinetic electron and fully kinetic ion (GeFi) particle-in-cell (PIC) simulation model with a realistic ion-to-electron mass ratio m _i/m _e. In contrast to the fully kinetic PIC simulation scheme, the fast electron cyclotron motion and plasma oscillations are systematically removed in the GeFi model, and hence one can employ the realistic m _i/m _e. The GeFi simulations are benchmarked against and show excellent agreement with both the fully kinetic PIC simulation and the analytical eigenmode theory. Our studies indicate that, for small wavenumbers, ky, along the current direction, the most unstable eigenmodes are peaked at the location wheremore »« less
Cited by 2
DOI: 10.1063/1.4954830

Full Text Available
A scalable, fully implicit algorithm for the reduced two-field low-β extended MHD model

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

Here, we demonstrate a scalable fully implicit algorithm for the two-field low-β extended MHD model. This reduced model describes plasma behavior in the presence of strong guide fields, and is of significant practical impact both in nature and in laboratory plasmas. The model displays strong hyperbolic behavior, as manifested by the presence of fast dispersive waves, which make a fully implicit treatment very challenging. In this study, we employ a Jacobian-free Newton–Krylov nonlinear solver, for which we propose a physics-based preconditioner that renders the linearized set of equations suitable for inversion with multigrid methods. As a result, the algorithm ismore »« less
Cited by 2
DOI: 10.1016/j.jcp.2016.09.007

Full Text Available

Similar Records