The nature of type Ia supernovae (SNIa) - thermonuclear explosions of white dwarf stars - is an open question in astrophysics. Virtually all existing theoretical models of normal, bright SNIa require the explosion to produce a detonation in order to consume all of stellar material, but the mechanism for the deflagration-to-detonation transition (DDT) remains unclear. We present a unified theory of turbulence-induced DDT that describes the mechanism and conditions for initiating detonation both in unconfined chemical and thermonuclear explosions. The model is validated using experiments with chemical flames and numerical simulations of thermonuclear flames. We use the developed theory to determine criteria for detonation initiation in the single degenerate Chandrasekhar-mass SNIa model, and show that DDT is almost inevitable at densities 10^7 - 10^8 g/cc.

(Abridged) A series of three-dimensional numerical simulations is used to study the intrinsic stability of high-speed turbulent flames. Calculations model the interaction of a fully-resolved premixed flame with a highly subsonic, statistically steady, homogeneous, isotropic turbulence. We consider a wide range of turbulent intensities and system sizes, corresponding to the Damk\"ohler numbers Da = 0.1-6.0. These calculations show that turbulent flames in the regimes considered are intrinsically unstable. In particular, we find three effects. 1) Turbulent flame speed develops pulsations with the observed peak-to-peak amplitude > 10 and a characteristic time scale close to a large-scale eddy turnover time. Such variability is caused by the interplay between turbulence, which continuously creates the flame surface, and highly intermittent flame collisions, which consume the flame surface. 2) Unstable burning results in the periodic pressure build-up and the formation of pressure waves or shocks, when the flame speed approaches or exceeds the speed of a Chapman-Jouguet deflagration. 3) Coupling of pressure gradients formed during pulsations with density gradients across the flame leads to the anisotropic amplification of turbulence inside the flame volume and flame acceleration. Such process, which is driven by the baroclinic term in the vorticity transport equation, is a reacting-flow analog of the mechanism underlying the Richtmyer-Meshkov instability. With the increase in turbulent intensity, the limit-cycle instability discussed here transitions to the regime described in our previous work, in which the growth of the flame speed becomes unbounded and produces a detonation.

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 study the dynamics and properties of a turbulent flame, formed in the presence of subsonic, high-speed, homogeneous, isotropic Kolmogorov-type turbulence in an unconfined system. Direct numerical simulations are performed with Athena-RFX, a massively parallel, fully compressible, high-order, dimensionally unsplit, reactive-flow code. A simplified reaction-diffusion model represents a stoichiometric H2-air mixture. The system being modeled represents turbulent combustion with the Damkohler number Da = 0.05 and with the turbulent velocity at the energy injection scale 30 times larger than the laminar flame speed. The simulations show that flame interaction with high-speed turbulence forms a steadily propagating turbulent flame with a flame brush width approximately twice the energy injection scale and a speed four times the laminar flame speed. A method for reconstructing the internal flame structure is described and used to show that the turbulent flame consists of tightly folded flamelets. The reaction zone structure of these is virtually identical to that of the planar laminar flame, while the preheat zone is broadened by approximately a factor of two. Consequently, the system evolution represents turbulent combustion in the thin-reaction zone regime. The turbulent cascade fails to penetrate the internal flame structure, and thus the action of small-scale turbulence is suppressed throughout most of the flame. Finally, our results suggest that for stoichiometric H2-air mixtures, any substantial flame broadening by the action of turbulence cannot be expected in all subsonic regimes.

(Abridged) Direct numerical simulations of the interaction of a premixed flame with driven, subsonic, homogeneous, isotropic, Kolmogorov-type turbulence in an unconfined system are used to study the mechanisms determining the turbulent flame speed, S_T, in the thin reaction zone regime. High intensity turbulence is considered with the r.m.s. velocity 35 times the laminar flame speed, S_L, resulting in the Damkohler number Da = 0.05. Here we show that: (1) The flame brush has a complex internal structure, in which the isosurfaces of higher fuel mass fractions are folded on progressively smaller scales. (2) Global properties of the turbulent flame are best represented by the structure of the region of peak reaction rate, which defines the flame surface. (3) In the thin reaction zone regime, S_T is predominantly determined by the increase of the flame surface area, A_T, caused by turbulence. (4) The observed increase of S_T relative to S_L exceeds the corresponding increase of A_T relative to the surface area of the planar laminar flame, on average, by ~14%, varying from only a few percent to ~30%. (5) This exaggerated response is the result of tight flame packing by turbulence, which causes frequent flame collisions and formation of regions of high flame curvature, or "cusps." (6) The local flame speed in the cusps substantially exceeds its laminar value, which results in a disproportionately large contribution of cusps to S_T compared with the flame surface area in them. (7) A criterion is established for transition to the regime significantly influenced by cusp formation. In particular, at Karlovitz numbers Ka > 20, flame collisions provide an important mechanism controlling S_T, in addition to the increase of A_T by large-scale motions and the potential enhancement of diffusive transport by small-scale turbulence.

