Deflagration-to-detonation transition (DDT) can occur in environments ranging from experimental and industrial systems to astrophysical thermonuclear (type Ia) supernovae explosions. Substantial progress has been made in explaining the nature of DDT in confined systems with walls, internal obstacles, or pre-existing shocks. It remains unclear, however, whether DDT can occur in unconfined media. Here we use direct numerical simulations (DNS) to show that for high enough turbulent intensities unconfined, subsonic, premixed, turbulent flames are inherently unstable to DDT. The associated mechanism, based on the nonsteady evolution of flames faster than the Chapman-Jouguet deflagrations, is qualitatively different from the traditionally suggested spontaneous reaction wave model, and thus does not require the formation of distributed flames. Critical turbulent flame speeds, predicted by this mechanism for the onset of DDT, are in agreement with DNS results.

We describe the implementation of the shearing box approximation for the study of the dynamics of accretion disks in the Athena magnetohydrodynamics (MHD) code. Second-order Crank-Nicholson time differencing is used for the Coriolis and tidal gravity source terms that appear in the momentum equation for accuracy and stability. We show this approach conserves energy for epicyclic oscillations in hydrodynamic flows to round-off error. In the energy equation, the tidal gravity source terms are differenced as the gradient of an effective potential in a way which guarantees that total energy (including the gravitational potential energy) is also conserved to round-off error. We introduce an orbital advection algorithm for MHD based on constrained transport to preserve the divergence-free constraint on the magnetic field. This algorithm removes the orbital velocity from the time step constraint, and makes the truncation error more uniform in radial position. Modifications to the shearing box boundary conditions applied at the radial boundaries are necessary to conserve the total vertical magnetic flux. In principle similar corrections are also required to conserve mass, momentum and energy, however in practice we find the orbital advection method conserves these quantities to better than 0.03% over hundreds of orbits. The algorithms have been applied to studies of the nonlinear regime of the MRI in very wide (up to 32 scale heights) horizontal domains.

A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher-order Godunov methods with the constrained transport (CT) technique to enforce the divergence-free constraint on the magnetic field. Discretization is based on cell-centered volume-averages for mass, momentum, and energy, and face-centered area-averages for the magnetic field. Novel features of the algorithm include (1) a consistent framework for computing the time- and edge-averaged electric fields used by CT to evolve the magnetic field from the time- and area-averaged Godunov fluxes, (2) the extension to MHD of spatial reconstruction schemes that involve a dimensionally-split time advance, and (3) the extension to MHD of two different dimensionally-unsplit integration methods. Implementation of the algorithm in both C and Fortran95 is detailed, including strategies for parallelization using domain decomposition. Results from a test suite which includes problems in one-, two-, and three-dimensions for both hydrodynamics and MHD are given, not only to demonstrate the fidelity of the algorithms, but also to enable comparisons to other methods. The source code is freely available for download on the web.

We present a single step, second-order accurate Godunov scheme for ideal MHD which is an extension of the method described by Gardiner & Stone (2005) to three dimensions. This algorithm combines the corner transport upwind (CTU) method of Colella for multidimensional integration, and the constrained transport (CT) algorithm for preserving the divergence-free constraint on the magnetic field. We describe the calculation of the PPM interface states for 3D ideal MHD which must include multidimensional ``MHD source terms'' and naturally respect the balance implicit in these terms by the ${\bf\nabla\cdot B}=0$ condition. We compare two different forms for the CTU integration algorithm which require either 6- or 12-solutions of the Riemann problem per cell per time-step, and present a detailed description of the 6-solve algorithm. Finally, we present solutions for test problems to demonstrate the accuracy and robustness of the algorithm.

We study the magnetic Rayleigh-Taylor instability in three dimensions, with focus on the nonlinear structure and evolution that results from different initial field configurations. We study strong fields in the sense that the critical wavelength l_c at which perturbations along the field are stable is a large fraction of the size of the computational domain. We consider magnetic fields which are initially parallel to the interface, but have a variety of configurations, including uniform everywhere, uniform in the light fluid only, and fields which change direction at the interface. Strong magnetic fields do not suppress instability, in fact by inhibiting secondary shear instabilities, they reduce mixing between the heavy and light fluid, and cause the rate of growth of bubbles and fingers to increase in comparison to hydrodynamics. Fields parallel to, but whose direction changes at, the interface produce long, isolated fingers separated by the critical wavelength l_c, which may be relevant to the morphology of the optical filaments in the Crab nebula.

