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.
The information necessary for application of the STARS-2P (shell theory automated for rotational structures-2(plasticity)) program is presented. This addition to the STARS system of programs retains the basic fractures characteristic of the system. For vol. 1, see N75-27419.
If an axion of mass between about 10(exp -3) eV and 1 eV exists, axion emission would have significantly affected the cooling of the nascent neutron star associated with SN 1987A. For an axion of mass less than about 10(exp -2) eV, axions produced deep inside the neutron star simply stream out; in a previous paper this case has been addressed. Remarkably, for an axion of mass greater than about 10(exp -2) eV axions would, like neutrinos, have a mean-free path that is smaller than the size of a neutron star, and thus would become 'trapped' and radiated from an axion sphere. In this paper the trapping regime is treated by using numerical models of the initial cooling of a hot neutron star that incorporate a leakage approximation scheme for axion-energy transport. The axion opacity is computed due to inverse nucleon-nucleon, axion bremsstrahlung, and numerical models are used to calculate the integrated axion luminosity, the temperature of the axion sphere, and the effect of axion emission on the neutrino bursts detected by the Kamiokande 2 (K2) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the axion mass, the stronger the trapping and the smaller the axion luminosity. The earlier estimate is confirmed and refined of the axion mass above which trapping is so strong that axion emission does not significantly affect the neutrino burst. Based upon the neutrino-burst duration--the most sensitive barometer of axion cooling--it is concluded that for an axion mass of greater than about 0.3 eV, axion emission would not have had a significant effect on the neutrino bursts detected by K2 and IMB. The present work, together with the previous work, strongly suggests that an axion with mass in the interval 10(exp -3) eV to 0.3 eV is excluded by SN 1987A.
Observations of gamma-2 Vel and zeta Pup are presented for the spectral region 5800-11,000 per cm, at resolutions of 2 and 4 per cm, respectively. This constitutes the first reported spectroscopy of a Wolf-Rayet star beyond 1.1 microns. The new spectral region contained numerous emission features for the WC8 component of gamma-2 Vel, but only the P-beta line in absorption for the O4f star zeta Pup. Line identifications are given for gamma-2 Vel and are shown to be consistent with visual spectroscopy. The species identified are H, He I-II, and C II-IV.
Time-independent spherical accretion by a neutron star is studied using general relativistic radiation hydrodynamics. Numerical integrations of the flow equations are presented. These show that when the luminosity is sufficiently close to (but below) the Eddington limit, the flow velocity increases with decreasing radius far from the neutron star, reaches a maximum at an intermediate radius, and decreases at small radii. A large fraction of the binding energy of the flow is transferred to the radiation through scattering before the flow strikes the surface of the neutron star. Following Miller's treatment of accretion at luminosites near the Eddington limit (which neglected general relativistic effects), analytic approximations for the decelerating phase of the flow's velocity profile are derived. The dependence of the solutions on the variable Eddington factor prescription chosen to close the radiation moment equations is also examined.
An algorithm searches a star catalog to identify guide stars within the field of view of a telescope or camera. The algorithm is fast: the number of computations needed to perform the search is approximately proportional to the logarithm of the number of stars in the catalog. The algorithm requires the prior organization of the star catalog into a hierarchy utilizing independent spherical coverings (see figure), such that each successively higher level contains fewer elements. In the lowest and most numerous level of the hierarchy, the elements are individual stars in the star catalog. The next higher level contains a spherical covering (a constellation of n points on a sphere that minimizes the maximum distance of any point on the sphere from the closest one of the n points), the next higher level contains a smaller spherical covering, and so forth, ending at the highest level, which contains one element representing the point of entry into the search structure. With necessary exceptions at the lowest and highest levels, each element at each level is labeled in terms of the element to which it is linked in the next higher level and the first element to which it is linked in the next lower level. Each element is also labeled in terms of (1) its coordinates on the celestial sphere and (2) the largest angular distance to any element in any lower level in the hierarchy. The elements at all levels of the hierarchy are numbered on a single list, such that the elements of each constellation at each level are numbered consecutively. The algorithm is recursive. The input required to start the algorithm comprises the coordinates of a point on the celestial sphere. Attention is then focused on individual elements of the hierarchy, starting from the topmost one, as follows: The angle between the input point and the element under consideration is calculated. If the calculated angle is larger than the sum of (1) the predetermined angle to the most distant element plus (2) the half field of view of the telescope, then no stars will be within the field of view and this recursive part of the algorithm is terminated.
A mesoscale numerical model of the Florida peninsula was formulated and applied to a dry, neutral atmosphere. The prospective use of the STAR-100 computer for the submesoscale model is discussed. The numerical model presented is tested under synoptically undisturbed conditions. Two cases, differing only in the direction of the prevailing geostrophic wind, are examined: a prevailing southwest wind and a prevailing southeast wind, both 6 m/sec at all levels initially.
The nature of star formation and Type Ia supernovae (SNIa) in galaxies in the field and in rich galaxy clusters are contrasted by juxtaposing the buildup of heavy metals in the universe inferred from observed star formation and supernovae rate histories with data on the evolution of Fe abundances in the intracluster medium (ICM). Models for the chemical evolution of Fe in these environments are constructed, subject to observational constraints, for this purpose. While models with a mean delay for SNIa of 3 Gyr and standard initial mass function (IMF) are fully consistent with observations in the field, cluster Fe enrichment immediately tracked a rapid, top-heavy phase of star formation - although transport of Fe into the ICM may have been more prolonged and star formation likely continued beyond redshift 1. The means of this prompt enrichment consisted of SNII yielding greater than or equal to 0.1 solar mass per explosion (if the SNIa rate normalization is scaled down from its value in the field according to the relative number of candidate progenitor stars in the 3 - 8 solar mass range) and/or SNIa with short delay times originating during the rapid star formation epoch. Star formation is greater than 3 times more efficient in rich clusters than in the field, mitigating the overcooling problem in numerical cluster simulations. Both the fraction of baryons cycled through stars, and the fraction of the total present-day stellar mass in the form of stellar remnants, are substantially greater in clusters than in the field.
Engineering programming information is presented for the STARS-2P (shell theory automated for rotational structures-2P (plasticity)) digital computer program, and FORTRAN 4 was used in writing the various subroutines. The execution of this program requires the use of thirteen temporary storage units. The program was initially written and debugged on the IBM 370-165 computer and converted to the UNIVAC 1108 computer, where it utilizes approximately 60,000 words of core. Only basic FORTRAN library routines are required by the program: sine, cosine, absolute value, and square root.
A special data debugging package called SAT-1P created for the STARS-2P computer program is described. The program was written exclusively in FORTRAN 4 for the IBM 370-165 computer, and then converted to the UNIVAC 1108.
The evolution of a neutron star undergoing a series of thermonuclear flashes in its accreted hydrogen-rich layer has been numerically followed to determine the effects of the history of the neutron star's thermal and compositional structure on the properties of the emitted X-ray bursts. Burst characteristics are studied for a range of mass accretion rates, CNO abundances in the accreted matter, and initial thermal states of the underlying neutron star core; the bursts exhibit erratic behavior for low CNO metal abundances and cool neutron star cores, with the burst recurrence time scales varying by 1-2 orders of magnitude. There is typically a continued presence of a substantial amount of unburnt hydrogen in the accreted layer throughout the series of the X-ray burst events. Convective mixing during the quiescent phase leads to the inward transport of helium to high densities and eventually to the initiation of the next outburst. The resulting bursts can be weak and, in such cases, are characterized by short recurrence time scales (1-2 hr), low peak luminosities (0.1-0.2 times the Eddington value), and low alpha-values (about 20).
Star formation and the creation of protostellar disks generally occur in a crowded environment. Nearby young stars and protostars can influence the disks of their closets neighbors by a combination of outflows and hard radiation. The central stars themselves can have a stellar wind and may produce sufficient UV and X-ray to ultimately destroy their surrounding disks. Here we describe the results of numerical simulations of the influence that an external UV source and a central star's wind can have on its circumstellar disk. The numerical method (axial symmetry assumed) is described elsewhere. We find that protostellar disks will be destroyed on a relatively short time scale (~ 10(sup 5)yr) unless they are well shielded from O-stars. Initially isotropic T-Tauri winds do not significantly influence their disks, but instead are focused toward the rotation axis by the disk wind from photoevaporation.
The physical processes determining the rate of star formation in protogalactic clouds and the resulting effects on the structure of the galaxies are examined with reference to analytical results and numerical simulations. For spheroidal systems, it is found that rapid star formation cannot be prevented once the density exceeds a typical value. The corresponding length scale is comparable to the length scales of spheroidal systems. The gas is converted into stars at a time scale of 2 Gyr or less, during which the total luminosity of the protogalaxy in ionizing photons is 10 exp 43-44 ergs/s, and the gas is enriched to solar abundance with a steep gradient. A simplistic analysis of the resulting stellar distribution finds that de Vaucouleurs profiles may naturally result. The Faber-Jackson relation is also found to follow naturally. The Fischer-Tully relation follows by an analysis similar to that used to derive the Faber-Jackson relation, provided that the peak rotational velocity is determined by the bulge mass.
Recent observations of X-ray pulsars at low luminosities allow, for the first time, to compare theoretical models for the emission from highly magnetized neutron star atmospheres at low mass accretion rates ( M ≲ 10(exp 15) g s(exp -1)) with the broadband X-ray data. The purpose of this paper is to investigate the spectral formation in the neutron star atmosphere at low M and to conduct a parameter study of physical properties of the emitting region. We obtain the structure of the static atmosphere, assuming that Coulomb collisions are the dominant deceleration process. The upper part of the atmosphere is strongly heated by the braking plasma, reaching temperatures of 30–40 keV, while its denser isothermal interior is much cooler (~2 keV). We numerically solve the polarized radiative transfer in the atmosphere with magnetic Compton scattering, free-free processes, and non-thermal cyclotron emission due to possible collisional excitations of electrons. The strongly polarized emitted spectrum has a double-hump shape that is observed in low-luminosity X-ray pulsars. A low-energy “thermal” component is dominated by extraordinary photons that can leave the atmosphere from deeper layers due to their long mean free path at soft energies. We find that a high-energy component is formed due to resonant Comptonization in the heated non-isothermal part of the atmosphere even in the absence of collisional excitations. The latter, however, affect the ratio of the two components. A strong cyclotron line originates from the optically thin, uppermost zone. A fit of the model to NuSTAR and Swift/XRT observations of GX3041 provides an accurate description of the data with reasonable parameters. The model can thus reproduce the characteristic double-hump spectrum observed in low-luminosity X-ray pulsars and provides insights into spectral formation.
Approximation techniques for use in numerical schemes for estimating spatially varying coefficients in continuum models such as those for Euler-Bernoulli beams are discussed. The techniques are based on quintic spline state approximations and cubic spline parameter approximations. Both theoretical and numerical results are presented. Previously announced in STAR as N83-28934
A homologous transformation is derived which permits the application of the numerical coronal models of Hammer from a star with solar mass and radius to other stars. This scaling requires a few approximations concerning the lower boundary conditions and the temperature dependence of the conductivity and emissivity. These approximations are discussed and found to be surprisingly mild. Therefore, the scaling of the coronal models to other stars is rather accurate; it is found to be particularly accurate for main-sequence stars. The transformation is used to derive an equation that gives the maximum temperature of open coronal regions as a function of stellar mass and radius, the coronal heating flux, and the characteristic damping length over which the corona is heated.
Numerical models of the development of star-formation bursts in collisional ring galaxies are presented. To extend the work of Appleton and Struck-Marcell (1987) target disks which have relatively high mean cloud mass and gas density are emphasized. In such cases, even relatively low mass intruder galaxies are capable of triggering intense star-formation bursts in the density waves. Although the bursts are very short-lived in any individual gas element, pressure effects stimulate neighboring gas elements to burst, which can result in a sustained enhancement in the net star-formation rate. The results are capable of explaining the high far-infrared fluxes observed in righ galaxies and provide clues to the development of starburst activity in other colliding galaxies.
A method is developed for determining the nuclear equation of state directly from a knowledge of the masses and radii of neutron stars. This analysis assumes only that equilibrium neutron-star matter has the stress-energy tensor of an isotropic fluid with a barotropic equation of state, and that general relativity describes a neutron star's internal gravitational field. We present numerical examples which illustrate how well this method will determine the equation of state when the appropriate observational data become available.
We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low-mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each simulation, adding realistic noise and observational limits. All three simulations show a decrease in the number of quiescent, isolated galaxies in the mass range M * = 10 9−10 M ꙩ , in broad agreement with observations. However, even after accounting for observational and selection biases, none of the simulations reproduce the observed absence of quiescent field galaxies below M * = 10 9 M ꙩ . We find that the low-mass quiescent populations selected via synthetic observations have consistent quenching timescales, despite an apparent variation in the late-time star formation histories. The effect of increased numerical resolution is not uniform across simulations and cannot fully mitigate the differences between the simulations and the observations. The framework presented here demonstrates a path toward more robust and accurate comparisons between theoretical simulations and galaxy survey observations, while the quenching threshold serves as a sensitive probe of feedback implementations.