Deflagration Model Research Articles

Three-dimensional pure deflagration models with nucleosynthesis and synthetic observables for Type Ia supernovae

We investigate whether pure deflagration models of Chandrasekhar-mass carbon-oxygen white dwarf stars can account for one or more subclass of the observed population of Type Ia supernova (SN Ia) explosions. We compute a set of 3D full-star hydrodynamic explosion models, in which the deflagration strength is parametrized using the multispot ignition approach. For each model, we calculate detailed nucleosynthesis yields in a post-processing step with a 384 nuclide nuclear network. We also compute synthetic observables with our 3D Monte Carlo radiative transfer code for comparison with observations. For weak and intermediate deflagration strengths (energy release E_nuc <~ 1.1 x 10^51 erg), we find that the explosion leaves behind a bound remnant enriched with 3 to 10 per cent (by mass) of deflagration ashes. However, we do not obtain the large kick velocities recently reported in the literature. We find that weak deflagrations with E_nuc ~ 0.5 x 10^51 erg fit well both the light curves and spectra of 2002cx-like SNe Ia, and models with even lower explosion energies could explain some of the fainter members of this subclass. By comparing our synthetic observables with the properties of SNe Ia, we can exclude the brightest, most vigorously ignited models as candidates for any observed class of SN Ia: their B - V colours deviate significantly from both normal and 2002cx-like SNe Ia and they are too bright to be candidates for other subclasses.

Modelling the interaction of thermonuclear supernova remnants with circumstellar structures: the case of Tycho's supernova remnant

The well-established Type Ia remnant of Tycho's supernova (SN 1572) reveals discrepant ambient medium density estimates based on either the measured dynamics or on the X-ray emission properties. This discrepancy can potentially be solved by assuming that the supernova remnant (SNR) shock initially moved through a stellar wind bubble, but is currently evolving in the uniform interstellar medium with a relatively low density. We investigate this scenario by combining hydrodynamical simulations of the wind-loss phase and the supernova remnant evolution with a coupled X-ray emission model, which includes non-equilibrium ionization. For the explosion models we use the well-known W7 deflagration model and the delayed detonation model that was previously shown to provide good fits to the X-ray emission of Tycho's SNR. Our simulations confirm that a uniform ambient density cannot simultaneously reproduce the dynamical and X-ray emission properties of Tycho. In contrast, models that considered that the remnant was evolving in a dense, but small, wind bubble reproduce reasonably well both the measured X-ray emission spectrum and the expansion parameter of Tycho's SNR. Finally, we discuss possible mass loss scenarios in the context of single- and double-degenerate models which possible could form such a small dense wind bubble.

An Eulerian–Lagrangian computational model for deflagration and detonation of high explosives

A solid phase explosives deflagration and detonation model capable of surface burning, convective bulk burning and detonation is formulated in the context of Eulerian–Lagrangian material mechanics. Well-validated combustion and detonation models, WSB and JWL++, are combined with two simple, experimentally indicated transition thresholds partitioning the three reaction regimes. Standard experiments are simulated, including the Aluminum Flyer Plate test, the Cylinder test, the rate stick test and the Steven test in order to validate the model. Cell and particle resolution dependence of simulation metrics are presented and global uncertainties assigned. Error quantification comparisons with experiments led to values generally below 7% (1σ). Finally, gas flow through porous media is implicated as the driving force behind the deflagration to detonation transition.

TYPE Ia SUPERNOVAE AS SITES OF THEp-PROCESS: TWO-DIMENSIONAL MODELS COUPLED TO NUCLEOSYNTHESIS

