BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

A. Almgren

A. Almgren contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.SRastro-ph.IMastro-ph.COphysics.comp-phhep-exhep-lathep-phhep-th

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2016arXiv

ASCR/HEP Exascale Requirements Review Report

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Carbon Deflagration in Type Ia Supernova: I. Centrally Ignited Models

A leading model for Type Ia supernovae (SNe Ia) begins with a white dwarf near the Chandrasekhar mass that ignites a degenerate thermonuclear runaway close to its center and explodes. In a series of papers, we shall explore the consequences of ignition at several locations within such dwarfs. Here we assume central ignition, which has been explored before, however, the problem is worth revisiting, if only to validate those previous studies and to further elucidate the relevant physics for future work. A perturbed sphere of hot iron ash with a radius of ~100 km is initialized at the middle of the star. The subsequent explosion is followed in several simulations using a thickened flame model in which the flame speed is either fixed --- within the range expected from turbulent combustion --- or based on the local turbulent intensity. Global results, including the explosion energy and bulk nucleosynthesis (e.g. 56Ni of 0.48--0.56 $\Msun$) turn out to be insensitive to this speed. In all completed runs, the energy released by the nuclear burning is adequate to unbind the star, but not enough to give the energy and brightness of typical SNe Ia. As found previously, the chemical stratification observed in typical events is not reproduced. These models produce a large amount of unburned carbon and oxygen in central low velocity regions, which is inconsistent with spectroscopic observations, and the intermediate mass elements and iron group elements are strongly mixed during the explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

CASTRO: A New Compressible Astrophysical Solver. III. Multigroup Radiation Hydrodynamics

We present a formulation for multigroup radiation hydrodynamics that is correct to order $O(v/c)$ using the comoving-frame approach and the flux-limited diffusion approximation. We describe a numerical algorithm for solving the system, implemented in the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically-rectangular variable-sized grids with simultaneous refinement in both space and time. In our multigroup radiation solver, the system is split into three parts, one part that couples the radiation and fluid in a hyperbolic subsystem, another part that advects the radiation in frequency space, and a parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem and the frequency space advection are solved explicitly with high-order Godunov schemes, whereas the parabolic part is solved implicitly with a first-order backward Euler method. Our multigroup radiation solver works for both neutrino and photon radiation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Hydrodynamic Origin of Neutron Star Kicks

We present results from a suite of axisymmetric, core-collapse supernova simulations in which hydrodynamic recoil from an asymmetric explosion produces large proto-neutron star (PNS) velocities. We use the adaptive-mesh refinement code CASTRO to self-consistently follow core-collapse, the formation of the PNS and its subsequent acceleration. We obtain recoil velocities of up to 620 km/s at ~1 s after bounce. These velocities are consistent with the observed distribution of pulsar kicks and with PNS velocities obtained in other theoretical calculations. Our PNSs are still accelerating at several hundred km/s at the end of our calculations, suggesting that even the highest velocity pulsars may be explained by hydrodynamic recoil in generic, core-collapse supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

CASTRO: A New Compressible Astrophysical Solver. II. Gray Radiation Hydrodynamics

We describe the development of a flux-limited gray radiation solver for the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of logically-rectangular variable-sized grids with simultaneous refinement in both space and time. The gray radiation solver is based on a mixed-frame formulation of radiation hydrodynamics. In our approach, the system is split into two parts, one part that couples the radiation and fluid in a hyperbolic subsystem, and another parabolic part that evolves radiation diffusion and source-sink terms. The hyperbolic subsystem is solved explicitly with a high-order Godunov scheme, whereas the parabolic part is solved implicitly with a first-order backward Euler method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Induced Rotation in 3D Simulations of Core Collapse Supernovae: Implications for Pulsar Spins

It has been suggested that the observed rotation periods of radio pulsars might be induced by a non-axisymmetric spiral-mode instability in the turbulent region behind the stalled supernova bounce shock, even if the progenitor core was not initially rotating. In this paper, using the three-dimensional AMR code CASTRO with a realistic progenitor and equation of state and a simple neutrino heating and cooling scheme, we present a numerical study of the evolution in 3D of the rotational profile of a supernova core from collapse, through bounce and shock stagnation, to delayed explosion. By the end of our simulation ($\sim$420 ms after core bounce), we do not witness significant spin up of the proto-neutron star core left behind. However, we do see the development before explosion of strong differential rotation in the turbulent gain region between the core and stalled shock. Shells in this region acquire high spin rates that reach $\sim$$150\,$ Hz, but this region contains too little mass and angular momentum to translate, even if left behind, into rapid rotation for the full neutron star. We find also that much of the induced angular momentum is likely to be ejected in the explosion, and moreover that even if the optimal amount of induced angular momentum is retained in the core, the resulting spin period is likely to be quite modest. Nevertheless, induced periods of seconds are possible.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

