Astrophysical explosions such as supernovae are fascinating events that require sophisticated algorithms and substantial computational power to model. Castro and MAESTROeX are nuclear astrophysics codes that simulate thermonuclear fusion in the context of supernovae and X-ray bursts. Examining these nuclear burning processes using high resolution simulations is critical for understanding how these astrophysical explosions occur. In this paper we describe the changes that have been made to these codes to transform them from standard MPI + OpenMP codes targeted at petascale CPU-based systems into a form compatible with the pre-exascale systems now online and the exascale systems coming soon. We then discuss what new science is possible to run on systems such as Summit and Perlmutter that could not have been achieved on the previous generation of supercomputers.

We investigate the structure of laterally-propagating flames through the highly-stratified burning layer in an X-ray burst. Two-dimensional hydrodynamics simulations of flame propagation are performed through a rotating plane-parallel atmosphere, exploring the structure of the flame. We discuss the approximations needed to capture the length and time scales at play in an X-ray burst and describe the flame acceleration observed. Our studies complement other multidimensional studies of burning in X-ray bursts.

Modelling long-time convective flows in the interiors of stars is extremely challenging using conventional compressible hydrodynamics codes due to the acoustic timestep limitation. Many of these flows are in the low Mach number regime, which allows us to exploit the relationship between acoustic and advective time scales to develop a more computationally efficient approach. MAESTROeX is an open source low Mach number stellar hydrodynamics code that allows much larger timesteps to be taken, therefore enabling systems to be modelled for much longer periods of time. This is particularly important for the problem of convection in the cores of rotating massive stars prior to core collapse. To fully capture the dynamics, it is necessary to model these systems in three dimensions at high resolution over many rotational periods. We present an overview of MAESTROeX's current capabilities, describe ongoing work to incorporate the effects of rotation and discuss how we are optimising the code to run on GPUs.

We describe recent developments to the Castro astrophysics simulation code, focusing on new features that enable our simulations of X-ray bursts. Two highlights of Castro's ongoing development are the new integration technique to couple hydrodynamics and reactions to high order and GPU offloading. We discuss how these features will help offset some of the computational expense in X-ray burst models.

The creation of "hybrid" white dwarfs, made of a C-O core within a O-Ne shell has been proposed, and studies indicate that ignition in the C-rich central region makes these viable progenitors for thermonuclear (type Ia) supernovae. Recent work found that the C-O core is mixed with the surrounding O-Ne as the white dwarf cools prior to accretion, which results in lower central C fractions in the massive progenitor than previously assumed. To further investigate the efficacy of hybrid white dwarfs as progenitors of thermonuclear supernovae, we performed simulations of thermonuclear supernovae from a new series of hybrid progenitors that include the effects of mixing during cooling. The progenitor white dwarf model was constructed with the one-dimensional stellar evolution code MESA and represented a star evolved through the phase of unstable interior mixing followed by accretion until it reached conditions for the ignition of carbon burning. This MESA model was then mapped to a two-dimensional initial condition for explosions simulated with FLASH. For comparison, similar simulations were performed for a traditional C-O progenitor white dwarf. By comparing the yields of the explosions, we find that, as with earlier studies, the lower C abundance in the hybrid progenitor compared to the traditional C-O progenitor leads to a lower average yield of 56Ni. Although the unmixed hybrid WD showed a similar decrement also in total iron group yield, the mixed case does not and produces a smaller fraction of iron group elements in the form of 56Ni. We attribute this to the higher central density required for ignition and the location, center or off-center, of deflagration ignition.

Thermonuclear (type Ia) supernovae are bright stellar explosions with the unique property that the light curves can be standardized, allowing them to be used as distance indicators for cosmological studies. Many fundamental questions bout these events remain, however. We provide a critique of our present understanding of these and present results of simulations assuming the single-degenerate progenitor model consisting of a white dwarf that has gained mass from a stellar companion. We present results from full three-dimensional simulations of convection with weak reactions comprising the A=23 Urca process in the progenitor white dwarf.

We discuss the challenges of modeling X-ray bursts in multi-dimensions, review the different calculations done to date, and discuss our new set of ongoing simulations. We also describe algorithmic improvements that may help in the future to offset some of the expense of these simulations, and describe what may be possible with exascale computing.

We present early results from a study addressing the question of how one treats the propagation of incertitude, that is, epistemic uncertainty, in input parameters in astrophysical simulations. As an example, we look at the propagation of incertitude in control parameters for stellar winds in MESA stellar evolution simulations. We apply two methods of incertitude propagation, the Cauchy Deviates method and the Quadratic Response Surface method, to quantify the output uncertainty in the final white dwarf mass given a range of values for wind parameters. The methodology we apply is applicable to the problem of propagating input incertitudes through any simulation code treated as a "black box," i.e. a code for which the algorithmic details are either inaccessible or prohibitively complicated. We have made the tools developed for this study freely available to the community.

pynucastro addresses two needs in the field of nuclear astrophysics: visual exploration of nuclear reaction rates or networks and automated code generation for integrating reaction network ODEs. pynucastro accomplishes this by interfacing with nuclear reaction rate parameterizations published by the JINA Reaclib project (Cyburt et al. 2010).

We describe the AMReX suite of astrophysics codes and their application to modeling problems in stellar astrophysics. Maestro is tuned to efficiently model subsonic convective flows while Castro models the highly compressible flows associated with stellar explosions. Both are built on the block-structured adaptive mesh refinement library AMReX. Together, these codes enable a thorough investigation of stellar phenomena, including Type Ia supernovae and X-ray bursts. We describe these science applications and the approach we are taking to make these codes performant on current and future many-core and GPU-based architectures.

Thermonuclear (Type Ia) supernovae are bright stellar explosions, the light curves of which can be calibrated to allow for use as "standard candles" for measuring cosmological distances. Contemporary research investigates how the brightness of an event may be influenced by properties of the progenitor system that follow from properties of the host galaxy such as composition and age. The goals are to better understand systematic effects and to assess the intrinsic scatter in the brightness, thereby reducing uncertainties in cosmological studies. We present the results from ensembles of simulations in the single-degenerate paradigm addressing the influence of age and metallicity on the brightness of an event and compare our results to observed variations of brightness that correlate with properties of the host galaxy. We also present results from "hybrid" progenitor models that incorporate recent advances in stellar evolution.

Motivated by recent results in stellar evolution that predict the existence of hybrid white dwarf (WD) stars with a C-O core inside an O-Ne shell, we simulate thermonuclear (Type Ia) supernovae from these hybrid progenitors. We use the FLASH code to perform multidimensional simulations in the deflagration to detonation transition (DDT) explosion paradigm. Our hybrid progenitor models were produced with the MESA stellar evolution code and include the effects of the Urca process, and we map the progenitor model to the FLASH grid. We performed a suite of DDT simulations over a range of ignition conditions consistent with the progenitor's thermal and convective structure assuming multiple ignition points. To compare the results from these hybrid WD stars to previous results from C-O white dwarfs, we construct a set of C-O WD models with similar properties and similarly simulate a suite of explosions. We find that despite significant variability within each suite, trends distinguishing the explosions are apparent in their $^{56}$Ni yields and the kinetic properties of the ejecta. We comment on the feasibility of these explosions as the source of some classes of observed subluminous events.