Source author record

David Arnett

David Arnett appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

13 published item(s)

preprint2016arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

Linking 1D Stellar Evolution to 3D Hydrodynamical Simulations

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.

preprint2014arXiv

Preparing for an Explosion: Hydrodynamic Instabilities and Turbulence in Presupernovae

Both observations and direct numerical simulations are discordant with predictions of conventional stellar evolution codes for the latest stages of a massive star's life prior to core collapse. We suggest that the problem lies in the treatment of turbulent convection in these codes, which ignores finite amplitude fluctuations in velocity and temperature, and their nonlinear interaction with nuclear burning. The hydrodynamic instabilities that may arise prompt us to discuss a number of far-reaching implications for the fates of massive stars. In particular, we explore connections to enhanced presupernova mass loss, unsteady nuclear burning and consequent eruptions, swelling of the stellar radius that may trigger violent interactions with a companion star, and potential modifications to the core structure that could dramatically impact calculations of the core-collapse mechanism itself. These modifications may be of fundamental importance to the interpretation of measured isotopic anomalies in meteorites, changing the predictions of both mixing and detailed nucleosynthesis, and of young core-collapse supernova remnants. They may also make possible the development of an early warning system for the detection of impending core collapse.

preprint2013arXiv

The Age and Stellar Parameters of the Procyon Binary System

The Procyon AB binary system (orbital period 40.838 years, a newly-refined determination), is near and bright enough that the component radii, effective temperatures, and luminosities are very well determined, although more than one possible solution to the masses has limited the claimed accuracy. Preliminary mass determinations for each component are available from HST imaging, supported by ground-based astrometry and an excellent Hipparcos parallax; we use these for our preferred solution for the binary system. Other values for the masses are also considered. We have employed the TYCHO stellar evolution code to match the radius and luminosity of the F5 IV-V primary star to determine the system's most likely age as 1.87 +/- 0.13 Gyr. Since prior studies of Procyon A found its abundance indistinguishable from solar, the solar composition of Asplund, Grevesse & Sauval (Z=0.014) is assumed for the HR Diagram flitting. An unsuccessful attempt to fit using the older solar abundance scale of Grevesse & Sauval (Z=0.019) is also reported. For Procyon B, eleven new sequences for the cooling of non-DA white dwarfs have been calculated, to investigate the dependence of the cooling age on (1) the mass, (2) the core composition, (3) the helium layer mass, and (4) heavy-element opacities in the helium envelope. Our calculations indicate a cooling age of 1.19+/-0.11 Gyr, which implies that the progenitor mass of Procyon B was 2.59(+0.44,-0.26) Msun. In a plot of initial vs final mass of white dwarfs in astrometric binaries or star clusters (all with age determinations), the Procyon B final mass lies several sigma below a straight line fit.

preprint2013arXiv

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

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.

preprint2012arXiv

Evidence for Type Ia Supernova Diversity from Ultraviolet Observations with the Hubble Space Telescope

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

preprint2010arXiv

Presupernova Structure of Massive Stars

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.

preprint2010arXiv

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

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.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

preprint2009arXiv

The Reddening Towards Cassiopeia A's Supernova: Constraining the 56Ni Yield

We present new reddening measurements towards the young supernova remnant Cassiopeia A, using two techniques not previously applied to this object. Our observations of the near-infrared [Fe II] 1.257 micron and 1.644 micron lines show the extinction to be highly variable across the remnant, increasing towards the west and the south, consistent with previous radio and X-ray observations. While the absolute value of AV as determined by the [Fe II] lines is uncertain due to conflicting calculations and observations of their intrinsic flux ratio, parts of the remnant without previous optical measurements show comparatively higher reddening. We find AV = 6.2 +/- 0.6 from the broadband shape of the infrared synchrotron emission of a knot within 13" of the expansion center. Given this reddening, the apparent faintness of the nascent supernova, and iron mass constraints from X-ray observations, we estimate an ejected mass of 56Ni of 0.058 - 0.16 MSun. Taken with gamma-ray observations of the 44Ti decay chain, this nickel mass is broadly consistent with the solar 44Ca/56Fe ratio.

preprint2008arXiv

Turbulent Convection in Stellar Interiors. II. The Velocity Field

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.