On the numerical stability of computations of stellar evolution
Stellar evolution computations numerical instability due to coupling between hydrostatic equilibrium and thermal processes
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Stellar evolution computations numerical instability due to coupling between hydrostatic equilibrium and thermal processes
Stellar evolution, discussing mass loss to interstellar medium, galactic distances determination and supernova explosions
Stellar evolutionary phases from gravitational contraction to exhaustion of nuclear fuel
Aspects of normal stellar evolution are discussed together with evolution near the main sequence, stellar evolution from main sequence to white dwarf or carbon ignition, the structure of massive main-sequence stars, and problems of stellar stability and stellar pulsation. Other subjects considered include variable stars, white dwarfs, close binaries, novae, early supernova luminosity, neutron stars, the photometry of field horizontal-branch stars, and stellar opacity. Transport mechanisms in stars are examined together with thermonuclear reactions and nucleosynthesis, the instability problem in nuclear burning shells, stellar coalescence, and intense magnetic fields in astrophysics. Individual items are announced in this issue.
Stellar evolution and evolutionary sequences, use of high speed computer methods and application to age determination of stars
PreHayashi phase of stellar evolution, discussing protostellar disks /stellisks/ formation from flattened fragments of collapsed interstellar cloud under turbulent viscosity
A study of stellar chromospheres based on the internal structure of particular stars is presented. Used are complex flow diagrams of the linkage paths between mass loss, angular momentum loss, magnetic field from the turbulent dynamo and its relations to differential rotations and the convection zone, and stellar evolution.
The stellar-evolution program developed by Sweigart (1973) was used to compute main sequence turnoff, red giant sequences, and horizontal branch sequences for many values of helium abundance and heavy-element abundance. The immediate aim has been to produce a self-consistent set of evolutionary sequences for wide ranges in the composition and mass. The long-term aim involves application to cluster HR diagrams and the integrated properties of galaxies. Turnoff sequences, red giant sequences, and horizontal branch sequences are discussed.
Stellar evolution model compared with observation to determine interior structure of real stars
The present conference discusses theoretical and observational views of star formation, spectroscopic constraints on the evolution of massive stars, very low mass stars and brown dwarfs, asteroseismology, globular clusters as tests of stellar evolution, observational tests of stellar evolution, and mass loss from cool evolved giant stars. Also discussed are white dwarfs and hot subdwarfs, neutron stars and black holes, supernovae from single stars, close binaries with evolved components, accretion disks in interacting binaries, supernovae in binary systems, stellar evolution and galactic chemical evolution, and interacting binaries containing compact components.
Disks of gas accreting onto supermassive black holes may host numerous stellar-mass objects, formed within the disk or captured from a nuclear star cluster. We present a simplified model of stellar evolution in these dense environments, which exhibits exceptional agreement with full stellar evolution calculations at a minuscule fraction of the cost. Although the model presented here is limited to stars burning hydrogen in their cores, it is sufficient to determine the evolutionary fate of disk-embedded stars: whether they proceed to later stages of nuclear burning and leave behind a compact remnant, reach a quasi-steady state where mass loss and accretion balance one another, or whether accretion proceeds faster than stellar structure can adjust, causing a runaway. We highlight how various disk parameters and phenomena such as gap opening affect stellar evolution outcomes. We also highlight how our model can accommodate time-varying conditions, such as those experienced by a star on an eccentric orbit, and can couple to N -body integrations. This model will enable more detailed studies of stellar populations and their interaction with accretion disks than have previously been possible.
Explore the source record for details and available documents.
A large number of stellar evolution models with (Fe/H) = -2.3 and -3.3 have been calculated in order to determine the primordial Li-7 abundance and to test current stellar evolution models by a comparison to the extensive database of accurate Li abundances in extremely metal-poor halo stars observed by Thorburn (1994). Standard models with gray atmospheres do a very good job of fitting the observed Li abundances in stars hotter than approximately 5600 K. They predict a primordial. Li-7 abundance of log N(Li) = 2.24 +/- 0.03. Models which include microscopic diffusion predict a downward curvature in the Li-7 destruction isochrones at hot temperatures which is not present in the observations. Thus, the observations clearly rule out models which include uninhibited microscopic diffusion of Li-7 from the surface of the star. Rotational mixing inhibits the microscopic diffusion and the (Fe/H) = -2.28 stellar models which include both diffusion and rotational mixing provide an excellent match to the mean trend in T(sub eff) which is present in the observations. Both the plateau stars and the heavily depleted cool stars are well fit by these models. The rotational mixing leads to considerable Li-7 depletion in these models and the primordial Li-7 abundance inferred from these models is log N(Li) = 3.08 +/- 0.1. However, the (Fe/H) = -3.28 isochrones reveal problems with the combined models. These isochrones predict a trend of decreasing log N(Li) with increasing T(sub eff) which is not present in the observations. Possible causes for this discrepancy are discussed.
Numerical stability in computations of stellar evolution taking into account hydrostatic equilibrium and thermal processes
Analytic models of stellar evolution, describing homogeneous and inhomogeneous stages, energy conservation and transport radiative and convective transfer mechanisms, etc
The development of appropriate computer programs has made it possible to conduct studies of stellar evolution which are more detailed and accurate than the investigations previously feasible. However, the use of such programs can also entail some serious drawbacks which are related to the time and expense required for the work. One approach for overcoming these drawbacks involves the employment of simplified stellar evolution codes which incorporate the essential physics of the problem of interest without attempting either great generality or maximal accuracy. Rappaport et al. (1982) have developed a simplified code to study the evolution of close binary stellar systems composed of a collapsed object and a low-mass secondary. The present investigation is concerned with a more general, but still simplified, technique for calculating the evolution of close binary systems with collapsed binaries and mass-losing secondaries.
In this paper a program is described for stellar evolution computations, based on the Henyey method. The program shall be used for very different phases of stellar evolution. Therefore, the main program for the solution of the time dependent basic equations is kept as general as possible, while all assumptions about special physical conditions in the star occur only in subprograms. Hence the range of applicability of the program can easily be extended. The system of difference equations and details of the integration method used in the program are given in two appendices. For a first test the program has been applied to the post-main sequence evolution of a star of 7 M; the results will appear in part II and III. Some technical details of these calculations are given and discussed here. The program turned out to be applicable to calculations of the evolution near the main sequence, of phases of core contraction and of nuclear burning in shell sources.