An Axisymmetric, Radiative Bow Shock Model with a Realistic Treatment of Ionization and Cooling
arXiv:astro-ph/9610134 · doi:10.1086/312987
Abstract
We have chosen a reduced set of 18 ionization rate equations (for ions of H, C, N, O, S and Ne), which allow us to obtain a moderately accurate estimate of the non-equilibrium radiative cooling function. We evaluate the accuracy of this approach by comparing our cooling function with previous calculations, computed with a more extended set of ions, for the case of gas that cools from a high temperature at constant density. We also compute steady, plane shock models, which we find to compare well with models calculated with much more detailed microphysics. Using our reduced set of rate equations, we present a simulation of a radiative bow shock formed by a supersonic, plane stream impinging on a rigid sphere. The parameters for the calculation are chosen as to approximately represent typical values found for Herbig-Haro objects, and to give a cooling distance to bow shock radius ratio of 1/10. This simulation is done with an adaptive grid code, which allows a reasonably high resolution (with greater than or equal to 25 points) of the standoff distance between the bow shock and the rigid obstacle. Contrary to the standard expectation, we find that the bow shock standoff distance in the stagnation region is considerably shorter than the cooling distance behind a plane-parallel shock. Also, the centrifugal pressure is found to play an important role in determining the structure of the recombin- ation region. This appears to partially invalidate previous bow shock models based on a ``quasi-1D'' approach, at least for the particular parameters chosen for the present simulation. Finally, we present tabulations of the cooling rates that have been used (for the different species), in order to facilitate the inclusion of this treatment in other gasdynamic codes.
28 pages, AASTeX v.4.0, 10 Postscript figures; To appear in ApJS. Also available at ftp://www.astro.su.se/pub/supernova/preprints.html#RML
Cited by in corpus (26)
- PLUTO: a Numerical Code for Computational Astrophysics
- Photoevaporation of Cosmological Minihalos during Reionization
- C^2-Ray: A new method for photon-conserving transport of ionizing radiation
- Cosmological Radiative Transfer Codes Comparison Project I: The Static Density Field Tests
- Are the Models for Type Ia Supernova Progenitors Consistent with the Properties of Supernova Remnants?
- Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests
- Jet/cloud collision, 3D gasdynamic simulations of HH 110
- Incorporating the molecular gas phase in galaxy-size numerical simulations: first applications in dwarf galaxies
- Models of the circumstellar medium of evolving, massive runaway stars moving through the Galactic plane
- Periodic bursts of Star Formation in Irregular Galaxies
- The X-ray Luminosities of HH Objects
- MHD Models of Axisymmetric Protostellar Jets
- Temporal evolution of magnetic molecular shocks I. Moving grid simulations
- Nonequilibrium ionization states and cooling rates of the photoionized enriched gas
- Simulating radiative astrophysical flows with the PLUTO code: A non-equilibrium, multi-species cooling function
- A 3-mode, Variable Velocity Jet Model for HH 34
- Stellar wind bubbles around WR and [WR] stars
- Radiative transfer of energetic photons: X-rays and helium ionization in C2-Ray
- Dynamics of H II regions around exiled O stars
- Hydrodynamical models of cometary HII regions
- Interacting Jets from Binary Protostars
- Mixing of materials in magnetised core-collapse supernova remnants
- Fate of clumps in damped Ly systems
- Simulating Radiative Magnetohydrodynamical Flows with AstroBEAR: Implementation and Applications of Non-equilibrium Cooling
- A multi-ion non-equilibrium solver for ionised astrophysical plasmas with arbitrary elemental abundances
- Radio polarization from runaway star bowshocks-I. The general case