A Hybrid Gyrokinetic Ion and Isothermal Electron Fluid Code for Astrophysical Plasma

Abstract

This paper describes a new code for simulating astrophysical plasmas that solves a hybrid model composed of gyrokinetic ions (GKI) and an isothermal electron fluid (ITEF) [A. Schekochihin et al., Astrophys. J. Suppl. \textbf{182}, 310 (2009)]. This model captures ion kinetic effects that are important near the ion gyro-radius scale while electron kinetic effects are ordered out by an electron-ion mass ratio expansion. The code is developed by incorporating the ITEF approximation into {\tt AstroGK}, an Eulerian $\delta f$ gyrokinetics code specialized to a slab geometry [R. Numata et al., J. Compute. Pays. \textbf{229}, 9347 (2010)]. The new code treats the linear terms in the ITEF equations implicitly while the nonlinear terms are treated explicitly. We show linear and nonlinear benchmark tests to prove the validity and applicability of the simulation code. Since the fast electron timescale is eliminated by the mass ratio expansion, the Courant--Friedrichs--Lewy condition is much less restrictive than in full gyrokinetic codes; the present hybrid code runs $\sim 2\sqrt{m_{\mathrm{i}}/m_{\mathrm{e}}}\sim 100$ times faster than {\tt AstroGK}\ with a single ion species and kinetic electrons where $m_{\mathrm{i}}/m_{\mathrm{e}}$ is the ion-electron mass ratio. The improvement of the computational time makes it feasible to execute ion scale gyrokinetic simulations with a high velocity space resolution and to run multiple simulations to determine the dependence of turbulent dynamics on parameters such as electron--ion temperature ratio and plasma beta.

References in corpus (15)

• Kinetic Simulations of Magnetized Turbulence in Astrophysical Plasmas
• Linearized model Fokker-Planck collision operators for gyrokinetic simulations. I. Theory
• Multiscale nature of the dissipation range in gyrokinetic simulations of Alfvénic turbulence
• Gyrokinetic turbulence: a nonlinear route to dissipation through phase space
• Magnetorotational Turbulence and Dynamo in a Collisionless Plasma
• Diagnosing collisionless energy transfer using field-particle correlations: gyrokinetic turbulence
• The Structure of Plasma Heating in Gyrokinetic Alfvénic Turbulence
• A Prospectus on Kinetic Heliophysics
• Kinetic microtearing modes and reconnecting modes in strongly magnetised slab plasmas
• Viriato: a Fourier-Hermite spectral code for strongly magnetised fluid-kinetic plasma dynamics
• Plasma turbulence at ion scales: a comparison between PIC and Eulerian hybrid-kinetic approaches
• Energy Dissipation and Landau Damping in Two- and Three-Dimensional Plasma Turbulence
• Verification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code
• Multiscale Gyrokinetics for Rotating Tokamak Plasmas II: Reduced Models for Electron Dynamics
• Gyrokinetic simulation of entropy cascade in two-dimensional electrostatic turbulence

Cited by in corpus (5)

• Thermal disequilibration of ions and electrons by collisionless plasma turbulence
• Constraints on ion vs. electron heating by plasma turbulence at low beta
• Ion versus electron heating in compressively driven astrophysical gyrokinetic turbulence
• Calliope: Pseudospectral shearing magnetohydrodynamics code with a pencil decomposition
• The Core-shift of Sagittarius A* as a Discriminant between Disk and Jet Emission Models with millimeter-VLBI