Search

The development of turbulent gas flows in the intra-cluster medium and in the core of a galaxy cluster is studied by means of adaptive mesh refinement (AMR) cosmological simulations. A series of six runs was performed, employing identical simulation parameters but different criteria for triggering the mesh refinement. In particular, two different AMR strategies were followed, based on the regional variability of control variables of the flow and on the overdensity of subclumps, respectively. We show that both approaches, albeit with different results, are useful to get an improved resolution of the turbulent flow in the ICM. The vorticity is used as a diagnostic for turbulence, showing that the turbulent flow is not highly volume-filling but has a large area-covering factor, in agreement with previous theoretical expectations. The measured turbulent velocity in the cluster core is larger than 200 km/s, and the level of turbulent pressure contribution to the cluster hydrostatic equilibrium is increased by using the improved AMR criteria., Comment: 13 pages, 14 figures, accepted for publication in MNRAS. Section 4.3.1 (convergence tests of the AMR criteria) and further minor changes added during the reviewing stage. Errors have been corrected in Table 3, but the conclusions are not affected

The problem of the resolution of turbulent flows in adaptive mesh refinement (AMR) simulations is investigated by means of 3D hydrodynamical simulations in an idealised setup, representing a moving subcluster during a merger event. AMR simulations performed with the usual refinement criteria based on local gradients of selected variables do not properly resolve the production of turbulence downstream of the cluster. Therefore we apply novel AMR criteria which are optimised to follow the evolution of a turbulent flow. We demonstrate that these criteria provide a better resolution of the flow past the subcluster, allowing us to follow the onset of the shear instability, the evolution of the turbulent wake and the subsequent back-reaction on the subcluster core morphology. We discuss some implications for the modelling of cluster cold fronts., Comment: 11 pages, 14 figures. Small changes to match the version accepted by MNRAS

The onset of the thermonuclear runaway in a Chandrasekhar-mass white dwarf, leading to the explosion as a type Ia supernova, is studied with hydrodynamical simulations. We investigate the evolution of temperature fluctuations (``bubbles'') in the WD's convective core by means of 2D numerical simulations. We show how the occurrence of the thermonuclear runaway depends on various bubble parameters. The relevance of the progenitor's composition for the ignition process is also discussed., Comment: 4 pages, 2 figures, proceedings of the Catania Workshop on Nuclear and Neutrino Astrophysics (WNNA 2007)

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 simulations' (LES). Such consistency requires to capture the onset of the turbulent cascade on resolved scales. This is achieved by computing the dynamical evolution on a 1024$^3$ 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., Comment: 8 pages, 5 figures, accepted for publication in the ApJ

Aims: We present the first full-star three-dimensional explosion simulations of thermonuclear supernovae including parameterized deflagration-to-detonation transitions that occur once the flame enters the distributed burning regime. Methods: Treating the propagation of both the deflagration and the detonation waves in a common front-tracking approach, the detonation is prevented from crossing ash regions. Results: Our criterion triggers the detonation wave at the outer edge of the deflagration flame and consequently it has to sweep around the complex structure and to compete with expansion. Despite the impeded detonation propagation, the obtained explosions show reasonable agreement with global quantities of observed type Ia supernovae. By igniting the flame in different numbers of kernels around the center of the exploding white dwarf, we set up three different models shifting the emphasis from the deflagration phase to the detonation phase. The resulting explosion energies and iron group element productions cover a large part of the diversity of type Ia supernovae. Conclusions: Flame-driven deflagration-to-detonation transitions, if hypothetical, remain a possibility deserving further investigation., Comment: 4 pages, 1 figure

In the context of explosion models for Type Ia Supernovae, we present one- and two-dimensional simulations of fully resolved detonation fronts in degenerate C+O White Dwarf matter including clumps of previously burned material. The ability of detonations to survive the passage through sheets of nuclear ashes is tested as a function of the width and composition of the ash region. We show that detonation fronts are quenched by microscopically thin obstacles with little sensitivity to the exact ash composition. Front-tracking models for detonations in macroscopic explosion simulations need to include this effect in order to predict the amount of unburned material in delayed detonation scenarios., Comment: 6 pages, 9 figures, uses isotope.sty, accepted for publication in A&A

The dynamics of the explosive burning process is highly sensitive to the flame speed model in numerical simulations of type Ia supernovae. Based upon the hypothesis that the effective flame speed is determined by the unresolved turbulent velocity fluctuations, we employ a new subgrid scale model which includes a localised treatment of the energy transfer through the turbulence cascade in combination with semi-statistical closures for the dissipation and non-local transport of turbulence energy. In addition, subgrid scale buoyancy effects are included. In the limit of negligible energy transfer and transport, the dynamical model reduces to the Sharp-Wheeler relation. According to our findings, the Sharp-Wheeler relation is insuffcient to account for the complicated turbulent dynamics of flames in thermonuclear supernovae. The application of a co-moving grid technique enables us to achieve very high spatial resolution in the burning region. Turbulence is produced mostly at the flame surface and in the interior ash regions. Consequently, there is a pronounced anisotropy in the vicinity of the flame fronts. The localised subgrid scale model predicts significantly enhanced energy generation and less unburnt carbon and oxygen at low velocities compared to earlier simulations., Comment: 13 pages, 10 figures, accepted for publication in Astron. Astrophys.; 3D visualisations not included; complete PDF version can be downloaded from http://www.astro.uni-wuerzburg.de/%7Eschmidt/Paper/SGSModel_II_AA.pdf

