Showing 1–50 of 51 results for author: Meakin, C

3D simulations of convective shell Neon-burning in a massive star

Authors: C. Georgy, F. Rizzuti, R. Hirschi, V. Varma, W. D. Arnett, C. Meakin, M. Mocak, A. StJ. Murphy, T. Rauscher

Abstract: The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever increasing computing power makes now possible to simulate in 3D part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulen… ▽ More The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever increasing computing power makes now possible to simulate in 3D part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same timescale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 17 pages, 20 figures, accepted for publication in MNRAS

3D stellar evolution: hydrodynamic simulations of a complete burning phase in a massive star

Authors: F. Rizzuti, R. Hirschi, W. D. Arnett, C. Georgy, C. Meakin, A. StJ. Murphy, T. Rauscher, V. Varma

Abstract: Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multi-dimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an… ▽ More Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multi-dimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M$_\odot$ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution, thus revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions and compact remnants. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 12 pages, 15 figures. Accepted for publication in MNRAS

Realistic 3D hydrodynamics simulations find significant turbulent entrainment in massive stars

Authors: F. Rizzuti, R. Hirschi, C. Georgy, W. D. Arnett, C. Meakin, A. StJ. Murphy

Abstract: Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multi-dimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multi-dimensional processes over short timescal… ▽ More Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multi-dimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multi-dimensional processes over short timescales, thanks to the recent advances in computing resources. Among these processes, turbulent entrainment leading to mixing across convective boundaries is one of the least understood and most impactful. Here we present the results from a set of hydrodynamics simulations of the neon-burning shell in a massive star, and interpret them in the framework of the turbulent entrainment law from geophysics. Our simulations differ from previous studies in their unprecedented degree of realism in reproducing the stellar environment. Importantly, the strong entrainment found in the simulations highlights the major flaws of the current implementation of convective boundary mixing in 1D stellar models. This study therefore calls for major revisions of how convective boundaries are modelled in 1D, and in particular the implementation of entrainment in these models. This will have important implications for supernova theory, nucleosynthesis, neutron stars and black holes physics. △ Less

Submitted 16 August, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 7 pages, 6 figures, accepted for publication in MNRAS

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.

Comments: 20 pages, 19 figures, new Fig. B.2, new clarifications, electronic materials updated, accepted for publication in A&A

Journal ref: A&A 659, A193 (2022)

Convective core entrainment in 1D main sequence stellar models

Authors: L. J. A. Scott, R. Hirschi, C. Georgy, W. D. Arnett, C. Meakin, E. A. Kaiser, S. Ekström, N. Yusof

Abstract: 3D hydrodynamics models of deep stellar convection exhibit turbulent entrainment at the convective-radiative boundary which follows the entrainment law, varying with boundary penetrability. We implement the entrainment law in the 1D Geneva stellar evolution code. We then calculate models between 1.5 and 60 M$_{\odot}$ at solar metallicity ($Z=0.014$) and compare them to previous generations of mod… ▽ More 3D hydrodynamics models of deep stellar convection exhibit turbulent entrainment at the convective-radiative boundary which follows the entrainment law, varying with boundary penetrability. We implement the entrainment law in the 1D Geneva stellar evolution code. We then calculate models between 1.5 and 60 M$_{\odot}$ at solar metallicity ($Z=0.014$) and compare them to previous generations of models and observations on the main sequence. The boundary penetrability, quantified by the bulk Richardson number, $Ri_{\mathrm{B}}$, varies with mass and to a smaller extent with time. The variation of $Ri_{\mathrm{B}}$ with mass is due to the mass dependence of typical convective velocities in the core and hence the luminosity of the star. The chemical gradient above the convective core dominates the variation of $Ri_{\mathrm{B}}$ with time. An entrainment law method can therefore explain the apparent mass dependence of convective boundary mixing through $Ri_{\mathrm{B}}$. New models including entrainment can better reproduce the mass dependence of the main sequence width using entrainment law parameters $A \sim 2 \times 10^{-4}$ and $n=1$. We compare these empirically constrained values to the results of 3D hydrodynamics simulations and discuss implications. △ Less

Submitted 11 March, 2021; v1 submitted 10 March, 2021; originally announced March 2021.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

Dependence of Convective Boundary Mixing on Boundary Properties and Turbulence Strength

Authors: A. Cristini, R. Hirschi, C. Meakin, D. Arnett, C. Georgy, I. Walkington

Abstract: Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with th… ▽ More Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with the PROMPI code. We performed a mean field analysis of the simulations within the Reynolds-averaged Navier-Stokes framework. Both the vertical RMS velocity within the convective region and the bulk Richardson number of the boundaries are found to scale with the driving luminosity as expected from theory. The positions of the convective boundaries were estimated through the composition profiles across them, and the strength of convective boundary mixing was determined by analysing the boundaries within the framework of the entrainment law. We find that the entrainment is approximately inversely proportional to the bulk Richardson number. Although the entrainment law does not encompass all the processes occurring at boundaries, our results support the use of the entrainment law to describe convective boundary mixing in 1D models, at least for the advanced phases. The next steps and challenges ahead are also discussed. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS

3D Simulations and MLT: II. Onsager's Ideal Turbulence

Authors: W. David Arnett, Raphael Hirschi, Simon W. Campbell, Miroslav Mocák, Cyril Georgy, Casey Meakin, Andrea Cristini, Laura J. A. Scott, Etienne A. Kaiser, Maxime Viallet

Abstract: We simulate stellar convection at high Reynolds number (Re$\lesssim$7000) with causal time stepping but no explicit viscosity. We use the 3D Euler equations with shock capturing (Colella & Woodward 1984). Anomalous dissipation of turbulent kinetic energy occurs as an emergent feature of advection ("Onsager damping"), caused by the moderate shocks which terminate the turbulent kinetic energy spectr… ▽ More We simulate stellar convection at high Reynolds number (Re$\lesssim$7000) with causal time stepping but no explicit viscosity. We use the 3D Euler equations with shock capturing (Colella & Woodward 1984). Anomalous dissipation of turbulent kinetic energy occurs as an emergent feature of advection ("Onsager damping"), caused by the moderate shocks which terminate the turbulent kinetic energy spectrum; see also (Perry 2021). In strongly stratified stellar convection the asymptotic limit for the global damping length of turbulent kinetic energy is $\ell_d \sim \langle u^3 \rangle /\langle ε\rangle$. This "dissipative anomaly" (Onsager 1949) fixes the value of the "mixing length parameter", $α= \ell_{\rm MLT}/H_P =\overline{\langleΓ_1\rangle}$, which is $\sim\, 5/3$ for complete ionization. The estimate is numerically robust, agrees to within 10% with estimates from stellar evolution with constant $α$. For weak stratification $\ell_d$ shrinks to the depth of a thin convective region. Our flows are filamentary, produce surfaces of separation at boundary layers, resolve the energy-containing eddies, and develop a turbulent cascade down to the grid scale which agrees with the $4096^3$ direct numerical simulation of Kaneda (2003). The cascade converges quickly, and satisfies a power-law velocity spectrum similar to Kolmogorov (1941). Our flows exhibit intermittency, anisotropy, and interactions between coherent structures, features missing from K41 theory. We derive a dissipation rate from Reynolds stresses which agrees with (i) our flows, (ii) experiment (Warhaft 2002), and (iii) high Re simulations of the Navier-Stokes equations (Iyer, et al. 2018). △ Less

Submitted 23 April, 2021; v1 submitted 10 October, 2018; originally announced October 2018.