We present the method for computation of fluid flows that are characterized by the large degree of expansion/contraction and in which the fluid velocity is dominated by the bulk component associated with the expansion/contraction and/or rotation of the flow. We consider the formulation of Euler equations of fluid dynamics in a homologously expanding/contracting and/or rotating reference frame. The frame motion is adjusted to minimize local fluid velocities. Such approach allows to accommodate very efficiently large degrees of change in the flow extent. Moreover, by excluding the contribution of the bulk flow to the total energy the method eliminates the high Mach number problem in the flows of interest. An important practical advantage of the method is that it can be easily implemented with virtually any implicit or explicit Eulerian hydrodynamic scheme and adaptive mesh refinement (AMR) strategy. We also consider in detail equation invariance and existence of conservative formulation of equations for special classes of expanding/contracting reference frames. Special emphasis is placed on extensive numerical testing of the method for a variety of reference frame motions, which are representative of the realistic applications of the method. We study accuracy, conservativity, and convergence properties of the method both in problems which are not its optimal applications as well as in systems in which the use of this method is maximally beneficial. Such detailed investigation of the numerical solution behavior is used to define the requirements that need to be considered in devising problem-specific fluid motion feedback mechanisms.

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.

We present a design for high energy density laboratory experiments studying the interaction of hypersonic shocks with a large number of inhomogeneities. These ``clumpy'' flows are relevant to a wide variety of astrophysical environments including the evolution of molecular clouds, outflows from young stars, Planetary Nebulae and Active Galactic Nuclei. The experiment consists of a strong shock (driven by a pulsed power machine or a high intensity laser) impinging on a region of randomly placed plastic rods. We discuss the goals of the specific design and how they are met by specific choices of target components. An adaptive mesh refinement hydrodynamic code is used to analyze the design and establish a predictive baseline for the experiments. The simulations confirm the effectiveness of the design in terms of articulating the differences between shocks propagating through smooth and clumpy environments. In particular, we find significant differences between the shock propagation speeds in a clumpy medium compared to a smooth one with the same average density. The simulation results are of general interest for foams in both inertial confinement fusion and laboratory astrophysics studies. Our results highlight the danger of using average properties of inhomogeneous astrophysical environments when comparing timescales for critical processes such as shock crossing and gravitational collapse times.

We consider the stability of an accretion disk wind to cloud formation when subject to a central radiation force. For a vertical launch velocity profile that is Keplerian or flatter and the presence of a significant radiation pressure, the wind flow streamlines cross in a conical layer. We argue that such regions are highly unstable, and are natural sites for supersonic turbulence and, consequently, density compressions. We suggest that combined with thermal instability these will all conspire to produce clouds. Such clouds can exist in dynamical equilibrium, constantly dissipating and reforming. As long as there is an inner truncation radius to the wind, our model emerges with a biconical structure similar to that inferred by Elvis (2000) for the broad line region (BLR) of active galactic nuclei (AGN). Our results may also apply to other disk-wind systems.

Many astrophysical flows occur in inhomogeneous (clumpy) media. We present results of a numerical study of steady, planar shocks interacting with a system of embedded cylindrical clouds. Our study uses a two-dimensional geometry. Our numerical code uses an adaptive mesh refinement allowing us to achieve sufficiently high resolution both at the largest and the smallest scales. We neglect any radiative losses, heat conduction, and gravitational forces. Detailed analysis of the simulations shows that interaction of embedded inhomogeneities with the shock/postshock wind depends primarily on the thickness of the cloud layer and arrangement of the clouds in the layer. The total cloud mass and the total number of individual clouds is not a significant factor. We define two classes of cloud distributions: thin and thick layers. We define the critical cloud separation along the direction of the flow and perpendicular to it distinguishing between the interacting and noninteracting regimes of cloud evolution. Finally we discuss mass-loading and mixing in such systems.

Many astrophysical flows occur in inhomogeneous media. The broad-line regions (BLR) of active galactic nuclei (AGNs) are one of the important examples where emission-line clouds interact with the outflow. We present results of a numerical study of the interaction of a steady, planar shock / supersonic postshock flow with a system of embedded cylindrical clouds in a two-dimensional geometry. Detailed analysis shows that the interaction of embedded inhomogeneities with the shock / postshock wind depends primarily on the thickness of the cloud layer and the arrangement of the clouds in the layer, as opposed to the total cloud mass and the total number of individual clouds. This allows us to define two classes of cloud distributions: thin and thick layers. We define the critical cloud separation along the direction of the flow and perpendicular to it. This definition allows us to distinguish between the interacting and noninteracting regimes of cloud evolution. Finally we discuss mass-loading in such systems.