Astronomical observations, analytical solutions and numerical simulations have provided the building blocks to formulate the current theory of young stellar object jets. Although each approach has made great progress independently, it is only during the last decade that significant efforts are being made to bring the separate pieces together. Building on previous work that combined analytical solutions and numerical simulations, we apply a sophisticated cooling function to incorporate optically thin energy losses in the dynamics. On the one hand, this allows a self-consistent treatment of the jet evolution and on the other, it provides the necessary data to generate synthetic emission maps. Firstly, analytical disk and stellar outflow solutions are properly combined to initialize numerical two-component jet models inside the computational box. Secondly, magneto-hydrodynamical simulations are performed in 2.5D, following properly the ionization and recombination of a maximum of $29$ ions. Finally, the outputs are post-processed to produce artificial observational data. The first two-component jet simulations, based on analytical models, that include ionization and optically thin radiation losses demonstrate promising results for modeling specific young stellar object outflows. The generation of synthetic emission maps provides the link to observations, as well as the necessary feedback for the further improvement of the available models.

We investigate supersonic, axisymmetric magnetohydrodynamic (MHD) jets with a time-dependent injection velocity by numerical simulations with the PLUTO code. Using a comprehensive set of parameters, we explore different jet configurations in the attempt to construct models that can be directly compared to observational data of microjets. In particular, we focus our attention on the emitting properties of traveling knots and construct, at the same time, accurate line intensity ratios and surface brightness maps. Direct comparison of the resulting brightness and line intensity ratios distributions with observational data of microjets shows that a closer match can be obtained only when the jet material is pre-ionized to some degree. A very likely source for a pre-ionized medium is photoionization by X-ray flux coming from the central object.

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. 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. 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. 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. 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.

Time-dependent cooling processes are of paramount importance in the evolution of astrophysical gaseous nebulae and, in particular, when radiative shocks are present. The present work introduces a necessary set of tools that can be used to model radiative astrophysical flows in the optically-thin plasma limit. We aim to provide reliable and accurate predictions of emission line ratios and radiative cooling losses in astrophysical simulations of shocked flows. Moreover, we discuss numerical implementation aspects to ease future improvements and implementation in other MHD numerical codes. The most important source of radiative cooling for our plasma conditions comes from the collisionally-excited line radiation. We evolve a chemical network, including 29 ion species, to compute the ionization balance in non-equilibrium conditions. After a series of validations and tests, typical astrophysical setups are simulated in 1D and 2D, employing both the present cooling model and a simplified one. The influence of the cooling model on structure morphologies can become important, especially for emission line diagnostic purposes.

We present a new numerical code, PLUTO, for the solution of hypersonic flows in 1, 2 and 3 spatial dimensions and different systems of coordinates. The code provides a multi-physics, multi-algorithm modular environment particularly oriented towards the treatment of astrophysical flows in presence of discontinuities. Different hydrodynamic modules and algorithms may be independently selected to properly describe Newtonian, relativistic, MHD or relativistic MHD fluids. The modular structure exploits a general framework for integrating a system of conservation laws, built on modern Godunov-type shock-capturing schemes. Although a plethora of numerical methods has been successfully developed over the past two decades, the vast majority shares a common discretization recipe, involving three general steps: a piecewise polynomial reconstruction followed by the solution of Riemann problems at zone interfaces and a final evolution stage. We have checked and validated the code against several benchmarks available in literature. Test problems in 1, 2 and 3 dimensions are discussed.