Comments: 31 pages, 6 figures. The time evolution and the fly-through movies may be found at \url{http://www.astro.keele.ac.uk/shyne/321D/convection-and-convective-boundary-mixing/visualisations}

3D Simulations and MLT: I. Renzini's Critique

Authors: W. David Arnett, Casey Meakin, Raphael Hirschi, Andrea Cristini, Cyril Georgy, Simon Campbell, Laura Scott, Etienne Kaiser

Abstract: Renzini (1987) wrote an influential critique of mixing-length theory (MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy. Implicit large eddy simulations (ILES)… ▽ More Renzini (1987) wrote an influential critique of mixing-length theory (MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy. Implicit large eddy simulations (ILES) automatically connect large scale turbulence to a Kolmogorov cascade below the grid scale, allowing turbulent boundary layers to remove singularities that appear in the theory. Interactions between coherent structures give multi-modal behavior, driving intermittency and fluctuations. Reynolds averaging (RA) allows us to abstract the essential features of this dynamical behavior of boundaries which are appropriate to stellar evolution, and consider how they relate static boundary conditions (Richardson, Schwarzschild or Ledoux). We clarify several questions concerning when and why MLT works, and does not work, using both analytical theory and 3D high resolution numerical simulations. The composition gradients and boundary layer structure which are produced by our simulations suggest a self-consistent approach to boundary layers, removing the need for ad hoc procedures for 'convective overshooting' and `semi-convection'. In a companion paper we quantify the adequacy of our numerical resolution, determine of the length scale of dissipation (the `mixing length') without astronomical calibration, quantify agreement with the four-fifths law of Kolmogorov for weak stratification, and extend MLT to deal with strong stratification. △ Less

Submitted 10 October, 2018; originally announced October 2018.

Turbulent Mixing and Nuclear Burning in Stellar Interiors

Authors: Miroslav Mocák, Casey Meakin, Simon Wattana Campbell, David Arnett

Abstract: The turbulent burning of nuclei is a common phenomenon in the evolution of stars. Here we examine a challenging case: the merging of the neon and oxygen burning shells in a 23 M$_{\odot}$ star. A previously unknown quasi-steady state is established by the interplay between mixing, turbulent transport, and nuclear burning. The resulting stellar structure has two burning shells within a single conve… ▽ More The turbulent burning of nuclei is a common phenomenon in the evolution of stars. Here we examine a challenging case: the merging of the neon and oxygen burning shells in a 23 M$_{\odot}$ star. A previously unknown quasi-steady state is established by the interplay between mixing, turbulent transport, and nuclear burning. The resulting stellar structure has two burning shells within a single convection zone. We find that the new neon burning layer covers an extended region of the convection zone, with the burning peak occurring substantially below where the Damköhler number first becomes equal to unity. These characteristics differ from those predicted by 1D stellar evolution models of similar ingestion events. We develop the mean-field turbulence equations that govern compositional evolution, and use them to interpret our data set. An important byproduct is a means to quantify sub-grid-scale effects intrinsic to the numerical hydrodynamic scheme. For implicit large eddy simulations, the analysis method is particularly powerful because it can reveal where and how simulated flows are modified by resolution, and provide straightforward physical interpretations of the effects of dissipation or induced transport. Focusing on the mean-field composition variance equations for our analysis, we recover a Kolmogorov rate of turbulent dissipation without it being imposed, in agreement with previous results which used the turbulent kinetic energy equation. △ Less

Submitted 20 August, 2018; originally announced August 2018.

Comments: 17 pages, 8 figures, Submitted to MNRAS

Progenitors of Core-Collapse Supernovae

Authors: Raphael Hirschi, David Arnett, Andrea Cristini, Cyril Georgy, Casey Meakin, Ian Walkington

Abstract: Massive stars have a strong impact on their surroundings, in particular when they produce a core-collapse supernova at the end of their evolution. In these proceedings, we review the general evolution of massive stars and their properties at collapse as well as the transition between massive and intermediate-mass stars. We also summarise the effects of metallicity and rotation. We then discuss som… ▽ More Massive stars have a strong impact on their surroundings, in particular when they produce a core-collapse supernova at the end of their evolution. In these proceedings, we review the general evolution of massive stars and their properties at collapse as well as the transition between massive and intermediate-mass stars. We also summarise the effects of metallicity and rotation. We then discuss some of the major uncertainties in the modelling of massive stars, with a particular emphasis on the treatment of convection in 1D stellar evolution codes. Finally, we present new 3D hydrodynamic simulations of convection in carbon burning and list key points to take from 3D hydrodynamic studies for the development of new prescriptions for convective boundary mixing in 1D stellar evolution codes. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 10 pages, 5 figures, "SN 1987A, 30 years later", Proceedings IAU Symposium No. 331, 2017; A. Marcowith, G. Dubner, A. Ray, A. Bykov, & M. Renaud, eds

The First 3D Simulations of Carbon Burning in a Massive Star

Authors: Andrea Cristini, Casey Meakin, Raphael Hirschi, David Arnett, Cyril Georgy, Maxime Viallet

Abstract: We present the first detailed three-dimensional hydrodynamic implicit large eddy simulations of turbulent convection for carbon burning. The simulations start with an initial radial profile mapped from a carbon burning shell within a 15 solar mass stellar evolution model. We considered 4 resolutions from 128^3 to 1024^3 zones. These simulations confirm that convective boundary mixing (CBM) occurs… ▽ More We present the first detailed three-dimensional hydrodynamic implicit large eddy simulations of turbulent convection for carbon burning. The simulations start with an initial radial profile mapped from a carbon burning shell within a 15 solar mass stellar evolution model. We considered 4 resolutions from 128^3 to 1024^3 zones. These simulations confirm that convective boundary mixing (CBM) occurs via turbulent entrainment as in the case of oxygen burning. The expansion of the boundary into the surrounding stable region and the entrainment rate are smaller at the bottom boundary because it is stiffer than the upper boundary. The results of this and similar studies call for improved CBM prescriptions in 1D stellar evolution models. △ Less

Submitted 19 December, 2017; originally announced December 2017.

Comments: 5 pages, 3 figures. Published in IAUS 329 on 28/07/17

Journal ref: Proceedings of the International Astronomical Union, 12(S329), 237-241 (2017)

3D Hydrodynamic Simulations of Carbon Burning in Massive Stars

Authors: Andrea Cristini, Casey Meakin, Raphael Hirschi, David Arnett, Cyril Georgy, Maxime Viallet

Abstract: We present the first detailed three-dimensional (3D) hydrodynamic implicit large eddy simulations of turbulent convection of carbon burning in massive stars. Simulations begin with radial profiles mapped from a carbon burning shell within a 15$\,\textrm{M}_\odot$ one-dimensional stellar evolution model. We consider models with $128^3$, $256^3$, $512^3$ and $1024^3$ zones. The turbulent flow proper… ▽ More We present the first detailed three-dimensional (3D) hydrodynamic implicit large eddy simulations of turbulent convection of carbon burning in massive stars. Simulations begin with radial profiles mapped from a carbon burning shell within a 15$\,\textrm{M}_\odot$ one-dimensional stellar evolution model. We consider models with $128^3$, $256^3$, $512^3$ and $1024^3$ zones. The turbulent flow properties of these carbon burning simulations are very similar to the oxygen burning case. We performed a mean field analysis of the kinetic energy budgets within the Reynolds-averaged Navier-Stokes framework. For the upper convective boundary region, we find that the numerical dissipation is insensitive to resolution for linear mesh resolutions above 512 grid points. For the stiffer, more stratified lower boundary, our highest resolution model still shows signs of decreasing sub-grid dissipation suggesting it is not yet numerically converged. We find that the widths of the upper and lower boundaries are roughly 30% and 10% of the local pressure scale heights, respectively. The shape of the boundaries is significantly different from those used in stellar evolution models. As in past oxygen-shell burning simulations, we observe entrainment at both boundaries in our carbon-shell burning simulations. In the large Péclet number regime found in the advanced phases, the entrainment rate is roughly inversely proportional to the bulk Richardson number, Ri$_{\rm B}$ ($\propto $Ri${\rm_B}^{-α}$, $0.5\lesssim α\lesssim 1.0$). We thus suggest the use of Ri$_{\rm B}$ as a means to take into account the results of 3D hydrodynamics simulations in new 1D prescriptions of convective boundary mixing. △ Less

Submitted 18 June, 2017; v1 submitted 17 October, 2016; originally announced October 2016.

Comments: 26 pages, 15 figures, accepted for publication in MNRAS, movie available at the following URL: http://www.astro.keele.ac.uk/shyne/321D/convection-and-convective- boundary-mixing/visualisations/very-high-resolution-movie-of-the-c-shell/view

Key Issues Review: Numerical studies of turbulence in stars

Authors: W. David Arnett, Casey Meakin

Abstract: The numerical simulation of turbulence in stars has led to a rich set of possibilities regarding stellar pulsations, asteroseismology, thermonuclear yields, and formation of neutron stars and black holes. The breaking of symmetry by turbulent flow grows in amplitude as collapse is approached, which insures that the conditions at the onset of collapse are not spherical. This lack of spherical symme… ▽ More The numerical simulation of turbulence in stars has led to a rich set of possibilities regarding stellar pulsations, asteroseismology, thermonuclear yields, and formation of neutron stars and black holes. The breaking of symmetry by turbulent flow grows in amplitude as collapse is approached, which insures that the conditions at the onset of collapse are not spherical. This lack of spherical symmetry has important implications for the mechanism of explosion and ejected nucleosynthesis products. Numerical resolution of several different types of three--dimensional (3D) stellar simulations are compared; it is suggested that core collapse simulations may be under-resolved. New physical effects which appear in 3D are summarized. Connections between simulations of progenitor explosion and observations of supernova remnants (SNR) are discussed. Present treatment of boundaries, for mixing regions during He--burning, requires revision. △ Less

Submitted 17 March, 2016; originally announced March 2016.

Comments: 8 pages, 1 figure, 1 table, submitted to Reports on Progress in Physics

Linking 1D Evolutionary to 3D Hydrodynamical Simulations of Massive Stars

Authors: Andréa Cristini, Casey Meakin, Raphael Hirschi, David Arnett, Cyril Georgy, Maxime Viallet

Abstract: Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective… ▽ More Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective boundaries in stellar models has been studied for decades but still lacks a definitive solution. This paper presents initial results of a study on convective boundary mixing (CBM) in massive stars. The 'stiffness' of a convective boundary can be quantified using the bulk Richardson number ($\textrm{Ri}_B$), the ratio of the potential energy for restoration of the boundary to the kinetic energy of turbulent eddies. A 'stiff' boundary ($\textrm{Ri}_B \sim 10^4$) will suppress CBM, whereas in the opposite case a 'soft' boundary ($\textrm{Ri}_B \sim 10$) will be more susceptible to CBM. One of the key results obtained so far is that lower convective boundaries (closer to the centre) of nuclear burning shells are 'stiffer' than the corresponding upper boundaries, implying limited CBM at lower shell boundaries. This is in agreement with 3D hydrodynamic simulations carried out by Meakin and Arnett [The Astrophysical Journal 667:448-475, 2007]. This result also has implications for new CBM prescriptions in massive stars as well as for nuclear burning flame front propagation in Super-Asymptotic Giant Branch stars and also the onset of novae. △ Less

Submitted 7 January, 2016; originally announced January 2016.

Comments: Accepted for publication (12/12/15) in the Physica Scripta focus issue on Turbulent Mixing and Beyond

Towards 21st Century Stellar Models: Star Clusters, Supercomputing, and Asteroseismology

Authors: S. W. Campbell, T. N. Constantino, V. D'Orazi, C. Meakin, D. Stello, J. Christensen-Dalsgaard, C. Kuehn, G. M. De Silva, W. D. Arnett, J. C. Lattanzio, B. T. MacLean

Abstract: Stellar models provide a vital basis for many aspects of astronomy and astrophysics. Recent advances in observational astronomy -- through asteroseismology, precision photometry, high-resolution spectroscopy, and large-scale surveys -- are placing stellar models under greater quantitative scrutiny than ever. The model limitations are being exposed and the next generation of stellar models is neede… ▽ More Stellar models provide a vital basis for many aspects of astronomy and astrophysics. Recent advances in observational astronomy -- through asteroseismology, precision photometry, high-resolution spectroscopy, and large-scale surveys -- are placing stellar models under greater quantitative scrutiny than ever. The model limitations are being exposed and the next generation of stellar models is needed as soon as possible. The current uncertainties in the models propagate to the later phases of stellar evolution, hindering our understanding of stellar populations and chemical evolution. Here we give a brief overview of the evolution, importance, and substantial uncertainties of core helium burning stars in particular and then briefly discuss a range of methods, both theoretical and observational, that we are using to advance the modelling. △ Less

Submitted 15 December, 2015; originally announced December 2015.

Comments: To appear in Astronomische Nachrichten, special issue "Reconstruction the Milky Way's History: Spectroscopic surveys, Asteroseismology and Chemo-dynamical models", Guest Editors C. Chiappini, J. Montalbán, and M. Steffen, AN 2016 (in press)

Toward a consistent use of overshooting parametrizations in 1D stellar evolution codes

Authors: Maxime Viallet, Casey Meakin, Vincent Prat, David Arnett

Abstract: Several parametrizations for overshooting in 1D stellar evolution calculations coexist in the literature. These parametrizations are used somewhat arbitrarily in stellar evolution codes, based on what works best for a given problem, or even for historical reasons related to the development of each code. We bring attention to the fact that these different parametrizations correspond to different ph… ▽ More Several parametrizations for overshooting in 1D stellar evolution calculations coexist in the literature. These parametrizations are used somewhat arbitrarily in stellar evolution codes, based on what works best for a given problem, or even for historical reasons related to the development of each code. We bring attention to the fact that these different parametrizations correspond to different physical regimes of overshooting, depending whether the effects of radiation are dominant, marginal, or negligible. Our analysis is based on previously published theoretical results, as well as multidimensional hydrodynamical simulations of stellar convection where the interaction between the convective region and a stably-stratified region is observed. Although the underlying hydrodynamical processes are the same, the outcome of the overshooting process is profoundly affected by radiative effects. Using a simple picture of the scales involved in the overshooting process, we show how three regimes are obtained, depending on the importance of radiative effects. These three regimes correspond to the different behaviors observed in hydrodynamical simulations so far, and to the three types of parametrizations used in 1D codes. We suggest that the existing parametrizations for overshooting should coexist in 1D stellar evolution codes, and should be applied consistently at convective boundaries depending on the local physical conditions. △ Less

Submitted 26 June, 2015; v1 submitted 9 June, 2015; originally announced June 2015.

Comments: 5 pages, 2 figures, to appear in A&A as a regular paper. Last version: language editing + typos in Eq. (6) & (9) corrected

Journal ref: A&A 580, A61 (2015)

Beyond Mixing-length Theory: a step toward 321D

Authors: W. David Arnett, Casey Meakin, Maxime Viallet, Simon W. Campbell, John Lattanzio, Miroslav Mocák

Abstract: We examine the physical basis for algorithms to replace mixing-length theory (MLT) in stellar evolutionary computations. Our 321D procedure is based on numerical solutions of the Navier-Stokes equations. These implicit large eddy simulations (ILES) are three-dimensional (3D), time-dependent, and turbulent, including the Kolmogorov cascade. We use the Reynolds-averaged Navier-Stokes (RANS) formulat… ▽ More We examine the physical basis for algorithms to replace mixing-length theory (MLT) in stellar evolutionary computations. Our 321D procedure is based on numerical solutions of the Navier-Stokes equations. These implicit large eddy simulations (ILES) are three-dimensional (3D), time-dependent, and turbulent, including the Kolmogorov cascade. We use the Reynolds-averaged Navier-Stokes (RANS) formulation to make concise the 3D simulation data, and use the 3D simulations to give closure for the RANS equations. We further analyze this data set with a simple analytical model, which is non-local and time-dependent, and which contains both MLT and the Lorenz convective roll as particular subsets of solutions. A characteristic length (the damping length) again emerges in the simulations; it is determined by an observed balance between (1) the large-scale driving, and (2) small-scale damping. The nature of mixing and convective boundaries is analyzed, including dynamic, thermal and compositional effects, and compared to a simple model. We find that (1) braking regions (boundary layers in which mixing occurs) automatically appear {\it beyond} the edges of convection as defined by the Schwarzschild criterion, (2) dynamic (non-local) terms imply a non-zero turbulent kinetic energy flux (unlike MLT), (3) the effects of composition gradients on flow can be comparable to thermal effects, and (4) convective boundaries in neutrino-cooled stages differ in nature from those in photon-cooled stages (different Péclet numbers). The algorithms are based upon ILES solutions to the Navier-Stokes equations, so that, unlike MLT, they do not require any calibration to astronomical systems in order to predict stellar properties. Implications for solar abundances, helioseismology, asteroseismology, nucleosynthesis yields, supernova progenitors and core collapse are indicated. △ Less

Submitted 9 June, 2015; v1 submitted 1 March, 2015; originally announced March 2015.

Comments: 22 pages, 4 figures, 2 tables; significantly re-written, critique of Pasetto, et al. model added, accepted for publication by ApJ

Linking 1D Stellar Evolution to 3D Hydrodynamical Simulations

Authors: Andrea Cristini, Raphael Hirschi, Cyril Georgy, Casey Meakin, David Arnett, Maxime Viallet

Abstract: In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness'. Models for a 15Mo star were computed. We f… ▽ More In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness'. Models for a 15Mo star were computed. We found that for convective shells above the core, the lower (in radius or mass) boundaries are `stiffer' according to the bulk Richardson number than the relative upper (Schwarzschild) boundaries. Thus, we expect reduced CBM at the lower boundaries in comparison to the upper. This has implications on flame front propagation and the onset of novae. △ Less

Submitted 28 October, 2014; originally announced October 2014.

Comments: 2 pages, 1 figure. To appear in proceedings of the IAU Symposium 307: New Windows on Massive Stars: Asteroseismology, Interferometry and Spectropolarimetry

Compressible Hydrodynamic Mean-Field Equations in Spherical Geometry and their Application to Turbulent Stellar Convection Data

Authors: Miroslav Mocák, Casey Meakin, Maxime Viallet, David Arnett

Abstract: We present a statistical analysis of turbulent convection in stars within our Reynolds-Averaged Navier Stokes (RANS) framework in spherical geometry which we derived from first principles. The primary results reported in this document include: (1) an extensive set of mean-field equations for compressible, multi-species hydrodynamics, and (2) corresponding mean-field data computed from various simu… ▽ More We present a statistical analysis of turbulent convection in stars within our Reynolds-Averaged Navier Stokes (RANS) framework in spherical geometry which we derived from first principles. The primary results reported in this document include: (1) an extensive set of mean-field equations for compressible, multi-species hydrodynamics, and (2) corresponding mean-field data computed from various simulation models. Some supplementary scale analysis data is also presented. The simulation data which is presented includes: (1) shell convection during oxygen burning in a 23 solar mass supernova progenitor, (2) envelope convection in a 5 solar mass red giant, (3) shell convection during the helium flash, and (4) a hydrogen injection flash in a 1.25 solar mass star. These simulations have been partially described previously in Meakin [2006], Meakin and Arnett [2007a,b, 2010], Arnett et al. [2009, 2010], Viallet et al. [2011, 2013a,b] and Mocak et al. [2009, 2011]. New data is also included in this document with several new domain and resolution configurations as well as some variations in the physical model such as convection zone depth and driving source term. The long term goal of this work is to aid in the development of more sophisticated models for treating hydrodynamic phenomena (e.g., turbulent convection) in the field of stellar evolution by providing a direct link between 3D simulation data and the mean fields which are modeled by 1D stellar evolution codes. As such, this data can be used to test previously proposed turbulence models found in the literature and sometimes used in stellar modeling. This data can also serve to test basic physical principles for model building and inspire new prescriptions for use in 1D evolution codes. △ Less

Submitted 21 January, 2014; originally announced January 2014.

Comments: 176 Pages, 115 Figures

Chaos and Turbulent Nucleosynthesis Prior to a Supernova Explosion

Authors: W. David Arnett, Casey Meakin, Maxime Viallet

Abstract: Three-dimensional (3D), time dependent numerical simulations, of flow of matter in stars, now have sufficient resolution to be fully turbulent. The late stages of the evolution of massive stars, leading up to core collapse to a neutron star (or black hole), and often to supernova explosion and nucleosynthesis, are strongly convective because of vigorous neutrino cooling and nuclear heating. Unlike… ▽ More Three-dimensional (3D), time dependent numerical simulations, of flow of matter in stars, now have sufficient resolution to be fully turbulent. The late stages of the evolution of massive stars, leading up to core collapse to a neutron star (or black hole), and often to supernova explosion and nucleosynthesis, are strongly convective because of vigorous neutrino cooling and nuclear heating. Unlike models based on current stellar evolutionary practice, these simulations show a chaotic dynamics characteristic of highly turbulent flow. Theoretical analysis of this flow, both in the Reynolds-averaged Navier-Stokes (RANS) framework and by simple dynamic models, show an encouraging consistency with the numerical results. It may now be possible to develop physically realistic and robust procedures for convection and mixing which (unlike 3D numerical simulation) may be applied throughout the long life times of stars. In addition, a new picture of the presupernova stages is emerging which is more dynamic and interesting (i.e., predictive of new and newly observed phenomena) than our previous one. △ Less

Submitted 28 January, 2014; v1 submitted 11 December, 2013; originally announced December 2013.

Comments: 11 pages, 2 figures, Submitted to AIP Advances: Stardust, added figures and modest rewriting

Turbulent convection in stellar interiors. III. Mean-field analysis and stratification effects

Authors: Maxime Viallet, Casey Meakin, David Arnett, Miroslav Mocak

Abstract: We present 3D implicit large eddy simulations (ILES) of the turbulent convection in the envelope of a 5 Msun red giant star and in the oxygen-burning shell of a 23 Msun supernova progenitor. The numerical models are analyzed in the framework of 1D Reynolds-Averaged Navier-Stokes (RANS) equations. The effects of pressure fluctuations are more important in the red giant model, owing to larger strati… ▽ More We present 3D implicit large eddy simulations (ILES) of the turbulent convection in the envelope of a 5 Msun red giant star and in the oxygen-burning shell of a 23 Msun supernova progenitor. The numerical models are analyzed in the framework of 1D Reynolds-Averaged Navier-Stokes (RANS) equations. The effects of pressure fluctuations are more important in the red giant model, owing to larger stratification of the convective zone. We show how this impacts different terms in the mean-field equations. We clarify the driving sources of kinetic energy, and show that the rate of turbulent dissipation is comparable to the convective luminosity. Although our flows have low Mach number and are nearly adiabatic, our analysis is general and can be applied to photospheric convection as well. The robustness of our analysis of turbulent convection is supported by the insensitivity of the mean-field balances to linear mesh resolution. We find robust results for the turbulent convection zone and the stable layers in the oxygen-burning shell model, and robust results everywhere in the red giant model, but the mean fields are not well converged in the narrow boundary regions (which contain steep gradients) in the oxygen-burning shell model. This last result illustrates the importance of unresolved physics at the convective boundary, which governs the mixing there. △ Less

Submitted 13 March, 2013; v1 submitted 27 December, 2012; originally announced December 2012.

Comments: 26 pages, 20 figures, Accepted for publication in ApJ

Hydrodynamic Studies of the Evolution of Recurrent, Symbiotic, and Dwarf Novae: The White Dwarf Components are Growing in Mass

Authors: S. Starrfield, F. X. Timmes, C. Iliadis, W. R. Hix, W. D. Arnett, C. Meakin, W. M. Sparks

Abstract: Symbiotic binaries are systems containing white dwarfs (WDs) and red giants. Symbiotic novae are those systems in which thermonuclear eruptions occur on the WD components. These are to be distinguished from events driven by accretion disk instabilities analogous to dwarf novae eruptions in cataclysmic variable outbursts. Another class of symbiotic systems are those in which the WD is extremely lum… ▽ More Symbiotic binaries are systems containing white dwarfs (WDs) and red giants. Symbiotic novae are those systems in which thermonuclear eruptions occur on the WD components. These are to be distinguished from events driven by accretion disk instabilities analogous to dwarf novae eruptions in cataclysmic variable outbursts. Another class of symbiotic systems are those in which the WD is extremely luminous and it seems likely that quiescent nuclear burning is ongoing on the accreting WD. A fundamental question is the secular evolution of the WD. Do the repeated outbursts or quiescent burning in these accreting systems cause the WD to gain or lose mass? If it is gaining mass, can it eventually reach the Chandrasekhar Limit and become a supernova (a SN Ia if it can hide the hydrogen and helium in the system)? In order to better understand these systems, we have begun a new study of the evolution of Thermonuclear Runaways (TNRs) in the accreted envelopes of WDs using a variety of initial WD masses, luminosities and mass accretion rates. We use our 1-D hydro code, NOVA, which includes the new convective algorithm of Arnett, Meakin and Young, the Hix and Thielemann nuclear reaction solver, the Iliadis reaction rate library, the Timmes equation of state, and the OPAL opacities. We assume a solar composition (Lodders abundance distribution) and do not allow any mixing of accreted material with core material. This assumption strongly influences our results. We report here (1) that the WD grows in mass for all simulations so that canonical `steady burning' does not occur, and (2) that only a small fraction of the accreted matter is ejected in some (but not all) simulations. We also find that the accreting systems, before thermonuclear runaway, are too cool to be seen in X-ray searches for SN Ia progenitors. △ Less

Submitted 26 November, 2012; originally announced November 2012.

Comments: 11 pages 10 figures in the 2011 Asiago Symbiotic Workshop

Journal ref: Baltic Astronomy, Vol 21, pages 76 - 87, 2012

General-Relativistic Simulations of Three-Dimensional Core-Collapse Supernovae

Authors: C. D. Ott, E. Abdikamalov, P. Moesta, R. Haas, S. Drasco, E. O'Connor, C. Reisswig, C. Meakin, E. Schnetter

Abstract: We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects… ▽ More We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27-solar-mass progenitor was studied in 2D by B. Mueller et al. (ApJ 761:72, 2012), who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability. △ Less

Submitted 25 March, 2013; v1 submitted 24 October, 2012; originally announced October 2012.

Comments: 24 emulateapj pages, 15 figures. Accepted for publication in the Astrophysical Journal

Theoretical Studies of Accretion of Matter onto White Dwarfs and the Single Degenerate Scenario for Supernovae of Type Ia

Authors: S. Starrfield, C. Iliadis, F. X. Timmes, W. R. Hix, W. D. Arnett, C. Meakin, W. M. Sparks

Abstract: We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and s… ▽ More We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and steady burning does not occur. For a broad range in WD mass (0.4 Solar masses to 1.35 Solar Masses), the maximum ejected material occurs for the 1.25 Solar Mass sequences and then decreases as the white dwarf mass decreases. Therefore, the white dwarfs are growing in mass as a consequence of the accretion of solar material and as long as there is no mixing of accreted material with core material. In contrast, a thermonuclear runaway in the accreted hydrogen-rich layers on the low luminosity WDs in close binary systems where mixing of core matter with accreted material has occurred is the outburst mechanism for Classical, Recurrent, and Symbiotic novae. The differences in characteristics of these systems is likely the WD mass and mass accretion rate. The high levels of enrichment of CN ejecta in elements ranging from carbon to sulfur confirm that there is dredge-up of matter from the core of the WD and enable them to contribute to the chemical enrichment of the interstellar medium. Therefore, studies of CNe can lead to an improved understanding of Galactic nucleosynthesis, some sources of pre-solar grains, and the Extragalactic distance scale. The characteristics of the outburst depend on the white dwarf mass, luminosity, mass accretion rate, and the chemical composition of both the accreting material and WD material. The properties of the outburst also depends on when, how, and if the accreted layers are mixed with the WD core and the mixing mechanism is still unknown. △ Less

Submitted 22 October, 2012; originally announced October 2012.

Comments: 25 Pages, Bulletin of the Astronomical Society of India (BASI) in press

The Detonation Mechanism of the Pulsationally-Assisted Gravitationally-Confined Detonation Model of Type Ia Supernovae

Authors: G. C. Jordan IV, C. Graziani, R. T. Fisher, D. M. Townsley, C. Meakin, K. Weide, L. B. Reid, J. Norris, R. Hudson, D. Q. Lamb

Abstract: We describe the detonation mechanism comprising the "Pulsationally Assisted" Gravitationally Confined Detonation (GCD) model of Type Ia supernovae SNe Ia. This model is analogous to the previous GCD model reported in Jordan et al.(2008); however, the chosen initial conditions produce a substantively different detonation mechanism, resulting from a larger energy release during the deflagration phas… ▽ More We describe the detonation mechanism comprising the "Pulsationally Assisted" Gravitationally Confined Detonation (GCD) model of Type Ia supernovae SNe Ia. This model is analogous to the previous GCD model reported in Jordan et al.(2008); however, the chosen initial conditions produce a substantively different detonation mechanism, resulting from a larger energy release during the deflagration phase. The resulting final kinetic energy and nickel-56 yields conform better to observational values than is the case for the "classical" GCD models. In the present class of models, the ignition of a deflagration phase leads to a rising, burning plume of ash. The ash breaks out of the surface of the white dwarf, flows laterally around the star, and converges on the collision region at the antipodal point from where it broke out. The amount of energy released during the deflagration phase is enough to cause the star to rapidly expand, so that when the ash reaches the antipodal point, the surface density is too low to initiate a detonation. Instead, as the ash flows into the collision region (while mixing with surface fuel), the star reaches its maximally expanded state and then contracts. The stellar contraction acts to increase the density of the star, including the density in the collision region. This both raises the temperature and density of the fuel-ash mixture in the collision region and ultimately leads to thermodynamic conditions that are necessary for the Zel'dovich gradient mechanism to produce a detonation. We demonstrate feasibility of this scenario with three 3-dimensional (3D), full star simulations of this model using the FLASH code. We characterized the simulations by the energy released during the deflagration phase, which ranged from 38% to 78% of the white dwarf's binding energy. We show that the necessary conditions for detonation are achieved in all three of the models. △ Less

Submitted 17 September, 2012; v1 submitted 17 February, 2012; originally announced February 2012.

Comments: 22 pages, 8 figures; ApJ

A new stellar mixing process operating below shell convection zones following off-center ignition

Authors: M. Mocák, Casey A. Meakin, E. Müller, L. Siess

Abstract: During most stages of stellar evolution the nuclear burning of lighter to heavier elements results in a radial composition profile which is stabilizing against buoyant acceleration, with light material residing above heavier material. However, under some circumstances, such as off-center ignition, the composition profile resulting from nuclear burning can be destabilizing, and characterized by an… ▽ More During most stages of stellar evolution the nuclear burning of lighter to heavier elements results in a radial composition profile which is stabilizing against buoyant acceleration, with light material residing above heavier material. However, under some circumstances, such as off-center ignition, the composition profile resulting from nuclear burning can be destabilizing, and characterized by an outwardly increasing mean molecular weight. The potential for instabilities under these circumstances, and the consequences that they may have on stellar structural evolution, remain largely unexplored. In this paper we study the development and evolution of instabilities associated with unstable composition gradients in regions which are initially stable according to linear Schwarzschild and Ledoux criteria. In particular, we explore the mixing taking place under various conditions with multi-dimensional hydrodynamic convection models based on stellar evolutionary calculations of the core helium flash in a 1.25 \Msun star, the core carbon flash in a 9.3\,\Msun star, and of oxygen shell burning in a star with a mass of 23\,\Msun. The results of our simulations reveal a mixing process associated with regions having outwardly increasing mean molecular weight that reside below convection zones. The mixing is not due to overshooting from the convection zone, nor is it due directly to thermohaline mixing which operates on a timescale several orders of magnitude larger than the simulated flows. Instead, the mixing appears to be due to the presence of a wave field induced in the stable layers residing beneath the convection zone which enhances the mixing rate by many orders of magnitude and allows a thermohaline type mixing process to operate on a dynamical, rather than thermal, timescale. We discuss our results in terms of related laboratory phenomena and associated theoretical developments. △ Less

Submitted 22 August, 2011; originally announced August 2011.

Comments: accepted for publication in Astrophysical Journal, 9 pages, 8 figures

A Global Turbulence Model for Neutrino-Driven Convection in Core-Collapse Supernovae

Authors: Jeremiah W. Murphy, Casey Meakin

Abstract: Simulations of core-collapse supernovae (CCSNe) result in successful explosions once the neutrino luminosity exceeds a critical curve, and recent simulations indicate that turbulence further enables explosion by reducing this critical neutrino luminosity. We propose a theoretical framework to derive this result and take the first steps by deriving the governing mean-field equations. Using Reynolds… ▽ More Simulations of core-collapse supernovae (CCSNe) result in successful explosions once the neutrino luminosity exceeds a critical curve, and recent simulations indicate that turbulence further enables explosion by reducing this critical neutrino luminosity. We propose a theoretical framework to derive this result and take the first steps by deriving the governing mean-field equations. Using Reynolds decomposition, we decompose flow variables into background and turbulent flows and derive self-consistent averaged equations for their evolution. As basic requirements for the CCSN problem, these equations naturally incorporate steady-state accretion, neutrino heating and cooling, non-zero entropy gradients, and turbulence terms associated with buoyant driving, redistribution, and dissipation. Furthermore, analysis of two-dimensional (2D) CCSN simulations validate these Reynolds-averaged equations, and we show that the physics of turbulence entirely accounts for the differences between 1D and 2D CCSN simulations. As a prelude to deriving the reduction in the critical luminosity, we identify the turbulent terms that most influence the conditions for explosion. Generically, turbulence equations require closure models, but these closure models depend upon the macroscopic properties of the flow. To derive a closure model that is appropriate for CCSNe, we cull the literature for relevant closure models and compare each with 2D simulations. These models employ local closure approximations and fail to reproduce the global properties of neutrino-driven turbulence. Motivated by the generic failure of these local models, we propose an original model for turbulence which incorporates global properties of the flow. This global model accurately reproduces the turbulence profiles and evolution of 2D CCSN simulations. △ Less

Submitted 13 October, 2011; v1 submitted 27 June, 2011; originally announced June 2011.

Comments: 20 pages in emulateapj, 14 figures, accepted by ApJ

Asymmetry and the Nucleosynthetic Signature of Nearly Edge-Lit Detonation in White Dwarf Cores

Authors: David A. Chamulak, Casey A. Meakin, Ivo R. Seitenzahl, James W. Truran

Abstract: Most of the leading explosion scenarios for Type Ia supernovae involve the nuclear incineration of a white dwarf star through a detonation wave. Several scenarios have been proposed as to how this detonation may actually occur, but the exact mechanism and environment in which it takes place remain unknown. We explore the effects of an off-center initiated detonation on the spatial distribution of… ▽ More Most of the leading explosion scenarios for Type Ia supernovae involve the nuclear incineration of a white dwarf star through a detonation wave. Several scenarios have been proposed as to how this detonation may actually occur, but the exact mechanism and environment in which it takes place remain unknown. We explore the effects of an off-center initiated detonation on the spatial distribution of the nucleosynthetic yield products in a toy model -- a pre-expanded near Chandrasekhar-mass white dwarf. We find that a single-point near edge-lit detonation results in asymmetries in the density and thermal profiles, notably the expansion timescale, throughout the supernova ejecta. We demonstrate that this asymmetry of the thermodynamic trajectories should be common to off-center detonations where a small amount of the star is burned prior to detonation. The sensitivity of the yields on the expansion timescale results in an asymmetric distribution of the elements synthesized as reaction products. We tabulate the shift in the center of mass of the various elements produced in our model supernova and find an odd-even pattern for elements past silicon. Our calculations show that off-center single-point detonations in carbon-oxygen white dwarfs are marked by significant composition asymmetries in their remnants which bear potentially observable signatures in both velocity and coordinate space, including an elemental nickel mass fraction which varies by a factor of two to three from one side of the remnant to the other. △ Less

Submitted 16 September, 2011; v1 submitted 15 March, 2011; originally announced March 2011.

Comments: 7 pages, 7 figures, accepted for publication in the Astrophysical Journal

Towards Realistic Progenitors of Core-Collapse Supernovae

Authors: W. David Arnett, Casey Meakin

Abstract: Two-dimensional (2D) hydrodynamical simulations of progenitor evolution of a 23 solar mass star, close to core collapse (about 1 hour, in 1D), with simultaneously active C, Ne, O, and Si burning shells, are presented and contrasted to existing 1D models (which are forced to be quasi-static). Pronounced asymmetries, and strong dynamical interactions between shells are seen in 2D. Although instigate… ▽ More Two-dimensional (2D) hydrodynamical simulations of progenitor evolution of a 23 solar mass star, close to core collapse (about 1 hour, in 1D), with simultaneously active C, Ne, O, and Si burning shells, are presented and contrasted to existing 1D models (which are forced to be quasi-static). Pronounced asymmetries, and strong dynamical interactions between shells are seen in 2D. Although instigated by turbulence, the dynamic behavior proceeds to sufficiently large amplitudes that it couples to the nuclear burning. Dramatic growth of low order modes is seen, as well as large deviations from spherical symmetry in the burning shells. The vigorous dynamics is more violent than that seen in earlier burning stages in the 3D simulations of a single cell in the oxygen burning shell, or in 2D simulations not including an active Si shell. Linear perturbative analysis does not capture the chaotic behavior of turbulence (e.g., strange attractors such as that discovered by Lorenz), and therefore badly underestimates the vigor of the instability. The limitations of 1D and 2D models are discussed in detail. The 2D models, although flawed geometrically, represent a more realistic treatment of the relevant dynamics than existing 1D models, and present a dramatically different view of the stages of evolution prior to collapse. Implications for interpretation of SN1987A, abundances in young supernova remnants, pre-collapse outbursts, progenitor structure, neutron star kicks, and fallback are outlined. While 2D simulations provide new qualitative insight, fully 3D simulations are needed for a quantitative understanding of this stage of stellar evolution. The necessary properties of such simulations are delineated. △ Less

Submitted 28 January, 2011; originally announced January 2011.

Comments: 26 pages, 1 table, 4 figures

Turbulent Cells in Stars: I. Fluctuations in Kinetic Energy and Luminosity

Authors: W. David Arnett, Casey Meakin

Abstract: Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett, 2007b) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a general instability of the convective cell, requiring only a localized source of heating or cooling. Such fluctuations are shown to be suppressed in simulations of stellar evolution which use mixing-len… ▽ More Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett, 2007b) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a general instability of the convective cell, requiring only a localized source of heating or cooling. Such fluctuations are shown to be suppressed in simulations of stellar evolution which use mixing-length theory (MLT). Quantitatively similar behavior occurs in the model of a convective roll (cell) of Lorenz (1963), which is known to have a strange attractor that gives rise to chaotic fluctuations in time of velocity and, as we show, luminosity. Study of simulations suggests that the behavior of a Lorenz convective roll may resemble that of a cell in convective flow. We examine some implications of this simplest approximation, and suggest paths for improvement. Using the Lorenz model as representative of a convective cell, a multiple-cell model of a convective layer gives total luminosity fluctuations which are suggestive of irregular variables (red giants and supergiants (Schwarzschild 1975)), and of the long secondary period feature in semi-regular AGB variables (Stothers 2010, Wood, Olivier and Kawaler 2004). This "tau-mechanism" is a new source for stellar variability, which is inherently non-linear (unseen in linear stability analysis), and one closely related to intermittency in turbulence. It was already implicit in the 3D global simulations of Woodward, Porter and Jacobs (2003). This fluctuating behavior is seen in extended 2D simulations of CNeOSi burning shells (Arnett and Meakin 2011b), and may cause instability which leads to eruptions in progenitors of core collapse supernovae PRIOR to collapse. △ Less

Submitted 30 July, 2011; v1 submitted 8 December, 2010; originally announced December 2010.

Comments: 30 pages, 13 figures

Time-dependent Turbulence in Stars

Authors: W. David Arnett, Casey Meakin

Abstract: Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett 2007) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a global instability in the convective region, which has been suppressed in calculations of stellar evolution which use mixing-length theory (MLT). Quantitatively similar behavior occurs in the model of a c… ▽ More Three-dimensional (3D) hydrodynamic simulations of shell oxygen burning (Meakin and Arnett 2007) exhibit bursty, recurrent fluctuations in turbulent kinetic energy. These are shown to be due to a global instability in the convective region, which has been suppressed in calculations of stellar evolution which use mixing-length theory (MLT). Quantitatively similar behavior occurs in the model of a convective roll (cell) of Lorenz (1963), which is known to have a strange attractor that gives rise to random fluctuations in time.An extension of the Lorenz model, which includes Kolmogorov damping and nuclear burning, is shown to exhibit bursty, recurrent fluctuations like those seen in the 3D simulations. A simple model of a convective layer (composed of multiple Lorenz cells) gives luminosity fluctuations which are suggestive of irregular variables (red giants and supergiants, Schwarzschild 1975). Apparent inconsistencies between Arnett, Meakin, and Young (2009) and Nordlund, Stein, and Asplund (2009) on the nature of convective driving have been resolved, and are discussed. △ Less

