Engineering Papers⌕ Search

Engineering topics

Arnett, David

Publications and source records attributed to Arnett, David.

Thermonuclear runaways in nova outbursts. 2: Effect of strong, instantaneous, local fluctuations

In an attempt to understand the manner in which nova outbursts are initiated on the surface of a white dwarf, we investigate the effects fluctuations have on the evolution of a thermonuclear runaway. Fluctuations in temperature density, or the composition of material in the burning shell may arise due to the chaotic flow field generated by convection when it occurs, or by the accretion process itself. With the aid of two-dimensional reactive flow calculations, we consider cases where a strong fluctutation in temperature arises during the early, quiescent accretion phase or during the later, more dynamic, explosion phase. In all cases we find that an instantaneous, local temperature fluctuation causes the affected material to become Rayleigh-Taylor unstable. The rapid rise and subsequent expansion of matter immediately cools the hot blob, which prevents the lateral propagation of burning. This suggests that local temperature fluctuations do not play a significant role in directly initiating the runaway, especially during the early stages. However, they may provide an efficient mechanism of mixing core material into the envelope (thereby pre-enriching the fuel for subsequent episodes of explosive hydrogen burning) and of mixing substantial amounts of the radioactive nucleus N-13 into the surface layers, making novae potential gamma-ray sources. This suggests that it is the global not the local, evolution of the core-envelope interface to high temperatures which dominates the development of the runaway. We also present a possible new scenario for the initiation of nova outbursts based on our results.

Shankar, Anurag↗

Convection, nucleosynthesis, and core collapse

We use a piecewise parabolic method hydrodynamics code (PROMETHEUS) to study convective burning in two dimensions in an oxygen shell prior to core collapse. Significant mixing beyond convective boundaries determined by mixing-length theory brings fuel (C-12) into the convective regon, causing hot spots of nuclear burning. Plumes dominate the velocity structure. Finite perturbations arise in a region in which O-16 will be explosively burned to Ni-56 when the star explodes; the resulting instabilities and mixing are likely to distribute Ni-56 throughout the supernova envelope. Inhomogeneities in Y(sub e) may be large enough to affect core collapse and will affect explosive nucleosynthesis. The nature of convective burning is dramatically different from that assumed in one-dimensional simulations; quantitative estimates of nucleosynthetic yields, core masses, and the approach to core collapse will be affected.

Bazan, Grant↗

Oxgen-burning hydrodynamics. 1: Steady shell burning

With new hydrodynamic techniques, the relatively fast evolutionary stages of a star prior to core collapse may be explicitly computed in two spatial dimensions, with a treatment of the microphysics (e.g., nuclear reactions, equation of state, neutrino cooling) which is comparable to typical one-dimensional simulations. The nature of shell oxygen burning in a massive star, prior to core collapse, is used as a first example; it is of particular interest because it is (1) the region in which Ni-56 will be produced by the supernova shock, (2) the region of the 'mass cut', which will separate the collapsed core from the ejected mantle, (3) the site of much of the explosive nucleosynthesis, and (4) a suggested source of symmetry breaking to drive mixing instabilities which were observed in SN 1987A. The nature of the shell burning affects the size of the core which will collapse. The method is illustrated on this test case, and the character of the convection is examined.

Arnett, David↗

The delayed-detonation model of a type Ia supernovae. 1: The deflagration phase

The nature of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is investigated by using two-dimensional numerical hydrodynamics simulations. A new algorithm is used to treat the deflagration front. Assuming that it propagates locally at the laminar flame speed, the deflagration is insufficient to unbind the star. Expansion shuts of the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. The burning front does become wrinkled, but the wavelength of the instability is much larger than the computational grid size and is resolved; this is consistent with previous analysis. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. This paper deals with the deflagration phase, from the onset of burning, through expansion and quenching of the flame, to the first contraction.

Arnett, David↗

The delayed-detonation model of Type Ia supernovae. 2: The detonation phase

