Researcher profile

A. S. Almgren

A. S. Almgren contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
4topics
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

12 published item(s)

preprint2015arXiv

Comparisons of Two- and Three-Dimensional Convection in Type I X-ray Bursts

We perform the first detailed three-dimensional simulation of low Mach number convection preceding thermonuclear ignition in a mixed H/He X-ray burst. Our simulations include a moderate-sized, approximate network that captures hydrogen and helium burning up through rp-process breakout. We look in detail at the difference between two- and three-dimensional convective fields, including the details of the turbulent convection.

preprint2014arXiv

Cosmological Fluid Mechanics with Adaptively Refined Large Eddy Simulations

We investigate turbulence generated by cosmological structure formation by means of large eddy simulations using adaptive mesh refinement. In contrast to the widely used implicit large eddy simulations, which resolve a limited range of length scales and treat the effect of turbulent velocity fluctuations below the grid scale solely by numerical dissipation, we apply a subgrid-scale model for the numerically unresolved fraction of the turbulence energy. For simulations with adaptive mesh refinement, we utilize a new methodology that allows us to adjust the scale-dependent energy variables in such a way that the sum of resolved and unresolved energies is globally conserved. We test our approach in simulations of randomly forced turbulence, a gravitationally bound cloud in a wind, and the Santa Barbara cluster. To treat inhomogeneous turbulence, we introduce an adaptive Kalman filtering technique that separates turbulent velocity fluctuations on resolved length scales from the non-turbulent bulk flow. From the magnitude of the fluctuating component and the subgrid-scale turbulence energy, a total turbulent velocity dispersion of several 100 km/s is obtained for the Santa Barbara cluster, while the low-density gas outside the accretion shocks is nearly devoid of turbulence. The energy flux through the turbulent cascade and the dissipation rate predicted by the subgrid-scale model correspond to dynamical time scales around 5 Gyr, independent of numerical resolution.

preprint2014arXiv

Influence of adaptive mesh refinement and the hydro solver on shear-induced mass stripping in a minor-merger scenario

We compare two different codes for simulations of cosmological structure formation to investigate the sensitivity of hydrodynamical instabilities to numerics, in particular, the hydro solver and the application of adaptive mesh refinement (AMR). As a simple test problem, we consider an initially spherical gas cloud in a wind, which is an idealized model for the merger of a subcluster or galaxy with a big cluster. Based on an entropy criterion, we calculate the mass stripping from the subcluster as a function of time. Moreover, the turbulent velocity field is analyzed with a multi-scale filtering technique. We find remarkable differences between the commonly used PPM solver with directional splitting in the Enzo code and an unsplit variant of PPM in the Nyx code, which demonstrates that different codes can converge to systematically different solutions even when using uniform grids. For the test case of an unbound cloud, AMR simulations reproduce uniform-grid results for the mass stripping quite well, although the flow realizations can differ substantially. If the cloud is bound by a static gravitational potential, however, we find strong sensitivity to spurious fluctuations which are induced at the cutoff radius of the potential and amplified by the bow shock. This gives rise to substantial deviations between uniform-grid and AMR runs performed with Enzo, while the mass stripping in Nyx simulations of the subcluster is nearly independent of numerical resolution and AMR. Although many factors related to numerics are involved, our study indicates that unsplit solvers with advanced flux limiters help to reduce grid effects and to keep numerical noise under control, which is important for hydrodynamical instabilities and turbulent flows.

preprint2014arXiv

Multidimensional Modeling of Type I X-ray Bursts. II. Two-Dimensional Convection in a Mixed H/He Accretor