Submitted 27 October, 2010; originally announced October 2010.

Comments: 8 pages, 2 figures, IAU Symposium 271 "Astrophysical Dynamics: From Galaxies to Stars", Nice, FR, 2010

Spectra of Type Ia Supernovae from Double Degenerate Mergers

Authors: Chris L. Fryer, Ashley J. Ruiter, Krzysztof Belczynski, Peter J. Brown, Filomena Bufano, Steven Diehl, Christopher J. Fontes, Lucille H. Frey, Stephen T. Holland, Aimee L. Hungerford, Stefan Immler, Paolo Mazzali, Casey Meakin, Peter A. Milne, Cody Raskin, Francis X. Timmes

Abstract: The merger of two white dwarfs (a.k.a. double degenerate merger) has often been cited as a potential progenitor of type Ia supernovae. Here we combine population synthesis, merger and explosion models with radiation-hydrodynamics light-curve models to study the implications of such a progenitor scenario on the observed type Ia supernova population. Our standard model, assuming double degenerate me… ▽ More The merger of two white dwarfs (a.k.a. double degenerate merger) has often been cited as a potential progenitor of type Ia supernovae. Here we combine population synthesis, merger and explosion models with radiation-hydrodynamics light-curve models to study the implications of such a progenitor scenario on the observed type Ia supernova population. Our standard model, assuming double degenerate mergers do produce thermonuclear explosions, produces supernova light-curves that are broader than the observed type Ia sample. In addition, we discuss how the shock breakout and spectral features of these double degenerate progenitors will differ from the canonical bare Chandrasekhar-massed explosion models. We conclude with a discussion of how one might reconcile these differences with current observations. △ Less

