Source author record

W. David Arnett

W. 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

14works
4topics
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

14 published item(s)

preprint2020arXiv

Biographical Memoir of Stirling Colgate

Stirling Colgate was a remarkably imaginative physicist, an independent thinker with a wide breadth of interests and contagious enthusiasm, a born leader with enduring drive to attack fundamental problems in science. Among his many achievements, he founded the quantitative theory of stellar collapse and supernova explosions, and introduced numerical simulation into the astrophysical toolbox. He brought strong physical intuition to both theory and experiment, in the sciences of nuclear weapons, magnetic and inertial fusion, as well as astrophysics.

preprint2020arXiv

Relative Importance of Convective Uncertainties in Massive Stars

In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called `overshoot', namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary the amount of CBM. We calculate each grid with the initial masses $15$, $20$ and $25\,\rm{M}_\odot$. We present the stellar structure of the models during the hydrogen and helium burning phases. In the latter, we examine the impact on the nucleosynthesis. We find a broadening of the main-sequence with more CBM, which is more in agreement with observations. Furthermore during the core hydrogen burning phase there is a convergence of the convective boundary location due to CBM. The uncertainties of the intermediate convective zone remove this convergence. The behaviour of this convective zone strongly affects the surface evolution of the model, i.e. how fast it evolves red-wards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the $^{12}$C to $^{16}$O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduce and we find differences of up to $70\%$ for the core masses and the total mass of the star.

preprint2018arXiv

Relative Importance of Convective Uncertainties

Convection plays a key role in the evolution of stars due to energy transport and mixing of composition. Despite its importance, this process is still not well understood. One longstanding conundrum in all 1D stellar evolution codes is the treatment of convective boundaries. In this study we compare two convective uncertainties, the boundary location (Ledoux versus Schwarzschild) and the amount of extra mixing, and their impact on the early evolution of massive stars. With increasing convective boundary mixing (CBM), we find a convergence of the two different boundary locations, a decreasing blue to red super giant ratio and a reduced importance of semiconvection.

preprint2016arXiv

Convective Properties of Rotating Two-Dimensional Core-Collapse Supernova Progenitors

We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in one-dimensional (1D) stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions.

preprint2015arXiv

Beyond Mixing-length Theory: a step toward 321D

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.

preprint2015arXiv

Exploring stellar evolution models of sdB stars using MESA

Stellar evolution calculations have had great success reproducing the observed atmospheric properties of different classes of stars. Recent detections of g-mode pulsations in evolved He burning stars allow a rare comparison of their internal structure with stellar models. Asteroseismology of subdwarf B stars suggests convective cores of $0.22-0.28\,M_\odot$, $\gtrsim 45\,\%$ of the total stellar mass. Previous studies found significantly smaller convective core masses ($\lesssim 0.19\,M_\odot$) at a comparable evolutionary stage. We evolved stellar models with MESA (Modules for Experiments in Stellar Astrophysics) to explore how well the interior structure inferred from asteroseismology can be reproduced by standard algorithms. Our qualitative evolutionary paths, position in the $\log g-T_{\rm{eff}}$ diagram and model timescales are consistent with previous results. SdB masses from our full evolutionary sequences fall within the range of the empirical sdB mass distribution, but are nearly always lower than the median. Using standard MLT with atomic diffusion we find convective core masses of $\sim 0.17-0.18 M_\odot$, averaged over the entire sdB lifetime. We can increase the convective core sizes to be as large as those inferred from asteroseismology, but only for extreme values of the overshoot parameter (overshoot gives numerically unstable and physically unrealistic behavior at the boundary). High resolution three-dimensional (3D) simulations of turbulent convection in stars suggest that the Schwarzschild criterion for convective mixing sytematically {\it underestimates the actual extent of mixing because a boundary layer forms.} Accounting for this would decrease the errors in both sdB total and convective core masses.

preprint2015arXiv

The Three Dimensional Evolution to Core Collapse of a Massive Star

We present the first three dimensional (3D) simulation of the final minutes of iron core growth in a massive star, up to and including the point of core gravitational instability and collapse. We self-consistently capture the development of strong convection driven by violent Si burning in the shell surrounding the iron core. This convective burning builds the iron core to its critical (Chandrasekhar) mass and collapse ensues, driven by electron capture and photodisintegration. The non-spherical structure and motion (turbulent fluctuations) generated by 3D convection is substantial at the point of collapse. We examine the impact of such physically-realistic 3D initial conditions on the core-collapse supernova mechanism using 3D simulations including multispecies neutrino leakage. We conclude that non-spherical progenitor structure should not be ignored, and has a significant and favorable impact on the likelihood for neutrino-driven explosions.

preprint2014arXiv

Chaos and Turbulent Nucleosynthesis Prior to a Supernova Explosion

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.

preprint2012arXiv

The Ring Nebula Around the Blue Supergiant SBW1: Pre-Explosion Snapshot of a SN 1987A Twin

SBW1 is a B supergiant surrounded by a ring nebula that is a nearby twin of SN 1987A. We present images and spectra of SBW1 obtained with HST, Spitzer, and Gemini South. HST images of SBW1 do not exhibit long Rayleigh-Taylor fingers, which are presumed to cause the hotspots in the SN1987A ring, but instead show a geometrically thin clumpy ring. The radial mass distribution and size scales of inhomogeneities in SBW1's ring closely resemble those in the SN1987A ring, but the more complete disk expected to reside at the base of the RT fingers is absent. This structure may explain why portions of the SN1987A ring between the hotspots have not yet brightened after more than 15 years. The model we suggest does not require a fast wind colliding with a previous red supergiant wind. More surprisingly, images of SBW1 also reveal diffuse emission filling the interior of the ring in H-alpha and thermal-IR emission; 190K dust dominates the 8-20 micron luminosity (but contains only 1e-5Msun). Cooler (85K dust resides in the equatorial ring itself (dust mass of 5e-3Msun). Diffuse emission extends inward to 1 arcsecond from the central star, where a paucity of emission suggests an inner hole excavated by the wind. We propose that diffuse emission inside the ring arises from an ionized flow of material photoevaporated from the dense ring, and it prevents the supergiant wind from advancing in the equator. This inner emission could correspond to a structure hypothesized to reside around SN1987A that was never directly detected. We suggest that photoionization can play an important dynamical role in shaping the ring nebula, and we speculate that this might help explain the origin of the polar rings around SN1987A. abridged.

preprint2011arXiv

Towards Realistic Progenitors of Core-Collapse Supernovae

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.

preprint2011arXiv

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

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.

preprint2010arXiv

Time-dependent Turbulence in Stars

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.

preprint2009arXiv

Convection Theory and Sub-photospheric Stratification

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}.

preprint2009arXiv

Turbulent Mixing in Stars: Theoretical Hurdles

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.