We present a one-equation subgrid scale model that evolves the turbulence energy corresponding to unresolved velocity fluctuations in large eddy simulations. The model is derived in the context of the Germano consistent decomposition of the hydrodynamical equations. The eddy-viscosity closure for the rate of energy transfer from resolved toward subgrid scales is localised by means of a dynamical procedure for the computation of the closure parameter. Therefore, the subgrid scale model applies to arbitrary flow geometry and evolution. For the treatment of microscopic viscous dissipation a semi-statistical approach is used, and the gradient-diffusion hypothesis is adopted for turbulent transport. A priori tests of the localised eddy-viscosity closure and the gradient-diffusion closure are made by analysing data from direct numerical simulations. As an a posteriori testing case, the large eddy simulation of thermonuclear combustion in forced isotropic turbulence is discussed. We intend the formulation of the subgrid scale model in this paper as a basis for more advanced applications in numerical simulations of complex astrophysical phenomena involving turbulence., Comment: 18 pages, 10 figures, accepted for publication in Astron. Astrophys.; 3D visualisations not included; complete PDF version can be downloaded from http://www.astro.uni-wuerzburg.de/%7Eschmidt/Paper/SGSModel_I_AA.pdf

In the framework of the Chandrasekhar-mass deflagration model for Type Ia supernovae (SNe Ia), a persisting free parameter is the initial morphology of the flame front, which is linked to the ignition process in the progenitor white dwarf. Previous analytical models indicate that the thermal runaway is driven by temperature perturbations (''bubbles'') that develop in the white dwarf's convective core. In order to probe the conditions at ignition (diameters, temperatures and evolutionary timescales), we have performed hydrodynamical 2D simulations of buoyant bubbles in white dwarf interiors. Our results show that fragmentation occurring during the bubble rise affects the outcome of the bubble evolution. Possible implications for the ignition process of SNe Ia are discussed., Comment: 12 pages, 11 figures, submitted to A&A

We present a systematic survey of the capabilities of type Ia supernova explosion models starting from a number of flame seeds distributed around the center of the white dwarf star. To this end we greatly improved the resolution of the numerical simulations in the initial stages. This novel numerical approach facilitates a detailed study of multi-spot ignition scenarios with up to hundreds of ignition sparks. Two-dimensional simulations are shown to be inappropriate to study the effects of initial flame configurations. Based on a set of three-dimensional models, we conclude that multi-spot ignition scenarios may improve type Ia supernova models towards better agreement with observations. The achievable effect reaches a maximum at a limited number of flame ignition kernels as shown by the numerical models and corroborated by a simple dimensional analysis., Comment: 14 pages, 12 figures with reduced resolution to meet astro-ph file size restriction, full-resolution version available from http://www.mpa-garching.mpg.de/~fritz/publications/astro-ph/multispot.pdf submitted to A&A

We apply an ad hoc model for dynamical ignition in three-dimensional numerical simulations of thermonuclear supernovae assuming pure deflagrations. The model makes use of the statistical description of temperature fluctuations in the pre-supernova core proposed by Wunsch & Woosley (2004). Randomness in time is implemented by means of a Poisson process. We are able to vary the explosion energy and nucleosynthesis depending on the free parameter of the model which controls the rapidity of the ignition process. However, beyond a certain threshold, the strength of the explosion saturates and the outcome appears to be robust with respect to number of ignitions. In the most energetic explosions, we find about 0.75 solar masses of iron group elements. Other than in simulations with simultaneous multi-spot ignition, the amount of unburned carbon and oxygen at radial velocities of a few 1000 km/s tends to be reduced for an ever increasing number of ignition events and, accordingly, more pronounced layering results., Comment: 7 pages, 6 figures, accepted for publication in Astron. Astrophys.; PDF version with full resolution figures available from http://www.astro.uni-wuerzburg.de/~schmidt/Paper/StochIgnt_AA.pdf

