Volume 507, Number 2, November IV 2009
|Page(s)||581 - 588|
|Published online||17 September 2009|
Time-dependent MHD shocks and line intensity ratios in the HH 30 jet: a focus on cooling function and numerical resolution
RCAPA – Department of Physics, University of Bucharest, CP MG-11, RO-077125, Bucureşti-Mãgurele, Romania
2 Dipartimento di Fisica Generale dell'Università, via Pietro Giuria 1, 10125 Torino, Italy e-mail: firstname.lastname@example.org; [tesilean;massaglia;mignone]@ph.unito.it
3 INAF – Osservatorio Astronomico di Torino, Viale Osservatorio 20, 10025 Pino Torinese, Italy e-mail: email@example.com
4 INAF – Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy e-mail: firstname.lastname@example.org
Accepted: 18 August 2009
Context. The coupling between time-dependent, multidimensional MHD numerical codes and radiative line emission is of utmost importance in the studies of the interplay between dynamical and radiative processes in many astrophysical environments, with particular interest for problems involving radiative shocks. There is a widespread consensus that line emitting knots observed in Herbig-Haro jets can be interpreted as radiative shocks. Velocity perturbations at the jet base steepen into shocks to emit the observed spectra. To derive the observable characteristics of the emitted spectra, such as line intensity ratios, one has to study physical processes that involve the solution of the MHD equations coupled with radiative cooling in non-equilibrium conditions.
Aims. In this paper we address two different aspects relevant to the time-dependent calculations of the line intensity ratios of forbidden transitions, resulting from the excitation by planar, time-dependent radiative shocks traveling in a stratified medium. The first one concerns the impact of the radiation and ionization processes included in the cooling model, and the second one the effects of the numerical grid resolution.
Methods. Dealing with both dynamical and radiative processes in the same numerical scheme means to treat phenomena characterized by different time and length scales. This may be especially arduous and computationally expensive when discontinuities are involved, such as in the case of shocks. Adaptive mesh refinement (AMR) methods have been introduced in order to alleviate these difficulties. In this paper we apply the AMR methodology to the treatment of radiating shocks and show how this method is able to vastly reduce the integration time.
Results. The technique is applied to the knots of the HH 30 jet to obtain the observed line intensity ratios and derive the physical parameters, such as density, temperature and ionization fraction. We consider the impact of two different cooling functions and different grid resolutions on the results.
Conclusions. We conclude that the use of different cooling routines has effects on results whose weight depends upon the line ratio considered. Moreover, we find the minimum numerical resolution of the simulation grid behind the shock to achieve convergence in the results. This is crucial for the forthcoming 2D calculations of radiative shocks.
Key words: ISM: jets and outflows / ISM: Herbig-Haro objects / magnetohydrodynamics (MHD) / shock waves / methods: numerical
© ESO, 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.