skip to main content


Title: octo-tiger : a new, 3D hydrodynamic code for stellar mergers that uses hpx parallelization
ABSTRACT octo-tiger is an astrophysics code to simulate the evolution of self-gravitating and rotating systems of arbitrary geometry based on the fast multipole method, using adaptive mesh refinement. octo-tiger is currently optimized to simulate the merger of well-resolved stars that can be approximated by barotropic structures, such as white dwarfs (WDs) or main-sequence stars. The gravity solver conserves angular momentum to machine precision, thanks to a ‘correction’ algorithm. This code uses hpx parallelization, allowing the overlap of work and communication and leading to excellent scaling properties, allowing for the computation of large problems in reasonable wall-clock times. In this paper, we investigate the code performance and precision by running benchmarking tests. These include simple problems, such as the Sod shock tube, as well as sophisticated, full, WD binary simulations. Results are compared to analytical solutions, when known, and to other grid-based codes such as flash. We also compute the interaction between two WDs from the early mass transfer through to the merger and compare with past simulations of similar systems. We measure octo-tiger’s scaling properties up to a core count of ∼80 000, showing excellent performance for large problems. Finally, we outline the current and planned areas of development aimed at tackling a number of physical phenomena connected to observations of transients.  more » « less
Award ID(s):
1814967
NSF-PAR ID:
10293152
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
504
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
5345 to 5382
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    As the earliest relics of star formation episodes of the Universe, the most massive galaxies are the key to our understanding of the stellar population, cosmic structure, and supermassive black hole (SMBH) evolution. However, the details of their formation histories remain uncertain. We address these problems by planning a large survey sample of 101 ultramassive galaxies (z ≤ 0.3, |δ + 24°| < 45°, |b| > 8°), including 76 per cent ellipticals, 17 per cent lenticulars, and 7 per cent spirals brighter than MK ≤ −27 mag (stellar mass 2 × 1012 ≲ M⋆ ≲ 5 × 1012 M⊙) with ELT/HARMONI. Our sample comprises diverse galaxy environments ranging from isolated to dense-cluster galaxies. The primary goals of the project are to (1) explore the stellar dynamics inside galaxy nuclei and weigh SMBHs, (2) constrain the black hole scaling relations at the highest mass, and (3) probe the late-time assembly of these most massive galaxies through the stellar population and kinematical gradients. We describe the survey, discuss the distinct demographics and environmental properties of the sample, and simulate their HARMONI Iz-, Iz + J-, and H + K-band observations by combining the inferred stellar-mass models from Pan-STARRS observations, an assumed synthetic spectrum of stars, and SMBHs with masses estimated based on different black hole scaling relations. Our simulations produce excellent state-of-the-art integral field spectrography and stellar kinematics (ΔVrms ≲ 1.5 per cent) in a relatively short exposure time. We use these stellar kinematics in combination with the Jeans anisotropic model to reconstruct the SMBH mass and its error using a Markov chain Monte Carlo simulation. Thus, these simulations and modellings can be benchmarks to evaluate the instrument models and pipelines dedicated to HARMONI to exploit the unprecedented capabilities of ELT.

     
    more » « less
  2. ABSTRACT

    The inner structure of core helium burning (CHeB) stars remains uncertain due to the yet unknown nature of mixing at the boundary of their cores. Large convective cores beyond a bare Schwarzschild model are favoured both from theoretical arguments and from asteroseismological constraints. However, the exact nature of this extra mixing, and in particular the possible presence of semiconvective layers, is still debated. In this work, we approach this problem through a new avenue by performing the first full-sphere 3D hydrodynamics simulations of the interiors of CHeB stars. We use the ppmstar explicit gas dynamics code to simulate the inner 0.45$\, {\rm M}_{\odot }$ of a 3 M⊙ CHeB star. Simulations are performed using different Cartesian grid resolutions (7683, 11523, and 17283) and heating rates. We use two different initial states, one based on mesas's predictive mixing scheme (which significantly extends the core beyond the Schwarzschild boundary) and one based on the convective premixing approach (which exhibits a semiconvective interface). The general behaviour of the flow in the convective core and in the stable envelope (where internal gravity waves are observed) is consistent with our recent simulations of core convection in massive main-sequence stars, and so are the various luminosity scaling relations. The semiconvective layers are dominated by strong internal gravity waves that do not produce measurable species mixing, but overshooting motions from the convective core gradually homogenize the semiconvective interface. This process can possibly completely erase the semiconvective layers, which would imply that CHeB stars do not harbour a semiconvection zone.

     
    more » « less
  3. Abstract

    White dwarfs (WDs) offer unrealized potential in solving two problems in astrophysics: stellar age accuracy and precision. WD cooling ages can be inferred from surface temperatures and radii, which can be constrained with precision by high-quality photometry and parallaxes. Accurate and precise Gaia parallaxes along with photometric surveys provide information to derive cooling and total ages for vast numbers of WDs. Here we analyze 1372 WDs found in wide binaries with main-sequence (MS) companions and report on the cooling and total age precision attainable in these WD+MS systems. The total age of a WD can be further constrained if its original metallicity is known because the MS lifetime depends on metallicity at fixed mass, yet metallicity is unavailable via spectroscopy of the WD. We show that incorporating spectroscopic metallicity constraints from 38 wide binary MS companions substantially decreases internal uncertainties in WD total ages compared to a uniform constraint. Averaged over the 38 stars in our sample, the total (internal) age uncertainty improves from 21.04% to 16.77% when incorporating the spectroscopic constraint. Higher mass WDs yield better total age precision; for eight WDs with zero-age MS masses ≥2.0M, the mean uncertainty in total ages improves from 8.61% to 4.54% when incorporating spectroscopic metallicities. We find that it is often possible to achieve 5% total age precision for WDs with progenitor masses above 2.0Mif parallaxes with ≤1% precision and Pan-STARRSg,r, andiphotometry with ≤0.01 mag precision are available.

     
    more » « less
  4. Octo-Tiger is a code for modeling three-dimensional self-gravitating astrophysical fluids. It was particularly designed for the study of dynamical mass transfer between interacting binary stars. Octo-Tiger is parallelized for distributed systems using the asynchronous many-task runtime system, the C++ standard library for parallelism and concurrency (HPX) and utilizes CUDA for its gravity solver. Recently, we have remodeled Octo-Tiger’s hydro solver to use a three-dimensional reconstruction scheme. In addition, we have ported the hydro solver to GPU using CUDA kernels. We present scaling results for the new hydro kernels on ORNL’s Summit machine using a Sedov-Taylor blast wave problem. We also compare Octo-Tiger’s new hydro scheme with its old hydro scheme, using a rotating star as a test problem. 
    more » « less
  5. Contagion dynamics on networks are used to study many problems, including disease and virus epidemics, incarceration, obesity, protests and rebellions, needle sharing in drug use, and hurricane and other natural disaster events. Simulators to study these problems range from smaller-scale serial codes to large-scale distributed systems. In recent years, Python-based simulation systems have been built. In this work, we describe a new Python-based agent-based simulator called CSonNet. It differs from codes such as Epidemics on Networks in that it performs discrete time simulations based on the graph dynamical systems formalism. CSonNet is a parallel code; it implements concurrency through an embarrassingly parallel approach of running multiple simulation instances on a user-specified number of forked processes. It has a modeling framework whereby agent models are composed using a set of pre-defined state transition rules. We provide strong-scaling performance results and case studies to illustrate its features. 
    more » « less