A Nonlinear General Implicit Runge-Kutta Integrator for Multifluid Charged Dust and Gas

Abstract

Charged dust in magnetized gas is subject to aerodynamic drag and the Lorentz force, both of which can be difficult to model numerically in a fluid simulation where their characteristic timescales are shorter than the timestep set by the Courant condition. We present a two-stage general implicit Runge-Kutta (GIRK) source integrator for multifluid charged dust coupled to an ideal magnetohydrodynamic (MHD) gas, implemented in the open source GPU-accelerated quokka code. Drag and Lorentz forces are…