We explore SNe Ia as p-process sites in the framework of two-dimensional SN Ia delayed detonation and pure deflagration models. The WD precursor is assumed to have reached the Chandrasekhar mass in a binary system by mass accretion from a giant/main sequence companion. We use enhanced s-seed distributions, obtained from a sequence of thermal pulse instabilities both in the AGB phase and in the accreted material. We apply the tracer-particle method to reconstruct the nucleosynthesis by the thermal histories of Lagrangian particles, passively advected in the hydrodynamic calculations. For each particle we follow the explosive nucleosynthesis with a detailed network for all isotopes up to 209Bi. We find that SNe Ia can produce a large amount of p-nuclei, both the light p-nuclei below A=120 and the heavy-p nuclei, at quite flat average production factors, tightly related to the s-process seed distribution. For the first time, we find a stellar source able to produce both, light and heavy p-nuclei almost at the same level as 56Fe, including the very debated neutron magic 92,94Mo and 96,98Ru. We also find that there is an important contribution from p-process nucleosynthesis to the s-only nuclei 80Kr, 86Sr, to the neutron magic 90Zr, and to the neutron-rich 96Zr. Finally, we investigate the metallicity effect on p-process. Starting with different s-process seed distributions, for two metallicities Z = 0.02 and Z = 0.001, running SNe Ia models with different initial composition, we estimate that SNe Ia can contribute to, at least, 50% of the solar p-process composition.

Abundance stratification in Type Ia supernovae - III. The normal SN 2003du

The element abundance distributions in the ejecta of Type Ia supernova (SN) is studied by modelling a time series of optical spectra of SN 2003du until ~1 year after the explosion. Since SN 2003du is a very normal Type Ia SN both photometrically and spectroscopically, the abundance distribution derived for it can be considered as representative of normal Type Ia SNe. We find that the innermost layers are dominated by stable Fe-group elements, with a total mass of ~ 0.2 Msun, which are synthesized through electron capture. Above the core of stable elements there are thick 56Ni-rich layers. The total mass of 56Ni is 0.65 Msun. The Si- and S-rich layers are located above the 56Ni-rich layers. The dominant element in the outermost layers (M_r > 1.1 Msun, v > 13000 km/s) is O, with a small amount of Si. Little unburned C remains, with an upper limit of 0.016 Msun. The element distributions in the ejecta are moderately mixed, but not fully mixed as seen in three-dimensional deflagration models.

A versatile approach to explosion vent design that considers the influence of flame acceleration from congestion

AbstractCalculations for explosion venting in flammable vapor systems are often done using established methods documented in NFPA 68. It is understood that these equations may under predict vent area where enclosures are congested or the environment in the enclosure is initially turbulent. For new designs, it is helpful for design teams to have the flexibility to choose from a number of potential venting configurations. This article reviews a methodology that was developed and implemented for several projects where computational fluid dynamics was used to evaluate enclosure design parameters for a range of vent configurations. The results of the modeling simulations were then fitted to a lumped parameter model for deflagration venting developed by Molkov et al. The end result is an equation that can be used by the design team as needed to finalize the design minimizing the use of computer modeling resources. © 2010 American Institute of Chemical Engineers Process Saf Prog, 2010

NUCLEOSYNTHESIS IN TWO-DIMENSIONAL DELAYED DETONATION MODELS OF TYPE Ia SUPERNOVA EXPLOSIONS

The nucleosynthetic characteristics of various explosion mechanisms of Type Ia supernovae (SNe Ia) is explored based on three two-dimensional explosion simulations representing extreme cases: a pure turbulent deflagration, a delayed detonation following an approximately spherical ignition of the initial deflagration, and a delayed detonation arising from a highly asymmetric deflagration ignition. Apart from this initial condition, the deflagration stage is treated in a parameter-free approach. The detonation is initiated when the turbulent burning enters the distributed burning regime. This occurs at densities around $10^{7}$ g cm$^{-3}$ -- relatively low as compared to existing nucleosynthesis studies for one-dimensional spherically symmetric models. The burning in these multidimensional models is different from that in one-dimensional simulations as the detonation wave propagates both into unburned material in the high density region near the center of a white dwarf and into the low density region near the surface. Thus, the resulting yield is a mixture of different explosive burning products, from carbon-burning products at low densities to complete silicon-burning products at the highest densities, as well as electron-capture products synthesized at the deflagration stage. In contrast to the deflagration model, the delayed detonations produce a characteristic layered structure and the yields largely satisfy constraints from Galactic chemical evolution. In the asymmetric delayed detonation model, the region filled with electron capture species (e.g., $^{58}$Ni, $^{54}$Fe) is within a shell, showing a large off-set, above the bulk of $^{56}$Ni distribution, while species produced by the detonation are distributed more spherically (abridged).