We present a method for simulating the evolution of HII regions driven by point sources of ionizing radiation in magnetohydrodynamic media, implemented in the three-dimensional Athena MHD code. We compare simulations using our algorithm to analytic solutions and show that the method passes rigorous tests of accuracy and convergence. The tests reveal several conditions that an ionizing radiation-hydrodynamic code must satisfy to reproduce analytic solutions. As a demonstration of our new method, we present the first three-dimensional, global simulation of an HII region expanding into a magnetized gas. The simulation shows that magnetic fields suppress sweeping up of gas perpendicular to magnetic field lines, leading to small density contrasts and extremely weak shocks at the leading edge of the HII region's expanding shell.

We study the nonlinear evolution of the magnetic Rayleigh-Taylor instability using three-dimensional MHD simulations. We consider the idealized case of two inviscid, perfectly conducting fluids of constant density separated by a contact discontinuity perpendicular to the effective gravity g, with a uniform magnetic field B parallel to the interface. Modes parallel to the field with wavelengths smaller than l_c = [B B/(d_h - d_l) g] are suppressed (where d_h and d_l are the densities of the heavy and light fluids respectively), whereas modes perpendicular to B are unaffected. We study strong fields with l_c varying between 0.01 and 0.36 of the horizontal extent of the computational domain. Even a weak field produces tension forces on small scales that are significant enough to reduce shear (as measured by the distribution of the amplitude of vorticity), which in turn reduces the mixing between fluids, and increases the rate at which bubbles and finger are displaced from the interface compared to the purely hydrodynamic case. For strong fields, the highly anisotropic nature of unstable modes produces ropes and filaments. However, at late time flow along field lines produces large scale bubbles. The kinetic and magnetic energies transverse to gravity remain in rough equipartition and increase as t^4 at early times. The growth deviates from this form once the magnetic energy in the vertical field becomes larger than the energy in the initial field. We comment on the implications of our results to Z-pinch experiments, and a variety of astrophysical systems.

Using fully three-dimensional hydrodynamic simulations, we investigate the effect of the Coriolis force on the hydrodynamic and observable properties of colliding wind binary systems. To make the calculations tractable, we assume adiabatic, constant velocity winds. The neglect of radiative driving, gravitational deceleration, and cooling limit the application of our models to real systems. However, these assumptions allow us to isolate the effect of the Coriolis force, and by simplifying the calculations, allow us to use a higher resolution (up to 640^3) and to conduct a larger survey of parameter space. We study the dynamics of collidng winds with equal mass loss rates and velocities emanating from equal mass stars on circular orbits, with a range of values for the ratio of the wind to orbital velocity. We also study the dynamics of winds from stars on elliptical orbits and with unequal strength winds. Orbital motion of the stars sweeps the shocked wind gas into an Archimedean spiral, with asymmetric shock strengths and therefore unequal postshock temperatures and densities in the leading and trailing edges of the spiral. We observe the Kelvin-Helmholtz instability at the contact surface between the shocked winds in systems with orbital motion even when the winds are identical. The change in shock strengths caused by orbital motion increases the volume of X-ray emitting post-shock gas with T > 0.59 keV by 63% for a typical system as the ratio of wind velocity to orbital velocity decreases to V_w/V_o = 2.5. This causes increased free-free emission from systems with shorter orbital periods and an altered time-dependence of the wind attenuation. We comment on the importance of the effects of orbital motion on the observable properties of colliding wind binaries.

We carry out three-dimensional, high resolution (up to $1024^2\times 256$) hydrodynamic simulations of the evolution of vortices in vertically unstratified Keplerian disks using the shearing sheet approximation. The transient amplification of incompressible, linear amplitude leading waves (which has been proposed as a possible route to nonlinear hydrodynamical turbulence in disks) is used as one test of our algorithms; our methods accurately capture the predicted amplification, converges at second-order, and is free from aliasing. Waves expected to reach nonlinear amplitude at peak amplification become unstable to Kelvin-Helmholtz modes when $\mid W_{\rm max}\mid\gtrsim \Omega$ (where $W_{\rm max}$ is the local maximum of vorticity and $\Omega$ the angular velocity). We study the evolution of a power-law distribution of vorticity consistent with Kolmogorov turbulence; in two-dimensions long-lived vortices emerge and decay slowly, similar to previous studies. In three-dimensions, however, vortices are unstable to bending modes, leading to rapid decay. Only vortices with a length to width ratio smaller than one survive; in three-dimensions the residual kinetic energy and shear stress is at least one order of magnitude smaller than in two-dimensions. No evidence for sustained hydrodynamical turbulence and transport is observed in three-dimensions. Instead, at late times the residual transport is determined by the amplitude of slowly decaying, large-scale vortices (with horizontal extent comparable to the scale height of the disk), with additional contributions from nearly incompressible inertial waves possible. Evaluating the role that large-scale vortices play in astrophysical accretion disks will require understanding the mechanisms that generate and destroy them.