We study the dynamics of thermonuclear flames propagating in fuel stirred by stochastic forcing. The fuel consists of carbon and oxygen in a state which is encountered in white dwarfs close to the Chandrasekhar limit. The level set method is applied to represent the flame fronts numerically. The computational domain for the numerical simulations is cubic, and periodic boundary conditions are imposed. The goal is the development of a suitable flame speed model for the small-scale dynamics of turbulent deflagration in thermonuclear supernovae. Because the burning process in a supernova explosion is transient and spatially inhomogeneous, the localised determination of subgrid scale closure parameters is essential. We formulate a semi-localised model based on the dynamical equation for the subgrid scale turbulence energy $k_{\mathrm{sgs}}$. The turbulent flame speed $s_{\mathrm{t}}$ is of the order $\sqrt{2k_{\mathrm{sgs}}}$. In particular, the subgrid scale model features a dynamic procedure for the calculation of the turbulent energy transfer from resolved toward subgrid scales, which has been successfully applied to combustion problems in engineering. The options of either including or suppressing inverse energy transfer in the turbulence production term are compared. In combination with the piece-wise parabolic method for the hydrodynamics, our results favour the latter option. Moreover, different choices for the constant of proportionality in the asymptotic flame speed relation, $s_{\mathrm{t}}\propto\sqrt{2k_{\mathrm{sgs}}}$, are investigated., Comment: 36 pages, 11 figures (separately in PNG format), accepted for publication in Combust. Theory Modelling; PDF version with figures included can be downloaded from http://www.astro.uni-wuerzburg.de/~schmidt/Paper/TurbDefl_CTM.pdf

We show that a flame tracking/capturing scheme originally developed for deflagration fronts can be used to model thermonuclear detonations in multidimensional explosion simulations of type Ia supernovae. After testing the accuracy of the front model, we present a set of two-dimensional simulations of delayed detonations with a physically motivated off-center deflagration-detonation-transition point. Furthermore, we demonstrate the ability of the front model to reproduce the full range of possible interactions of the detonation with clumps of burned material. This feature is crucial for assessing the viability of the delayed detonation scenario., Comment: 7 pages, accepted by A&A

The piece-wise parabolic method (PPM) is applied to simulations of forced isotropic turbulence with Mach numbers $\sim 0.1... 1$. The equation of state is dominated by the Fermi pressure of an electron-degenerate fluid. The dissipation in these simulations is of purely numerical origin. For the dimensionless mean rate of dissipation, we find values in agreement with known results from mostly incompressible turbulence simulations. The calculation of a Smagorinsky length corresponding to the rate of numerical dissipation supports the notion of the PPM supplying an implicit subgrid scale model. In the turbulence energy spectra of various flow realisations, we find the so-called bottleneck phenomenon, i.e., a flattening of the spectrum function near the wavenumber of maximal dissipation. The shape of the bottleneck peak in the compensated spectrum functions is comparable to what is found in turbulence simulations with hyperviscosity. Although the bottleneck effect reduces the range of nearly inertial length scales considerably, we are able to estimate the value of the Kolmogorov constant. For steady turbulence with a balance between energy injection and dissipation, it appears that $C\approx 1.7$. However, a smaller value is found in the case of transonic turbulence with a large fraction of compressive components in the driving force. Moreover, we discuss length scales related to the dissipation, in particular, an effective numerical length scale $\Delta_{\mathrm{eff}}$, which can be regarded as the characteristic smoothing length of the implicit filter associated with the PPM., Comment: 23 pages, 7 figures. Revised version accepted by Comp. Fluids. Not all figures included due to size restriction. Complete PDF available at http://www.astro.uni-wuerzburg.de/%7Eschmidt/Paper/NumDiss_CF.pdf

We present numerical calculations of the parameters $C_{\nu}$, $C_{\epsilon}$ and $C_{\kappa}$ associated with the common closures for turbulence production, dissipation and diffusion. In the case of homogeneous and isotropic turbulence, these parameters are expected to be statistically scale-invariant within the inertial subrange. In order to scrutinise this conjecture, we utilised a generalisation of the Germano filtering formalism, which is applicable to compressible flows as well. The filtering of data obtained from three-dimensional direct numerical simulations of forced isotropic turbulence with Mach numbers in the range $\sim 0.1...1$ then yielded values of the closure parameters associated with different length scales. The results indicate that the closure parameters are nearly universal for subsonic or moderately transonic flows, although the resolution of $432^{3}$ grid cells in our simulations is not quite sufficient to clearly establish scale invariance. In addition, it was found that the customary assumption of a kinetic Prandtl number of about unity for the gradient-diffusion closure is flawed due to the misalignment between turbulent flux and the gradient of the turbulence energy. Nevertheless, sound correlation can be achieved if the flux magnitude rather than the flux vector is locally matched. This conclusion is particularly useful for the family of subgrid scale models based on the turbulence energy equation. Furthermore, the parameter of production $C_{\nu}$ was computed in the fashion of dynamical procedures. Thereby, superior agreement between modelled and explicitly evaluated turbulence stresses in comparison to the eddy-viscosity closure with constant $C_{\nu}$ was verified., Comment: 17 pages, 4 figures; submitted to J. Fluid Mech. PDF file with full resolution figures available at http://www.astro.uni-wuerzburg.de/~schmidt/SimlTurb_JFM.pdf

