Engineering Papers⌕ Search

Engineering topics

Arnett, W. David

Publications and source records attributed to Arnett, W. David.

14 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Pre-nebular Light Curves of SNe I

We compare analytic predictions of supernova light curves with recent high-quality data from SN2011fe (Ia), KSN2011b (Ia), and the Palomar Transient Factory and the La Silla-QUEST variability survey (LSQ) (Ia). Because of the steady, fast cadence of observations, KSN2011b provides unique new information on SNe Ia: the smoothness of the light curve, which is consistent with significant large-scale mixing during the explosion, possibly due to 3D effects (e.g., Rayleigh–Taylor instabilities), and provides support for a slowly varying leakage (mean opacity). For a more complex light curve (SN2008D, SN Ib), we separate the luminosity due to multiple causes and indicate the possibility of a radioactive plume. The early rise in luminosity is shown to be affected by the opacity (leakage rate) for thermal and non-thermal radiation. A general derivation of Arnett’s rule again shows that it depends upon all processes heating the plasma, not just radioactive ones, so that SNe Ia will differ from SNe Ibc if the latter have multiple heating processes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Sirius System and Its Astrophysical Puzzles: Hubble Space Telescope and Ground-based Astrometry

Sirius, the seventh-nearest stellar system, is a visual binary containing the metallic-line A1 V star Sirius A, the brightest star in the sky, orbited in a 50.13 year period by Sirius B, the brightest and nearest white dwarf (WD). Using images obtained over nearly two decades with the Hubble Space Telescope (HST), along with photographic observations covering almost 20 years and nearly 2300 historical measurements dating back to the 19th century, we determine precise orbital elements for the visual binary. Combined with the parallax and the motion of the A component, these elements yield dynamical masses of 2.063±0.023 M{sub ⊙} and 1.018±0.011 M{sub ⊙} for Sirius A and B, respectively. Our precise HST astrometry rules out third bodies orbiting either star in the system, down to masses of ∼15–25 M{sub Jup}. The location of Sirius B in the Hertzsprung–Russell diagram is in excellent agreement with theoretical cooling tracks for WDs of its dynamical mass, and implies a cooling age of ∼126 Myr. The position of Sirius B on the mass–radius plane is also consistent with WD theory, assuming a carbon–oxygen core. Including the pre-WD evolutionary timescale of the assumed progenitor, the total age of Sirius B is about 228 ± 10 Myr. We calculated evolutionary tracks for stars with the dynamical mass of Sirius A, using two independent codes. We find it necessary to assume a slightly subsolar metallicity, of about 0.85 Z{sub ⊙}, to fit its location on the luminosity–radius plane. The age of Sirius A based on these models is about 237–247 Myr, with uncertainties of ±15 Myr, consistent with that of the WD companion. We discuss astrophysical puzzles presented by the Sirius system, including the probability that the two stars must have interacted in the past, even though there is no direct evidence for this and the orbital eccentricity remains high.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Synergies between asteroseismology and three-dimensional simulations of stellar turbulence