In these proceedings we present recent efforts to understand the energetics of magnetohydrodynamic (MHD) turbulence driven by the magneto-rotational instability (MRI). These studies are carried out in the local (shearing box) approximation using the Athena simulation code. Athena is a higher order Godunov algorithm based on the piecewise parabolic method (PPM), the corner transport upwind (CTU) integration algorithm, and the constrained transport (CT) algorithm for evolving the magnetic field. This algorithm is particularly suited for these studies owing to the conservation properties of a Godunov scheme and the particular implementation of the shearing box source terms used here. We present a variety of calculations which may be compared directly to previously published results and discuss them in some detail. The only significant discrepancy found between the results presented here and in the published literature involves the turbulent heating rate. We observe the presence of recurrent channel solutions in calculations involving a mean vertical magnetic field and the associated time lag between the energy injection and thermalization rate. We also present the results of a shearing box calculation which includes an optically thin radiative term with a cooling rate selected to match the turbulent heating rate. Some properties which we find uniformly present in all of the calculations presented here are compressible fluctuations, spiral waves and weak shocks. It is found that these compressible modes dominate the temporal fluctuations in the probability distribution functions for most of the thermodynamic variables; only the specific entropy is relatively immune to their effects.

We describe a new numerical scheme for MHD which combines a higher order Godunov method (PPM) with Constrained Transport. The results from a selection of multidimensional test problems are presented. The complete test suite used to validate the method, as well as implementations of the algorithm in both F90 and C, are available from the web. A fully three-dimensional version of the algorithm has been developed, and is being applied to a variety of astrophysical problems including the decay of supersonic MHD turbulence, the nonlinear evolution of the MHD Rayleigh-Taylor instability, and the saturation of the magnetorotational instability in the shearing box. Our new simulations of the MRI represent the first time that a higher-order Godunov scheme has been applied to this problem, providing a quantitative check on the accuracy of previous results computed with ZEUS; the latter are found to be reliable.

We describe a single step, second-order accurate Godunov scheme for ideal MHD based on combining the piecewise parabolic method (PPM) for performing spatial reconstruction, the corner transport upwind (CTU) method of Colella for multidimensional integration, and the constrained transport (CT) algorithm for preserving the divergence-free constraint on the magnetic field. We adopt the most compact form of CT, which requires the field be represented by area-averages at cell faces. We demonstrate that the fluxes of the area-averaged field used by CT can be made consistent with the fluxes of the volume-averaged field returned by a Riemann solver if they obey certain simple relationships. We use these relationships to derive new algorithms for constructing the CT fluxes at grid cell corners which reduce exactly to the equivalent one-dimensional solver for plane-parallel, grid-aligned flow. We show that the PPM reconstruction algorithm must include multidimensional terms for MHD, and we describe a number of important extensions that must be made to CTU in order for it to be used for MHD with CT. We present the results of a variety of test problems to demonstrate the method is accurate and robust.

We present two-dimensional, cylindrically symmetric simulations of hydrodynamic and magnetohydrodynamic (MHD) wide-angle winds interacting with a collapsing environment. These simulations have direct relevance to young stellar objects (YSOs). The results may also be of use in the study of collimated outflows from proto-planetary and planetary nebulae. We study a range of wind configurations consistent with asymptotic MHD wind collimation. The degree of collimation is parameterized by the ratio of the wind density at the pole to that of the equator. We find that a toroidal magnetic field can have a significant influence on the resulting outflow, giving rise to a very dense, jet-like flow in the post-shock region. The properties of the flow in this region are similar to the asymptotic state of a collimated MHD wind. We conclude that wide-angle MHD winds are quite likely capable of driving molecular outflows. Due to difficulty in treating MHD winds ab-initio in simulations we choose magnetic field strengths in the wind consistent slow magnetic rotators. While MHD launched winds will be in the fast rotator regime we discuss how our results, which rely on toroidal pinch effects, will hold for stronger field strengths.

