We present results of a series of magnetohydrodynamic (MHD) and hydro- dynamic (HD) 2.5D simulations of the morphology of outflows driven by nested wide-angle winds - i.e. winds which eminate from a central star as well as from an orbiting accretion disk. While our results are broadly relevent to nested wind systems we have tuned the parameters of the simulations to touch on issues in both Young Stellar Objects and Planetary Nebula studies. In particular our studies connect to open issues in the early evolution of Planetary Nebulae. We find that nested MHD winds exhibit marked morphological differences from the single MHD wind case along both dimensions of the flow. Nested HD winds on the other hand give rise mainly to geometric distortions of an outflow that is topologically similar to the flow arising from a single stellar HD wind. Our MHD results are insensitive to changes in ambient temperature between ionized and un-ionized circumstellar environments. The results are sensitive to the relative mass-loss rates, and to the relative speeds of the stellar and disk winds. We also present synthetic emission maps of both nested MHD and HD simulations. We find that nested MHD winds show knots of emission appearing on-axis that do not appear in the HD case.

A description is given of the algorithms implemented in the AstroBEAR adaptive mesh refinement code for ideal magnetohydrodynamics. The code provides several high resolution, shock capturing schemes which are constructed to maintain conserved quantities of the flow in a finite volume sense. Divergence free magnetic field topologies are maintained to machine precision by collating the components of the magnetic field on a cell-interface staggered grid and utilizing the constrained transport approach for integrating the induction equations. The maintenance of magnetic field topologies on adaptive grids is achieved using prolongation and restriction operators which preserve the divergence and curl of the magnetic field across co-located grids of different resolution. The robustness and correctness of the code is demonstrated by comparing the numerical solution of various tests with analytical solutions or previously published numerical solutions obtained by other codes.

Herbig-Haro (HH) jets are commonly thought of as homogeneous beams of plasma traveling at hypersonic velocities. Structure within jet beams is often attributed to periodic or ``pulsed'' variations of conditions at the jet source. Simulations based on this scenario result in knots extending across the jet diameter. Observations and recent high energy density laboratory experiments shed new light on structures below this scale and indicate they may be important for understanding the fundamentals of jet dynamics. In this paper we offer an alternative to ``pulsed'' models of protostellar jets. Using direct numerical simulations we explore the possibility that jets are chains of sub-radial clumps propagating through a moving inter-clump medium. Our models explore an idealization of this scenario by injecting small ($r<r_{jet}$), dense ($\rho>\rho_{jet}$) spheres embedded in an otherwise smooth inter-clump jet flow. The spheres are initialized with velocities differing from the jet velocity by $\sim15$%. We find the consequences of shifting from homogeneous to heterogeneous flows are significant as clumps interact with each other and with the inter-clump medium in a variety of ways. Structures which mimic what is expected from pulsed-jet models can form, as can previously unseen ``sub-radial'' behaviors including backward facing bow shocks and off-axis working surfaces. While these small-scale structures have not been seen before in simulation studies, they are found in high resolution jet observations. We discuss implications of our simulations for the interpretation of protostellar jets with regard to characterization of knots by a ``lifetime'' or ``velocity history'' approach as well as linking observed structures with central engines which produce the jets.

Many proto-planetary nebulae (PPN) appear as narrow collimated structures sometimes showing multiple, roughly aligned lobes. In addition, many PPN flows have been shown to have short acceleration times. In this paper we explore whether jet or ``bullet'' (a massive clump) models fit the observations of individual collimated lobes adequately by comparing simulations of both radiatively cooled (stable) jets and bullets. We find that the clump model is somewhat favored over jets because (1) it leads to greater collimation of outflows (2) it accounts better and more naturally for ring-like structures observed in the PPN CRL 618, and (3) it is more successful in reproducing the Hubble-flow character of observed kinematics in some PPN. In addition, bullets naturally account for observed multipolar flows, since the likely MHD launch mechanisms required to drive outflows make multiple non-aligned jets unlikely. Thus we argue that PPN outflows may be driven by explosive MHD launch mechanisms such as those discussed in the context of supernovae (SNe) and gamma-ray bursts(GRB).

We have performed 2D hydrodynamic simulations of a pulsed astrophysical jet propagating through a medium that is populated with spherical inhomogeneities, or "clumps," which are smaller than the jet radius. The clumps are seen to affect the jet in several ways, such as impeding jet propagation and deflecting the jet off-axis. While there has been some debate as to the prevalence of these types of condensations in the ISM or in molecular clouds, the exploration of this region of parameter space nonetheless both shows the potential for these clumps to disrupt astrophysical jets and yields results which recover aspects of recent observations of Herbig-Haro objects. We find that the propagation of the jet and the vorticity induced in the clump/ambient medium correlate well with a "dynamic filling function" $f_d$ across all the simulations.

We present two dimensional cylindrically symmetric hydrodynamic simulations and synthetic emission maps of a stellar wind propagating into an infalling, rotating environment. The resulting outflow morphology, collimation and stability observed in these simulations have relevance to the study of young stellar objects, Herbig-Haro jets and molecular outflows. Our code follows hydrogen gas with molecular, atomic and ionic components tracking the associated time dependent molecular chemistry and ionization dynamics with radiative cooling appropriate for a dense molecular gas. We present tests of the code as well as new simulations which indicate the presence of instabilities in the wind-blown bubble's swept-up shell.

Many astrophysical flows occur in inhomogeneous media. We briefly discuss some general properties of the adiabatic and radiative inhomogeneous systems and discuss the relevance of those properties to the planetary nebulae systems. We then focus on radiative hypersonic bullets and the applicability of this model to planetary and protoplanetary systems such as CRL 618, NGC 6543, Hen 3-1475.

We present the results of a numerical study focusing on the propagation of a hypersonic bullet subject to radiative cooling. Our goal is to explore the feasibility of such a model for the formation of ``strings'' observed in the Eta Carinae Homunculus nebula. Our simulations were performed in cylindrical symmetry with the adaptive mesh refinement code AstroBEAR. The radiative cooling of the system was followed using the cooling curve by Dalgarno & McCray (1972). In this letter we discuss the evolution and overall morphology of the system as well as key kinematic properties. We find that radiative bullets can produce structures with properties similar to those of the Eta Carinae strings, i.e. high length-to-width ratios and Hubble-type flows in the form of a linear velocity increase from the base of the wake to the bullet head. These features, along with the appearance of periodic ``ring-like'' structures, may also make this model applicable to other astrophysical systems such as planetary nebulae, e.g. CRL 618 and NGC 6543, young stellar objects, etc.