• ### Progenitors of Core-Collapse Supernovae(1802.07503)

Feb. 21, 2018 astro-ph.SR
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.
• ### A. G. W. Cameron 1925-2005, Biographical Memoir, National Academy of Sciences(1708.05429)

Alastair Graham Walker Cameron was an astrophysicist and planetary scientist of broad interests and exceptional originality. A founder of the field of nuclear astrophysics, he developed the theoretical understanding of the chemical elementsâ origins and made pioneering connections between the abundances of elements in meteorites to advance the theory that the Moon originated from a giant impact with the young Earth by an object at least the size of Mars. Cameron was an early and persistent exploiter of computer technology in the theoretical study of complex astronomical systemsâincluding nuclear reactions in supernovae, the structure of neutron stars, and planetary collisions.
• ### 3D Hydrodynamic Simulations of Carbon Burning in Massive Stars(1610.05173)

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\'eclet 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}^{-\alpha}$, $0.5\lesssim \alpha \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.
• ### Linking 1D Evolutionary to 3D Hydrodynamical Simulations of Massive Stars(1601.01572)

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.
• ### Toward a consistent use of overshooting parametrizations in 1D stellar evolution codes(1506.03100)

June 26, 2015 astro-ph.SR
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.
• ### Preparing for an Explosion: Hydrodynamic Instabilities and Turbulence in Presupernovae(1307.5035)

Feb. 5, 2014 astro-ph.SR, astro-ph.HE
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.
• ### The Age and Stellar Parameters of the Procyon Binary System(1305.0587)

May 2, 2013 astro-ph.SR
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.
• ### Turbulent convection in stellar interiors. III. Mean-field analysis and stratification effects(1212.6365)

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.
• 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.
• ### Presupernova Structure of Massive Stars(1006.0513)

June 2, 2010 astro-ph.SR
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.
• ### Some Properties of the Kinetic Energy Flux and Dissipation in Turbulent Stellar Convection Zones(1001.4138)

Jan. 23, 2010 astro-ph.SR
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.
• 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.
• ### Turbulent Convection in Stellar Interiors. II. The Velocity Field(0809.1625)

Nov. 25, 2008 astro-ph
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.
• ### Turbulent Convection in Stellar Interiors. I. Hydrodynamic Simulation(astro-ph/0611315)

