The Thermonuclear Explosion Of Chandrasekhar Mass White Dwarfs

The flame born in the deep interior of a white dwarf that becomes a Type Ia supernova is subject to several instabilities. We briefly review these instabilities and the corresponding flame acceleration. We discuss the conditions necessary for each of the currently proposed explosion mechanisms and the attendant uncertainties. A grid of critical masses for detonation in the range $10^7$ - $2 \times 10^9$ g cm$^{-3}$ is calculated and its sensitivity to composition explored. Prompt detonations are physically improbable and appear unlikely on observational grounds. Simple deflagrations require some means of boosting the flame speed beyond what currently exists in the literature. ``Active turbulent combustion'' and multi-point ignition are presented as two plausible ways of doing this. A deflagration that moves at the ``Sharp-Wheeler'' speed, $0.1 g_{\rm eff} t$, is calculated in one dimension and shows that a healthy explosion is possible in a simple deflagration if the front moves with the speed of the fastest floating bubbles. The relevance of the transition to the ``distributed burning regime'' is discussed for delayed detonations. No model emerges without difficulties, but detonation in the distributed regime is plausible, will produce intermediate mass elements, and warrants further study.Comment: 28 pages, 4 figures included, uses aaspp4.sty. Submitted to Ap

Emergence of pulled fronts in fermionic microscopic particle models

We study the emergence and dynamics of pulled fronts described by the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation in the microscopic reaction-diffusion process A + A A$ on the lattice when only a particle is allowed per site. To this end we identify the parameter that controls the strength of internal fluctuations in this model, namely, the number of particles per correlated volume. When internal fluctuations are suppressed, we explictly see the matching between the deterministic FKPP description and the microscopic particle model.Comment: 4 pages, 4 figures. Accepted for publication in Phys. Rev. E as a Rapid Communicatio

Conditional-moment Closure with Differential Diffusion for Soot Evolution in Fire

The conditional-moment closure (CMC) equation for the evolution of a large Lewis number scalar, soot, is derived starting from the joint probability density function (pdf) equation for the gas-phase mixture fraction, ξ g , and the soot mass fraction, Y s . Unlike previous approaches starting with the joint pdf, the residual terms that result from the typical closure models were retained. A new formulation of the one-dimensional turbulence (ODT) model suitable for spatially evolving flows with buoyant acceleration and radiative transport in participating media was employed to carry out simulations of a prototypical ethene fire. The resulting ODT evolution of ξ g and Y s was used to assess the significance of various terms in the CMC equation including the residual correlations. The terms involving differential diffusion are found to be important along with the soot source terms and the large-scale evolution of both ξ g and Y s . Of particular importance in the regions in mixture-fraction space around the soot production and consumption is a residual term, not previously identified, related to the correlation between the differential diffusion and Y s . This term results in a diffusion-like behavior of Y s in the mixture fraction coordinate that has an apparent Lewis number near unity. In scenarios where the large Lewis number component is a non-negligible component of the mixture fraction (i.e., large soot loading), it is found easier to employ a mixture fraction neglecting this component. Such a mixture-fraction variable has a chemical source term, but this appears easier to model than the differential diffusion and dissipation terms that result when the large Lewis number component is retained in the mixture-fraction definition

Shift in the velocity of a front due to a cut-off

We consider the effect of a small cut-off epsilon on the velocity of a traveling wave in one dimension. Simulations done over more than ten orders of magnitude as well as a simple theoretical argument indicate that the effect of the cut-off epsilon is to select a single velocity which converges when epsilon tends to 0 to the one predicted by the marginal stability argument. For small epsilon, the shift in velocity has the form K(log epsilon)^(-2) and our prediction for the constant K agrees very well with the results of our simulations. A very similar logarithmic shift appears in more complicated situations, in particular in finite size effects of some microscopic stochastic systems. Our theoretical approach can also be extended to give a simple way of deriving the shift in position due to initial conditions in the Fisher-Kolmogorov or similar equations.Comment: 12 pages, 3 figure

Fronts in randomly advected and heterogeneous media and nonuniversality of Burgers turbulence: Theory and numerics

A recently established mathematical equivalence--between weakly perturbed Huygens fronts (e.g., flames in weak turbulence or geometrical-optics wave fronts in slightly nonuniform media) and the inviscid limit of white-noise-driven Burgers turbulence--motivates theoretical and numerical estimates of Burgers-turbulence properties for specific types of white-in-time forcing. Existing mathematical relations between Burgers turbulence and the statistical mechanics of directed polymers, allowing use of the replica method, are exploited to obtain systematic upper bounds on the Burgers energy density, corresponding to the ground-state binding energy of the directed polymer and the speedup of the Huygens front. The results are complementary to previous studies of both Burgers turbulence and directed polymers, which have focused on universal scaling properties instead of forcing-dependent parameters. The upper-bound formula can be heuristically understood in terms of renormalization of a different kind from that previously used in combustion models, and also shows that the burning velocity of an idealized turbulent flame does not diverge with increasing Reynolds number at fixed turbulence intensity, a conclusion that applies even to strong turbulence. Numerical simulations of the one-dimensional inviscid Burgers equation using a Lagrangian finite-element method confirm that the theoretical upper bounds are sharp within about 15% for various forcing spectra (corresponding to various two-dimensional random media). These computations provide a new quantitative test of the replica method. The inferred nonuniversality (spectrum dependence) of the front speedup is of direct importance for combustion modeling.Comment: 20 pages, 2 figures, REVTeX 4. Moved some details to appendices, added figure on numerical metho

Can Deflagration-Detonation-Transitions occur in Type Ia Supernovae?

The mechanism for deflagration-detonation-transition (DDT) by turbulent preconditioning, suggested to explain the possible occurrence of delayed detonations in Type Ia supernova explosions, is argued to be conceptually inconsistent. It relies crucially on diffusive heat losses of the burned material on macroscopic scales. Regardless of the amplitude of turbulent velocity fluctuations, the typical gradient scale for temperature fluctuations is shown to be the laminar flame width or smaller, rather than the factor of thousand more required for a DDT. Furthermore, thermonuclear flames cannot be fully quenched in regions much larger than the laminar flame width as a consequence of their simple ``chemistry''. Possible alternative explosion scenarios are briefly discussed.Comment: 8 pages, uses aastex; added references. Accepted by ApJ Letter

On the small-scale stability of thermonuclear flames in Type Ia supernovae

We present a numerical model which allows us to investigate thermonuclear flames in Type Ia supernova explosions. The model is based on a finite-volume explicit hydrodynamics solver employing PPM. Using the level-set technique combined with in-cell reconstruction and flux-splitting schemes we are able to describe the flame in the discontinuity approximation. We apply our implementation to flame propagation in Chandrasekhar-mass Type Ia supernova models. In particular we concentrate on intermediate scales between the flame width and the Gibson-scale, where the burning front is subject to the Landau-Darrieus instability. We are able to reproduce the theoretical prediction on the growth rates of perturbations in the linear regime and observe the stabilization of the flame in a cellular shape. The increase of the mean burning velocity due to the enlarged flame surface is measured. Results of our simulation are in agreement with semianalytical studies.Comment: 9 pages, 7 figures, Uses AASTEX, emulateapj5.sty, onecolfloat.sty. Replaced with accepted version (ApJ), Figures 1 and 3 are ne