MHD-kinetic hybrid code based on structure-preserving finite elements with particles-in-cell

JOURNAL OF COMPUTATIONAL PHYSICS

卷:433

Article

Holderied, Florian^1,2  Possanner, Stefan^1,3  Wang, Xin¹ 

[1] Max Planck Inst Plasma Phys, Boltzmannstr 2, D-85748 Garching, Germany

[2] Tech Univ Munich, Dept Phys, James Franck Str 1, D-85748 Garching, Germany

[3] Tech Univ Munich, Dept Math, Boltzmannstr 3, D-85748 Garching, Germany

关键词: Plasma simulation;    Finite element exterior calculus;    Particle-in-cell;    Structure-preserving;   

DOI  :  10.1016/j.jcp.2021.110143

来源: Elsevier

【 摘 要 】

We present a STRUcture-Preserving HYbrid code - STRUPHY - for the simulation of magneto-hydrodynamic (MHD) waves interacting with a population of energetic particles far from thermal equilibrium (kinetic species). The implemented model features linear, ideal MHD equations in curved, three-dimensional space, coupled nonlinearly to the fullorbit Vlasov equations via a current coupling scheme. The algorithm is based on finite element exterior calculus (FEEC) for MHD and particle-in-cell (PIC) methods for the kinetic part; it provably conserves mass, energy, and the divergence-free magnetic field, irrespective of metric (= space curvature), mesh parameters and chosen order of the scheme. These properties enable reliable long-time simulations of energetic particle physics in complex geometries, covering the whole range of MHD waves. In STRUPHY, the finite element spaces are built from tensor products of univariate B-splines on the logical cuboid and can be made high-order by increasing the polynomial degree. Time-stepping is based on a skew-symmetric splitting with implicit sub-steps, mitigating CFL conditions from fast magneto-acoustic waves. High-order time splitting schemes can be used in this regard. (C) 2021 Elsevier Inc. All rights reserved.

【 授权许可】

Free