-
Practical Effects of Integrating Temperature with Strang Split Reactions
Authors:
M. Zingale,
M. P. Katz,
D. E. Willcox,
A. Harpole
Abstract:
For astrophysical reacting flows, operator splitting is commonly used to couple hydrodynamics and reactions. Each process operates independent of one another, but by staggering the updates in a symmetric fashion (via Strang splitting) second order accuracy in time can be achieved. However, approximations are often made to the reacting system, including the choice of whether or not to integrate tem…
▽ More
For astrophysical reacting flows, operator splitting is commonly used to couple hydrodynamics and reactions. Each process operates independent of one another, but by staggering the updates in a symmetric fashion (via Strang splitting) second order accuracy in time can be achieved. However, approximations are often made to the reacting system, including the choice of whether or not to integrate temperature with the species. Here we demonstrate through a simple convergence test that integrating an energy equation together with reactions achieves the best convergence when modeling reactive flows with Strang splitting. Additionally, second order convergence cannot be achieved without integrating an energy or temperature equation.
△ Less
Submitted 23 March, 2021;
originally announced March 2021.
-
Dynamics of Laterally Propagating Flames in X-ray Bursts. II. Realistic Burning & Rotation
Authors:
A. Harpole,
N. M. Ford,
K. Eiden,
M. Zingale,
D. E. Willcox,
Y. Cavecchi,
M. P. Katz
Abstract:
We continue to investigate two-dimensional laterally propagating flames in type I X-ray bursts using fully compressible hydrodynamics simulations. In the current study we relax previous approximations where we artificially boosted the flames. We now use more physically realistic reaction rates, thermal conductivities, and rotation rates, exploring the effects of neutron star rotation rate and ther…
▽ More
We continue to investigate two-dimensional laterally propagating flames in type I X-ray bursts using fully compressible hydrodynamics simulations. In the current study we relax previous approximations where we artificially boosted the flames. We now use more physically realistic reaction rates, thermal conductivities, and rotation rates, exploring the effects of neutron star rotation rate and thermal structure on the flame. We find that at lower rotation rates the flame becomes harder to ignite, whereas at higher rotation rates the nuclear burning is enhanced by increased confinement from the Coriolis force and the flame propagates steadily. At higher crustal temperatures, the flame moves more quickly and accelerates as it propagates through the atmosphere. If the temperature is too high, instead of a flame propagating across the surface the entire atmosphere burns steadily. All of the software used for these simulations is freely available.
△ Less
Submitted 26 March, 2021; v1 submitted 29 January, 2021;
originally announced February 2021.
-
Preparing Nuclear Astrophysics for Exascale
Authors:
Max P. Katz,
Ann Almgren,
Maria Barrios Sazo,
Kiran Eiden,
Kevin Gott,
Alice Harpole,
Jean M. Sexton,
Don E. Willcox,
Weiqun Zhang,
Michael Zingale
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 July, 2020;
originally announced July 2020.
-
Dynamics of Laterally Propagating Flames in X-ray Bursts. I. Burning Front Structure
Authors:
Kiran Eiden,
Michael Zingale,
Alice Harpole,
Donald Willcox,
Yuri Cavecchi,
Max P. Katz
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 10 December, 2019;
originally announced December 2019.
-
Modelling low Mach number stellar hydrodynamics with MAESTROeX
Authors:
A. Harpole,
D. Fan,
M. P. Katz,
A. J. Nonaka,
D. E. Willcox,
M. Zingale
Abstract:
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 so…
▽ More
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.
△ Less
Submitted 28 October, 2019;
originally announced October 2019.
-
The Castro AMR Simulation Code: Current and Future Developments
Authors:
M. Zingale,
A. S. Almgren,
M. Barrios Sazo,
J. B. Bell,
K. Eiden,
A. Harpole,
M. P. Katz,
A. J. Nonaka,
D. E. Willcox,
W. Zhang
Abstract:
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 model…
▽ More
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.
△ Less
Submitted 28 October, 2019;
originally announced October 2019.
-
MAESTROeX: A Massively Parallel Low Mach Number Astrophysical Solver
Authors:
Duoming Fan,
Andrew Nonaka,
Ann S. Almgren,
Alice Harpole,
Michael Zingale
Abstract:
We present MAESTROeX, a massively parallel solver for low Mach number astrophysical flows. The underlying low Mach number equation set allows for efficient, long-time integration for highly subsonic flows compared to compressible approaches. MAESTROeX is suitable for modeling full spherical stars as well as well as planar simulations of dynamics within localized regions of a star, and can robustly…
▽ More
We present MAESTROeX, a massively parallel solver for low Mach number astrophysical flows. The underlying low Mach number equation set allows for efficient, long-time integration for highly subsonic flows compared to compressible approaches. MAESTROeX is suitable for modeling full spherical stars as well as well as planar simulations of dynamics within localized regions of a star, and can robustly handle several orders of magnitude of density and pressure stratification. Previously, we have described the development of the predecessor of MAESTROeX, called MAESTRO, in a series of papers. Here, we present a new, greatly simplified temporal integration scheme that retains the same order of accuracy as our previous approaches. We also explore the use of alternative spatial mapping of the one-dimensional base state onto the full Cartesian grid. The code leverages the new AMReX software framework for block-structured adaptive mesh refinement (AMR) applications, allowing for scalability to large fractions of leadership-class machines. Using our previous studies on the convective phase of single-degenerate progenitor models of Type Ia supernovae as a guide, we characterize the performance of the code and validate the new algorithmic features. Like MAESTRO, MAESTROeX is fully open source.
△ Less
Submitted 9 August, 2019;
originally announced August 2019.
-
pyro: a framework for hydrodynamics explorations and prototyping
Authors:
Alice Harpole,
Michael Zingale,
Ian Hawke,
Taher Chegini
Abstract:
pyro is a Python-based simulation framework designed for ease of implementation and exploration of hydrodynamics methods. It is built in a object-oriented fashion, allowing for the reuse of the core components and fast prototyping of new methods.
pyro is a Python-based simulation framework designed for ease of implementation and exploration of hydrodynamics methods. It is built in a object-oriented fashion, allowing for the reuse of the core components and fast prototyping of new methods.
△ Less
Submitted 22 February, 2019;
originally announced February 2019.
-
Catching Element Formation In The Act
Authors:
Chris L. Fryer,
Frank Timmes,
Aimee L. Hungerford,
Aaron Couture,
Fred Adams,
Wako Aoki,
Almudena Arcones,
David Arnett,
Katie Auchettl,
Melina Avila,
Carles Badenes,
Eddie Baron,
Andreas Bauswein,
John Beacom,
Jeff Blackmon,
Stephane Blondin,
Peter Bloser,
Steve Boggs,
Alan Boss,
Terri Brandt,
Eduardo Bravo,
Ed Brown,
Peter Brown,
Steve Bruenn. Carl Budtz-Jorgensen,
Eric Burns
, et al. (194 additional authors not shown)
Abstract:
Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray…
▽ More
Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.
△ Less
Submitted 7 February, 2019;
originally announced February 2019.
-
Toward Resolved Simulations of Burning Fronts in Thermonuclear X-ray Bursts
Authors:
M. Zingale,
K. Eiden,
Y. Cavecchi,
A. Harpole,
J. B. Bell,
M. Chang,
I. Hawke,
M. P. Katz,
C. M. Malone,
A. J. Nonaka,
D. E. Willcox,
W. Zhang
Abstract:
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 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.
△ Less
Submitted 8 January, 2019; v1 submitted 10 October, 2018;
originally announced October 2018.
-
Effects of tangential velocity in the reactive relativistic Riemann problem
Authors:
A. Harpole,
I. Hawke
Abstract:
Type I X-ray bursts are thermonuclear burning events which occur on the surfaces of accreting neutron stars. Burning begins in a localised spot in the star's ocean layer before propagating across the entire surface as a deflagration. On the scale of the entire star, the burning front can be thought of as discontinuity. To model this, we investigated the reactive Riemann problem for relativistic de…
▽ More
Type I X-ray bursts are thermonuclear burning events which occur on the surfaces of accreting neutron stars. Burning begins in a localised spot in the star's ocean layer before propagating across the entire surface as a deflagration. On the scale of the entire star, the burning front can be thought of as discontinuity. To model this, we investigated the reactive Riemann problem for relativistic deflagrations and detonations and developed a numerical solver. Unlike for the Newtonian Riemann problem, where only the velocity perpendicular to the interface is relevant, in the relativistic case the tangential velocity becomes coupled through the Lorentz factor and can alter the waves present in the solution. We investigated whether a fast tangential velocity may be able to cause a deflagration wave to transition to a detonation. We found that such a transition is possible, but only for tangential velocities that are a significant fraction of the speed of light or for systems already on the verge of transitioning. Consequently, it is highly unlikely that this transition would occur for a burning front in a neutron star ocean without significant contributions from additional multidimensional effects.
△ Less
Submitted 2 September, 2019; v1 submitted 19 June, 2018;
originally announced June 2018.
