2023/04/04 by Yingzhe Li, Li, Yingzhe, Martin Campos Pinto +7 · 1 citation
Mathematics · Physics and Astronomy · #65M75 #FOS: Mathematics #FOS: Physical sciences #G.1.8 #Gas Dynamics and Kinetic Theory #Ionosphere and magnetosphere dynamics #Magnetic confinement fusion research #Numerical Analysis (math.NA) #Plasma Physics (physics.plasm-ph)
paper · pdf · doi:10.48550/arxiv.2304.01891
openalex publication_date 2023/04/04 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We explore the possibilities of applying structure-preserving numerical methods to a plasma hybrid model with kinetic ions and mass-less fluid electrons satisfying the quasi-neutrality relation. The numerical schemes are derived by finite element methods in the framework of finite element exterior calculus (FEEC) for field variables, particle-in-cell (PIC) methods for the Vlasov equation, and splitting methods in time based on an anti-symmetric bracket proposed. Conservation properties of energy, quasi-neutrality relation, positivity of density, and divergence-free property of the magnetic field are given irrespective of the used resolution and metric. Local quasi-interpolation is used for dealing with the current terms in order to make the proposed methods more efficient. The implementation has been done in the framework of the Python package STRUPHY [1], and has been verified by extensive numerical experiments.