arXiv preprint ·
A Nonlinear General Implicit Runge-Kutta Integrator for Multifluid Charged Dust and Gas
- Published
- Source
- arXiv preprint, posted before peer review
- Original
- Read it at arXiv
- Topics
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 advanced together in one implicit source update, avoiding errors in the terminal drift that can arise when the two forces are advanced in separate substeps. The timestep is compared with a drag-Lorentz timescale to select between resolved and stiff branches. The resolved branch uses Gauss-Legendre coefficients, while the stiff branch uses a coefficient set designed to recover the terminal drift set by the drag-Lorentz balance. We also allow the dust stopping time and charge-to-mass ratio to depend on the evolving state, evaluating them separately at the two implicit stage states and solving the resulting nonlinear stage equations by Picard iteration. The source update conserves total gas-dust momentum, and its computational cost scales linearly with the number of dust species. Tests reproduce the expected drag damping and terminal drift, confirm the predicted gyromotion amplitude and phase errors, demonstrate accurate treatment of state-dependent source coefficients, and recover the expected wave, shock, and instability dynamics in dusty MHD flows.
The authors’ abstract, as arXiv lists it. arXiv’s listings, abstracts included, are free for anyone to reuse.
Read it on ReadSearch
Each opens after you create an account, and brings you back to this paper.
More in Astrophysics of Galaxies
Sulphur-bearing molecule survey towards prestellar and starless cores in the Perseus molecular cloud
Hot Molecular Hydrogen in Galaxies at Redshift ~ 0.2-1.5 Revealed by JWST
A radio counterpart of the local little-red-dot J092359+415736
Large-scale Environmental Dependence of AGN Activity Revealed by JWST Pa$α$ Emitters
Mergers may not be the main channel for galaxy stellar mass growth at redshift $z\gtrsim5$: Predictions from cosmological hierarchical models and simulations
Get the few papers worth your time in Astrophysics of Galaxies every morning.
ReadSearch reads the new papers in your field each day and brings you the few worth reading, each with its key finding in a sentence.