Magnetohydrodynamic Models of Axisymmetric Protostellar Jets

Abstract

We present the results of a series of axisymmetric time-dependent magnetohydrodynamic (MHD) simulations of the propagation of cooling, overdense jets. Our numerical models are motivated by the properties of outÑows associated with young stellar objects. A variety of initial Ðeld strengths and con- Ðgurations are explored for both steady and time-variable (pulsed) jets. For the parameters of protostel- lar jets adopted here, even apparently weak magnetic Ðelds with strengths kG in the preshocked BZ60 jet beam can have a signiÐcant e†ect on the dynamics, for example, by altering the density, width, and fragmentation of thin shells formed by cooling gas. Strong toroidal Ðelds (º100 kG) with a radial proÐle that peaks near the surface of the jet result in the accumulation of dense shocked gas in a ““ nose cone ÏÏ at the head of jet. We suggest that this structure is unstable in three dimensions. A linear analysis pre- dicts that axisymmetric pinch modes of the MHD Kelvin-Helmholtz instability should grow only slowly for the highly supermagnetosonic jets studied here ; we Ðnd no evidence for them in our simulations. Some of our models appear unstable to current-driven pinch modes ; however, the resulting pressure and density variations induced in the jet beam are not large, making this mechanism an unlikely source of emission knots in the jet beam. In the case of pulsed jets, radial hoop stresses conÐne shocked jet material in the pulses to the axis, resulting in a higher density in the pulses in comparison to purely hydrodynamic models. In addition, if the magnetic Ðeld strength varies with radius, signiÐcant radial structure is produced in the pulses (the density is strongly axially peaked, for example) even if the density and velocity in the jet follow a constant ““ top-hat ÏÏ proÐle initially.