Many problems in stellar astrophysics feature flows at low Mach numbers. Conventional compressible hydrodynamics schemes frequently used in the field have been developed for the transonic regime and exhibit excessive numerical dissipation for these f
lows. While schemes were proposed that solve hydrodynamics strictly in the low Mach regime and thus restrict their applicability, we aim at developing a scheme that correctly operates in a wide range of Mach numbers. Based on an analysis of the asymptotic behavior of the Euler equations in the low Mach limit we propose a novel scheme that is able to maintain a low Mach number flow setup while retaining all effects of compressibility. This is achieved by a suitable modification of the well-known Roe solver. Numerical tests demonstrate the capability of this new scheme to reproduce slow flow structures even in moderate numerical resolution. Our scheme provides a promising approach to a consistent multidimensional hydrodynamical treatment of astrophysical low Mach number problems such as convection, instabilities, and mixing in stellar evolution.
A promising model for normal Type Ia supernova (SN Ia) explosions are delayed detonations of Chandrasekhar-mass white dwarfs, in which the burning starts out as a subsonic deflagration and turns at a later phase of the explosion into a supersonic det
onation. The mechanism of the underlying deflagration-to-detonation transition (DDT) is unknown in detail, but necessary conditions have been determined recently. The region of detonation initiation cannot be spatially resolved in multi-dimensional full-star simulations of the explosion. We develop a subgrid-scale (SGS) model for DDTs in thermonuclear supernova simulations that is consistent with the currently known constraints. The probability for a DDT to occur is calculated from the distribution of turbulent velocities measured on the grid scale in the vicinity of the flame and the fractal flame surface area that satisfies further physical constraints, such as fuel fraction and fuel density. The implementation of our DDT criterion provides a solid basis for simulations of thermonuclear supernova explosions in the delayed detonation scenario. It accounts for the currently known necessary conditions for the transition and avoids the inclusion of resolution-dependent quantities in the model. The functionality of our DDT criterion is demonstrated on the example of one three-dimensional thermonuclear supernova explosion simulation.
There are strong indications that the process of conversion of a neutron star into a strange quark star proceeds as a strong deflagration implying that in a few milliseconds almost the whole star is converted. Starting from the three-dimensional hydr
odynamic simulations of the combustion process which provide the temperature profiles inside the newly born strange star, we calculate for the first time the neutrino signal that is to be expected if such a conversion process takes place. The neutrino emission is characterized by a luminosity and a duration that is typical for the signal expected from protoneutron stars and represents therefore a powerful source of neutrinos which could be possibly directly detected in case of events occurring close to our Galaxy. We discuss moreover possible connections between the birth of strange stars and explosive phenomena such as supernovae and gamma-ray-bursts.
We compare models for Type Ia supernova (SN Ia) light curves and spectra with an extensive set of observations. The models come from a recent survey of 44 two-dimensional delayed-detonation models computed by Kasen, Roepke & Woosley (2009), each view
ed from multiple directions. The data include optical light curves of 251 SNe Ia and 2231 low-dispersion spectra from the Center for Astrophysics, plus data from the literature. The analysis uses standard techniques employed by observers, including MLCS2k2, SALT2, and SNooPy for light-curve analysis, and the Supernova Identification (SNID) code of Blondin & Tonry for spectroscopic comparisons to assess how well the models match the data. We show that the models that match observed spectra best lie systematically on the observed width-luminosity relation. Conversely, we reject six models with highly asymmetric ignition conditions and a large amount (>1 M_sun) of synthesized 56Ni that yield poor matches to observed SN Ia spectra. More subtle features of the comparison include the general difficulty of the models to match the U-band flux at early times, caused by a hot ionized ejecta that affect the subsequent redistribution of flux at longer wavelengths. We examine ways in which the asymptotic kinetic energy of the explosion affects both the predicted velocity and velocity gradient in the Si II and Ca II lines. Models with an asymmetric distribution of 56Ni are found to result in a larger variation of photometric and spectroscopic properties with viewing angle, regardless of the initial ignition setup. We discuss more generally whether highly anisotropic ignition conditions are ruled out by observations, and how detailed comparisons between models and observations involving both light curves and spectra can lead to a better understanding of SN Ia explosion mechanisms.
Delayed detonations of Chandrasekhar-mass white dwarfs (WDs) have been very successful in explaining the spectra, light curves, and the width-luminosity relation of spectroscopically normal Type Ia supernovae (SNe Ia). The ignition of the thermonucle
ar deflagration flame at the end of the convective carbon simmering phase in the core of the WD is still not well understood and much about the ignition kernel distribution remains unknown. Furthermore, the central density at the time of ignition depends on the still uncertain screened carbon fusion reaction rates, the accretion history and cooling time of the progenitor, and the composition. We present the results of twelve high-resolution three-dimensional delayed detonation SN Ia explosion simulations that employ a new criterion to trigger the deflagration to detonation transition (DDT). All simulations trigger our DDT criterion and the resulting delayed detonations unbind the star. We find a trend of increasing iron group element (IGE) production with increasing central density for bright, faint, and intermediate SNe. The total 56Ni yield, however, remains more or less constant, even though increased electron captures at high density result in a decreasing 56Ni mass fraction of the IGE material. We attribute this to an approximate balance of 56Ni producing and destroying effects. The deflagrations that were ignited at higher density initially have a faster growth rate of subgrid-scale turbulence. Hence, the effective flame speed increases faster, which triggers the DDT criterion earlier, at a time when the central density of the expanded star is higher. This leads to an overall increase of IGE production, which off-sets the percental reduction of 56Ni due to neutronization.
The theory of radiative transfer provides the link between the physical conditions in an astrophysical object and the observable radiation which it emits. Thus accurately modelling radiative transfer is often a necessary part of testing theoretical m
odels by comparison with observations. We describe a new radiative transfer code which employs Monte Carlo methods for the numerical simulation of radiation transport in expanding media. We discuss the application of this code to the calculation of synthetic spectra and light curves for a Type Ia supernova explosion model and describe the sensitivity of the results to certain approximations made in the simulations.
