AHKASH: a new Hybrid particle-in-cell code for simulations of astrophysical collisionless plasma

Abstract

We introduce $\text{A}$strophysical $\text{H}$ybrid-$\text{K}$inetic simulations with the $\text{flASH}$ code ($\text{AHKASH}$) -- a new Hybrid particle-in-cell (PIC) code developed within the framework of the multi-physics code $\text{FLASH}$. The new code uses a second-order accurate Boris integrator and a predictor-predictor-corrector algorithm for advancing the Hybrid-kinetic equations, using the constraint transport method to ensure that magnetic fields are divergence-free. The code supports various interpolation schemes between the particles and grid cells, with post-interpolation smoothing to reduce finite particle noise. We further implement a $\delta f$ method to study instabilities in weakly collisional plasmas. The new code is tested on standard physical problems such as the motion of charged particles in uniform and spatially varying magnetic fields, the propagation of Alfvén and whistler waves, and Landau damping of ion acoustic waves. We test different interpolation kernels and demonstrate the necessity of performing post-interpolation smoothing. We couple the $\text{TurbGen}$ turbulence driving module to the new Hybrid PIC code, allowing us to test the code on the highly complex physical problem of the turbulent dynamo. To investigate steady-state turbulence with a fixed sonic Mach number, it is important to maintain isothermal plasma conditions. Therefore, we introduce a novel cooling method for Hybrid PIC codes and provide tests and calibrations of this method to keep the plasma isothermal. We describe and test the `hybrid precision' method, which significantly reduces (by a factor $\sim 1.5$) the computational cost, without compromising the accuracy of the numerical solutions. Finally, we test the parallel scalability of the new code, showing excellent scaling up to 10,000~cores.

23 pages, 20 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

References in corpus (28)

• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• PLUTO: a Numerical Code for Computational Astrophysics
• Athena: A New Code for Astrophysical MHD
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• The Quiescent Intracluster Medium in the Core of the Perseus Cluster
• The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained
• Simulations of relativistic collisionless shocks: shock structure and particle acceleration
• Magnetic field amplification in turbulent astrophysical plasmas
• The sonic scale revealed by the world's largest supersonic turbulence simulation
• The density structure and star formation rate of non-isothermal polytropic turbulence
• Mach number dependence of the onset of dynamo action
• Predicted Impacts of Proton Temperature Anisotropy on Solar Wind Turbulence
• Hybrid-Kinetic Simulations of Ion Heating in Alfvénic Turbulence
• Constraining Gas Motions in the Intra-Cluster Medium
• Validity of the Taylor Hypothesis for Linear Kinetic Waves in the Weakly Collisional Solar Wind
• Plasma turbulence in the interstellar medium
• Turbulent dynamo in the two-phase interstellar medium
• Seed magnetic fields in turbulent small-scale dynamos
• Measuring sloshing, merging and feedback velocities in the Virgo cluster
• Hybrid-VPIC: an Open-Source Kinetic/Fluid Hybrid Particle-in-Cell Code
• Efficient highly-subsonic turbulent dynamo and growth of primordial magnetic fields
• Fundamental scales in the kinematic phase of the turbulent dynamo
• The small-scale dynamo in a multiphase supernova-driven medium
• The velocity structure of the Intracluster Medium of the Centaurus cluster
• Measuring the ICM velocity structure in the Ophiuchus cluster
• Electron-ion heating partition in imbalanced solar-wind turbulence
• Critical magnetic Reynolds number of the turbulent dynamo in collisionless plasmas
• Hybrid codes (massless electron fluid)