Type Ia supernovae have recently received considerable attention because it appears that they can be used as "standard candles" to measure cosmic distances out to billions of light years away from us. Observations of type Ia supernovae seem to indicate that we are living in a universe that started to accelerate its expansion when it was about half its present age. These conclusions rest primarily on phenomenological models which, however, lack proper theoretical understanding, mainly because the explosion process, initiated by thermonuclear fusion of carbon and oxygen into heavier elements, is difficult to simulate even on supercomputers. Here, we investigate a new way of modeling turbulent thermonuclear deflagration fronts in white dwarfs undergoing a type Ia supernova explosion. Our approach is based on a level set method which treats the front as a mathematical discontinuity and allows for full coupling between the front geometry and the flow field. New results of the method applied to the problem of type Ia supernovae are obtained. It is shown that in 2-D with high spatial resolution and a physically motivated subgrid scale model for the nuclear flames numerically "converged" results can be obtained, but for most initial conditions the stars do not explode. In contrast, simulations in 3-D, do give the desired explosions and many of their properties, such as the explosion energies, lightcurves and nucleosynthesis products, are in very good agreement with observed type Ia supernovae., Comment: 17 pages, 13 figures, some figures are reduced in resolution, full-resolution version of this text available at http://www.mpa-garching.mpg.de/mpa/publications/preprints/pp2004/MPA1621.pdf -- to appear in "Analysis and Numerics of Conservation Laws", ed. G. Warnecke (Springer-Verlag)

We investigate the interaction of thermonuclear flames in Type Ia supernova explosions with vortical flows by means of numerical simulations. In our study, we focus on small scales, where the flame propagation is no longer dominated by the turbulent cascade originating from large-scale effects. Here, the flame propagation proceeds in the cellular burning regime, resulting from a balance between the Landau-Darrieus instability and its nonlinear stabilization. The interaction of a cellularly stabilized flame front with a vortical fuel flow is explored applying a variety of fuel densities and strengths of the velocity fluctuations. We find that the vortical flow can break up the cellular flame structure if it is sufficiently strong. In this case the flame structure adapts to the imprinted flow field. The transition from the cellularly stabilized front to the flame structure dominated by vortices of the flow proceeds in a smooth way. The implications of the results of our simulations for Type Ia Supernova explosion models are discussed., Comment: 13 pages, 2 tables, 10 figures, resolution of Figs. 4 and 5 degraded due to archive file size restrictions, submitted to A&A

We present a numerical investigation of the cellular burning regime in Type Ia supernova explosions. This regime holds at small scales (i.e. below the Gibson scale), which are unresolved in large-scale Type Ia supernova simulations. The fundamental effects that dominate the flame evolution here are the Landau-Darrieus instability and its nonlinear stabilization, leading to a stabilization of the flame in a cellular shape. The flame propagation into quiescent fuel is investigated addressing the dependence of the simulation results on the specific parameters of the numerical setup. Furthermore, we investigate the flame stability at a range of fuel densities. This is directly connected to the questions of active turbulent combustion (a mechanism of flame destabilization and subsequent self-turbulization) and a deflagration-to-detonation transition of the flame. In our simulations we find no substantial destabilization of the flame when propagating into quiescent fuels of densities down to ~10^7 g/cm^3, corroborating fundamental assumptions of large-scale SN Ia explosion models. For these models, however, we suggest an increased lower cutoff for the flame propagation velocity to take the cellular burning regime into account., Comment: 12 pages, 2 tables, 10 figures, resolution of figures degraded due to archive file size restrictions, submitted to A&A

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 new

We present the results of three-dimensional hydrodynamical simulations of the subsonic thermonuclear burning phase in type Ia supernovae. The burning front model contains no adjustable parameters so that variations of the explosion outcome can be linked directly to changes in the initial conditions. In particular, we investigate the influence of the initial flame geometry on the explosion energy and find that it appears to be weaker than in 2D. Most importantly, our models predict global properties such as the produced nickel masses and ejecta velocities within their observed ranges without any fine tuning., Comment: 7 pages, 5 figures, accepted by A&A

As a simple model for unknown Planck scale physics, we assume that the quantum modes responsible for producing primordial curvature perturbations during inflation are placed in their instantaneous adiabatic vacuum when their proper momentum reaches a fixed high energy scale M. The resulting power spectrum is derived and presented in a form that exhibits the amplitude and frequency of the superimposed oscillations in terms of H/M and the slow roll parameter epsilon. The amplitude of the oscillations is proportional to the third power of H/M. We argue that these small oscillations give the lower bound of the modifications of the power spectrum if the notion of free mode propagation ceases to exist above the critical energy scale M., Comment: 10 pages; matches version accepted by PRD

