2018/03/19 by Adam Stanier, A. Stanier, Luis Chacòn +3 · 29 citations
Physics and Astronomy · #Algorithm #Computer science #Electron #Instability #Ion #Ionosphere and magnetosphere dynamics #Magnetic confinement fusion research #Mechanics #Particle-in-cell #Physics #Population #Quantum mechanics #Smoothing #Solar and Space Plasma Dynamics #Statistical physics #physics.comp-ph #physics.plasm-ph #physics.space-ph
paper · pdf · doi:10.1016/j.jcp.2018.09.038
published in Journal of Computational Physics 376, 597-616 (Elsevier BV)
arxiv created 2018/03/19 · openalex publication_date 2018/09/21 · arxiv updated 2018/11/14 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/05
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. Here, a novel, implicit, particle-in-cell based scheme for the hybrid model is derived for multi-dimensional electromagnetic problems with multiple ion species, which features global mass, momentum and energy conservation. The scheme includes sub-cycling and orbit averaging of 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.