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.

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

References in corpus (51)

• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• Astrophysical magnetic fields and nonlinear dynamo theory
• PLUTO: a Numerical Code for Computational Astrophysics
• Athena: A New Code for Astrophysical MHD
• Comparing the statistics of interstellar turbulence in simulations and observations: Solenoidal versus compressive turbulence forcing
• SMILEI: a collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Simulations of small-scale turbulent dynamo
• The Quiescent Intracluster Medium in the Core of the Perseus Cluster
• The multi-scale nature of the solar wind
• On the universality of supersonic turbulence
• Simulations of particle acceleration beyond the classical synchrotron burnoff limit in magnetic reconnection: An explanation of the Crab flares
• Mach Number Dependence of Turbulent Magnetic Field Amplification: Solenoidal versus Compressive Flows
• Firehose and Mirror Instabilities in a Collisionless Shearing Plasma
• The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained
• Turbulence, magnetic fields and plasma physics in clusters of galaxies
• Simulations of relativistic collisionless shocks: shock structure and particle acceleration
• A prescription for the turbulent heating of astrophysical plasmas
• Magnetic field amplification in turbulent astrophysical plasmas
• The sonic scale revealed by the world's largest supersonic turbulence simulation
• A thermally stable heating mechanism for the intracluster medium: turbulence, magnetic fields and plasma instabilities
• The density structure and star formation rate of non-isothermal polytropic turbulence
• Magnetic Field Amplification by Small-Scale Dynamo Action: Dependence on Turbulence Models and Reynolds and Prandtl Numbers
• 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
• Pegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics
• dHybrid: a massively parallel code for hybrid simulations of space plasmas
• A nonlinear theory of the parallel firehose and gyrothermal instabilities in a weakly collisional plasma
• On the Boris solver in particle-in-cell simulation
• Plasma turbulence in the interstellar medium
• The Heating of Test Particles in Numerical Simulations of Alfvenic Turbulence
• Turbulent dynamo in a collisionless plasma
• Turbulent dynamo in the two-phase interstellar medium
• Fluctuation Dynamo in a Collisionless, Weakly Magnetized Plasma
• Seed magnetic fields in turbulent small-scale dynamos
• SHARP: A Spatially Higher-order, Relativistic Particle-in-Cell Code
• Cold Fronts and Shocks Formed by Gas Streams in Galaxy Clusters
• 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
• Runko: Modern multiphysics toolbox for plasma simulations
• Hybrid codes (massless electron fluid)