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 compute numerical simulations of spherical collapse triggered by a slow increase in external pressure. We compare isothermal models to models including cooling with a simple but self-consistent treatment of the coupling between gas, grains and radiation field temperatures. The hydrostatic equilibrium appears to hold past the marginally stable state, until the collapse proceeds. The last hydrostatic state before collapse has a lower central gas temperature in the centre due to the enhanced coupling between gas, grains and radiation field. This results in slightly lower pressure gradients in the bulk of the envelope which is hence slightly more extended than in the isothermal case. Due to the sensitivity of the collapse on these initial conditions, protostellar infall velocities in the envelope turn out to be much slower in the case with cooling.Our models also compute the radiative transfer and a rather large chemical network coupled to gas dynamics. However, we note that the steady-state chemisorption of CO is sufficient to provide an accurate cooling function of the gas. This justifies the use of post-processing techniques to account for the abundance of observed molecules. Existing observations of infall signatures put very stringent constraints on the kinematics and temperature profile of the class 0 protostar IRAM 04191+1522. We show that isothermal models fail to account for the innermost slow infall motions observed, even with the most hydrostatic initial conditions. In contrast, models with cooling reproduce the general shape of the temperature profile inferred from observations and are in much better agreement with the infall signatures in the inner 3000 AU.

We investigate the influence of non equilibrium thermodynamics on cosmological structure formation. In this paper, we consider the collapse of planar perturbations usually called "Zel'dovich pancakes". We have developed for that purpose a new two fluids (gas and dark matter) hydrodynamical code, with three different thermodynamical species: electrons, ions and neutral particles (T_e\ne T_i \ne T_n). We describe in details the complex structure of accretion shock waves. We include several relevant processes for a low density, high temperature, collisional plasma such as non-equilibrium chemical reactions, cooling, shock heating, thermal energy equipartition between electrons, ions and neutral particles and electronic conduction. We find two different regions in the pancake structure: a thermal precursor ahead of the compression front and an equipartition wave after the compression front where electrons and ions temperatures differ significantly. This complex structure may have two interesting consequences: pre-heating of unshocked regions in the vicinity of massive X-ray clusters and ions and electrons temperatures differences in the outer regions of X-rays clusters.

We investigate the role of radiative cooling within the core of large X-ray clusters using multi-phase hydrodynamics. We developed for that purpose a spherically symmetric hydrodynamical code, coupled to a "fluid model" that describes accurately the dark matter component. Cooling is included using a self-consistent multi-phase approach, leading to "cooled gas" mass deposition throughout the flow. We simulate the collapse and the subsequent evolution of a Coma-like X-ray cluster, avoiding the well-known "cooling catastrophe". The total mass profile of our simulated cluster is very similar to the "universal" profile proposed by Navarro, Frenk & White (1995). More interestingly, we also obtain a quasi-isothermal temperature profile, which is a direct consequence of multi-phase cooling within such a potential well.

We simulate the formation of a large X-ray cluster using a fully 3D hydrodynamical code coupled to a Particle-Mesh scheme which models the dark matter component. We focus on a possible decoupling between electrons and ions temperatures. We then solve the energy transfer equations between electrons, ions and neutrals without assuming thermal equilibrium between the three gases (T_e <> T_i <> T_n). We solve self-consistently the chemical equations for an hydrogen/helium primordial plasma without assuming ionization-recombination equilibrium. We find that the electron temperature differs from the true dynamical temperature by 20% at the Virial radius of our simulated cluster. This could lead marginally to an underestimate of the total mass in the outer regions of large X-ray clusters.

We investigate the formation of spherical cosmological structures following both dark matter and gas components. We focus on the dynamical aspect of the collapse assuming an adiabatic, $\gamma = 5/3$, fully ionized primordial plasma. We use for that purpose a fully Lagrangian hydrodynamical code designed to describe highly compressible flows in spherical geometry. We investigate also a "fluid approach" to describe the mean physical quantities of the dark matter flow. We test its validity for a wide range of initial density contrast. We show that an homogeneous isentropic core forms in the gas distribution, surrounded by a self-similar hydrostatic halo, with much higher entropy generated by shock dissipation. We derive analytical expressions for the size, density and temperature of the core, as well as for the surrounding halo. We show that, unless very efficient heating processes occur in the intergalactic medium, we are unable to reproduce within adiabatic models the typical core sizes in X-ray clusters. We also show that, for dynamical reasons only, the gas distribution is naturally antibiased relative to the total mass distribution, without invoking any reheating processes. This could explain why the gas fraction increases with radius in very large X-ray clusters. As a preparation for the next study devoted to the thermodynamical aspect of the collapse, we investigate the initial entropy level required to solve the core problem in X-ray clusters.

We reconsider the collapse of cosmic structures in an Einstein-de Sitter Universe, using the self similar initial conditions of Fillmore & Goldreich (1984). We first derive a new approximation to describe the dark matter dynamics in spherical geometry, that we refer to the "fluid approach". This method enables us to recover the self-similarity solutions of Fillmore & Goldreich for dark matter. We derive also new self-similarity solutions for the gas. We thus compare directly gas and dark matter dynamics, focusing on the differences due to their different dimensionalities in velocity space. This work may have interesting consequences for gas and dark matter distributions in large galaxy clusters, allowing to explain why the total mass profile is always steeper than the X-ray gas profile. We discuss also the shape of the dark matter density profile found in N-body simulations in terms of a change of dimensionality in the dark matter velocity space. The stable clustering hypothesis has been finally considered in the light of this analytical approach.

We present and test a general-purpose code, called PPASPH, for evolving self-gravitating fluids in astrophysics, both with and without a collisionless component. In PPASPH, hydrodynamical properties are computed by using the SPH (Smoothed Particle Hydrodynamics) method while, unlike most previous implementations of SPH, gravitational forces are computed by a PP (Particle-Particle) approach. another important feature of this code is that hydrodynamics equations optionally include the correction terms appearing when $h(t,\br)$ is not constant. Our code has been implemented by using the data parallel programming model on CM5(CM). PPASPH has been applied to study the importance of adaptative smoothing correction terms on the entropy conservation. We confirm Hernquist's (1993) interpretation of the entropy violation observed in previous SPH simulations as a result of having neglected these terms. An improvement on the entropy conservation is not found by merely considering larger numbers of particles or different $N_S$ choices. The correct continuum description is only obtained if the $\bn h$ correction terms are included. Otherwise, the entropy conservation is always rather poor as compared to that found for the total energy.