This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.

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 2-D 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 bubble's 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.

Critically understanding the `standard candle'-like behavior of Type Ia supernovae requires understanding their explosion mechanism. One family of models for Type Ia Supernovae begins with a deflagration in a Carbon-Oxygen white dwarf which greatly accelerates through wrinkling and flame instabilities. While the planar speed and behavior of astrophysically-relevant flames is increasingly well understood, more complex behavior, such as the flame's response to stretch and curvature, has not been extensively explored in the astrophysical literature; this behavior can greatly enhance or suppress instabilities and local flame-wrinkling, which in turn can increase or decrease the bulk burning rate. In this paper, we explore the effects of curvature on both nuclear flames and simpler model flames to understand the effect of curvature on the flame structure and speed.

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.