We present first results of astrophysically relevant experiments where highly supersonic plasma jets are generated via conically convergent flows. The convergent flows are created by electrodynamic acceleration of plasma in a conical array of fine metallic wires (a modification of the wire array Z-pinch). Stagnation of plasma flow on the axis of symmetry forms a standing conical shock effectively collimating the flow in the axial direction. This scenario is essentially similar to that discussed by Canto\' ~and collaborators as a purely hydrodynamic mechanism for jet formation in astrophysical systems. Experiments using different materials (Al, Fe and W) show that a highly supersonic ($M\sim 20$), well-collimated jet is generated when the radiative cooling rate of the plasma is significant. We discuss scaling issues for the experiments and their potential use for numerical code verification. The experiments also may allow direct exploration of astrophysically relevant issues such as collimation, stability and jet-cloud interactions.

We discuss the propagation of radiative MHD jets and outflows focusing on outflows driven by magneto-centrifugal rotators. Our goal is to link the properties of the jets with the physics of the sources which produce them. We find that density and magnetic field stratification (with radius) in jets from magnetized rotators leads to new behavior including the development of a dense inner jet core and a low density collar. We also report on more general studies of ambipolar diffusion and field geometry in pulsed jets. Finally we describe a new work designed to study the effects of magnetized winds on circumstellar environments appropriate to YSOs and PNe.

Recent studies have focused on the the role of initially weak toroidal magnetic fields embedded in a stellar wind as the agent for collimation in planetary nebulae. In these models the wind is assumed to be permeated by a helical magnetic field in which the poloidal component falls off faster than the toroidal component. The collimation only occurs after the wind is shocked at large distances from the stellar source. In this paper we re-examine assumptions built into this ``Magnetized Wind Blown Bubble'' (MWBB) model. We show that a self-consistent study of the model leads to a large parameter regime where the wind is self-collimated before the shock wave is encountered. We also explore the relation between winds in the MWBB model and those which are produced via magneto-centrifugal processes. We conclude that a more detailed examination of the role of self-collimation is needed in the context of PNe studies.

Hypersonic, highly collimated, mass outflows ({\it jets}) are a ubiquitous phenomena in astrophysics. While the character of the jets differ, many exhibit some form of quasi-periodic clumping indicating the jet source is episodic or {\it pulsed}. The presence of pulsed jets in so many astrophysical contexts suggests a common formation mechanism. Such a process seems to have been found in {\it Magneto-centrifugal launching}, the combination of magnetic and centrifugal forces that occurs when a magnetized gaseous accretion disk orbits a central gravitating source. Observations of strong magnetic fields in jets are, however, rare or indirect. Thus the presence and effects of magnetic fields in YSO jets remains an unresolved issue of the highest importance. In this letter we focus on what should be expected of the structure of the fields in pulsed YSO jets. We show that combining velocity variability with an initial field configuration consistent with collimated, Magneto-centrifugally launched jets leads to a clear set of predictions concerning the geometry and relative strength of the magnetic field components in evolving YSO (and perhaps other) jets.

We present simulations of the propagation of magnetized jets. This work differs from previous studies in that the cross-sectional distributions of the jets's state variables are derived from analytical models for magneto-centrifugal launching. The source is a magnetized rotator whose properties are specfied as boundary conditions. The jets in these simulations are considerably more complex than the ``top-hat''constant density etc. profiles used in previous work. We find that density and magnetic field stratification (with radius) in the jet leads to new behavior including the separation of an inner jet core from a low density collar. We find this {\it jet within a jet} structure, along with the magnetic stresses, leads to propagation behaviors not observed in previous simulation studies. Our methodology allows us to compare MHD jets from different types of sources whose properties could ultimately be derived from the behavior of the propagating jets.

We address the issue of ambipolar diffusion in Herbig-Haro jets. The current consensus holds that these jets are launched and collimated via MHD forces. Observations have, however, shown that the jets can be mildly to weakly ionized. Beginning with a simple model for cylindrical equilibrium between neutral, plasma and magnetic pressures we calculate the characteristic time-scale for ambipolar diffusion. Our results show that a significant fraction of HH jets will have ambipolar diffusion time-scales equivalent to, or less than the dynamical time-scales. This implies that MHD equilibria established at the base of a HH jet may not be maintained as the jet propagates far from its source. For typical jet parameters one finds that the length scale where ambipolar diffusion should become significant corresponds to the typical size of large (parsec) scale jets. We discuss the significance of these results for the issue of magnetic fields in parsec-scale jets.