Type I X-ray Bursts (XRBs) are thermonuclear explosions of accreted material on the surfaces of a neutron stars in low mass X-ray binaries. Prior to the ignition of a subsonic burning front, runaway burning at the base of the accreted layer drives convection that mixes fuel and heavy-element ashes. In this second paper in a series, we explore the behavior of this low Mach number convection in mixed hydrogen/helium layers on the surface of a neutron star using two-dimensional simulations with the Maestro code. Maestro takes advantage of the highly subsonic flow field by filtering dynamically unimportant sound waves while retaining local compressibility effects, such as those due to stratification and energy release from nuclear reactions. In these preliminary calculations, we find that the rp-process approximate network creates a convective region that is split into two layers. While this splitting appears artificial due to the approximations of the network regarding nuclear flow out of the breakout reaction 18Ne(a,p)21Na, these calculations hint at further simplifications and improvements of the burning treatment for use in subsequent calculations in three dimensions for a future paper.

preprint2013arXiv

The Deflagration Stage of Chandrasekhar Mass Models For Type Ia Supernovae: I. Early Evolution

We present high-resolution, full-star simulations of the post-ignition phase of Type Ia supernovae using the compressible hydrodynamics code Castro. Initial conditions, including the turbulent velocity field and ignition site, are imported directly from a simulation of the last few hours of presupernova convection using a low Mach number code, Maestro. Adaptive mesh refinement allows the initial burning front to be modeled with an effective resolution of 36,864^3 zones (~136 m/zone). The initial rise and expansion of the deflagration front are tracked until burning reaches the star's edge and the role of the background turbulence on the flame is investigated. The effect of artificially moving the ignition location closer to the star's center is explored. The degree to which turbulence affects the burning front decreases with increasing ignition radius since the buoyancy force is stronger at larger radii. Even central ignition --- in the presence of a background convective flow field --- is rapidly carried off-center as the flame is carried by the flow field. We compare our results to analytic models for burning thermals, and find that they reproduce the general trends of the bubble's size and mass, but underpredict the amount of buoyant acceleration due to simplifying assumptions of the bubble's properties. Overall, we find that the amount of mass that burns prior to flame break out is small, consistent with a "gravitationally confined detonation" occurring at a later epoch, but additional burning will occur following breakout that may modify this conclusion.

preprint2012arXiv

Low Mach Number Modeling of Convection in Helium Shells on Sub-Chandrasekhar White Dwarfs. I. Methodology

We assess the robustness of a low Mach number hydrodynamics algorithm for modeling helium shell convection on the surface of a white dwarf in the context of the sub-Chandrasekhar model for Type Ia supernovae. We use the low Mach number stellar hydrodynamics code, MAESTRO, to perform three-dimensional, spatially-adaptive simulations of convection leading up to the point of the ignition of a burning front. We show that the low Mach number hydrodynamics model provides a robust description of the system.

preprint2011arXiv

From Convection to Explosion: End-to-End Simulation of Type Ia Supernovae

We present our end-to-end capability for computing the convective phase through the explosion phase of Type Ia supernovae. We compute the convective phase up to the time of ignition using our low Mach number code, MAESTRO, and the subsequent explosion phase using our compressible code, CASTRO. Both codes share the same BoxLib software framework and use finite-volume, block-structured adaptive mesh refinement (AMR) to enable high-resolution, three-dimensional full-star simulations that scale to 100,000+ cores. We present preliminary results from the first-ever simulations of convection preceding ignition using MAESTRO with AMR. We also demonstrate our ability to initialize a compressible simulation of the explosion phase in CASTRO using data obtained directly from MAESTRO just before ignition. Some care must be taken during this initialization procedure when interpreting the size and distribution of hot spots.

preprint2011arXiv

High-Resolution Simulations of Convection Preceding Ignition in Type Ia Supernovae Using Adaptive Mesh Refinement

We extend our previous three-dimensional, full-star simulations of the final hours of convection preceding ignition in Type Ia supernovae to higher resolution using the adaptive mesh refinement capability of our low Mach number code, MAESTRO. We report the statistics of the ignition of the first flame at an effective 4.34 km resolution, and general flow field properties at an effective 2.17 km resolution. We find that off-center ignition is likely, with radius of 50 km most favored and a likely range of 40 to 75 km. This is consistent with our previous coarser (8.68 km resolution) simulations, implying that we have achieved sufficient resolution in our determination of likely ignition radii. The dynamics of the last few hot spots preceding ignition suggest that a multiple ignition scenario is not likely. With improved resolution, we can more clearly see the general flow pattern in the convective region, characterized by a strong outward plume with a lower speed recirculation. We show that the convective core is turbulent with a Kolmogorov spectrum and has a lower turbulent intensity and larger integral length scale than previously thought (on the order of 16 km s$^{-1}$ and 200 km, respectively), and we discuss the potential consequences for the first flames.

