Side-entrainment in a jet embedded in a sidewind

Numerical simulations of HH jets never show side-entrainment of environmental material into the jet beam. This is because the bow shock associated with the jet head pushes the surrounding environment into a dense shell, which is never in direct contact with the sides of the jet beam. We present 3D simulations in which a side-streaming motion (representing the motion of the outflow source through the surrounding medium) pushes the post-bow shock shell into direct contact with the jet beam. This is a possible mechanism for modelling well collimated "molecular jets" as an atomic/ionic flow which entrains molecules initially present only in the surrounding environment.Comment: 8 pages, 12 figures, 1 table, accepted for publication in Ap

Numerical Simulations of HH 555

We present 3D gasdynamic simulations of the Herbig Haro object HH 555. HH 555 is a bipolar jet emerging from the tip of an elephant trunk entering the Pelican Nebula from the adjacent molecular cloud. Both beams of HH 555 are curved away from the center of the H II region. This indicates that they are being deflected by a side-wind probably coming from a star located inside the nebula or by the expansion of the nebula itself. HH 555 is most likely an irradiated jet emerging from a highly embedded protostar, which has not yet been detected. In our simulations we vary the incident photon flux, which in one of our models is equal to the flux coming from a star 1 pc away emitting 5x10^48 ionizing (i. e., with energies above the H Lyman limit) photons per second. An external, plane-parallel flow (a ``side-wind'') is coming from the same direction as the photoionizing flux. We have made four simulations, decreasing the photon flux by a factor of 10 in each simulation. We discuss the properties of the flow and we compute Halpha emission maps (integrated along lines of sight). We show that the level of the incident photon flux has an important influence on the shape and visibility of the jet. If the flux is very high, it causes a strong evaporation of the neutral clump, producing a photoevaporated wind traveling in the direction opposite to the incident flow. The interaction of the two flows creates a double shock ``working surface'' around the clump protecting it and the jet from the external flow. The jet only starts to curve when it penetrates through the working surface.Comment: 14 pages, 4 figures, accepted by Ap

The 3-Dimensional Structure of HH 32 from GMOS IFU Spetroscopy