3D simulations of Rayleigh-Taylor mixing in core-collapse SNe with CASTRO

We present multidimensional simulations of the post-explosion hydrodynamics in three different 15 solar mass supernova models with zero, 10^{-4} solar metallicity, and solar metallicities. We follow the growth of the Rayleigh-Taylor instability that mixes together the stellar layers in the wake of the explosion. Models are initialized with spherically symmetric explosions and perturbations are seeded by the grid. Calculations are performed in two-dimensional axisymmetric and three-dimensional Cartesian coordinates using the new Eulerian hydrodynamics code, CASTRO. We find as in previous work, that Rayleigh-Taylor perturbations initially grow faster in 3D than in 2D. As the Rayleigh-Taylor fingers interact with one another, mixing proceeds to a greater degree in 3D than in 2D, reducing the local Atwood number and slowing the growth rate of the instability in 3D relative to 2D. By the time mixing has stopped, the width of the mixed region is similar in 2D and 3D simulations provided the Rayleigh-Taylor fingers show significant interaction. Our results imply that 2D simulations of light curves and nucleosynthesis in supernovae (SNe) that die as red giants may capture the features of an initially spherically symmetric explosion in far less computational time than required by a full 3D simulation. However, capturing large departures from spherical symmetry requires a significantly perturbed explosion. Large scale asymmetries cannot develop through an inverse cascade of merging Rayleigh-Taylor structures; they must arise from asymmetries in the initial explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dimension as a Key to the Neutrino Mechanism of Core-Collapse Supernova Explosions

We explore the dependence on spatial dimension of the viability of the neutrino heating mechanism of core-collapse supernova explosions. We find that the tendency to explode is a monotonically increasing function of dimension, with 3D requiring $\sim$40$-$50\% lower driving neutrino luminosity than 1D and $\sim$15$-$25\% lower driving neutrino luminosity than 2D. Moreover, we find that the delay to explosion for a given neutrino luminosity is always shorter in 3D than 2D, sometimes by many hundreds of milliseconds. The magnitude of this dimensional effect is much larger than the purported magnitude of a variety of other effects, such as nuclear burning, inelastic scattering, or general relativity, which are sometimes invoked to bridge the gap between the current ambiguous and uncertain theoretical situation and the fact of robust supernova explosions. Since real supernovae occur in three dimensions, our finding may be an important step towards unraveling one of the most problematic puzzles in stellar astrophysics. In addition, even though in 3D we do see pre-explosion instabilities and blast asymmetries, unlike the situation in 2D, we do not see an obvious axially-symmetric dipolar shock oscillation. Rather, the free energy available to power instabilites seems to be shared by more and more degrees of freedom as the dimension increases. Hence, the strong dipolar axisymmetry seen in 2D and previously identified as a fundamental characteristic of the shock hydrodynamics may not survive in 3D as a prominent feature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

MAESTRO, CASTRO, and SEDONA -- Petascale Codes for Astrophysical Applications

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Nucleosynthetic Imprint of 15-40 Solar Mass Primordial Supernovae on Metal-Poor Stars

The inclusion of rotationally-induced mixing in stellar evolution can alter the structure and composition of presupernova stars. We survey the effects of progenitor rotation on nucleosynthetic yields in Population III and II supernovae using the new adaptive mesh refinement (AMR) code CASTRO. We examine spherical explosions in 15, 25 and 40 solar mass stars at Z = 0 and 10^-4 solar metallicity with three explosion energies and two rotation rates. Rotation in the Z = 0 models resulted in primary nitrogen production and a stronger hydrogen burning shell which led all models to die as red supergiants. On the other hand, the Z=10^-4 solar metallicity models that included rotation ended their lives as compact blue stars. Because of their extended structure, the hydrodynamics favors more mixing and less fallback in the metal free stars than the Z = 10^-4 models. As expected, higher energy explosions produce more enrichment and less fallback than do lower energy explosions, and less massive stars produce more enrichment and leave behind smaller remnants than do more massive stars. We compare our nucleosynthetic yields to the chemical abundances in the three most iron-poor stars yet found and reproduce the abundance pattern of one, HE 0557-4840, with a zero metallicity 15 solar mass, 2.4 x 10^51 erg supernova. A Salpeter IMF averaged integration of our yields for Z=0 models with explosion energies of 2.4x10^51 ergs or less is in good agreement with the abundances observed in larger samples of extremely metal-poor stars, provided 15 solar mass stars are included. Since the abundance patterns of extremely metal-poor stars likely arise from a representative sample of progenitors, our yields suggest that low-mass supernovae contributed the bulk of the metals to the early universe.

0 citations0 reviewsNo code linkedOpen access