-
Large-scale ordered magnetic fields generated in mergers of helium white dwarfs
Authors:
Rüdiger Pakmor,
Ingrid Pelisoli,
Stephen Justham,
Abinaya S. Rajamuthukumar,
Friedrich K. Röpke,
Fabian R. N. Schneider,
Selma E. de Mink,
Sebastian T. Ohlmann,
Philipp Podsiadlowski,
Javier Moran Fraile,
Marco Vetter,
Robert Andrassy
Abstract:
Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1\%$.
We present two high-resolution magnetohydrod…
▽ More
Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1\%$.
We present two high-resolution magnetohydrodynamical (MHD) simulations of the merger of two helium white dwarfs in a binary system with the same total mass of $0.6\,M_\odot$. We analyse one equal-mass merger with two $0.3\,M_\odot$ white dwarfs, and one unequal-mass merger with a $0.25\,M_\odot$ white dwarf and a $0.35\,M_\odot$ white dwarf. We simulate the inspiral, merger, and further evolution of the merger remnant for about $50$ rotations.
We find efficient magnetic field amplification in both mergers via a small-scale dynamo, reproducing previous results of stellar merger simulations. The magnetic field saturates at similar strength for both simulations.
We then identify a second phase of magnetic field amplification in both merger remnants that happens on a timescale of several tens of rotational periods of the merger remnant. This phase generates a large-scale ordered azimuthal field. We identify it as a large-scale dynamo driven by the magneto-rotational instability (MRI).
Finally, we suggest that in the unequal-mass merger remnant, helium burning will eventually start in a shell around a cold core. The convection zone this generates will coincide with the region that contains most of the magnetic energy, probably erasing the strong, ordered field. The equal-mass merger remnant instead will probably ignite burning in the center, retaining its ordered field. Therefore, the mass ratio of the initial merger could be the selecting factor that decides if a merger remnant will stay highly magnetised long after the merger.
△ Less
Submitted 2 July, 2024;
originally announced July 2024.
-
Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows
Authors:
G. Leidi,
R. Andrassy,
W. Barsukow,
J. Higl,
P. V. F. Edelmann,
F. K. Röpke
Abstract:
High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so t…
▽ More
High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. In this work, we use our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We consider Mach numbers in the range from $10^{-3}$ to $10^{-1}$ in a well-posed, 2D, Kelvin--Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. We find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a "low-dissipation" Riemann solver and a sextic reconstruction scheme, (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same, (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate, (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow, (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.
△ Less
Submitted 26 February, 2024;
originally announced February 2024.
-
Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell
Authors:
G. Leidi,
R. Andrassy,
J. Higl,
P. V. F. Edelmann,
F. K. Röpke
Abstract:
Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dyn…
▽ More
Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible Seven-League Hydro code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of $25\ M_\odot$. We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as $0.33$, and the average magnetic field strength is ${\sim}10^{10}\,\mathrm{G}$. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to $60\%$ of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ${\approx}\,20\%$ as compared to the purely hydrodynamic case.
△ Less
Submitted 29 September, 2023;
originally announced September 2023.
-
Towards a self-consistent model of the convective core boundary in upper-main-sequence stars
Authors:
R. Andrassy,
G. Leidi,
J. Higl,
P. V. F. Edelmann,
F. R. N. Schneider,
F. K. Roepke
Abstract:
There is strong observational evidence that convective cores of intermediate-mass and massive main-sequence stars are substantially larger than standard stellar-evolution models predict. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-time-scale process th…
▽ More
There is strong observational evidence that convective cores of intermediate-mass and massive main-sequence stars are substantially larger than standard stellar-evolution models predict. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-time-scale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number Seven-League Hydro (SLH) code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a $15$ M$_\odot$ zero-age main-sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor ranging from $10^3$ to $10^6$. After reaching thermal equilibrium, all of our models show a clear penetration layer. Its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. As the luminosity is decreased, the penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification. This structure corresponds to the adiabatic ,,step overshoot'' model often employed in stellar-evolution calculations. The thickness of the penetration layer slowly decreases with decreasing luminosity. Depending on how we extrapolate our 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from $0.09$ to $0.44$ pressure scale heights, which are broadly compatible with observations.
△ Less
Submitted 8 July, 2023;
originally announced July 2023.
-
3D hydrodynamic simulations of massive main-sequence stars. I. Dynamics and mixing of convection and internal gravity waves
Authors:
Falk Herwig,
Paul R. Woodward,
Huaqing Mao,
William R. Thompson,
Pavel Denissenkov,
Josh Lau,
Simon Blouin,
Robert Andrassy,
Adam Paul
Abstract:
We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications.…
▽ More
We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications. As in previous works, the turbulent mass entrainment rate extrapolated to nominal heating is unrealistically high ($1.58\times 10^{-4}\ \mathrm{M_\odot/yr}$), which is discussed in terms of the non-equilibrium response of the simulations to the initial stratification. We measure quantitatively the effect of mixing due to internal gravity waves excited by core convection interacting with the boundary in our simulations. The wave power spectral density as a function of frequency and wavelength agrees well with the GYRE eigenmode predictions based on the 1D spherically averaged radial profile. A diffusion coefficient profile that reproduces the spherically averaged abundance distribution evolution is determined for each simulation. Through a combination of eigenmode analysis and scaling relations it is shown that in the $N^2$-peak region, mixing is due to internal gravity waves and follows the scaling relation $D_\mathrm{IGW-hydro} \propto L^{4/3}$ over a $\gtrapprox 2\ \mathrm{dex}$ range of heating factors. Different extrapolations of the mixing efficiency down to nominal heating are discussed. If internal gravity wave mixing is due to thermally-enhanced shear mixing, an upper limit is $D_\mathrm{IGW} \lessapprox 2$ to $3\times 10^4\ \mathrm{cm^2/s}$ at nominal heating in the $N^2$-peak region above the convective core.
△ Less
Submitted 19 July, 2023; v1 submitted 9 March, 2023;
originally announced March 2023.
-
A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors
Authors:
G. Leidi,
C. Birke,
R. Andrassy,
J. Higl,
P. V. F. Edelmann,
G. Wiest,
C. Klingenberg,
F. K. Röpke
Abstract:
Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers (M_son < 0.01 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to b…
▽ More
Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers (M_son < 0.01 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant-Friedrichs-Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume Seven-League Hydro (SLH) code, and it makes use of a low-Mach version of the five-wave Harten-Lax-van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit-explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at M_son ~ 0.001 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids.
△ Less
Submitted 4 October, 2022;
originally announced October 2022.
-
Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
Authors:
R. Andrassy,
J. Higl,
H. Mao,
M. Mocák,
D. G. Vlaykov,
W. D. Arnett,
I. Baraffe,
S. W. Campbell,
T. Constantino,
P. V. F. Edelmann,
T. Goffrey,
T. Guillet,
F. Herwig,
R. Hirschi,
L. Horst,
G. Leidi,
C. Meakin,
J. Pratt,
F. Rizzuti,
F. K. Roepke,
P. Woodward
Abstract:
Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility o…
▽ More
Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.
△ Less
Submitted 26 January, 2022; v1 submitted 1 November, 2021;
originally announced November 2021.
-
Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star
Authors:
L. Horst,
R. Hirschi,
P. V. F. Edelmann,
R. Andrassy,
F. K. Roepke
Abstract:
Context. Multidimensional hydrodynamic simulations of convection in stellar interiors are numerically challenging, especially for flows at low Mach numbers.
Methods. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional si…
▽ More
Context. Multidimensional hydrodynamic simulations of convection in stellar interiors are numerically challenging, especially for flows at low Mach numbers.
Methods. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M$_\odot$ star model. The results obtained with the low-Mach AUSM$^{+}$-up solver were compared to results when using its non low-Mach variant AUSM$_\mathrm{B}^{+}$-up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from $180 \times 90^2$ to $810 \times 540^2$ cells and the nuclear energy release was boosted by factors of $3 \times 10^3$, $1 \times 10^4$, and $3 \times 10^4$ to study the dependence of the results on these parameters.
Results. The boosted energy input results in convection at Mach numbers in the range of $10^{-2}$ to $10^{-3}$. Standard mixing-length theory (MLT) predicts convective velocities of about $1.6 \times 10^{-4}$ if no boosting is applied. Simulations with AUSM$^{+}$-up show a Kolmogorov-like inertial range in the kinetic energy spectrum that extends further toward smaller scales compared with its non low-Mach variant. The kinetic energy dissipation of the AUSM$^{+}$-up solver already converges at a lower resolution compared to AUSM$^{+}_{\mathrm{B}}$ -up. The extracted entrainment rates at the boundaries of the convection zone are well represented by the bulk Richardson entrainment law and the corresponding fitting parameters are in agreement with published results for carbon shell burning.
△ Less
Submitted 5 July, 2021;
originally announced July 2021.
-
Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers
Authors:
P. V. F. Edelmann,
L. Horst,
J. P. Berberich,
R. Andrassy,
J. Higl,
G. Leidi,
C. Klingenberg,
F. K. Roepke
Abstract:
Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable s…
▽ More
Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable size can lead to spurious flows orders of magnitude faster than the physical flow. Well-balanced numerical schemes can deal with this problem. Three such schemes were applied in the implicit, finite-volume Seven-League Hydro (SLH) code in combination with a low-Mach-number numerical flux function. We compare how the schemes perform in four numerical experiments addressing some of the challenges imposed by typical problems in stellar hydrodynamics. We find that the $α$-$β$ and deviation well-balancing methods can accurately maintain hydrostatic solutions provided that gravitational potential energy is included in the total energy balance. They accurately conserve minuscule entropy fluctuations advected in an isentropic stratification, which enables the methods to reproduce the expected scaling of convective flow speed with the heating rate. The deviation method also substantially increases accuracy of maintaining stationary orbital motions in a Keplerian disk on long timescales. The Cargo-LeRoux method fares substantially worse in our tests, although its simplicity may still offer some merits in certain situations. Overall, we find the well-balanced treatment of gravity in combination with low Mach number flux functions essential to reproducing correct physical solutions to challenging stellar slow-flow problems on affordable collocated grids.
△ Less
Submitted 8 July, 2021; v1 submitted 25 February, 2021;
originally announced February 2021.
-
Fully compressible simulations of waves and core convection in main-sequence stars
Authors:
L. Horst,
P. V. F. Edelmann,
R. Andrassy,
F. K. Roepke,
D. M. Bowman,
C. Aerts,
R. P. Ratnasingam
Abstract:
Context. Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and rather high artificial viscosity and thermal diffusivity are needed for numerical stability. Direct comparison with obse…
▽ More
Context. Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and rather high artificial viscosity and thermal diffusivity are needed for numerical stability. Direct comparison with observations is thus limited. Aims. We explore the capabilities of our compressible multidimensional hydrodynamics code SLH to simulate stellar oscillations. Methods. We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic $3\mathrm{M}_\odot$ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results. Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with observed spectra of massive stars. Conclusion. The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows the coupling of gravity and pressure waves to be studied too.
△ Less
Submitted 14 September, 2020; v1 submitted 4 June, 2020;
originally announced June 2020.
-
3D1D hydro-nucleosynthesis simulations. I. Advective-reactive post-processing method and its application to H ingestion into He-shell flash convection in rapidly accreting white dwarfs
Authors:
David Stephens,
Falk Herwig,
Paul Woodward,
Pavel Denissenkov,
Robert Andrassy,
Huaqing Mao
Abstract:
We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective-reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with $\mathrm{[Fe/H]} = -2.6$, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simula…
▽ More
We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective-reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with $\mathrm{[Fe/H]} = -2.6$, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simulations. A simpler approach uses diffusion coefficients calculated from the same simulations. All 3D simulations include the energy feedback of the $^{12}$C(p,$γ$)$^{13}$N reaction from the H entrainment. Global oscillations of shell H ingestion in two of the RAWD simulations cause bursts of entrainment of H and non-radial hydrodynamic feedback. With the same nuclear network as in the 3D simulations, the 1D advective two-stream model reproduces the rate and location of the H burning within the He shell closely matching the 3D simulation predictions, as well as qualitatively displaying the asymmetry of the $X_{\mathrm{H}}$ profiles between the up- and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the up- and downstream. For example, $^{89}$Kr and $^{90}$Kr with half-lives of 3.18 min and 32.3 s differ by a factor 2-10 in the two streams. In this particular application the diffusion approach provides globally the same abundance distribution as the advective two-stream mixing model. The resulting i-process yields are in excellent agreement with observations of the exemplary CEMP-r/s star CS31062-050.
△ Less
Submitted 17 February, 2021; v1 submitted 29 January, 2020;
originally announced January 2020.
-
Simulating 3-D Stellar Hydrodynamics using PPM and PPB Multifluid Gas Dynamics on CPU and CPU+GPU Nodes
Authors:
Paul R. Woodward,
Pei-Hung Lin,
Huaqing Mao,
Robert Andrassy,
Falk Herwig
Abstract:
The special computational challenges of simulating 3-D hydrodynamics in deep stellar interiors are discussed, and numerical algorithmic responses described. Results of recent simulations carried out at scale on the NSF's Blue Waters machine at the University of Illinois are presented, with a special focus on the computational challenges they address. Prospects for future work using GPU-accelerated…
▽ More
The special computational challenges of simulating 3-D hydrodynamics in deep stellar interiors are discussed, and numerical algorithmic responses described. Results of recent simulations carried out at scale on the NSF's Blue Waters machine at the University of Illinois are presented, with a special focus on the computational challenges they address. Prospects for future work using GPU-accelerated nodes such as those on the DoE's new Summit machine at Oak Ridge National Laboratory are described, with a focus on numerical algorithmic accommodations that we believe will be necessary.
△ Less
Submitted 31 October, 2018;
originally announced October 2018.
-
H-He Shell Interactions and Nucleosynthesis in Massive Population III Stars
Authors:
Ondrea Clarkson,
Falk Herwig,
Robert Andrassy,
Marco Pignatari,
Paul Woodward,
Huaqing Mao
Abstract:
We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M…
▽ More
We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M$_\odot$. Finally, we describe our progress in simulating a Pop III double convection zone in the PPMStar hydrodynamics code.
△ Less
Submitted 29 October, 2018;
originally announced October 2018.
-
The i-process yields of rapidly-accreting white dwarfs from multicycle He-shell flash stellar evolution models with mixing parameterizations from 3D hydrodynamics simulations
Authors:
Pavel Denissenkov,
Falk Herwig,
Paul Woodward,
Robert Andrassy,
Marco Pignatari,
Samuel Jones
Abstract:
We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for $0.7\leq M_\mathrm{RAWD}/M_\odot\leq 0.75$ and [Fe/H] ranging from $0$ to $-2.6$. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash in…
▽ More
We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for $0.7\leq M_\mathrm{RAWD}/M_\odot\leq 0.75$ and [Fe/H] ranging from $0$ to $-2.6$. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash ingests protons from the accreted H-rich surface layer, which results in maximum neutron densities $N_\mathrm{n,max}\approx 10^{13}$-$10^{15}\ \mathrm{cm}^{-3}$. The H-ingestion rate and the convective boundary mixing (CBM) parameter $f_\mathrm{top}$ adopted in the one-dimensional nucleosynthesis and stellar evolution models are constrained through 3D hydrodynamic simulations. The mass ingestion rate and, for the first time, the scaling laws for the CBM parameter $f_\mathrm{top}$ have been determined from 3D hydrodynamic simulations. We confirm our previous result that the high-metallicity RAWDs have a low mass retention efficiency ($η< 10\%$). A new result is that RAWDs with [Fe/H]$< -2$ have $η> 20\%$, therefore their masses may reach the Chandrasekhar limit and they may eventually explode as SNeIa. This result and the good fits of the i-process yields from the metal-poor RAWDs to the observed chemical composition of the CEMP-r/s stars suggest that some of the present-day CEMP-r/s stars could be former distant members of triple systems, orbiting close binary systems with RAWDs that may have later exploded as SNeIa.
△ Less
Submitted 9 July, 2019; v1 submitted 10 September, 2018;
originally announced September 2018.
-
3D hydrodynamic simulations of C ingestion into a convective O shell
Authors:
Robert Andrassy,
Falk Herwig,
Paul Woodward,
Christian Ritter
Abstract:
Interactions between convective shells in evolved massive stars have been linked to supernova impostors, to the production of the odd-Z elements Cl, K, and Sc, and they might also help generate the large-scale asphericities that are known to facilitate shock revival in supernova explosion models. We investigate the process of ingestion of C-shell material into a convective O-burning shell, includi…
▽ More
Interactions between convective shells in evolved massive stars have been linked to supernova impostors, to the production of the odd-Z elements Cl, K, and Sc, and they might also help generate the large-scale asphericities that are known to facilitate shock revival in supernova explosion models. We investigate the process of ingestion of C-shell material into a convective O-burning shell, including the hydrodynamic feedback from the nuclear burning of the ingested material. Our 3D hydrodynamic simulations span almost 3 dex in the total luminosity $L_\mathrm{tot}$. All but one of the simulations reach a quasi-stationary state with the entrainment rate and convective velocity proportional to $L_\mathrm{tot}$ and $L_\mathrm{tot}^{1/3}$, respectively. Carbon burning provides $14\,$--$\,33\%$ of the total luminosity, depending on the set of reactions considered. Equivalent simulations done on $768^3$ and $1152^3$ grids are in excellent quantitative agreement. The flow is dominated by a few large-scale convective cells. An instability leading to large-scale oscillations with Mach numbers in excess of $0.2$ develops in an experimental run with the energy yield from C burning increased by a factor of 10. This run represents most closely the conditions expected in a violent O-C shell merger, which is a potential production site for odd-Z elements such as K and Sc and which may seed asymmetries in the supernova progenitor. 1D simulations may underestimate the energy generation from the burning of ingested material by as much as a factor two owing to their missing the effect of clumpiness of entrained material on the nuclear reaction rate.
△ Less
Submitted 21 October, 2019; v1 submitted 12 August, 2018;
originally announced August 2018.
-
Cyberhubs: Virtual Research Environments for Astronomy
Authors:
Falk Herwig,
Robert Andrassy,
Nic Annau,
Ondrea Clarkson,
Benoit Cote,
Aaron D'Sa,
Sam Jones,
Belaid Moa,
Jericho O'Connell,
David Porter,
Christian Ritter,
Paul Woodward
Abstract:
Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to med…
▽ More
Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to medium-sized collaborations based on the Jupyter and Docker technology, that allows web-browser enabled, remote, interactive analytic access to shared data. It offers an initial step to address these challenges. The features and deployment steps of the system are described, as well as the requirements collection through an account of the different approaches to data structuring, handling and available analytic tools for the NuGrid and PPMstar collaborations. NuGrid is an international collaboration that creates stellar evolution and explosion physics and nucleosynthesis simulation data. The PPMstar collaboration performs large-scale 3D stellar hydrodynamics simulation of interior convection in the late phases of stellar evolution. Examples of science that is presently performed on cyberhubs, in the areas 3D stellar hydrodynamic simulations, stellar evolution and nucleosynthesis and Galactic chemical evolution, are presented.
△ Less
Submitted 6 February, 2018;
originally announced February 2018.
-
Convective-reactive nucleosynthesis of K, Sc, Cl and p-process isotopes in O-C shell mergers
Authors:
C. Ritter,
R. Andrassy,
B. Côté,
F. Herwig,
P. R. Woodward,
M. Pignatari,
S. Jones
Abstract:
We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1…
▽ More
We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1D nucleosynthesis simulations which show that such convective-reactive ingestion events can be a production site for P, Cl, K and Sc. An entrainment rate of $10^{-3}{{\rm M}_\odot}$/s features overproduction factors $OP_\mathrm{s} \approx 7$. Full O-C shell mergers in our 1D stellar evolution massive star models have overproduction factors $OP_\mathrm{m}>1 \mathrm{dex}$ but for such cases 3D hydrodynamic simulations suggest deviations from spherical symmetry. p-process species can be produced with overproduction factors of $OP_\mathrm{m}>1 \mathrm{dex}$, e.g. for $^{130,132}Ba$. Using the uncertain prediction of the $15{{\rm M}_\odot}$, $Z=0.02$ massive star model ($OP_\mathrm{m} \approx 15$) as representative for merger or entrainment convective-reactive events involving O- and C-burning shells, and assume that such events occur in more than 50% of all stars, our chemical evolution models reproduce the observed Galactic trends of the odd-Z elements.
△ Less
Submitted 25 September, 2017; v1 submitted 19 April, 2017;
originally announced April 2017.
-
Idealised hydrodynamic simulations of turbulent oxygen-burning shell convection in 4π geometry
Authors:
Samuel Jones,
Robert Andrassy,
Stou Sandalski,
Austin Davis,
Paul Woodward,
Falk Herwig
Abstract:
This work investigates the properties of convection in stars with particular emphasis on entrainment across the upper convective boundary (CB). Idealised simulations of turbulent convection in the O-burning shell of a massive star are performed in $4π$ geometry on $768^3$ and $1536^3$ grids, driven by a representative heating rate. A heating series is also performed on the $768^3$ grid. The…
▽ More
This work investigates the properties of convection in stars with particular emphasis on entrainment across the upper convective boundary (CB). Idealised simulations of turbulent convection in the O-burning shell of a massive star are performed in $4π$ geometry on $768^3$ and $1536^3$ grids, driven by a representative heating rate. A heating series is also performed on the $768^3$ grid. The $1536^3$ simulation exhibits an entrainment rate at the upper CB of $1.33\times10^{-6}~M_\odot~\mathrm{s}^{-1}$. The $768^3$ simulation with the same heating rate agrees within 17 per cent. The entrainment rate at the upper convective boundary is found to scale linearly with the driving luminosity and with the cube of the shear velocity at the upper boundary, while the radial RMS fluid velocity scales with the cube root of the driving luminosity, as expected. The mixing is analysed in a 1D diffusion framework, resulting in a simple model for CB mixing. The analysis confirms previous findings that limiting the MLT mixing length to the distance to the CB in 1D simulations better represents the spherically-averaged radial velocity profiles from the 3D simulations and provides an improved determination of the reference diffusion coefficient $D_0$ for the exponential diffusion CB mixing model in 1D. From the 3D simulation data we adopt as the convective boundary the location of the maximum gradient in the horizontal velocity component which has $2σ$ spatial fluctuations of $\approx0.17 H_P$ . The exponentially decaying diffusion CB mixing model with $f = 0.03$ reproduces the spherically-averaged 3D abundance profiles.
△ Less
Submitted 2 November, 2016; v1 submitted 12 May, 2016;
originally announced May 2016.
-
Convective settling in main sequence stars: Li and Be depletion
Authors:
R. Andrássy,
H. C. Spruit
Abstract:
The process of convective settling is based on the assumption that a small fraction of the low-entropy downflows sink from the photosphere down to the bottom of the star's envelope convection zone retaining a substantial entropy contrast. We have previously shown that this process could explain the slow Li depletion observed in the Sun. We construct a parametric model of convective settling to inv…
▽ More
The process of convective settling is based on the assumption that a small fraction of the low-entropy downflows sink from the photosphere down to the bottom of the star's envelope convection zone retaining a substantial entropy contrast. We have previously shown that this process could explain the slow Li depletion observed in the Sun. We construct a parametric model of convective settling to investigate the dependence of Li and Be depletion on stellar mass and age. Our model is generally in good agreement with the Li abundances measured in open clusters and solar twins, although it seems to underestimate the Li depletion in the first ~1 Gyr. The model is also compatible with the Be abundances measured in a sample of field stars.
△ Less
Submitted 15 June, 2015;
originally announced June 2015.
-
Overshooting by differential heating
Authors:
R. Andrássy,
H. C. Spruit
Abstract:
On the long nuclear time scale of stellar main-sequence evolution, even weak mixing processes can become relevant for redistributing chemical species in a star. We investigate a process of "differential heating," which occurs when a temperature fluctuation propagates by radiative diffusion from the boundary of a convection zone into the adjacent radiative zone. The resulting perturbation of the hy…
▽ More
On the long nuclear time scale of stellar main-sequence evolution, even weak mixing processes can become relevant for redistributing chemical species in a star. We investigate a process of "differential heating," which occurs when a temperature fluctuation propagates by radiative diffusion from the boundary of a convection zone into the adjacent radiative zone. The resulting perturbation of the hydrostatic equilibrium causes a flow that extends some distance from the convection zone. We study a simplified differential-heating problem with a static temperature fluctuation imposed on a solid boundary. The astrophysically relevant limit of a high Reynolds number and a low Péclet number (high thermal diffusivity) turns out to be interestingly non-intuitive. We derive a set of scaling relations for the stationary differential heating flow. A numerical method adapted to a high dynamic range in flow amplitude needed to detect weak flows is presented. Our two-dimensional simulations show that the flow reaches a stationary state and confirm the analytic scaling relations. These imply that the flow speed drops abruptly to a negligible value at a finite height above the source of heating. We approximate the mixing rate due to the differential heating flow in a star by a height-dependent diffusion coefficient and show that this mixing extends about $4\%$ of the pressure scale height above the convective core of a $10\,M_\odot$ zero-age main sequence star.
△ Less
Submitted 19 February, 2015;
originally announced February 2015.
-
Overshooting by convective settling
Authors:
R. Andrássy,
H. C. Spruit
Abstract:
We study a process of slow mixing in stars with convective envelopes, which is driven by the settling of cool downward plumes below the base of the convection zone. If a small fraction (of order $10^{-7}$) of the material cooled at the surface retains a significant entropy deficit while descending in plumes, it can reach the depth where lithium burning takes place. The model calculates the thermal…
▽ More
We study a process of slow mixing in stars with convective envelopes, which is driven by the settling of cool downward plumes below the base of the convection zone. If a small fraction (of order $10^{-7}$) of the material cooled at the surface retains a significant entropy deficit while descending in plumes, it can reach the depth where lithium burning takes place. The model calculates the thermal response and mixing below the convection zone due to the settling process, assuming that the plumes arrive at the base of the convection zone with a broad range of entropy contrasts. We obtain a good fit to the observed lithium depletion in the Sun by assuming that the settling mass flux is distributed with respect to the entropy contrast as a power law with a slope around -2. We find convective settling to have a negligible influence on the stratification below the convection zone, although mixing induced by it could modify the gradient of helium concentration.
△ Less
Submitted 30 October, 2013;
originally announced October 2013.