Simulating radiative magnetohydrodynamical flows withastrobear: implementation and applications of non-equilibrium cooling

Abstract

Radiative cooling plays a crucial role in the dynamics of many astrophysical flows, and is particularly important in the dense shocked gas within Herbig-Haro (HH) objects and stellar jets. Simulating cooling processes accurately is necessary to compare numerical simulations with existing and planned observations of HH objects, such as those from the Hubble Space Telescope and the James Webb Space Telescope. In this paper we discuss a new, non-equilibrium cooling scheme we have implemented into the 3-D magnetohydrodynamic (MHD) code AstroBEAR. The new cooling function includes ionization, recombination, and excitation of all the important atomic species that cool below 10000 K. We tested the routine by comparing its predictions with those from the well-tested 1-D Cox-Raymond shock code (Raymond 1979). The results show thatAstroBEAR accurately tracks the ionization fraction, temperature, and other MHD variables for all low-velocity (.90 km/s) magnetized radiative shock waves. The new routine allows us to predict synthetic emission maps in all the bright forbidden and permitted lines observed in stellar jets, including H{\alpha}, [NII], [OI], and [SII]. We present an example as to how these synthetic maps facilitate a direct comparison with narrowband images of HH objects.