We present a numerical model to study the behavior of thermonuclear flames in the discontinuity approximation. This model is applied to investigate the Landau-Darrieus instability under conditions found in Type Ia supernova explosions of Chandrasekhar mass white dwarfs. This is a first step to explore the flame microstructure in these events. The model reproduces Landau's linearized stability analysis in early stages of the flame evolution and the stabilization in a cellular flame structure in the nonlinear stage., Comment: 6 pages, 3 figures, to appear in Proc. 11th Workshop on Nuclear Astrophysics, Ringberg Castle, Germany, 2002

Type Ia supernovae have become an indispensable tool for studying the expansion history of the universe, yet our understanding of the explosion mechanism is still incomplete. We describe the variety of discussed scenarios, sketch the most relevant physics, and report recent advances in multidimensional simulations of Chandrasekhar mass white dwarf explosions., Comment: Based on an invited review at the XIII Rencontres de Blois - Frontiers of the Universe, June 17-23,2001

Following up on earlier work on this topic (Reinecke et al. 1999, A&A 347, pp. 724 and 739), we present an improved set of numerical models for simulations of white dwarfs exploding as Type Ia supernovae (SNe Ia). Two-dimensional simulations were used to test the reliability and numerical robustness of these algorithms; the results indicate that integral quantities like the total energy release are insensitive to changes of the grid resolution (above a certain threshold), which was not the case for our former code. The models were further enhanced to allow fully three-dimensional simulations of SNe Ia. A direct comparison of a 2D and a 3D calculation with identical initial conditions shows that the explosion is considerably more energetic in three dimensions; this is most likely caused by the assumption of axisymmetry in 2D, which inhibits the growth of flame instabilities in the azimuthal direction and thereby decreases the flame surface., Comment: 8 pages, 6 figures, submitted to A&A

It is shown that varying speed of light cosmology follows from a string-inspired minimal length uncertainty relation. Due to the reduction of the available phase space volume per quantum mode at short wavelengths, the equation of state of ultrarelativistic particles stiffens at very high densities. This causes a stronger than usual deceleration of the scale factor which competes with a higher than usual propagation speed of the particles. Various measures for the effective propagation speed are analyzed: the group and phase velocity in the high energy tail, the thermal average of the group and phase velocity, and the speed of sound. Of these three groups, only the first provides a possible solution to the cosmological horizon problem., Comment: 5 pages, 2 figures

In this article we discuss the first simulations of two- and three-dimensional Type Ia supernovae with an improved hydrodynamics code. After describing the various enhancements, the obtained results are compared to those of earlier code versions, observational data and the findings of other researchers in this field., Comment: 7 pages, 4 figures

It has been pointed out that the perturbation spectrum predicted by inflation may be sensitive to a natural ultraviolet cutoff, thus potentially providing an experimentally accessible window to aspects of Planck scale physics. A priori, a natural ultraviolet cutoff could take any form, but a fairly general classification of possible Planck scale cutoffs has been given. One of those categorized cutoffs, also appearing in various studies of quantum gravity and string theory, has recently been implemented into the standard inflationary scenario. Here, we continue this approach by investigating its effects on the predicted perturbation spectrum. We find that the size of the effect depends sensitively on the scale separation between cutoff and horizon during inflation., Comment: 6 pages; matches version accepted by PRD

We investigate the insensitivity of the predictions of inflationary models with respect to modifications of Planck energy physics. The modification we consider consists in replacing the usual dispersion relation by nonlinear ones. This way of addressing the problem has recently recieved attention and contradictory results were found. Our main result is to show that the adiabaticity of the mode propagation and the separation of two scales of interest, the Planck scale and the cosmological horizon scale, are sufficient conditions for the predictions to be unchanged. We then show that almost all models satisfy the first condition if the second is met. Therefore the introduction of a nonlinear dispersion is unlikely to have any discernable effects on the power spectrum of cosmological perturbations., Comment: 4 pages; v.2: minor modifications, added references

Type Ia supernovae (SN Ia) are generally believed to be the result of the thermonuclear disruption of Chandrasekhar-mass carbon-oxygen white dwarfs, mainly because such thermonuclear explosions can account for the right amount of nickel, which is needed to explain the light curves and the late-time spectra, and the abundances of intermediate-mass nuclei which dominate the spectra near maximum light. Because of their enormous brightness and apparent homogeneity SN Ia have become an important tool to measure cosmological parameters. In this article the present understanding of the physics of thermonuclear explosions is reviewed. In particular, we focus our attention on subsonic ("deflagration") fronts, i.e. we investigate fronts propagating by heat diffusion and convection rather than by compression. Models based upon this mode of nuclear burning have been applied very successfully to the SN Ia problem, and are able to reproduce many of their observed features remarkably well. However, the models also indicate that SN Ia may differ considerably from each other, which is of importance if they are to be used as standard candles., Comment: 9 pages, 3 figures, conference proceedings