TURBULENCE IN A THREE-DIMENSIONAL DEFLAGRATION MODEL FOR TYPE Ia SUPERNOVAE. II. INTERMITTENCY AND THE DEFLAGRATION-TO-DETONATION TRANSITION PROBABILITY

The delayed detonation model describes the observational properties of the majority of type Ia supernovae very well. Using numerical data from a three-dimensional deflagration model for type Ia supernovae, the intermittency of the turbulent velocity field and its implications on the probability of a deflagration-to-detonation (DDT) transition are investigated. From structure functions of the turbulent velocity fluctuations, we determine intermittency parameters based on the log-normal and the log-Poisson models. On the other hand, the analysis of the turbulent velocity fluctuations in the vicinity of the flame front by Roepke suggests a much higher probability of large velocity fluctuations on the grid scale in comparison to the log-normal intermittency model. Following Pan et al., we computed probability density functions for a DDT for the different distributions. Assuming that a DDT can occur in the stirred flame regime, as proposed by Woosley et al., the log-normal model would imply a delayed detonation between 0.7 and 0.8 seconds after the beginning of the deflagration phase for the multi-spot ignition scenario used in the simulation. However, the probability drops to virtually zero if a DDT is further constrained by the requirement that the turbulent velocity fluctuations reach about 500 km/s. Under this condition, delayed detonations are only possible if the distribution of the velocity fluctuations is not log-normal. From our calculations follows that the distribution obtained by Roepke allow for multiple DDTs around 0.8 seconds after ignition at a transition density close to 1x10^7 g/cm^3.

TURBULENCE IN A THREE-DIMENSIONAL DEFLAGRATION MODEL FOR TYPE Ia SUPERNOVAE. I. SCALING PROPERTIES

We analyze the statistical properties of the turbulent velocity field in the deflagration model for Type Ia supernovae. In particular, we consider the question of whether turbulence is isotropic and consistent with the Kolmogorov theory at small length scales. Using numerical data from a high-resolution simulation of a thermonuclear supernova explosion, spectra of the turbulence energy and velocity structure functions are computed. We show that the turbulent velocity field is isotropic at small length scales and follows a scaling law that is consistent with the Kolmogorov theory until most of the nuclear fuel is burned. At length scales greater than a certain characteristic scale, turbulence becomes anisotropic. Here, the radial velocity fluctuations follow the scaling law of the Rayleigh-Taylor instability, whereas the angular component still obeys Kolmogorov scaling. In the late phase of the explosion, this characteristic scale drops below the numerical resolution of the simulation. The analysis confirms that a subgrid-scale model for the unresolved turbulence energy is required for the consistent calculation of the flame speed in deflagration models of Type Ia supernovae, and that the assumption of isotropy on these scales is appropriate.

Thermonuclear supernovae: a multi-scale astrophysical problem challenging numerical simulations and visualization

The numerical modeling of type Ia supernovae is a demanding astrophysical task. Relevant physical processes take place on vastly different length- and timescales. This multi-scale character of the object poses challenges to the numerical approaches. We discuss an implementation that accounts for these problems by employing a large eddy simulation (LES) strategy for treating turbulence effects and a level-set technique to represent the thin thermonuclear flames. It is demonstrated that this approach works efficiently in simulations of the deflagration model and the delayed detonation model of type Ia supernovae. The resulting data reflect the multi-scale nature of the problem. Therefore, visualization has to be tackled with special techniques. We describe an approach that enables the interactive exploration of large datasets on commodity hardware. To this end, out-of-core methods are employed and the rendering of the data is achieved by a hybrid particle-based and texture-based volume-rendering technique.

A three-dimensional picture of the delayed-detonation model of type Ia supernovae

