Viriato: a Fourier-Hermite spectral code for strongly magnetised fluid-kinetic plasma dynamics

Abstract

We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihin et al., Astrophys. J. Suppl. 182:310 (2009)]. Two main applications of these equations are magnetised (Alfvenic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge-Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor-corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag-Tang-type decaying turbulence, both in fluid and kinetic regimes.

42 pages, 15 figures, submitted to J. Comp. Phys

References in corpus (8)

• Instability of current sheets and formation of plasmoid chains
• Kinetic Simulations of Magnetized Turbulence in Astrophysical Plasmas
• Numerical simulations of possible finite time singularities in the incompressible Euler equations: comparison of numerical methods
• Fourier-Hermite spectral representation for the Vlasov-Poisson system in the weakly collisional limit
• Kinetic microtearing modes and reconnecting modes in strongly magnetised slab plasmas
• An iterative semi-implicit scheme with robust damping
• Gyrokinetic and kinetic particle-in-cell simulations of guide-field reconnection. I: Macroscopic effects of the electron flows
• Linear Collisionless Landau Damping in Hilbert Space

Cited by in corpus (31)

• Improved quantum algorithms for linear and nonlinear differential equations
• Fully kinetic versus reduced-kinetic modelling of collisionless plasma turbulence
• A Prospectus on Kinetic Heliophysics
• Laguerre-Hermite Pseudo-Spectral Velocity Formulation of Gyrokinetics
• Constraints on ion vs. electron heating by plasma turbulence at low beta
• Multi-scale dynamics of magnetic flux tubes and inverse magnetic energy transfer
• A gyrokinetic model for the plasma periphery of tokamak devices
• A Drift-Kinetic Analytical Model for SOL Plasma Dynamics at Arbitrary Collisionality
• Magnetic island merger as a mechanism for inverse magnetic energy transfer
• Electromagnetic instabilities and plasma turbulence driven by electron-temperature gradient
• The multi-dimensional Hermite-discontinuous Galerkin method for the Vlasov-Maxwell equations
• Theory of the Drift-Wave Instability at Arbitrary Collisionality
• Dynamic Phase Alignment in Inertial Alfven Turbulence
• Gyrokinetic simulations of plasma turbulence in a Z-pinch using a moment based approach and advanced collision operators
• Linear Theory of Electron-Plasma Waves at Arbitrary Collisionality
• Characterizing magnetic reconnection regions using Gaussian mixture models on particle velocity distributions
• Plasmoid instability in the semi-collisional regime
• A Hybrid Gyrokinetic Ion and Isothermal Electron Fluid Code for Astrophysical Plasma
• Physics-based adaptivity of a spectral method for the Vlasov-Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space
• Phase space dynamics of unmagnetized plasmas: collisionless and collisional regimes
• Calliope: Pseudospectral shearing magnetohydrodynamics code with a pencil decomposition
• Four-dimensional equations for the study of electromagnetic plasma turbulence in a drift kinetic limit
• The effects of finite electron inertia on helicity-barrier-mediated turbulence
• Machine-learning Closure for Vlasov-Poisson Dynamics in Fourier-Hermite Space
• A Parallel-Kinetic-Perpendicular-Moment Model for Magnetized Plasmas
• Turbulent heating in collisionless low-beta plasmas: imbalance, Landau damping, and electron-ion energy partition
• Four-dimensional drift-kinetic model for scrape-off layer plasmas
• Symmetries of a reduced fluid-gyrokinetic system
• The gyrokinetic field invariant and electromagnetic temperature-gradient instabilities in `good-curvature' plasmas
• Transition from weak to strong turbulence in magnetized plasmas
• The Machine Learning Approach to Moment Closure Relations for Plasma: A Review