We present new high resolution spectroscopic observations of the Herbig-Haro object HH 32 from System Verification observations made with the GMOS IFU at Gemini North Observatory. The 3D spectral data covers a 8''.7 x 5''.85 spatial field and 4820 - 7040 Angstrom spectral region centered on the HH~32 A knot complex. We show the position-dependent line profiles and radial velocity channel maps of the Halpha line, as well as line ratio velocity channel maps of [OIII]5007/Halpha, [OI]6300/Halpha, [NII]6583/Halpha, [SII](6716+6730)/Halpha and [SII]6716/6730. We find that the line emission and the line ratios vary significantly on spatial scales of ~1'' and over velocities of ~50 km/s. A ``3/2-D'' bow shock model is qualitatively successful at reproducing the general features of the radial velocity channel maps, but it does not show the same complexity as the data and it fails to reproduce the line ratios in our high spatial resolution maps. The observations of HH 32 A show two or three superimposed bow shocks with separations of ~3'', which we interpret as evidence of a line of sight superposition of two or three working surfaces located along the redshifted body of the HH 32 outflow.Comment: Accepted for Publication in the Astronomical Journal (January 2004

A model of Mira's cometary head/tail entering the Local Bubble

We model the cometary structure around Mira as the interaction of an AGB wind from Mira A, and a streaming environment. Our simulations introduce the following new element: we assume that after 200 kyr of evolution in a dense environment Mira entered the Local Bubble (low density coronal gas). As Mira enters the bubble, the head of the comet expands quite rapidly, while the tail remains well collimated for a 100 kyr timescale. The result is a broad-head/narrow-tail structure that resembles the observed morphology of Mira's comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail.Comment: 12 pages, 8 figures (4 in color). Accepted for publication in The Astrophysical Journa

Emission lines from rotating proto-stellar jets with variable velocity profiles. I. Three-dimensional numerical simulation of the non-magnetic case

Using the Yguazu-a three-dimensional hydrodynamic code, we have computed a set of numerical simulations of heavy, supersonic, radiatively cooling jets including variabilities in both the ejection direction (precession) and the jet velocity (intermittence). In order to investigate the effects of jet rotation on the shape of the line profiles, we also introduce an initial toroidal rotation velocity profile, in agreement with some recent observational evidence found in jets from T Tauri stars which seems to support the presence of a rotation velocity pattern inside the jet beam, near the jet production region. Since the Yguazu-a code includes an atomic/ionic network, we are able to compute the emission coefficients for several emission lines, and we generate line profiles for the H, [O I]6300, [S II]6716 and [N II]6548 lines. Using initial parameters that are suitable for the DG Tau microjet, we show that the computed radial velocity shift for the medium-velocity component of the line profile as a function of distance from the jet axis is strikingly similar for rotating and non-rotating jet models. These findings lead us to put forward some caveats on the interpretation of the observed radial velocity distribution from a few outflows from young stellar objects, and we claim that these data should not be directly used as a doubtless confirmation of the magnetocentrifugal wind acceleration models.Comment: 15 pages, 8 figures. Accepted to publication in Astronomy and Astrophysic

A latitude-dependent wind model for Mira's cometary head

We present a 3D numerical simulation of the recently discovered cometary structure produced as Mira travels through the galactic ISM. In our simulation, we consider that Mira ejects a steady, latitude-dependent wind, which interacts with a homogeneous, streaming environment. The axisymmetry of the problem is broken by the lack of alignment between the direction of the relative motion of the environment and the polar axis of the latitude-dependent wind. With this model, we are able to produce a cometary head with a ``double bow shock'' which agrees well with the structure of the head of Mira's comet. We therefore conclude that a time-dependence in the ejected wind is not required for reproducing the observed double bow shock.Comment: 4 pages, 4 figures, accepted for publication in ApJ

The precession of the giant HH34 outflow: a possible jet deceleration mechanism

The giant jets represent a fundamental trace of the historical evolution of the outflow activity over timescales which are comparable to the accretion time of the outflow sources in their main protostellar phase. The study of such huge jets provides the possibility of retrieving important elements related to the life of the outflow sources. In this paper, we study the role of precession (combined with jet velocity-variability and the resulting enhanced interaction with the surrounding environment) as a deceleration mechanism for giant jets using a numerical approach. We obtain predictions of H alpha intensity maps and position-velocity diagrams from 3D simulations of the giant HH 34 jet (including an appropriate ejection velocity time-variability and a precession of the outflow axis), and we compare them with previously published observations of this object. Our simulations represent a step forward from previous numerical studies of HH objects, in that the use of a 7-level, binary adaptive grid has allowed us to compute models which appropiately cover all relevant scales of a giant jet, from the ~ 100 AU jet radius close to the source to the ~ 1 pc length of the outflow. A good qualitative and quantitative agreement is found between the model predictions and the observations. Moreover, we show that a critical parameter for obtaining a better or worse agreement with the observations is the ratio rho_j/rho_a between the jet and the environmental densities. The implications of this result in the context of the current star formation models are discussed (ABRIDGED).Comment: 19 pages, 8 eps figs.,uses aaspp4; accepted by the Ap

Interaction of Infall and Winds in Young Stellar Objects

The interaction of a stellar or disk wind with a collapsing environment holds promise for explaining a variety of outflow phenomena observed around young stars. In this paper we present the first simulations of these interactions. The focus here is on exploring how ram pressure balance between wind and ambient gas and post-shock cooling affects the shape of the resulting outflows. In our models we explore the role of ram pressure and cooling by holding the wind speed constant and adjusting the ratio of the inflow mass flux to the wind mass flux (Mdot_a/Mdot_w) Assuming non-spherical cloud collapse, we find that relatively strong winds can carve out wide, conical outflow cavities and that relatively weak winds can be strongly collimated into jet-like structures. If the winds become weak enough, they can be cut off entirely by the infalling environment. We identify discrepancies between results from standard snowplow models and those presented here that have important implications for molecular outflows. We also present mass vs. velocity curves for comparison with observations.Comment: 35 pages, 11 figures (PNG and EPS

Time-dependent ejection velocity model for the outflow of Hen 3--1475

We present 2D axisymmetric and 3D numerical simulations of the proto-planetary nebula Hen 3-1475, which is characterized by a remarkably highly collimated optical jet, formed by a string of shock-excited knots along the axis of the nebula. It has recently been suggested that the kinematical and morphological properties of the Hen 3-1475 jet could be the result of an ejection variability of the central source (Riera et al. 2003). The observations suggest a periodic variability of the ejection velocity superimposed on a smoothly increasing ejection velocity ramp. From our numerical simulations, we have obtained intensity maps (for different optical emission lines) and position-velocity diagrams, in order to make a direct comparison with the HST observations of this object. Our numerical study allows us to conclude that a model of a precessing jet with a time-dependent ejection velocity, which is propagating into an ISM previously perturbed by an AGB wind, can succesfully explain both the morphological and the kinematical characteristics of this proto-planetary nebula.Comment: Astronomy and Astrophysics (accepted) (8 figures