Deflagration models poorly explain the observed diversity of SNIa. Current multidimensional simulations of SNIa predict a significant amount of, so far unobserved, carbon and oxygen moving at low velocities. It has been proposed that these drawbacks can be resolved if there is a sudden jump to a detonation (delayed detonation), but this kind of models has been explored mainly in one dimension. Here we present new three-dimensional delayed detonation models in which the deflagraton-to-detonation transition (DDT) takes place in conditions like those favored by one-dimensional models. We have used a SPH code adapted to SNIa with algorithms devised to handle subsonic as well as supersonic combustion fronts. The starting point was a C-O white dwarf of 1.38 solar masses. When the average density on the flame surface reached 2-3x10^7 g/cm^3 a detonation was launched. The detonation wave processed more than 0.3 solar masses of carbon and oxygen, emptying the central regions of the ejecta of unburned fuel and raising its kinetic energy close to the fiducial 10^51 ergs expected from a healthy Type Ia supernova. The final amount of 56Ni synthesized also was in the correct range. However, the mass of carbon and oxygen ejected is still too high. The three-dimensional delayed detonation models explored here show an improvement over pure deflagration models, but they still fail to coincide with basic observational constraints. However, there are many aspects of the model that are still poorly known (geometry of flame ignition, mechanism of DDT, properties of detonation waves traversing a mixture of fuel and ashes). Therefore, it will be worth pursuing its exploration to see if a good SNIa model based on the three-dimensional delayed detonation scenario can be obtained.

A Three‐Dimensional Deflagration Model for Type Ia Supernovae Compared with Observations

A simulation of the thermonuclear explosion of a Chandrasekhar-mass C+O white dwarf, the most popular scenario of a Type Ia supernova (SN Ia), is presented. The underlying modeling is pursued in a self-consistent way, treating the combustion wave as a turbulent deflagration using well tested methods developed for laboratory combustion and based on the concept of large-eddy (LESs). Such consistency requires to capture the onset of the turbulent cascade on resolved scales. This is achieved by computing the dynamical evolution on a 10243 moving grid, which resulted in the best-resolved three-dimensional SN Ia simulation carried out thus far, reaching the limits of what can be done on present supercomputers. Consequently, the model has no free parameters other than the initial conditions at the onset of the explosion, and therefore it has considerable predictive power. Our main objective is to determine to which extent such a simulation can account for the observations of normal SNe Ia. Guided by previous simulations with less resolution and a less sophisticated flame model, initial conditions were chosen that yield a reasonably strong explosion and a sufficient amount of radioactive nickel for a bright display. We show that observables are indeed matched to a reasonable degree. In particular, good agreement is found with the light curves of normal SNe Ia. Moreover, the model reproduces the general features of the abundance stratification as inferred from the analysis of spectra. This indicates that it captures the main features of the explosion mechanism of SNe Ia. However, we also show that even a seemingly best-choice pure deflagration model has shortcomings that indicate the need for a different mode of nuclear burning at late times, perhaps the transition to a detonation at low density.

Signatures of Delayed Detonation, Asymmetry, and Electron Capture in the Mid‐Infrared Spectra of Supernovae 2003hv and 2005df

We present the first mid-infrared (5.2-15.2 micron) spectra of Type Ia supernovae (SNe 2003hv and 2005df). The detection of Ni emission in SN 2005df 135 days after the explosion provides direct observational evidence of high-density nuclear burning forming a significant amount of stable Ni in a Type Ia supernova. The observed emission line profiles in the SN 2005df spectrum indicate a chemically stratified ejecta structure. The SN 2005df Ar lines also exhibit a two-pronged emission profile implying that the Ar emission deviates significantly from spherical symmetry. The spectrum of SN 2003hv also shows signs of asymmetry, exhibiting blueshifted [Co III] which matches the blueshift of [Fe II] lines in nearly coeval NIR spectra. Finally, local thermodynamic equilibrium abundance estimates for the yield of radioactive 56Ni give ~0.5 Msun, for SN 2003hv, but only ~0.13-0.22 Msun for the apparently subluminous SN 2005df, supporting the notion that the luminosity of SNe Ia is primarily a function of the radioactive 56Ni yield. The chemically stratified ejecta structure observed in SN 2005df matches the predictions of delayed-detonation (DD) models, but is entirely incompatible with current three-dimensional deflagration models. Furthermore the degree that this layering persists to the innermost regions of the supernova is difficult to explain even in a DD scenario, where the innermost ejecta are still the product of deflagration burning. Thus, while these results are roughly consistent with a delayed detonation, it is clear that a key piece of physics is still missing from our understanding of the earliest phases of SN Ia explosions.

SN 2005hj: Evidence for Two Classes of Normal‐Bright SNe Ia and Implications for Cosmology