Submitted 4 July, 2010; originally announced July 2010.

Comments: 29 pages, 12 figures, comments welcome

Report number: LA-UR 10-01202

Presupernova Structure of Massive Stars

Authors: Casey A. Meakin, Tuguldur Sukhbold, David Arnett

Abstract: Issues concerning the structure and evolution of core collapse progenitor stars are discussed with an emphasis on interior evolution. We describe a program designed to investigate the transport and mixing processes associated with stellar turbulence, arguably the greatest source of uncertainty in progenitor structure, besides mass loss, at the time of core collapse. An effort to use precision obse… ▽ More Issues concerning the structure and evolution of core collapse progenitor stars are discussed with an emphasis on interior evolution. We describe a program designed to investigate the transport and mixing processes associated with stellar turbulence, arguably the greatest source of uncertainty in progenitor structure, besides mass loss, at the time of core collapse. An effort to use precision observations of stellar parameters to constrain theoretical modeling is also described. △ Less

Submitted 2 June, 2010; originally announced June 2010.

Comments: Proceedings for invited talk at High Energy Density Laboratory Astrophysics conference, Caltech, March 2010. Special issue of Astrophysics and Space Science, submitted for peer review: 7 pages, 3 figures

Some Properties of the Kinetic Energy Flux and Dissipation in Turbulent Stellar Convection Zones

