1,321 research outputs found
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
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
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
A 3-mode, Variable Velocity Jet Model for HH 34
Variable ejection velocity jet models can qualitatively explain the
appearance of successive working surfaces in Herbig-Haro (HH) jets. This paper
presents an attempt to explore which features of the HH-34 jet can indeed be
reproduced by such a model. From previously published data on this object, we
find evidence for the existence of a 3-mode ejection velocity variability, and
then explore the implications of such a variability. From simple, analytic
considerations it is possible to show that the longer period modes produce a
modulation on the shorter period modes, resulting in the formation of
``trains'' of multiple knots. The knots observed close to the source of HH-34
could correspond to such a structure. Finally, a numerical simulation with the
ejection velocity variability deduced from the HH-34 data is computed. This
numerical simulation shows a quite remarkable resemblance with the observed
properties of the HH-34 jet.Comment: 28 pages LaTex, 10 postscript figure
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
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
Filaments in Galactic Winds Driven by Young Stellar Clusters
The starburst galaxy M82 shows a system of H-emitting filaments which
extend to each side of the galactic disk. We model these filaments as the
result of the interaction between the winds from a distribution of Super
Stellar Clusters (SSCs). We first derive the condition necessary for producing
a radiative interaction between the cluster winds (a condition which is met by
the SSC distribution of M82). We then compute 3D simulations for SSC wind
distributions which satisfy the condition for a radiative interaction, and also
for distributions which do not satisfy this condition. We find that the highly
radiative models, that result from the interaction of high metallicity cluster
winds, produce a structure of H emitting filaments, which qualitatively
agrees with the observations of the M82, while the non-radiative SSC wind
interaction models do not produce filamentary structures. Therefore, our
criterion for radiative interactions (which depends on the mass loss rate and
the terminal velocity of the SSC winds, and the mean separation between SSCs)
can be used to predict whether or not an observed galaxy should have associated
H emitting filaments.Comment: 10 pages, 6 Figures. ApJ Accepted, August 7, 200
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
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
- …