HET optical spectra covering the evolution from about 6 days before to about 5 weeks after maximum light and the ROTSE-IIIb unfiltered light curve of the Type Ia Supernova SN 2005hj are presented. The host galaxy shows H II region lines at redshift of z = 0.0574, which puts the peak unfiltered absolute magnitude at a somewhat overluminous -19.6. The spectra show weak and narrow Si II lines. These profiles do not change in width or depth for 10 days beginning around maximum light, indicating a constant expansion velocity of ≈10,600 km s-1. We analyzed the observations based on detailed radiation dynamical models. Delayed detonation and deflagration models that have been used to explain the majority of SNe Ia do not predict a long velocity plateau in the Si II minimum with an unvarying line profile. Pulsating delayed detonations and merger scenarios form shell-like density structures with properties mostly related to the mass of the shell, Mshell. We discuss how these models may explain the observed Si II line evolution; however, these models are based on spherical calculations and other possibilities may exist. SN 2005hj is consistent with respect to the onset, duration, and velocity of the plateau, the peak luminosity and, within the uncertainties, with the intrinsic colors for models with Mshell = 0.2 M☉. Our analysis suggests a distinct class of events hidden within the Branch-normal SNe Ia. The predicted relations between observables are confirmed, they may provide a way to separate these two groups. We discuss the implications for cosmological studies.

Constraining Deflagration Models of Type Ia Supernovae through Intermediate‐Mass Elements

The physical structure of a nuclear flame is a basic ingredient of the theory of Type Ia supernovae (SNIa). Assuming an exponential density reduction with several characteristic times we have followed the evolution of a planar nuclear flame in an expanding background from an initial density 6.6 10^7 g/cm3 down to 2 10^6 g/cm3. The total amount of synthesized intermediate-mass elements (IME), from silicon to calcium, was monitored during the calculation. We have made use of the computed mass fractions, X_IME, of these elements to give an estimation of the total amount of IME synthesized during the deflagration of a massive white dwarf. Using X_IME and adopting the usual hypothesis that turbulence decouples the effective burning velocity from the laminar flame speed, so that the relevant flame speed is actually the turbulent speed on the integral length-scale, we have built a simple geometrical approach to model the region where IME are thought to be produced. It turns out that a healthy production of IME involves the combination of not too short expansion times, t_c > 0.2 s, and high turbulent intensities. According to our results it could be difficult to produce much more than 0.2 solar masses of intermediate-mass elements within the deflagrative paradigma. The calculations also suggest that the mass of IME scales with the mass of Fe-peak elements, making it difficult to conciliate energetic explosions with low ejected nickel masses, as in the well observed SN1991bg or in SN1998de. Thus a large production of Si-peak elements, especially in combination with a low or a moderate production of iron, could be better addressed by either the delayed detonation route in standard Chandrasekhar-mass models or, perhaps, by the off-center helium detonation in the sub Chandrasekhar-mass scenario.

Detonating Failed Deflagration Model of Thermonuclear Supernovae. I. Explosion Dynamics

We present a detonating failed deflagration model of Type Ia supernovae. In this model, the thermonuclear explosion of a massive white dwarf follows an off-center deflagration. We conduct a survey of asymmetric ignition configurations initiated at various distances from the stellar center. In all cases studied, we find that only a small amount of stellar fuel is consumed during deflagration phase, no explosion is obtained, and the released energy is mostly wasted on expanding the progenitor. Products of the failed deflagration quickly reach the stellar surface, polluting and strongly disturbing it. These disturbances eventually evolve into small and isolated shock-dominated regions which are rich in fuel. We consider these regions as seeds capable of forming self-sustained detonations that, ultimately, result in the thermonuclear supernova explosion. Preliminary nucleosynthesis results indicate the model supernova ejecta are typically composed of about 0.1-0.25 Msun of silicon group elements, 0.9-1.2 Msun of iron group elements, and are essentially carbon-free. The ejecta have a composite morphology, are chemically stratified, and display a modest amount of intrinsic asymmetry. The innermost layers are slightly egg-shaped with the axis ratio ~1.2-1.3 and dominated by the products of silicon burning. This central region is surrounded by a shell of silicon-group elements. The outermost layers of ejecta are highly inhomogeneous and contain products of incomplete oxygen burning with only small admixture of unburned stellar material. The explosion energies are ~1.3-1.5 10^51 erg.