Authors: Casey Meakin, David Arnett

Abstract: We investigate simulated turbulent flow within thermally driven stellar convection zones. Different driving sources are studied, including cooling at the top of the convectively unstable region, as occurs in surface convection zones; and heating at the base by nuclear burning. The transport of enthalpy and kinetic energy, and the distribution of turbulent kinetic energy dissipation are studied.… ▽ More We investigate simulated turbulent flow within thermally driven stellar convection zones. Different driving sources are studied, including cooling at the top of the convectively unstable region, as occurs in surface convection zones; and heating at the base by nuclear burning. The transport of enthalpy and kinetic energy, and the distribution of turbulent kinetic energy dissipation are studied. We emphasize the importance of global constraints on shaping the quasi-steady flow characteristics, and present an analysis of turbulent convection which is posed as a boundary value problem that can be easily incorporated into standard stellar evolution codes for deep, efficient convection. Direct comparison is made between the theoretical analysis and the simulated flow and very good agreement is found. Some common assumptions traditionally used to treat quasi-steady turbulent flow in stellar models are briefly discussed. The importance and proper treatment of convective boundaries are indicated. △ Less

Submitted 23 January, 2010; originally announced January 2010.

Comments: Refereed Conference Proceeding, HELAS Workshop on 'Synergies between solar and stellar modelling', Rome, June 2009, to be published in Astrophys. Space Sci