Turbulent mixing of chemical elements by convection has fundamental effects on the evolution of stars. The standard algorithm at present, mixing-length theory (MLT), is intrinsically local, and must be supplemented by extensions with adjustable parameters. As a step toward reducing this arbitrariness, we compare asteroseismically inferred internal structures of two Kepler slowly pulsating B stars (SPBs; M∼3.25M{sub ⊙}) to predictions of 321D turbulence theory, based upon well-resolved, truly turbulent three-dimensional simulations that include boundary physics absent from MLT. We find promising agreement between the steepness and shapes of the theoretically predicted composition profile outside the convective region in 3D simulations and in asteroseismically constrained composition profiles in the best 1D models of the two SPBs. The structure and motion of the boundary layer, and the generation of waves, are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

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 20 M {sub ⊙} star. Using the Modules for Experiments in Stellar Astrophysics we construct one-dimensional (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 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 the 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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

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

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

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 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 mass and collapse ensues, driven by electron capture and photodisintegration. The non-spherical structure and motion generated by 3D convection is substantial at the point of collapse, with convective speeds of several hundreds of km s{sup −1}. We examine the impact of such physically realistic 3D initial conditions on the core-collapse supernova mechanism using 3D simulations including multispecies neutrino leakage and find that the enhanced post-shock turbulence resulting from 3D progenitor structure aids successful explosions. We conclude that non-spherical progenitor structure should not be ignored, and should have a significant and favorable impact on the likelihood for neutrino-driven explosions. In order to make simulating the 3D collapse of an iron core feasible, we were forced to make approximations to the nuclear network making this effort only a first step toward accurate, self-consistent 3D stellar evolution models of the end states of massive stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

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 (sdB) stars suggests convective cores of 0.22–0.28 M{sub ⊙}, ≳45% of the total stellar mass. Previous studies found significantly smaller convective core masses (≲0.19 M{sub ⊙}) at a comparable evolutionary stage. We evolved stellar models with Modules for Experiments in Stellar Astrophysics (MESA) to explore how well the interior structures inferred from asteroseismology can be reproduced by standard algorithms. Our qualitative evolutionary paths, position in the logg−T{sub eff} diagram, and model timescales are consistent with previous results. The 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 ∼0.17–0.18 M{sub ⊙}, 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 simulations of turbulent convection in stars suggest that the Schwarzschild criterion for convective mixing systematically 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.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2012

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 A. Significant diversity is seen in the near-maximum-light spectra (approx.2000-3500 A) 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 the uvw lIF250W filter are found to correlate with the B-band light-curve shape parameter .(Delta)m15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., approx. 0.4 mag versus approx. 0.2 mag for those with 0.8 <.(Delta)m15(B) < 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 approx. 0,9 mag and approx. 2.0 mag in the uvwl1F250W 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

Wang, Xiaofeng↗
NASA NTRS1989

The late behavior of supernova 1987A. I - The light curve. II - Gamma-ray transparency of the ejecta

Observations of the late (t = 20-1500 days) bolometric light curve and the gamma-lines and X-rays from supernova 1987A are compared to theoretical models. It is found that 0.073 + or - 0.015 solar masses of freshly synthesized Ni-56 must be present to fit the bolometric light curve. The results place limits on the luminosity and presumed period of the newly formed pulsar/neutron star. In the second half of the paper, the problem of computing the luminosities in gamma-ray lines and in X-rays from supernova 1987A is addressed. High-energy observations suggest the development of large-scale clumping and bubbling of radioactive material in the ejecta. A model is proposed with a hydrogen envelope mass of about 7 solar masses, homologous scale expansion velocities of about 3000 km/s, and an approximately uniform mass distribution.

Arnett, W. David↗
NASA NTRS1989

On relative supernova rates and nucleosynthesis roles

It is shown that the Ni-56-Fe-56 observed in SN 1987A argues that core collapse supernovae may be responsible for more than 50 percent of the iron in the galaxy. Furthermore it is argued that the time averaged rate of thermonuclear driven Type I supernovae may be at least an order of magnitude lower than the average rate of core collapse supernovae. The present low rate of Type II supernovae (below their time averaged rate of approx. 1/10 yr) is either because the past rate was much higher because many core collapse supernovae are dim like SN 1987A. However, even in this latter case they are only an order of magnitude dimmer that normal Type II's due to the contribution of Ni-56 decay to the light curve.

Arnett, W. David↗
NASA NTRS1989

Supernova 1987A

Present understanding of SN 1987A is reviewed. The results of early optical and UV observations of the SN are summarized, and the reasons for the blue color of the progenitor star are theoretically examined. The physical processes involved in the explosion are discussed, and the neutrino burst, the light curve, and high-energy emissions are studied in detail. The spectral evaluation of SN 1987A is addressed, and issues for future study are briefly considered.

Arnett, W. David↗
NASA NTRS1988

On relative supernova rates and nucleosynthesis roles

It is shown that the Ni-56-Fe-56 observed in SN 1987A argues that core collapse supernovae may be responsible for more that 50 percent of the iron in the galaxy. Furthermore it is argued that the time averaged rate of thermonuclear driven Type I supernovae may be at least an order of magnitude lower than the average rate of core collapse supernovae. The present low rate of Type II supernovae (below their time averaged rate of approx. 1/10 yr) is either because the past rate was much higher because many core collapse supernovae are dim like SN 1987A. However, even in this latter case they are only an order of magnitude dimmer that normal Type II's due to the contribution of Ni-56 decay to the light curve.

Arnett, W. David↗