  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.
