Verification of a hybrid gyrokinetic model using the advanced semi-Lagrange code ssV
arXiv:2505.10195 · doi:10.1016/j.cpc.2025.109980
Abstract
The super simple Vlasov (ssV) code was developed to study instabilities, turbulence, and reconnection in weakly magnetized plasmas, such as the solar wind in the dissipation range and the edge of fusion plasmas. The ssV code overcomes the limitations of standard gyrokinetic theory by using a hybrid model that incorporates fully kinetic ions and gyrokinetic electrons. This hybrid gyrokinetic model enables accurate modeling in regimes characterized by steep gradients and high-frequency dynamics. To achieve this, ssV implements a set of semi-Lagrangian numerical schemes, including Positive Flux Conservative (PFC), Flux Conservative fifth-order (FCV), FCV with Umeda limiters, and a Semi-Lagrangian Monotonicity-Preserving fifth-order scheme (SLMP5). Benchmark problems such as Landau damping, ion-acoustic waves, ion Bernstein waves, and kinetic Alfven waves were employed to evaluate the schemes. The SLMP5 scheme consistently delivered the best overall accuracy and numerical stability performance. The code also addresses well-known electromagnetic gyrokinetic simulation issues, such as the Ampere cancellation problem, using carefully chosen velocity-space resolutions and accurate integral evaluation.
References in corpus (16)
- Astrophysical gyrokinetics: kinetic and fluid turbulent cascades in magnetized weakly collisional plasmas
- A Model of Turbulence in Magnetized Plasmas: Implications for the Dissipation Range in the Solar Wind
- Gyrokinetic Simulations of Solar Wind Turbulence from Ion to Electron Scales
- Recent progress in astrophysical plasma turbulence from solar wind observations
- Current Sheets and Collisionless Damping in Kinetic Plasma Turbulence
- Multiscale nature of the dissipation range in gyrokinetic simulations of Alfvénic turbulence
- Pegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics
- Fully kinetic versus reduced-kinetic modelling of collisionless plasma turbulence
- Electromagnetic full- gyrokinetics in the tokamak edge with discontinuous Galerkin methods
- Darwin-Vlasov Simulations of magnetized Plasmas
- Multidimensional VlasovPoisson Simulations with High-order Monotonicity- and Positivity-preserving Schemes
- ViDA: a Vlasov-DArwin solver for plasma physics at electron scales
- Coupled Vlasov and two-fluid codes on GPUs
- A new hybrid code (CHIEF) implementing the inertial electron fluid equation without approximation
- A performance portable implementation of the semi-Lagrangian algorithm in six dimensions
- Discovering exact, gauge-invariant, local energy-momentum conservation laws for the electromagnetic gyrokinetic system by high-order field theory on heterogeneous manifolds