Because calibrated light curves of Type Ia supernovae have become a major tool to determine the local expansion rate of the Universe and also its geometrical structure, considerable attention has been given to models of these events over the past couple of years. There are good reasons to believe that perhaps most Type Ia supernovae are the explosions of white dwarfs that have approached the Chandrasekhar mass, M_ch ~ 1.39 M_sun, and are disrupted by thermonuclear fusion of carbon and oxygen. However, the mechanism whereby such accreting carbon-oxygen white dwarfs explode continues to be uncertain. Recent progress in modeling Type Ia supernovae as well as several of the still open questions are addressed in this review. Although the main emphasis will be on studies of the explosion mechanism itself and on the related physical processes, including the physics of turbulent nuclear combustion in degenerate stars, we also discuss observational constraints., Comment: 38 pages, 4 figures, Annual Review of Astronomy and Astrophysics, in press

Type Ia supernovae (SN Ia) are generally believed to be the result of the thermonuclear disruption of Chandrasekhar-mass carbon-oxygen white dwarfs, mainly because such thermonuclear explosions can account for the right amount of Ni-56, which is needed to explain the light curves and the late-time spectra, and the abundances of intermediate-mass nuclei which dominate the spectra near maximum light. Because of their enormous brightness and apparent homogeneity SN Ia have become an important tool to measure cosmological parameters. In this article the present understanding of the physics of thermonuclear explosions is reviewed. In particular, we focus our attention on subsonic (``deflagration'') fronts, i.e. we investigate fronts propagating by heat diffusion and convection rather than by compression. Models based upon this mode of nuclear burning have been applied very successfully to the SN Ia problem, and are able to reproduce many of their observed features remarkably well. However, the models also indicate that SN Ia may differ considerably from each other, which is of importance if they are to be used as standard candles., Comment: 11 pages, 4 figures. To appear in Proc. 10th Ann. Astrophys. Conf. "Cosmic Explosions", Univ. of Maryland 1999, eds. S.S. Holt and W.W. Zhang

The implementation of a Planck-scale high frequency and short wavelength cutoff in quantum theories on expanding backgrounds may have potentially nontrivial implications, such as the breaking of local Lorentz invariance and the existence of a yet unknown mechanism for the creation of vacuum modes. In scenarios where inflation begins close to the cutoff scale, these effects could have observable consequences as trans-Planckian modes are redshifted to cosmological scales. In close analogy with similar studies of Hawking radiation, a simple theory of a minimally coupled scalar field in de Sitter space is studied, with a high frequency cutoff imposed by a nonlinear dispersion relation. Under certain conditions the model predicts deviations from the standard inflationary scenario. We also comment on the difficulties in generalizing fluid models of Hawking radiation to cosmological space-times., Comment: 8 pages

The burning regimes encountered by laminar deflagrations and ZND detonations propagating through helium-rich compositions in the presence of buoyancy-driven turbulence are analyzed. Particular attention is given to models of X-ray bursts which start with a thermonuclear runaway on the surface of a neutron star, and the thin shell helium instability of intermediate-mass stars. In the X-ray burst case, turbulent deflagrations propagating in the lateral or radial directions encounter a transition from the distributed regime to the flamlet regime at a density of 10^8 g cm^{-3}. In the radial direction, the purely laminar deflagration width is larger than the pressure scale height for densities smaller than 10^6 g cm^{-3}. Self-sustained laminar deflagrations travelling in the radial direction cannot exist below this density. Similarily, the planar ZND detonation width becomes larger than the pressure scale height at 10^7 g cm^{-3}, suggesting that a steady-state, self-sustained detonations cannot come into existance in the radial direction. In the thin helium shell case, turbulent deflagrations travelling in the lateral or radial directions encounter the distributed regime at densities below 10^7 g cm^{-3}, and the flamelet regime at larger densities. In the radial direction, the purely laminar deflagration width is larger than the pressure scale height for densities smaller than 10^4 g cm^{-3}, indicating that steady-state laminar deflagrations cannot form below this density. The planar ZND detonation width becomes larger than the pressure scale height at 5 10^4 g cm^{-3}, suggesting that steady-state, self-sustained detonations cannot come into existance in the radial direction., Comment: 13 pages, 2 figures. ApJ, in press

We present a new way of modeling turbulent thermonuclear deflagration fronts in Chandrasekhar-mass white dwarfs, consisting of carbon and oxygen, undergoing a type Ia supernova explosion. Our approach is a front capturing/tracking hybrid scheme, based on a level set method, which treats the front as a mathematical discontinuity and allows for full coupling between the front geometry and the flow field. First results of the method applied to the problem of type Ia supernovae are discussed. It will be shown that even in 2-D and even with a physically motivated sub-grid model numerically ``converged'' results are difficult to obtain., Comment: 13 pages, 5 figures. To appear in Proc. ``Two decades of numerical astrophysics'', Port Douglas, Australia 1999

