Abstract:We introduce $\texttt{A}$strophysical $\texttt{H}$ybrid-$\texttt{K}$inetic simulations with the $\texttt{flASH}$ code ($\texttt{AHKASH}$) -- a new Hybrid particle-in-cell (PIC) code developed within the framework of the multi-physics code $\texttt{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 $\texttt{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 $\sim1.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.

Coupling multi-fluid dynamics equipped with Landau closures to the particle-in-cell method

Two-way coupling of a global Hall magnetohydrodynamics model with a local implicit particle-in-cell model

PolyPIC: The Polymorphic-Particle-in-Cell Method for Fluid-Kinetic Coupling

High-order Implicit Particle-in-cell Method for Plasma Simulations at Solid Densities.

Mitigation of Numerical Cerenkov Radiation and Instability Using a Hybrid Finite Difference-Fft Maxwell Solver and a Local Charge Conserving Current Deposit

Hybrid fluid–particle modeling of shock-driven hydrodynamic instabilities in a plasma

MAGNETOHYDRODYNAMIC-PARTICLE-IN-CELL METHOD FOR COUPLING COSMIC RAYS WITH A THERMAL PLASMA: APPLICATION TO NON-RELATIVISTIC SHOCKS

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

Full-PIC Code Validation and Comparison against Fluid Models on Plasma Plume Expansions IEPC-2017-230

The magnetohydrodynamic-particle-in-cell module in athena plus plus : implementation and code tests

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

Canonical Symplectic Particle-In-Cell Method For Long-Term Large-Scale Simulations Of The Vlasov-Maxwell Equations

Magnetohydrodynamic with Adaptively Embedded Particle-in-Cell model: MHD-AEPIC

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

QuickPIC: a highly efficient fully parallelized PIC code for plasma-based acceleration

Coupling of Fluid and Particle-in-Cell Simulations of Ambipolar Plasma Thrusters

Implicit electrostatic particle-in-cell/Monte Carlo simulation for the magnetized plasma: Algorithms and application in gas-inductive breakdown*

Hybrid-VPIC: an Open-Source Kinetic/Fluid Hybrid Particle-in-Cell Code

Explicit Structure-Preserving Geometric Particle-in-Cell Algorithm in Curvilinear Orthogonal Coordinate Systems and Its Applications to Whole-Device 6D Kinetic Simulations of Tokamak Physics

Implementation of a Hybrid Particle Code with a PIC Description in R–z and a Gridless Description in Φ into OSIRIS

Controlling the Numerical Cerenkov Instability in PIC simulations using a customized finite difference Maxwell solver and a local FFT based current correction