Turbulent Mixing in Stars: Theoretical Hurdles

Authors: W. David Arnett, Casey Meakin

Abstract: A program is outlined, and first results described, in which fully three-dimensional, time dependent simulations of hydrodynamic turbulence are used as a basis for theoretical investigation of the physics of turbulence in stars. The inadequacy of the treatment of turbulent convection as a diffusive process is discussed. A generalization to rotation and magnetohydrodynamics is indicated, as are c… ▽ More A program is outlined, and first results described, in which fully three-dimensional, time dependent simulations of hydrodynamic turbulence are used as a basis for theoretical investigation of the physics of turbulence in stars. The inadequacy of the treatment of turbulent convection as a diffusive process is discussed. A generalization to rotation and magnetohydrodynamics is indicated, as are connection to simulations of 3D stellar atmospheres. △ Less

Submitted 15 December, 2009; originally announced December 2009.

Comments: 5 pages, 1 figure, IAU Symposium 265, 2009

Convection Theory and Sub-photospheric Stratification

Authors: W. David Arnett, Casey Meakin, Patrick A. Young

Abstract: As a step toward a complete theoretical integration of 3D compressible hydrodynamic simulations into stellar evolution, convection at the surface and sub-surface layers of the Sun is re-examined, from a restricted point of view, in the language of mixing-length theory (MLT) . Requiring that MLT use a hydrodynamically realistic dissipation length gives a new constraint on solar models. While the… ▽ More As a step toward a complete theoretical integration of 3D compressible hydrodynamic simulations into stellar evolution, convection at the surface and sub-surface layers of the Sun is re-examined, from a restricted point of view, in the language of mixing-length theory (MLT) . Requiring that MLT use a hydrodynamically realistic dissipation length gives a new constraint on solar models. While the stellar structure which results is similar to that obtained by YREC and Garching models, the theoretical picture differs. A new quantitative connection is made between macro-turbulence, micro-turbulence, and the convective velocity scale at the photosphere, which has finite values. The "geometric parameter" in MLT is found to correspond more reasonably with the size of the strong downward plumes which drive convection (Stein and Nordlund 1998), and thus has a physical interpretation even in MLT. Use of 3D simulations of both adiabatic convection and stellar atmospheres will allow the determination of the dissipation length and the geometric parameter (i.e., the entropy jump), with no astronomical calibration. A physically realistic treatment of convection in stellar evolution will require additional modifications beyond MLT, including effects of kinetic energy flux, entrainment (the most dramatic difference from MLT found by Meakin and Arnett 2007), rotation, and magnetic fields (Balbus 2009}. △ Less

Submitted 5 October, 2009; originally announced October 2009.

Comments: 16 pages, 5 figures, submitted to ApJ

Journal ref: Astrophys.J.710:1619-1626,2010

Initiation of the detonation in the gravitationally confined detonation model of Type Ia supernovae

Authors: Ivo R. Seitenzahl, Casey A. Meakin, Don Q. Lamb, James W. Truran

Abstract: We study the initiation of the detonation in the gravitationally confined detonation (GCD) model of Type Ia supernovae (SNe Ia). Initiation of the detonation occurs spontaneously in a region where the length scale of the temperature gradient extending from a flow (in which carbon burning is already occurring) into unburned fuel is commensurate to the range of critical length scales which have be… ▽ More We study the initiation of the detonation in the gravitationally confined detonation (GCD) model of Type Ia supernovae (SNe Ia). Initiation of the detonation occurs spontaneously in a region where the length scale of the temperature gradient extending from a flow (in which carbon burning is already occurring) into unburned fuel is commensurate to the range of critical length scales which have been derived from 1D simulations that resolve the initiation of a detonation. By increasing the maximum resolution in a truncated cone that encompasses this region, beginning somewhat before initiation of the detonation occurs, we successfully simulate in situ the first gradient-initiated detonation in a whole-star simulation. The detonation emerges when a compression wave overruns a pocket of fuel situated in a Kelvin-Helmholtz cusp at the leading edge of the inwardly directed jet of burning carbon. The compression wave pre-conditions the temperature in the fuel in such a way that the Zel'dovich gradient mechanism can operate and a detonation ensues. We explore the dependence of the length scale of the temperature gradient on spatial resolution and discuss the implications for the robustness of this detonation mechanism. We find that the time and the location at which initiation of the detonation occurs varies with resolution. In particular, initiation of a detonation had not yet occurred in our highest resolution simulation by the time we ended the simulation because of the computational demand it required. We suggest that the turbulent shear layer surrounding the inwardly directed jet provides the most favorable physical conditions, and therefore the most likely location, for initiation of a detonation in the GCD model. △ Less

Submitted 19 May, 2009; originally announced May 2009.

Comments: 28 pages, 12 figures, 1 table, accepted to ApJ

Journal ref: Astrophys.J.700:642-653,2009

Spontaneous Initiation of Detonations in White Dwarf Environments: Determination of Critical Sizes

Authors: Ivo R. Seitenzahl, Casey a. Meakin, Dean M. Townsley, Don Q. Lamb, James W. Truran

Abstract: Some explosion models for Type Ia supernovae (SN Ia), such as the gravitationally confined detonation (GCD) or the double detonation sub-Chandrasekhar (DDSC) models, rely on the spontaneous initiation of a detonation in the degenerate C/O material of a white dwarf. The length scales pertinent to the initiation of the detonation are notoriously unresolved in multi-dimensional stellar simulations,… ▽ More Some explosion models for Type Ia supernovae (SN Ia), such as the gravitationally confined detonation (GCD) or the double detonation sub-Chandrasekhar (DDSC) models, rely on the spontaneous initiation of a detonation in the degenerate C/O material of a white dwarf. The length scales pertinent to the initiation of the detonation are notoriously unresolved in multi-dimensional stellar simulations, prompting the use of results of 1D simulations at higher resolution, such as the ones performed for this work, as guidelines for deciding whether or not conditions reached in the higher dimensional full star simulations successfully would lead to the onset of a detonation. Spontaneous initiation relies on the existence of a suitable gradient in self-ignition (induction) times of the fuel, which we set up with a spatially localized non-uniformity of temperature -- a hot spot. We determine the critical (smallest) sizes of such hot spots that still marginally result in a detonation in white dwarf matter by integrating the reactive Euler equations with the hydrodynamics code FLASH. We quantify the dependences of the critical sizes of such hot spots on composition, background temperature, peak temperature, geometry, and functional form of the temperature disturbance, many of which were hitherto largely unexplored in the literature. We discuss the implications of our results in the context of modeling of SNe Ia. △ Less