We present a statistical model which shows the influence of turbulence on a thermonuclear flame propagating in C+O white dwarf matter. Based on a Monte Carlo description of turbulence, it provides a method for investigating the physics in the so-called distributed burning regime. Using this method we perform numerical simulations of turbulent flames and show that in this particular regime the flamelet model for the turbulent flame velocity loses its validity. In fact, at high turbulent intensities burning in the distributed regime can lead to a deceleration of the turbulent flame and thus induces a competing process to turbulent effects that cause a higher flame speed. It is also shown that in dense C+O matter turbulent heat transport is described adequately by the Peclet number, Pe, rather than by the Reynolds number, which means that flame propagation is decoupled from small-scale turbulence. Finally, at the onset of our results we argue that the available turbulent energy in an exploding C+O white dwarf is probably too low in order to make a deflagration to detonation transition possible., Comment: 27 pages, 7 figures, accepted for publication in the ApJ

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 Letters

Microscopic turbulence-flame interactions of thermonuclear fusion flames occuring in Type Ia Supernovae were studied by means of incompressible direct numerical simulations with a highly simplified flame description. The flame is treated as a single diffusive scalar field with a nonlinear source term. It is characterized by its Prandtl number, Pr << 1, and laminar flame speed, S_L. We find that if S_L ~ u', where u' is the rms amplitude of turbulent velocity fluctuations, the local flame propagation speed does not significantly deviate from S_L even in the presence of velocity fluctuations on scales below the laminar flame thickness. This result is interpreted in the context of subgrid-scale modeling of supernova explosions and the mechanism for deflagration-detonation-transitions., Comment: 8 pages, 6 figures, accepted by Astrophys. J

We present a numerical investigation of the gravitational collapse of horizon-size density fluctuations to primordial black holes (PBHs) during the radiation-dominated phase of the Early Universe. The collapse dynamics of three different families of initial perturbation shapes, imposed at the time of horizon crossing, is computed. The perturbation threshold for black hole formation, needed for estimations of the cosmological PBH mass function, is found to be $\delta_{\rm c} \approx 0.7$ rather than the generally employed $\delta_{\rm c} \approx 1/3$, if $\delta$ is defined as $\Delta M/\mh$, the relative excess mass within the initial horizon volume. In order to study the accretion onto the newly formed black holes, we use a numerical scheme that allows us to follow the evolution for long times after formation of the event horizon. In general, small black holes (compared to the horizon mass at the onset of the collapse) give rise to a fluid bounce that effectively shuts off accretion onto the black hole, while large ones do not. In both cases, the growth of the black hole mass owing to accretion is insignificant. Furthermore, the scaling of black hole mass with distance from the formation threshold, known to occur in near-critical gravitational collapse, is demonstrated to apply to primordial black hole formation., Comment: 10 pages, 8 figures, revtex style, submitted to PRD

Primordial black holes (PBHs) may form in the early universe when pre-existing adiabatic density fluctuations enter into the cosmological horizon and recollapse. It has been suggested that PBH formation may be facilitated when fluctuations enter into the horizon during a strongly first-order phase transition which proceeds in approximate equilibrium. We employ general-relativistic hydrodynamics numerical simulations in order to follow the collapse of density fluctuations during first-order phase transitions. We find that during late stages of the collapse fluctuations separate into two regimes, an inner part existing exclusively in the high-energy density phase with energy density $\epsilon_{\rm h}$, surrounded by an outer part which exists exclusively in the low-energy density phase with energy density $\epsilon_{\rm h}-L$, where $L$ is the latent heat of the transition. We confirm that the fluctuation density threshold $\delta\epsilon /\epsilon$ required for the formation of PBHs during first-order transitions decreases with increasing $L$ and falls below that for PBH formation during ordinary radiation dominated epochs. Our results imply that, in case PBHs form at all in the early universe, their mass spectrum is likely dominated by the approximate horizon masses during epochs when the universe undergoes phase transitions., Comment: 8 pages, 4 figures, revtex style, submitted to PRD

We present a new way of modeling deflagration fronts in reactive fluids, the main emphasis being on turbulent thermonuclear deflagration fronts in white dwarfs undergoing a Type Ia supernova explosion. Our approach is based on a level set method which treats the front as a mathematical discontinuity and allows full coupling between the front geometry and the flow field. With only minor modifications, this method can also be applied to describe contact discontinuities. Two different implementations are described and their physically correct behaviour for simple testcases is shown. First results of the method applied to the concrete problems of Type Ia supernovae and chemical hydrogen combustion are briefly discussed; a more extensive analysis of our astrophysical simulations is given in (Reinecke et al. 1998, MPA Green Report 1122b)., Comment: 11 pages, 13 figures, accepted by A&A, corrected and extended according to referee's comments