Constraints on the Physics of Type Ia Supernovae from the X‐Ray Spectrum of the Tycho Supernova Remnant

In this paper we use high-quality X-ray observations from XMM-Newton and Chandra to gain new insights into the explosion that originated Tycho's supernova 433 yr ago. We perform a detailed comparison between the ejecta emission from the spatially integrated X-ray spectrum of the supernova remnant and current models for Type Ia supernova explosions. We use a grid of synthetic X-ray spectra based on hydrodynamic models of the evolution of the supernova remnant and nonequilibrium ionization calculations for the state of the shocked plasma. We find that the fundamental properties of the X-ray emission in Tycho are well reproduced by a one-dimensional delayed detonation model with a kinetic energy of ~1.2 ? 1051 ergs. All the other paradigms for Type Ia explosions that we have tested, including one-dimensional deflagrations, pulsating delayed detonations, and sub-Chandrasekhar explosions, as well as deflagration models calculated in three dimensions, fail to provide a good approximation of the observed ejecta emission. Our results require that the supernova ejecta retain some degree of chemical stratification, with Fe peak elements interior to intermediate-mass elements. This strongly suggests that a supersonic burning front (i.e., a detonation) must be involved at some stage in the physics of Type Ia supernova explosions.

Spectral Modeling of SNe Ia Near Maximum Light: Probing the Characteristics of Hydrodynamical Models

We have performed detailed non-local thermodynamic equilibrium (NLTE) spectral synthesis modeling of two types of one-dimensional hydrodynamical models: the very highly parameterized deflagration model W7, and two delayed-detonation models. We find that, overall, both models do about equally well at fitting well-observed SNe Ia near maximum light. However, the Si II ?6150 feature of W7 is systematically too fast, whereas for the delayed-detonation models it is also somewhat too fast but significantly better than that of W7. We find that a parameterized mixed model does the best job of reproducing the Si II ?6150 line near maximum light, and we study the differences in the models that lead to better fits to normal SNe Ia. We discuss what is required of a hydrodynamical model to fit the spectra of observed SNe Ia near maximum light.

Theoretical light curves for deflagration models of type Ia supernova

Aims. We present synthetic bolometric and broad-band UBVRI light curves of SNe Ia for four selected 3D deflagration models of thermonuclear supernovae. Methods. The light curves are computed with the 1D hydro code stella, which models (multi-group time-dependent) nonequilibrium radiative transfer inside SN ejecta. Angle-averaged results from 3D hydrodynamical explosion simulations with the composition determined in a nucleosynthetic postprocessing step served as the input to the radiative transfer model. Results. The predicted model UBV light curves do agree reasonably well with the observed ones for SNe Ia in the range of low to normal luminosities, although the underlying hydrodynamical explosion models produced only a modest amount of radioactive 56 Ni (i.e. ∼0.24–0.42 M� ) and relatively low kinetic energy in the explosion (less than 0.7 × 10 51 erg). The evolution of predicted B and ∞

Type Ia supernova diversity in three-dimensional models

The use of type Ia supernovae as distance indicators for cosmology has initiated a search for theoretical arguments supporting the empirical calibration methods applied. To this end, as a first step, a sound understanding of the origin of the observed diversity in type Ia supernova properties is needed. Here we present a first systematic study of effects resulting from changing some physical parameters of three-dimensional deflagration models of thermonuclear supernovae. In our study we vary the progenitor's carbon-to-oxygen ratio and its central density prior to ignition because both properties are not well determined by stellar evolution theory and they may change from supernova to supernova. Next we compute for these explosion models the nucleosynthesis yields in a post-processing step. This, in addition, allows us to study variations in the progenitor's metallicity by means of different $^{22}$Ne mass fractions in the initial composition. We find that the progenitor's carbon-to-oxygen ratio and its central density affect the energy release of the models and thus the expansion velocity of the supernova. Moreover, we find that changing the metallicity and the central density changes the production of radioactive $^{56}$Ni and thus affects the luminosity. In contrast, the carbon-to-oxygen ratio has little effect on the $^{56}$Ni production. Implications of the found variations of the explosion energy and the produced $^{56}$Ni mass for the type Ia supernova diversity are discussed.