Open AccessA numerical scheme for multifluid magnetohydrodynamics
Author(s) -
Falle S. A. E. G.
Publication year - 2003
Publication title -
monthly notices of the royal astronomical society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.058
H-Index - 383
eISSN - 1365-2966
pISSN - 0035-8711
DOI - 10.1046/j.1365-8711.2003.06908.x
Subject(s) - physics , ambipolar diffusion , magnetohydrodynamics , inertia , shock (circulatory) , classical mechanics , magnetic field , flow (mathematics) , molecular cloud , limit (mathematics) , statistical physics , mechanics , scheme (mathematics) , astrophysics , plasma , quantum mechanics , mathematical analysis , stars , mathematics , medicine
This paper describes a numerical scheme for multifluid hydrodynamics in the limit of small mass densities of the charged particles. The inertia of the charged particles can then be neglected, which makes it possible to write an evolution equation for the magnetic field that can be solved using an implicit scheme. This avoids the severe restriction on the stable timestep that would otherwise arise at high resolution, or when the Hall effect is large. Numerical tests show that the scheme can accurately model steady multifluid shock structures both with and without subshocks. Although the emphasis is on shocks in molecular clouds, a multidimensional version of this code could be applied to any astrophysical flow in which ambipolar diffusion or the Hall effect, or both, play a significant role.