Submitted 23 January, 2009; originally announced January 2009.

Comments: 43 pages, 12 figures, 12 tables

Journal ref: Astrophys.J.696:515-527,2009

Turbulent Convection in Stellar Interiors. II. The Velocity Field

Authors: David Arnett, Casey Meakin, P. A. Young

Abstract: We analyze stellar convection with the aid of 3D hydrodynamic simulations, introducing the turbulent cascade into our theoretical analysis. We devise closures of the Reynolds-decomposed mean field equations by simple physical modeling of the simulations (we relate temperature and density fluctuations via coefficients); the procedure (CABS, Convection Algorithms Based on Simulations) is terrestri… ▽ More We analyze stellar convection with the aid of 3D hydrodynamic simulations, introducing the turbulent cascade into our theoretical analysis. We devise closures of the Reynolds-decomposed mean field equations by simple physical modeling of the simulations (we relate temperature and density fluctuations via coefficients); the procedure (CABS, Convection Algorithms Based on Simulations) is terrestrially testable and is amenable to systematic improvement. We develop a turbulent kinetic energy equation which contains both nonlocal and time dependent terms, and is appropriate if the convective transit time is shorter than the evolutionary time scale. The interpretation of mixing-length theory (MLT) as generally used in astrophysics is incorrect; MLT forces the mixing length to be an imposed constant. Direct tests show that the damping associated with the flow is that suggested by Kolmogorov. The eddy size is approximately the depth of the convection zone, and this dissipation length corresponds to the "mixing length". New terms involving local heating by turbulent dissipation should appear in the stellar evolution equations. The enthalpy flux ("convective luminosity") is directly connected to the buoyant acceleration, and hence the velocity scale. MLT tends to systematically underestimate this velocity scale. Quantitative comparison with a variety of 3D simulations reveals a previously recognized consistency. Examples of application to stellar evolution will be presented in subsequent papers in this series. △ Less

Submitted 25 November, 2008; v1 submitted 9 September, 2008; originally announced September 2008.

Comments: 47 pages, 7 figures, accepted by ApJ

Journal ref: Astrophys.J.690:1715-1729,2009

Study of the Detonation Phase in the Gravitationally Confined Detonation Model of Type Ia Supernovae

Authors: Casey A. Meakin, Ivo Seitenzahl, Dean Townsley, George C. Jordan IV, James Truran, Don Lamb

Abstract: We study the gravitationally confined detonation (GCD) model of Type Ia supernovae through the detonation phase and into homologous expansion. In the GCD model, a detonation is triggered by the surface flow due to single point, off-center flame ignition in carbon-oxygen white dwarfs. The simulations are unique in terms of the degree to which non-idealized physics is used to treat the reactive fl… ▽ More We study the gravitationally confined detonation (GCD) model of Type Ia supernovae through the detonation phase and into homologous expansion. In the GCD model, a detonation is triggered by the surface flow due to single point, off-center flame ignition in carbon-oxygen white dwarfs. The simulations are unique in terms of the degree to which non-idealized physics is used to treat the reactive flow, including weak reaction rates and a time dependent treatment of material in nuclear statistical equilibrium (NSE). Careful attention is paid to accurately calculating the final composition of material which is burned to NSE and frozen out in the rapid expansion following the passage of a detonation wave over the high density core of the white dwarf; and an efficient method for nucleosynthesis post-processing is developed which obviates the need for costly network calculations along tracer particle thermodynamic trajectories. Observational diagnostics are presented for the explosion models, including abundance stratifications and integrated yields. We find that for all of the ignition conditions studied here, a self regulating process comprised of neutronization and stellar expansion results in final \iso{Ni}{56} masses of $\sim$1.1\msun. But, more energetic models result in larger total NSE and stable Fe peak yields. The total yield of intermediate mass elements is $\sim0.1$\msun and the explosion energies are all around 1.5$\times10^{51}$ ergs. The explosion models are briefly compared to the inferred properties of recent Type Ia supernova observations. The potential for surface detonation models to produce lower luminosity (lower \iso{Ni}{56} mass) supernovae is discussed. △ Less

Submitted 30 June, 2008; originally announced June 2008.

Comments: 43 pages, 4 tables, 20 figures -- submitted to ApJ

Journal ref: Astrophys.J.693:1188-1208,2009

Hydrodynamic Processes in Massive Stars

Authors: Casey A. Meakin

Abstract: The hydrodynamic processes operating within stellar interiors are far richer than represented by the best stellar evolution model available. Although it is now widely understood, through astrophysical simulation and relevant terrestrial experiment, that many of the basic assumptions which underlie our treatments of stellar evolution are flawed, we lack a suitable, comprehensive replacement. This… ▽ More The hydrodynamic processes operating within stellar interiors are far richer than represented by the best stellar evolution model available. Although it is now widely understood, through astrophysical simulation and relevant terrestrial experiment, that many of the basic assumptions which underlie our treatments of stellar evolution are flawed, we lack a suitable, comprehensive replacement. This is due to a deficiency in our fundamental understanding of the transport and mixing properties of a turbulent, reactive, magnetized plasma; a deficiency in knowledge which stems from the richness and variety of solutions which characterize the inherently non-linear set of governing equations. The exponential increase in availability of computing resources, however, is ushering in a new era of understanding complex hydrodynamic flows; and although this field is still in its formative stages, the sophistication already achieved is leading to a dramatic paradigm shift in how we model astrophysical fluid dynamics. We highlight here some recent results from a series of multi-dimensional stellar interior calculations which are part of a program designed to improve our one-dimensional treatment of massive star evolution and stellar evolution in general. △ Less

Submitted 27 June, 2008; originally announced June 2008.

Comments: 10 pages, 4 figures, IAUS 252 Conference Proceeding (Sanya) - "The Art of Modeling Stars in the 21st Century"

A Two-Dimensional MagnetoHydrodynamics Scheme for General Unstructured Grids

Authors: Eli Livne, Luc Dessart, Adam Burrows, Casey A. Meakin

Abstract: We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which already includes radiation-hydrodynamics in various approximations and can be used with arbitrarily moving meshes (ALE). The MHD scheme, which consists of cell… ▽ More We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which already includes radiation-hydrodynamics in various approximations and can be used with arbitrarily moving meshes (ALE). The MHD scheme, which consists of cell-centered magnetic field variables, preserves the nodal finite difference representation of $div(\bB)$ by construction, and therefore any initially divergence-free field remains divergence-free through the simulation. In this paper, we describe the new scheme in detail and present comparisons of VULCAN/2D results with those of the code ZEUS/2D for several one-dimensional and two-dimensional test problems. The code now enables two-dimensional simulations of the collapse and explosion of the rotating, magnetic cores of massive stars. Moreover, it can be used to simulate the very wide variety of astrophysical problems for which multi-D radiation-magnetohydrodynamics (RMHD) is relevant. △ Less

Submitted 18 February, 2007; originally announced February 2007.

Comments: 22 pages, including 11 figures; Accepted to the Astrophysical Journal. Higher resolution figures available at http://zenith.as.arizona.edu/~burrows/mhd-code/

Anelastic and Compressible Simulations of Stellar Oxygen Burning

Authors: Casey A. Meakin, David Arnett

Abstract: In this paper we compare fully compressible (Meakin & Arnett 2006a,b) and anelastic (Kuhlen, Woosley, & Glatzmaier 2003) simulations of stellar oxygen shell burning. It is found that the two models are in agreement in terms of the velocity scale (v_c ~ 1e7 cm/s) and thermodynamic fluctuation amplitudes (e.g., rho'/<rho> ~ 2e-3) in the convective flow. Large fluctuations (~11%) arise in the compr… ▽ More In this paper we compare fully compressible (Meakin & Arnett 2006a,b) and anelastic (Kuhlen, Woosley, & Glatzmaier 2003) simulations of stellar oxygen shell burning. It is found that the two models are in agreement in terms of the velocity scale (v_c ~ 1e7 cm/s) and thermodynamic fluctuation amplitudes (e.g., rho'/<rho> ~ 2e-3) in the convective flow. Large fluctuations (~11%) arise in the compressible model, localized to the convective boundaries, and are due to internal waves excited in stable layers. Fluctuations on the several percent level are also present in the compressible model due to composition inhomogeneities from ongoing entrainment events at the convective boundaries. Comparable fluctuations (with amplitudes greater than ~1%) are absent in the anelastic simulation because they are due to physics not included in that model. We derive an analytic estimate for the expected density fluctuation amplitudes at convective boundaries by assuming that the pressure fluctuations due to internal waves at the boundary, p_w', balance the ram pressure of the convective motions, rho*v_c2. The predicted amplitudes agree well with the simulation data. The good agreement between the anelastic and the compressible solution within the convection zone and the agreement between the stable layer dynamics and analytic solutions to the non-radial wave equation indicate that the compressible hydrodynamic techniques used are robust for the simulated stellar convection model, even at the low Mach n umbers found M~0.01. △ Less

Submitted 9 November, 2006; originally announced November 2006.

Comments: 9 pages(emulateapj), 5 figures, 1 table. Recommended for publication in ApJ with minor revisions. Comments will be warmly welcomed. For high resolution version and related preprints please goto http://celsus.as.arizona.edu/preprints/

Journal ref: Astrophys.J.665:690-697,2007

Turbulent Convection in Stellar Interiors. I. Hydrodynamic Simulation

Authors: Casey A. Meakin, David Arnett