preprint2010arXiv

CASTRO: A New Compressible Astrophysical Solver. I. Hydrodynamics and Self-Gravity

We present a new code, CASTRO, that solves the multicomponent compressible hydrodynamic equations for astrophysical flows including self-gravity, nuclear reactions and radiation. CASTRO uses an Eulerian grid and incorporates adaptive mesh refinement (AMR). Our approach to AMR uses a nested hierarchy of logically-rectangular grids with simultaneous refinement in both space and time. The radiation component of CASTRO will be described in detail in the next paper, Part II, of this series.

preprint2010arXiv

MAESTRO: An Adaptive Low Mach Number Hydrodynamics Algorithm for Stellar Flows

Many astrophysical phenomena are highly subsonic, requiring specialized numerical methods suitable for long-time integration. In a series of earlier papers we described the development of MAESTRO, a low Mach number stellar hydrodynamics code that can be used to simulate long-time, low-speed flows that would be prohibitively expensive to model using traditional compressible codes. MAESTRO is based on an equation set derived using low Mach number asymptotics; this equation set does not explicitly track acoustic waves and thus allows a significant increase in the time step. MAESTRO is suitable for two- and three-dimensional local atmospheric flows as well as three-dimensional full-star flows. Here, we continue the development of MAESTRO by incorporating adaptive mesh refinement (AMR). The primary difference between MAESTRO and other structured grid AMR approaches for incompressible and low Mach number flows is the presence of the time-dependent base state, whose evolution is coupled to the evolution of the full solution. We also describe how to incorporate the expansion of the base state for full-star flows, which involves a novel mapping technique between the one-dimensional base state and the Cartesian grid, as well as a number of overall improvements to the algorithm. We examine the efficiency and accuracy of our adaptive code, and demonstrate that it is suitable for further study of our initial scientific application, the convective phase of Type Ia supernovae.

preprint2010arXiv

Multidimensional Modeling of Type I X-ray Bursts. I. Two-Dimensional Convection Prior to the Outburst of a Pure Helium Accretor

We present multidimensional simulations of the early convective phase preceding ignition in a Type I X-ray burst using the low Mach number hydrodynamics code, MAESTRO. A low Mach number approach is necessary in order to perform long-time integration required to study such phenomena. Using MAESTRO, we are able to capture the expansion of the atmosphere due to large-scale heating while capturing local compressibility effects such as those due to reactions and thermal diffusion. We also discuss the preparation of one-dimensional initial models and the subsequent mapping into our multidimensional framework. Our method of initial model generation differs from that used in previous multidimensional studies, which evolved a system through multiple bursts in one dimension before mapping onto a multidimensional grid. In our multidimensional simulations, we find that the resolution necessary to properly resolve the burning layer is an order of magnitude greater than that used in the earlier studies mentioned above. We characterize the convective patterns that form and discuss their resulting influence on the state of the convective region, which is important in modeling the outburst itself.

preprint2009arXiv

Low Mach Number Modeling of Type Ia Supernovae. IV. White Dwarf Convection

We present the first three-dimensional, full-star simulations of convection in a white dwarf preceding a Type Ia supernova, specifically the last few hours before ignition. For these long-time calculations we use our low Mach number hydrodynamics code, MAESTRO, which we have further developed to treat spherical stars centered in a three-dimensional Cartesian geometry. The main change required is a procedure to map the one-dimensional radial base state to and from the Cartesian grid. Our models recover the dipole structure of the flow seen in previous calculations, but our long-time integration shows that the orientation of the dipole changes with time. Furthermore, we show the development of gravity waves in the outer, stable portion of the star. Finally, we evolve several calculations to the point of ignition and discuss the range of ignition radii.