Performing high-resolution, high-fidelity, three-dimensional simulations of Type Ia supernovae (SNe Ia) requires not only algorithms that accurately represent the correct physics, but also codes that effectively harness the resources of the most powerful supercomputers. We are developing a suite of codes that provide the capability to perform end-to-end simulations of SNe Ia, from the early convective phase leading up to ignition to the explosion phase in which deflagration/detonation waves explode the star to the computation of the light curves resulting from the explosion. In this paper we discuss these codes with an emphasis on the techniques needed to scale them to petascale architectures. We also demonstrate our ability to map data from a low Mach number formulation to a compressible solver.

Type Ia supernovae (SNe Ia) are the largest thermonuclear explosions in the Universe. Their light output can be seen across great distances and has led to the discovery that the expansion rate of the Universe is accelerating. Despite the significance of SNe Ia, there are still a large number of uncertainties in current theoretical models. Computational modeling offers the promise to help answer the outstanding questions. However, even with today's supercomputers, such calculations are extremely challenging because of the wide range of length and time scales. In this paper, we discuss several new algorithms for simulations of SNe Ia and demonstrate some of their successes.

The convective period leading up to a Type Ia supernova (SN Ia) explosion is characterized by very low Mach number flows, requiring hydrodynamical methods well-suited to long-time integration. We continue the development of the low Mach number equation set for stellar scale flows by incorporating the effects of heat release due to external sources. Low Mach number hydrodynamics equations with a time-dependent background state are derived, and a numerical method based on the approximate projection formalism is presented. We demonstrate through validation with a fully compressible hydrodynamics code that this low Mach number model accurately captures the expansion of the stellar atmosphere as well as the local dynamics due to external heat sources. This algorithm provides the basis for an efficient simulation tool for studying the ignition of SNe Ia.

We introduce a low Mach number equation set for the large-scale numerical simulation of carbon-oxygen white dwarfs experiencing a thermonuclear deflagration. Since most of the interesting physics in a Type Ia supernova transpires at Mach numbers from 0.01 to 0.1, such an approach enables both a considerable increase in accuracy and savings in computer time compared with frequently used compressible codes. Our equation set is derived from the fully compressible equations using low Mach number asymptotics, but without any restriction on the size of perturbations in density or temperature. Comparisons with simulations that use the fully compressible equations validate the low Mach number model in regimes where both are applicable. Comparisons to simulations based on the more traditional anelastic approximation also demonstrate the agreement of these models in the regime for which the anelastic approximation is valid. For low Mach number flows with potentially finite amplitude variations in density and temperature, the low Mach number model overcomes the limitations of each of the more traditional models and can serve as the basis for an accurate and efficient simulation tool.

We extend a low Mach number hydrodynamics method developed for terrestrial combustion, to the study of thermonuclear flames in Type Ia supernovae. We discuss the differences between 2-D and 3-D Rayleigh-Taylor unstable flame simulations, and give detailed diagnostics on the turbulence, showing that the kinetic energy power spectrum obeys Bolgiano-Obukhov statistics in 2-D, but Kolmogorov statistics in 3-D. Preliminary results from 3-D reacting bubble calculations are shown, and their implications for ignition are discussed.

Flame instabilities play a dominant role in accelerating the burning front to a large fraction of the speed of sound in a Type Ia supernova. We present a three-dimensional numerical simulation of a Rayleigh-Taylor unstable carbon flame, following its evolution through the transition to turbulence. A low Mach number hydrodynamics method is used, freeing us from the harsh timestep restrictions imposed by sound waves. We fully resolve the thermal structure of the flame and its reaction zone, eliminating the need for a flame model. A single density is considered, $1.5\times 10^7 \gcc$, and half carbon\slash half oxygen fuel--conditions under which the flame propagated in the flamelet regime in our related two-dimensional study. We compare to a corresponding two-dimensional simulation, and show that while fire-polishing keeps the small features suppressed in two dimensions, turbulence wrinkles the flame on far smaller scales in the three-dimensional case, suggesting that the transition to the distributed burning regime occurs at higher densities in three dimensions. Detailed turbulence diagnostics are provided. We show that the turbulence follows a Kolmogorov spectrum and is highly anisotropic on the large scales, with a much larger integral scale in the direction of gravity. Furthermore, we demonstrate that it becomes more isotropic as it cascades down to small scales. Based on the turbulent statistics and the flame properties of our simulation, we compute the Gibson scale. We show the progress of the turbulent flame through a classic combustion regime diagram, indicating that the flame just enters the distributed burning regime near the end of our simulation.

Planar flames are intrinsically unstable in open domains due to the thermal expansion across the burning front--the Landau-Darrieus instability. This instability leads to wrinkling and growth of the flame surface, and corresponding acceleration of the flame, until it is stabilized by cusp formation. We look at the Landau-Darrieus instability for C/O thermonuclear flames at conditions relevant to the late stages of a Type Ia supernova explosion. Two-dimensional direct numerical simulations of both single-mode and multi-mode perturbations using a low Mach number hydrodynamics code are presented. We show the effect of the instability on the flame speed as a function of both the density and domain size, demonstrate the existence of the small scale cutoff to the growth of the instability, and look for the proposed breakdown of the non-linear stabilization at low densities. The effects of curvature on the flame as quantified through measurements of the growth rate and computation of the corresponding Markstein number. While accelerations of a few percent are observed, they are too small to have any direct outcome on the supernova explosion.

A Type Ia supernova explosion likely begins as a nuclear runaway near the center of a carbon-oxygen white dwarf. The outward propagating flame is unstable to the Landau-Darrieus, Rayleigh-Taylor, and Kelvin-Helmholtz instabilities, which serve to accelerate it to a large fraction of the speed of sound. We investigate the Rayleigh-Taylor unstable flame at the transition from the flamelet regime to the distributed-burning regime, around densities of $10^7$ g/cc, through detailed, fully resolved simulations. A low Mach number, adaptive mesh hydrodynamics code is used to achieve the necessary resolution and long time scales. As the density is varied, we see a fundamental change in the character of the burning--at the low end of the density range the Rayleigh-Taylor instability dominates the burning, whereas at the high end the burning suppresses the instability. In all cases, significant acceleration of the flame is observed, limited only by the size of the domain we are able to study. We discuss the implications of these results on the potential for a deflagration to detonation transition.