Nov. 9, 2006 astro-ph
(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/\sigma_H$, is well described by a power law dependence on the bulk Richardson number $Ri_B = \Delta b L/\sigma_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.
• ### Anelastic and Compressible Simulations of Stellar Oxygen Burning(astro-ph/0611317)

Nov. 9, 2006 astro-ph
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.
• ### Active Carbon and Oxygen Shell Burning Hydrodynamics(astro-ph/0601348)

Jan. 16, 2006 astro-ph
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.
• ### Constraints on the Progenitor of Cassiopeia A(astro-ph/0511806)

Nov. 29, 2005 astro-ph
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.
• ### A Model For the Formation of High Density Clumps in Proto-Planetary Nebulae(astro-ph/0309545)

Sept. 19, 2003 astro-ph
The detection of NaCl at large radii in the Egg Nebula, CRL 2688, requires densities of 10^7 - 10^8 cm^-3 in a thick shell of r ~ a few X 10^17 cm. To explain these results, a mechanism is needed for producing high densities at a considerable distance from the central star. In two dimensional simulations of the interaction of the fast wind with an ambient medium, the material becomes thermally unstable. The resulting clumps can achieve the requisite conditions for NaCl excitation. We present 2D models with simple physics as proof-of-principle calculations to show that the clumping behavior is robust. Clumping is a natural outcome of cooling in the colliding wind model and comparable to that inferred from observations.
• ### Stellar Hydrodynamics in Radiative Regions(astro-ph/0306237)

June 11, 2003 astro-ph
We present an analysis of the response of a radiative region to waves generated by a convective region of the star; this wave treatment of the classical problem of overshooting'' gives extra mixing relative to the treatment traditionally used in stellar evolutionary codes. The interface between convectively stable and unstable regions is dynamic and nonspherical, so that the nonturbulent material is driven into motion, even in the absence of penetrative overshoot.'' These motions may be described by the theory of nonspherical stellar pulsations, and are related to motion measured by helioseismology. Multi-dimensional numerical simulations of convective flow show puzzling features which we explain by this simplified physical model. Gravity waves generated at the interface are dissipated, resulting in slow circulation and mixing seen outside the formal convection zone. The approach may be extended to deal with rotation and composition gradients. Tests of this description in the stellar evolution code TYCHO produce carbon stars on the asymptotic giant branch (AGB), an isochrone age for the Hyades and three young clusters with lithium depletion ages from brown dwarfs, and lithium and beryllium depletion consistent with observations of the Hyades and Pleiades, all without tuning parameters. The insight into the different contributions of rotational and hydrodynamic mixing processes could have important implications for realistic simulation of supernovae and other questions in stellar evolution.
• ### Observational Tests and Predictive Stellar Evolution(astro-ph/0103390)

March 23, 2001 astro-ph
We compare eighteen binary systems with precisely determined radii and masses from 23 to 1.1 M_sol, and stellar evolution models produced with our newly revised code TYCHO. Overshooting'' and rotational mixing were suppressed in order to establish a baseline for isolating these and other hydrodynamic effects. Acceptable coeval fits are found for sixteen pairs without optimizing for heavy element or helium abundance. The precision of these tests is limited by the accuracies of the observed effective temperatures. High dispersion spectra and detailed atmospheric modeling should give more accurate effective temperatures and heavy element abundances. PV Cas, a peculiar early A system, EK Cep B, a known post-T Tauri star, and RS Cha, a member of a young OB association, are matched by pre-main sequence models. Predicted mass loss agrees with upper limits from IUE for CW Cep A and B. Relatively poor fits are obtained for binaries having at least one component in the mass range 1.7 < M/M_sol <2.6, whose evolution is sensitive to mixing. These discrepancies are robust and consistent with additional mixing in real stars. The predicted apsidal motion implies that massive star models are systematically less centrally condensed than the real stars. If these effects are due to overshooting, then the overshooting parameter alpha_OV increases with stellar mass. The apsidal motion constants are controlled by radiative opacity under conditions close to those directly measured in laser experiments, making this test more stringent than possible before.
• ### The Role of Mixing in Astrophysics(astro-ph/9909031)

Sept. 1, 1999 astro-ph
The role of hydrodynamic mixing in astrophysics is reviewed, emphasizing connections with laser physics experiments and inertial confinement fusion (ICF). Computer technology now allows two dimensional (2D) simulations, with complex microphysics, of stellar hydrodynamics and evolutionary sequences, and holds the promise for 3D. Careful validation of astrophysical methods, by laboratory experiment, by critical comparison of numerical and analytical methods, and by observation are necessary for the development of simulation methods with reliable predictive capability. Recent and surprising results from isotopic patterns in presolar grains, 2D hydrodynamic simulations of stellar evolution, and laser tests and computer simulations of Richtmeyer-Meshkov and Rayleigh-Taylor instabilities are discussed, and related to stellar evolution and supernovae.
• ### Universality in SNIae and the Phillips Relation(astro-ph/9908169)

Aug. 16, 1999 astro-ph
The use of supernovae of Type Ia for the determination of accurate distances rests upon the empirical Phillips relation, in which the brightest events are the broadest in time. Implications of new data upon the homogeneity of light curves under the operation of a stretch in time, of the parabolic luminosity increase at the earliest times, and of the time from explosion to maximum light are discussed. The early luminosity is in excellent agreement with predictions (Arnett 1982), and the lack of prominent higher modes of diffusion constrain progenitor and explosion models. Difficulties with reproducing the observed rise time are restricted to radiative transfer models (e.g., Hoeflich and Khokhlov), and probably due to an overestimate of thermal photon escape due to inadequate line lists. Because of the strong dependence of luminosity on 56Ni mass, some simple models can give a Phillips relation of the correct sense.
• ### Explosive Nucleosynthesis: Prospects(astro-ph/9908171)

Aug. 16, 1999 astro-ph
Explosive nucleosynthesis is a combination of the nuclear physics of thermonuclear reactions, and the hydrodynamics of the plasma in which the reactions occur. It depends upon the initial conditions---the stellar evolution up to the explosive instability, and the nature of the explosion mechanism. Some key issues for explosive nucleosynthesis are the interaction of burning with hydrodynamics, the degree of microscopic mixing in convective zones, and the breaking of spherical symmetry by convection and rotation. Recent experiments on high intensity lasers provides new opportunities for laboratory testing of astrophysical hydrodynamic codes. Implications of SN1987A, SN1998bw (GRB980425?), and eta Carina are discussed, as well as the formation of black holes or neutron stars.
• ### Two-Dimensional Hydrodynamics of Pre-Core Collapse: Oxygen Shell Burning(astro-ph/9702239)

Feb. 27, 1997 astro-ph
By direct hydrodynamic simulation, using the Piecewise Parabolic Method (PPM) code PROMETHEUS, we study the properties of a convective oxygen burning shell in a SN 1987A progenitor star prior to collapse. The convection is too heterogeneous and dynamic to be well approximated by one-dimensional diffusion-like algorithms which have previously been used for this epoch. Qualitatively new phenomena are seen. The simulations are two-dimensional, with good resolution in radius and angle, and use a large (90-degree) slice centered at the equator. The microphysics and the initial model were carefully treated. Many of the qualitative features of previous multi-dimensional simulations of convection are seen, including large kinetic and acoustic energy fluxes, which are not accounted for by mixing length theory. Small but significant amounts of carbon-12 are mixed non-uniformly into the oxygen burning convection zone, resulting in hot spots of nuclear energy production which are more than an order of magnitude more energetic than the oxygen flame itself. Density perturbations (up to 8%) occur at the `edges' of the convective zone and are the result of gravity waves generated by interaction of penetrating flows into the stable region. Perturbations of temperature and electron fraction at the base of the convective zone are of sufficient magnitude to create angular inhomogeneities in explosive nucleosynthesis products, and need to be included in quantitative estimates of yields. Combined with the plume-like velocity structure arising from convection, the perturbations will contribute to the mixing of nickel-56 throughout supernovae envelopes. Runs of different resolution, and angular extent, were performed to test the robustness of these