Alan Clark Calder

A comparative study of the turbulent Rayleigh–Taylor instability using high-resolution three-dimensional numerical simulations: The Alpha-Group collaboration

The turbulent Rayleigh–Taylor instability is investigated in the limit of strong mode-coupling using a variety of high-resolution, multimode, three dimensional numerical simulations (NS). The perturbations are initialized with only short wavelength modes so that the self-similar evolution (i.e., bubble diameter Db∝amplitude hb) occurs solely by the nonlinear coupling (merger) of saturated modes. After an initial transient, it is found that hb∼αbAgt2, where A=Atwood number, g=acceleration, and t=time. The NS yield Db∼hb/3 in agreement with experiment but the simulation value αb∼0.025±0.003 is smaller than the experimental value αb∼0.057±0.008. By analyzing the dominant bubbles, it is found that the small value of αb can be attributed to a density dilution due to fine-scale mixing in our NS without interface reconstruction (IR) or an equivalent entrainment in our NS with IR. This may be characteristic of the mode coupling limit studied here and the associated αb may represent a lower bound that is insensiti...

On validating an astrophysical simulation code

We present a case study of validating an astrophysical simulation code. Our study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code for studying the compressible, reactive flows found in many astrophysical environments. We describe the astrophysics problems of interest and the challenges associated with simulating these problems. We describe methodology and discuss solutions to difficulties encountered in verification and validation. We describe verification tests regularly administered to the code, present the results of new verification tests, and outline a method for testing general equations of state. We present the results of two validation tests in which we compared simulations to experimental data. The first is of a laser-driven shock propagating through a multilayer target, a configuration subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported against the force of gravity by a light fluid. Our simulations of the multilayer target experiments showed good agreement with the experimental results, but our simulations of the Rayleigh-Taylor instability did not agree well with the experimental results. We discuss our findings and present results of additional simulations undertaken to further investigate the Rayleigh-Taylor instability.

Type Ia supernova explosion: Gravitationally confined detonation

We present a new mechanism for Type Ia supernova explosions in massive white dwarfs. The scenario follows from relaxing assumptions of symmetry and involves a detonation born near the stellar surface. The explosion begins with an essentially central ignition of a deflagration that results in the formation of a buoyancy-driven bubble of hot material that reaches the stellar surface at supersonic speeds. The bubble breakout laterally accelerates fuel-rich outer stellar layers. This material, confined by gravity to the white dwarf, races along the stellar surface and is focused at the location opposite to the point of the bubble breakout. These streams of nuclear fuel carry enough mass and energy to trigger a detonation just above the stellar surface that will incinerate the white dwarf and result in an energetic explosion. The stellar expansion following the deflagration redistributes mass in a way that ensures production of intermediate-mass and iron group elements with ejecta having a strongly layered structure and a mild amount of asymmetry following from the early deflagration phase. This asymmetry, combined with the amount of stellar expansion determined by details of the evolution (principally the energetics of deflagration, timing of detonation, and structure of the progenitor), can be expected to create a family of mildly diverse Type Ia supernova explosions.

Three-Dimensional Simulations of the Deflagration Phase of the Gravitationally Confined Detonation Model of Type Ia Supernovae

We report the results of a series of three-dimensional (3D) simulations of the deflagration phase of the gravitationally confined detonation mechanism for Type Ia supernovae. In this mechanism, ignition occurs at one or several off-center points, resulting in a burning bubble of hot ash that rises rapidly, breaks through the surface of the star, and collides at a point opposite the breakout on the stellar surface. We find that detonation conditions are robustly reached in our 3D simulations for a range of initial conditions and resolutions. Detonation conditions are achieved as the result of an inwardly directed jet that is produced by the compression of unburnt surface material when the surface flow collides with itself. A high-velocity outwardly directed jet is also produced. The initial conditions explored in this paper lead to conditions at detonation that can be expected to produce large amounts of 56Ni and small amounts of intermediate-mass elements. These particular simulations are therefore relevant only to high-luminosity Type Ia supernovae. Recent observations of Type Ia supernovae imply a compositional structure that is qualitatively consistent with that expected from these simulations.

Morphology of Rising Hydrodynamic and Magnetohydrodynamic Bubbles from Numerical Simulations

Recent Chandra and XMM-Newton observations of galaxy cluster cooling flows have revealed X-ray emission voids of up to 30 kpc in size that have been identified with buoyant, magnetized bubbles. Motivated by these observations, we have investigated the behavior of rising bubbles in stratified atmospheres using the FLASH adaptive-mesh simulation code. We present results from two-dimensional simulations with and without the effects of magnetic fields and with varying bubble sizes and background stratifications. We find purely hydrodynamic bubbles to be unstable; a dynamically important magnetic field is required to maintain a bubbles integrity. This suggests that, even absent thermal conduction, for bubbles to be persistent enough to be regularly observed, they must be supported in large part by magnetic fields. Thermal conduction unmitigated by magnetic fields can dissipate the bubbles even faster. We also observe that the bubbles leave a tail as they rise; the structure of these tails can indicate the history of the dynamics of the rising bubble.