First results of simulations are presented which compute the dynamical evolution of a Chandrasekhar-mass white dwarf, consisting of equal amounts of carbon and oxygen, from the onset of violent thermonuclear burning, by means of a new two-dimensional numerical code. Since in the interior of such a massive white dwarf nuclear burning progresses on microscopic scales as a sharp discontinuity, a so-called flamelet, which cannot be resolved by any numerical scheme, and since on macroscopic scales the burning front propagates due to turbulence, we make an attempt to model both effects explicitly in the framework of a finite-volume hydrodynamics code. Turbulence is included by a sub-grid model, following the spirit of large eddy simulations, and the well-localized burning front is treated by means of a level set, which allows us to compute the geometrical structure of the front more accurately than with previous methods. The only free parameters of our simulations are the location and the amount of nuclear fuel that is ignited as an initial perturbation. We find that models in which explosive carbon burning is ignited at the center remain bound by the time the front reaches low densities, where we stopped the computations because our description of combustion is no longer applicable. In contrast, off-center ignition models give rise to explosions which, however, are still too weak for typical Type Ia supernovae. Possible reasons for this rather disappointing result are discussed., Comment: 10 pages, 10 figures, accepted by A&A, corrected and extended according to referee's comments

First results of a numerical investigation of primordial black hole formation in the radiation dominated phase of the Early Universe are presented. The simulations follow the gravitational collapse of three different families of high-amplitude density fluctuations imposed at the time of horizon crossing. The threshold for black hole formation, \delta_{c} \approx 0.7, is found to be nearly identical for all perturbation families if the control parameter, \delta, is chosen as the total excess mass within the initial horizon volume. Furthermore, we demonstrate that the scaling of black hole mass with distance from the formation threshold, known to occur in near-critical gravitational collapse, applies to primordial black hole formation., Comment: 4 pages, 1 figure included, uses prcrc.sty. To appear in the proceedings of Dark Matter '98, Los Angeles (ed. D. Cline). Also available at http://bigwhirl.uchicago.edu/jcn/pub_pbh.html

In the framework of the Chandrasekhar mass white dwarf model for Type Ia supernovae, various stages of the explosion are described in terms of the burning regimes of the thermonuclear flame front. In the early flamelet regime following the ``smoldering'' phase prior to the explosion, the flame is sufficiently thin and fast to remain laminar on small scales. As the white dwarf density declines, the thermal flame structure becomes subject to penetration by turbulent eddies, and it enters the ``distributed burning'' regime. A specific control parameter for this transition is proposed. Furthermore, we outline an argument for the coincidence of the transition between burning regimes with the onset of a deflagration-detonation-transition (DDT) in the late phase of the explosion., Comment: 5 pages, to appear in the proceedings of the Second Oak Ridge Symposium on Atomic and Nuclear Astrophysics. See also http://bigwhirl.uchicago.edu/jcn/pub_snIa.html

The recent discovery of critical phenomena arising in gravitational collapse near the threshold of black hole formation is used to estimate the initial mass function of primordial black holes (PBHs). It is argued that the universal scaling relation between black hole mass and initial perturbation found for a variety of collapsing space-times also applies to PBH formation, indicating the possibility of the formation of PBHs with masses much smaller than one horizon mass. Owing to the natural fine-tuning of initial conditions by the exponential decline of the probability distribution for primordial density fluctuations, sub-horizon mass PBHs are expected to form at all epochs. This result suggests that the constraints on the primordial fluctuation spectrum based on the abundance of PBHs at different mass scales may have to be revisited., Comment: 4 pages, uses revtex, also available at http://bigwhirl.uchicago.edu/jcn/pub_pbh.html . To appear in Phys. Rev. Lett

Gibson scaling and related properties of flame-surface geometry in turbulent premixed combustion are demonstrated using a novel computational model, Deterministic Turbulent Mixing (DTM). In DTM, turbulent advection is represented by a sequence of maps applied to the computational domain. The structure of the mapping sequence incorporates pertinent scaling properties of the turbulent cascade. Here, combustion in Kolmogorov turbulence (kinetic-energy cascade) and in Bolgiano-Obukhov convective turbulence (potential-energy cascade) is simulated. Implications with regard to chemical flames and astrophysical (thermonuclear) flames are noted., Comment: 16 pages, 5 ps figures included, accepted by Combust. Sci. Tech., also available at http://www.mpa-garching.mpg.de/~jcn/pub_flame.html

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 ApJ

A series of two dimensional numerical simulations of explosive nuclear burning is presented for white dwarfs near the Chandraskhar mass. We assume that the burning begins as a slow deflagration front at or near the center of the star, and continues until the density in the burning regions has declined to about 10$^7$ g cm$^{-3}$, where the flame is essentially extinguished. We employ a novel numerical representation of the turbulent flame brush based upon ideas previously developed for modelling laboratory combustion and explore in some detail the sensitivity of the outcome to the manner in which burning is initiated. In particular, we simulate 1) a centrally ignited deflagration, 2) off-center ignition at a single ``point", and 3) simultaneous off-center ignition at five ``points". We find that the amount of $^{56}$Ni that is produced and other observable properties depend sensitively upon how the fuel is ignited., Comment: 14 pages, 9 figures included (200kB), uses aas2pp4.sty. Accepted by ApJ