Abstract:The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run times and often to set-ups of less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conserving multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison with linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic time scales.

Use of a Hybrid Code for Global-Scale Plasma Simulation

Electromagnetic Formulation of Global Gyrokinetic Particle Simulation in Toroidal Geometry

Coupled MHD -- Hybrid Simulations of Space Plasmas

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

On the Development of Global Plasma-Ionosphere FDTD Algorithms for Electromagnetic Calculations in the Earth-ionosphere System

Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration

Maximus: a Hybrid Particle-in-Cell Code for Microscopic Modeling of Collisionless Plasmas

A new hybrid code (CHIEF) implementing the inertial electron fluid equation without approximation

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

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

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

A Kinetic-magnetohydrodynamic Model with Adaptive Mesh Refinement for Modeling Heliosphere Neutral-plasma Interaction

Hybrid particle-spectral method for kinetic plasma simulations

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

Three‐dimensional hybrid gyrokinetic‐magnetohydrodynamics simulation

Hybrid GRMHD and Force-Free Simulations of Black Hole Accretion

Simulating the Photospheric to Coronal Plasma Using Magnetohydrodynamic Characteristics. I. Data-driven Boundary Conditions

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

The muphyII Code: Multiphysics Plasma Simulation on Large HPC Systems

Importance of accurate consideration of the electron inertia in hybrid-kinetic simulations of collisionless plasma turbulence: 1. The 2D limit