Mapping Initial Hydrostatic Models in Godunov Codes

We look in detail at the process of mapping an astrophysical initial model from a stellar evolution code onto the computational grid of an explicit, Godunov-type code while maintaining hydrostatic equilibrium. This mapping process is common in astrophysical simulations, when it is necessary to follow short-timescale dynamics after a period of long-timescale buildup. We look at the effects of spatial resolution, boundary conditions, the treatment of the gravitational source terms in the hydrodynamics solver, and the initialization process itself. We conclude with a summary detailing the mapping process that yields the lowest ambient velocities in the mapped model.

Flame Evolution During Type Ia Supernovae and the Deflagration Phase in the Gravitationally Confined Detonation Scenario

We develop an improved method for tracking the nuclear flame during the deflagration phase of a Type Ia supernova and apply it in a study of the variation in outcomes expected from the gravitationally confined detonation (GCD) paradigm. A simplified three-stage burning model and a nonstatic ash state are integrated with an artificially thickened advection-diffusion-reaction (ADR) flame front in order to provide an accurate but highly efficient representation of the energy release and electron capture in and after the unresolvable flame. We demonstrate that neither our ADR nor our energy release methods generate significant acoustic noise, as has been a problem with previous ADR-based schemes. We proceed to model aspects of the deflagration, particularly the role of buoyancy of the hot ash, and find that our methods are reasonably well behaved with respect to numerical resolution. We show that if a detonation occurs in material swept up by the material ejected by the first rising bubble but gravitationally confined to the white dwarf (WD) surface (the GCD paradigm), the density structure of the WD at detonation is systematically correlated with the distance of the deflagration ignition point from the center of the star. Coupled to a suitably stochastic ignition process, this correlation may provide a plausible explanation for the variety of nickel masses seen in Type Ia supernovae.

Capturing the Fire: Flame Energetics and Neutronization for Type Ia Supernova Simulations

We develop and calibrate a realistic model flame for hydrodynamic simulations of deflagrations in white dwarf (Type Ia) supernovae. Our flame model builds on the advection-diffusion-reaction model of Khokhlov and includes electron screening and Coulomb corrections to the equation of state in a self-consistent way. We calibrate this model flame—its energetics and timescales for energy release and neutronization—with self-heating reaction network calculations that include both these Coulomb effects and up-to-date weak interactions. The burned material evolves postflame due to both weak interactions and hydrodynamic changes in density and temperature. We develop a scheme to follow the evolution, including neutronization, of the NSE state subsequent to the passage of the flame front. As a result, our model flame is suitable for deflagration simulations over a wide range of initial central densities and can track the temperature and electron fraction of the burned material through the explosion and into the expansion of the ejecta.

High-Performance Reactive Fluid Flow Simulations Using Adaptive Mesh Refinement on Thousands of Processors

We present simulations and performance results of nuclear burning fronts in super- novae on the largest domain and at the finest spatial resolution studied to date. These simulations were performed on the Intel ASCI-Red machine at Sandia National Laboratories using FLASH, a code developed at the Center for Astrophysical Thermonuclear Flashes at the University of Chicago. FLASH is a modular, adaptive mesh, parallel simulation code capable of handling compressible, reactive fluid flows in astrophysical environments. FLASH is written primarily in Fortran 90, uses the Message-Passing Interface library for inter-processor communication and portability, and employs the PARAMESH package to manage a block-structured adaptive mesh that places blocks only where resolution is required and tracks rapidly changing flow features, such as detonation fronts, with ease. We describe the key algorithms and their implementation as well as the optimizations required to achieve sustained performance of 238 GFLOPS on 6420 processors of ASCI-Red in 64 bit arithmetic.

On heavy element enrichment in classical novae

Many classical nova ejecta are enriched in CNO and Ne. Rosner and coworkers recently suggested that the enrichment might originate in the resonant interaction between large-scale shear flows in the accreted H/He envelope and gravity waves at the interface between the envelope and the underlying C/O white dwarf (WD). The shear flow amplifies the waves, which eventually form cusps and break. This wave breaking injects a spray of C/O into the superincumbent H/He. Using two-dimensional simulations, we formulate a quantitative expression for the amount of C/O per unit area that can be entrained, at saturation, into the H/He. The fraction of the envelope that is enriched depends on the horizontal distribution of shear velocity and the density contrast between the C/O WD and the H/He layer but is roughly independent of the vertical shape of the shear profile. Using this parameterization for the mixed mass, we then perform several one-dimensional Lagrangian calculations of an accreting WD envelope and consider two scenarios: that the wave breaking and mixing is driven by the convective flows and that the mixing occurs prior to the onset of convection. In the absence of enrichment prior to ignition, the base of the convective zone, as calculated from mixing-length theory with the Ledoux instability criterion, does not reach the C/O interface. As a result, there is no additional mixing, and the runaway is slow. In contrast, the formation of a mixed layer during the accretion of H/He, prior to ignition, causes a more violent runaway. The envelope can be enriched by 25% of C/O by mass (consistent with that observed in some ejecta) for shear velocities, over the surface, with Mach numbers 0.4.