The investigation, by use of two-dimensional numerical hydrodynamics simulations, of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is continued. Previously we found that the deflagration is insufficient to unbind the star. Expansion shuts off the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. The burning was allowed to develop into a detonation in these nonspherical models. The detonation grows toward spherical symmetry at late times. At these densities (rho approx. 10(exp 7) to 10(exp 8) g cm(exp -3)), either Ni-56 or nuclei of the Si-Ca group are the dominant products of the burning. The bulk yields are sensitive to the density of the star when the transition to detonation occurs. The relevance of the abundances, velocities, mixing, and total energy release to the theory and interpretation of Type Ia supernovae is discussed.

Arnett, David↗

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian↗

A theory of nonlocal mixing-length convection. I - The moment formalism

A flexible and potentially powerful theory of convection, based on the mixing length picture, is developed to make unbiased self-consistent predictions about overshooting and other complicated phenomena in convection. The basic formalism is set up, and the method's power is demonstrated by showing that a simplified version of the theory reproduces all the standard results of local convection. The second-order equations of the theory are considered in the limit of a steady state and vanishing third moments, and it is shown that they reproduce all the standard results of local mixing-length convection. There is a particular value of the superadiabatic gradient, below which the only possible steady state of a fluid is nonconvecting. Above this critical value, a fluid is convectively unstable. Two distinct regimes of convection, which are identified as efficient and inefficient convection, are determined.

Grossman, Scott A.↗

The ring around SN1987A

Stars in the 9-40 solar mass range play a prominent role in the hydrodynamical and chemical evolution of galaxies. Their stellar winds and supernova explosions are believed to create the hot component of the interstellar medium (ISM). In some galactic disks, the kiloparsec sized super bubbles formed around clusters of massive stars may blow out of the disk plane and release hot, metal enriched gas into the galaxy's halo. Additionally, the expanding shock front of a super bubble in the disk may trigger additional star formation. Furthermore, similar processes probably drive the galactic winds associated with star burst nuclei that enrich the intracluster and intergalactic mediums. Nonetheless, the explosion of a blue super giant in the Large Magellanic Cloud (LMC), SN1987A, illuminated the incompleteness of our understanding of massive stars. Evolutionary models of massive stars do not synthesize the observed super giant populations in either the Milky Way or LMC. Our modeling of the formation of SN1987A's ring will improve our knowledge of both the post-main-sequence evolution of massive stars and their coupling to the ISM in galaxies.

Martin, Crystal L.↗

Thermonuclear runaways in nova outbursts

Results of exploratory, two-dimensional numerical calculations of a local thermonuclear runaway on the surface of a white dwarf are reported. It is found that the energy released by the runaway can induce a significant amount of vorticity near the burning region. Such mass motions account naturally for mixing of core matter into the envelope during the explosion. A new mechanism for the lateral spread of nuclear burning is also discussed.

Shankar, Anurag↗

Massive star evolution and SN 1987A

The evolution of massive stars through hydrogen and helium burning is addressed. A set of stellar evolutionary sequences for mass/solar mass of 15, 20, and 25, and metallicity of 0.002, 0.005, 0.007, 0.010, and 0.20 are presented; semiconvection is restricted to operating slower than the local thermal time scale. Using these sequences, simple models of the massive star content of the LMC are found to agree moderately well with the new observational data of Fitzpatrick and Garmany (1990). LMC supergiants were detected only in their postmain-sequence phases, so that 5-10 times more massive stars are there but not identified as such. It is argued that SN 1987A exhibits the normal evolution of a single star of about 20 solar mases having LMC abundances. Despite the variety of envelope behavior, the structure of the core at collapse is rather similar for the stars of a given mass. Variations due to different rates of mass loss are likely to be larger than those due to composition.

Arnett, David↗

Instabilities in SN1987A and other supernovae

While giving a remarkably good description of many aspects of SN1987A, calculations assuming spherical symmetry have a number of flaws. Many of these problems naturally disappear in the development of a 2-D calculation, as nonspherical instabilities grow. This nonspherical behavior, made evident in SN1987A, has implications for other types of supernovae.

Arnett, David↗