Abstract: (Abridged) We describe the results of three-dimensional (3D) numerical simulations designed to study turbulent convection in the stellar interiors, and compare them to stellar mixing-length theory (MLT). Simulations in 2D are significantly different from 3D, both in terms of flow morphology and velocity amplitude. Convective mixing regions are better predicted using a [dynamic boundary condition… ▽ More (Abridged) We describe the results of three-dimensional (3D) numerical simulations designed to study turbulent convection in the stellar interiors, and compare them to stellar mixing-length theory (MLT). Simulations in 2D are significantly different from 3D, both in terms of flow morphology and velocity amplitude. Convective mixing regions are better predicted using a [dynamic boundary condition] based on the bulk Richardson number than by purely local, static criteria like Schwarzschild or Ledoux. MLT gives a good description of the velocity scale and temperature gradient for a mixing length of $\sim 1.1 H_p$ for shell convection, however there are other important effects that it does not capture near boundaries. Convective "overshooting" is best described as an elastic response by the convective boundary, rather than ballistic penetration of the stable layers by turbulent eddies. We find that the rate at which material entrainment proceeds at the boundaries is consistent with analogous laboratory experiments as well as simulation and observation of terrestrial atmospheric mixing. In particular, the normalized entrainment rate E=$u_E/σ_H$, is well described by a power law dependence on the bulk Richardson number $Ri_B = Δb L/σ_H^2$ for the conditions studied, $20\lesssim Ri_B \lesssim 420$. We find $E = A Ri_B^{-n}$, with best fit values, $\log A = 0.027 \pm 0.38$, and $n = 1.05 \pm 0.21$. We discuss the applicability of these results to stellar evolution calculations. △ Less

Submitted 9 November, 2006; originally announced November 2006.

Comments: 32 pages(emulateapj), 26 figures, 6 tables. Submitted to ApJ. Comments will be warmly welcomed. For version with full resolution figures and related pre-prints goto http://celsus.as.arizona.edu/preprints/

Journal ref: Astrophys.J.667:448-475,2007

Active Carbon and Oxygen Shell Burning Hydrodynamics

Authors: Casey Meakin, David Arnett

Abstract: We have simulated 2.5$\times10^3$ s of the late evolution of a $23 \rm M_\odot$ star with full hydrodynamic behavior. We present the first simulations of a multiple-shell burning epoch, including the concurrent evolution and interaction of an oxygen and carbon burning shell. In addition, we have evolved a 3D model of the oxygen burning shell to sufficiently long times (300 s) to begin to assess… ▽ More We have simulated 2.5$\times10^3$ s of the late evolution of a $23 \rm M_\odot$ star with full hydrodynamic behavior. We present the first simulations of a multiple-shell burning epoch, including the concurrent evolution and interaction of an oxygen and carbon burning shell. In addition, we have evolved a 3D model of the oxygen burning shell to sufficiently long times (300 s) to begin to assess the adequacy of the 2D approximation. We summarize striking new results: (1) strong interactions occur between active carbon and oxygen burning shells, (2) hydrodynamic wave motions in nonconvective regions, generated at the convective-radiative boundaries, are energetically important in both 2D and 3D with important consequences for compositional mixing, and (3) a spectrum of mixed p- and g-modes are unambiguously identified with corresponding adiabatic waves in these computational domains. We find that 2D convective motions are exaggerated relative to 3D because of vortex instability in 3D. We discuss the implications for supernova progenitor evolution and symmetry breaking in core collapse. △ Less

Submitted 16 January, 2006; originally announced January 2006.

Comments: 5 pages, 4 figures in emulateapj format. Accepted for publication in ApJ Letters. High resolution figure version available at http://spinach.as.arizona.edu

Journal ref: Astrophys.J. 637 (2006) L53-L56

Constraints on the Progenitor of Cassiopeia A

Authors: Patrick A. Young, Chris L. Fryer, Aimee Hungerford, David Arnett, Gabriel Rockefeller, F. X. Timmes, Benedict Voit, Casey Meakin, Kristoffer A. Eriksen

Abstract: We compare a suite of 3D explosion calculations and stellar models incorporating advanced physics with observational constraints on the progenitor of Cassiopeia A. We consider binary and single stars from 16 to 40 solar masses with a range of explosion energies and geometries. The parameter space allowed by observations of nitrogen rich high velocity ejecta, ejecta mass, compact remnant mass, an… ▽ More We compare a suite of 3D explosion calculations and stellar models incorporating advanced physics with observational constraints on the progenitor of Cassiopeia A. We consider binary and single stars from 16 to 40 solar masses with a range of explosion energies and geometries. The parameter space allowed by observations of nitrogen rich high velocity ejecta, ejecta mass, compact remnant mass, and 44Ti and 56Ni abundances individually and as an ensemble is considered. A progenitor of 15-25 solar masses which loses its hydrogen envelope to a binary interaction and undergoes an energetic explosion can match all the observational constraints. △ Less

Submitted 29 November, 2005; originally announced November 2005.

Comments: 31 pages, 4 figures, accepted for publication to the Astrophysical Journal

Report number: LA-UR-05-6933

Journal ref: Astrophys.J.640:891-900,2006

Young stars and dust in AFGL437: NICMOS/HST polarimetric imaging of an outflow source

Authors: Casey A. Meakin, Dean C. Hines, Roger I. Thompson

Abstract: We present near infrared broad band and polarimetric images of the compact star forming cluster AFGL437 obtained with the NICMOS instrument aboard HST. Our high resolution images reveal a well collimated bipolar reflection nebulosity in the cluster and allow us to identify WK34 as the illuminating source. The scattered light in the bipolar nebulosity centered on this source is very highly polari… ▽ More We present near infrared broad band and polarimetric images of the compact star forming cluster AFGL437 obtained with the NICMOS instrument aboard HST. Our high resolution images reveal a well collimated bipolar reflection nebulosity in the cluster and allow us to identify WK34 as the illuminating source. The scattered light in the bipolar nebulosity centered on this source is very highly polarized (up to 79%). Such high levels of polarization implies a distribution of dust grains lacking large grains, contrary to the usual dust models of dark clouds. We discuss the geometry of the dust distribution giving rise to the bipolar reflection nebulosity and make mass estimates for the underlying scattering material. We find that the most likely inclination of the bipolar nebulosity, south lobe inclined towards Earth, is consistent with the inclination of the large scale CO molecular outflow associated with the cluster, strengthening the identification of WK34 as the source powering it. △ Less

Submitted 1 August, 2005; originally announced August 2005.

Comments: 26 pages, 10 figues. Accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.634:1146-1154,2005

The Impact of Hydrodynamic Mixing on Supernova Progenitors

Authors: Patrick A. Young, Casey Meakin, David Arnett, Chris L. Fryer

Abstract: Recent multidimensional hydrodynamic simulations have demonstrated the importance of hydrodynamic motions in the convective boundary and radiative regions of stars to transport of energy, momentum, and composition. The impact of these processes increases with stellar mass. Stellar models which approximate this physics have been tested on several classes of observational problems. In this paper w… ▽ More Recent multidimensional hydrodynamic simulations have demonstrated the importance of hydrodynamic motions in the convective boundary and radiative regions of stars to transport of energy, momentum, and composition. The impact of these processes increases with stellar mass. Stellar models which approximate this physics have been tested on several classes of observational problems. In this paper we examine the implications of the improved treatment on supernova progenitors. The improved models predict substantially different interior structures. We present pre-supernova conditions and simple explosion calculations from stellar models with and without the improved mixing treatment at 23 solar masses. The results differ substantially. △ Less

Submitted 5 July, 2005; originally announced July 2005.

Comments: 12 pages, 2 figures, accepted for publication in the Astrophysical Journal Letters

Report number: LA-UR-05-3961

Journal ref: Astrophys.J.629:L101-L104,2005

Topics in Core-Collapse Supernova Theory

Authors: Adam Burrows, Christian D. Ott, Casey Meakin

Abstract: There are many interesting topics at the intersection of physics and astrophysics we call Supernova Theory. A small subset of them include the origin of pulsar kicks, gravitational radiation signatures of core bounce, and the possible roles of neutrinos and rotation in the mechanism of explosion. In this brief communication we summarize various recent ideas and calculations that bear on these th… ▽ More There are many interesting topics at the intersection of physics and astrophysics we call Supernova Theory. A small subset of them include the origin of pulsar kicks, gravitational radiation signatures of core bounce, and the possible roles of neutrinos and rotation in the mechanism of explosion. In this brief communication we summarize various recent ideas and calculations that bear on these themes. △ Less

Submitted 24 September, 2003; originally announced September 2003.

Comments: 10 pages, 2 color figures, to be published in the proceedings of ``3-D Signatures in Stellar Explosions: A Workshop honoring J. Craig Wheeler's 60th birthday," held June 10-13, 2003, Austin, Texas, USA

Revealing the Photodissociation Region: Hubble Space Telescope/NICMOS Imaging of NGC 7027

Authors: William B. Latter, Aditya Dayal, John H. Bieging, Casey Meakin, Joseph L. Hora, Douglas M. Kelly, A. G. G. M. Tielens

Abstract: We report results from a Hubble Space Telescope (HST) and Near-Infrared Camera and Multiobject Spectrometer (NICMOS) program to study the distribution of hot neutral (molecular hydrogen) and ionized circumstellar material in the young planetary nebulae NGC 7027. HST/NICMOS provided very high spatial resolution imaging in line and continuum emission, and the stability and large dynamic range need… ▽ More We report results from a Hubble Space Telescope (HST) and Near-Infrared Camera and Multiobject Spectrometer (NICMOS) program to study the distribution of hot neutral (molecular hydrogen) and ionized circumstellar material in the young planetary nebulae NGC 7027. HST/NICMOS provided very high spatial resolution imaging in line and continuum emission, and the stability and large dynamic range needed for investigating detailed structures in the circumstellar material. We present dramatic new images of NGC 7027 that have led to a new understanding of the structure in this important planetary nebula. The central star is clearly revealed, providing near-infrared fluxes that are used to directly determine the stellar temperature very accurately (T = 198,000 K). It is found that the photodissociation layer as revealed by near-infrared molecular hydrogen emission is very thin, and is biconical in shape. The interface region is structured and filamentary, suggesting the existence of hydrodynamic instabilities. We discuss evidence for the presence of one or more highly collimated, off-axis jets that might be present in NGC 7027. The evolutionary future of NGC 7027 is discussed. △ Less

Submitted 24 February, 2000; originally announced February 2000.

Comments: Astrophysical Journal, in press