Numerical experiments with a disk of stars
Computer model for dynamical evolution of isolated disks of stars with initial velocity dispersions corresponding to Toomre criterion, noting instability relative to large scale modes
SEARCH · Engineering Papers
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.
Computer model for dynamical evolution of isolated disks of stars with initial velocity dispersions corresponding to Toomre criterion, noting instability relative to large scale modes
Homogeneous models of luminous hydrogen-and-helium burning stars were constructed on the basis of Carson's new radiative opacities. These opacities exhibit a large 'bump' at moderate temperatures and low densities as a result of the ultimate ionization of the CNO group of elements and they induce in the envelopes of the more massive stars a strong local convection zone, a high central condensation and pulsational instability by means of the kappa-mechanism above a certain stellar mass. This critical mass for Population I hydrogen-burning stars is probably greater than 45 solar masses for the fundamental radial mode, with the overtones being more stable, while for the homogeneous helium-burning stars, the critical mass lies in the range 2-6 solar masses for all modes up to at least the third overtone. Convection is alternately treated by means of the mixing length theory and the assumption of strict adiabaticity.
Theoretical models developed to explain the emission of hot stars are examined, summarizing and analyzing the results of ground-based and space observations obtained at different wavelengths. The physical processes responsible for emission in the radio, IR, visible, UV, and X-ray bands are reviewed, and particular attention is given to observations of OB-star winds, B-emission stars, compact stellar cores, and hot stars in binary systems. The importance of the UV and EUV spectral regions and the need for a space observatory to replace the IUE when it becomes inoperative are stressed.
A model is developed to investigate the effect of wavelength dependence of magnetically generated broadband linear polarization (BLP) in late-type stars as a function of line blanketing, observation filters, and other variables. Linear polarization is assumed to be a single 'average' line and the line opacity, magnetic area filling factor, and continuum intensity are incorporated to establish an equation for disk-generated BLP. The method is investigated by computing special models for the flare star BY Dra and the sun which are compared to previous calculations. The model is shown to agree with the calculations of BLP, and the magnetic BLP dependent on wavelengths varies significantly from the simple power law dependences on wavelength. The model sets lower limits to the area occupied by magnetic regions, and an upper limit to the expected BLP also results from the formulation.
The high-density interior of a neutron star is expected to contain superconducting protons and superfluid neutrons. Theoretical estimates suggest that the protons will form a type II superconductor in which the stellar magnetic field is carried by flux tubes. The strong interaction between the flux tubes and the neutron rotational vortices could lead to strong ''pinning'', i.e., vortex motion could be impeded. This has important implications especially for pulsar glitch models as it would lead to a large part of the vorticity of the star being decoupled from the ''normal'' component to which the electromagnetic emission is locked. In this Letter, we explore the consequences of strong pinning in the core on the ''snowplow'' model for pulsar glitches, making use of realistic equations of state and relativistic background models for the neutron star. We find that, in general, a large fraction of the pinned vorticity in the core is not compatible with observations of giant glitches in the Vela pulsar. Thus, the conclusion is that either most of the core is in a type I superconducting state or the interaction between vortices and flux tubes is weaker than previously assumed.
The purpose of this study is to create a STAR-CCM+ model of a Belowground Vertical Dry Cask Simulator (BVDCS) at Sandia National Laboratories (SNL) and validate the model with SNL’s experimental results. The BVDCS consists of a single boiling water reactor assembly fitted with electric heaters encompassed by a containment vessel and shell to represent a belowground spent nuclear fuel (SNF) dry storage system. Blowers are located near the inlet and outlet of the BVDCS to simulate crosswind conditions. In addition to the experimental results, the STAR-CCM+ model developed for this study is compared with a previous computational fluid dynamics (CFD) model in a different software program, which is used as a software-to-software benchmark. The experimental results provide a dataset to compare the STAR-CCM+ model results for a variety of different conditions. The main objective is to validate and improve STAR-CCM+ CFD models for spent nuclear fuel storage systems with explicitly modeled external environments and “wind driven” crossflows. These CFD models aide in the study of external particle deposition in spent nuclear fuel storage systems, which is important to predicting the significance of chloride induced stress corrosion cracking (CISCC). In addition to experimental comparison, a sensitivity analysis study is performed using the STAR-CCM+ model. The sensitivity analysis provides a quantitative assessment of the sensitivity of various parameters. This helps provide information on various parameters that are of particular importance to constructing a model representative of real life systems. The STAR-CCM+ model compared well to the experimental results showing similar responses to changes in cross wind flow, and a number of parameters are identified for model improvement.
We present a grid of atmospheric models and synthetic spectral energy distributions (SEDs) for late-type dwarfs and giants of solar and 1/3 solar metallicity with many opacity sources computed in self-consistent non-local thermodynamic equilibrium (NLTE), and compare them to the LTE grid of Short and Hauschildt (Paper I). We describe, for the first time, how the NLTE treatment affects the thermal equilibrium of the atmospheric structure (T({tau}) relation) and the SED as a finely sampled function of T{sub eff}, log g, and [A/H] among solar metallicity and mildly metal-poor red giants. We compare the computed SEDs to the library of observed spectrophotometry described in Paper I across the entire visible band, and in the blue and red regions of the spectrum separately. We find that for the giants of both metallicities, the NLTE models yield best-fit T{sub eff} values that are 30-90 K lower than those provided by LTE models, while providing greater consistency between log g values, and, for Arcturus, T{sub eff} values, fitted separately to the blue and red spectral regions. There is marginal evidence that NLTE models give more consistent best-fit T{sub eff} values between the red and blue bands for earlier spectral classes among the solar metallicity GK giants than they do for the later classes, but no model fits the blue-band spectrum well for any class. For the two dwarf spectral classes that we are able to study, the effect of NLTE on derived parameters is less significant. We compare our derived T{sub eff} values to several other spectroscopic and photometric T{sub eff} calibrations for red giants, including one that is less model dependent based on the infrared flux method (IRFM). We find that the NLTE models provide slightly better agreement to the IRFM calibration among the warmer stars in our sample, while giving approximately the same level of agreement for the cooler stars.
In the present paper we investigate non-perturbatively and self-consistently the structure of neutron stars in R-squared gravity by simultaneously solving the interior and exterior problem. The mass-radius relations are obtained for several equations of state and for wide range of the R-squared gravity parameter a. Even though the deviation from general relativity for nonzero values of a can be large, they are still comparable with the variations due to different modern realistic equations of state. That is why the current observations of the neutron star masses and radii alone can not put constraints on the value of the parameter a. We also compare our results with those obtained within the perturbative method and we discuss the differences between them.
Seismic observations have led to doubts or ambiguities concerning the opacity calculations used in stellar physics. Here, we concentrate on the iron-group opacity peak, due to iron, nickel, and chromium, located around T = 200,000 K for densities from 10{sup −8} to 10{sup −4} g cm{sup −3}, which creates some convective layers in stellar radiative envelopes for masses between 3 and 18 M{sub ⊙}. These conditions were extensively studied in the 1980s. More recently, inconsistencies between OP and OPAL opacity calculations have complicated the interpretation of seismic observations as the iron-group opacity peak excites acoustic and gravity modes in SPB, β Cephei, and sdB stars. We investigate the reliability of the theoretical opacity calculations using the modern opacity codes ATOMIC and SCO-RCG. We show their temperature and density dependence for conditions that are achievable in the laboratory and equivalent to astrophysical conditions. We also compare new theoretical opacity spectra with OP spectra and quantify how different approximations impact the Rosseland mean calculations.This detailed study estimates new ATOMIC and SCO-RCG Rosseland mean values for astrophysical conditions which we compare to OP values. Some puzzling questions are still under investigation for iron, but we find a strong increase in the Rosseland mean nickel opacity of a factor between 2 and 6 compared to OP. This appears to be due to the use of extrapolated atomic data for the Ni opacity within the OP calculations. A study on chromium is also shown.
Pulsational stability of very massive stars initially composed of pure H, discussing effects of nonconstant scattering opacity
We consider the structure and emission spectra of nonthermal accretion disks around both strongly and weakly magnetized neutron stars. Such disks may be dissipating their gravitational binding energy and transferring their angular momentum via semicontinuous magnetic reconnections. We consider specifically the structure of the disk-stellar magnetospheric boundary where magnetic pressure balances the disk pressure. We consider energy dissipation via reconnection of the stellar field and small-scale disk turbulent fields of opposite polarity. Constraints on the disk emission spectrum are discussed.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Astrophysical models of primordial star formation require rate constants for three-body recombination as input. The current status of these rates for H{sub 2} due to collisions with H is far from satisfactory, with published rate constants showing orders of magnitude disagreement at the temperatures relevant for H{sub 2} formation in primordial gas. This Letter presents an independent calculation of this recombination rate constant as a function of temperature. An analytic expression is provided for the rate constant which should be more reliable than ones currently being used in astrophysical models.
Abstract In this paper, we derive correlations between core-collapse supernova observables and progenitor core structures that emerge from our suite of 20 state-of-the-art 3D core-collapse supernova simulations carried to late times. This is the largest such collection of 3D supernova models ever generated and allows one to witness and derive testable patterns that might otherwise be obscured when studying one or a few models in isolation. From this panoramic perspective, we have discovered correlations between explosion energy, neutron star gravitational birth masses, 56 Ni and α -rich freezeout yields, and pulsar kicks and theoretically important correlations with the compactness parameter of progenitor structure. We find a correlation between explosion energy and progenitor mantle binding energy, suggesting that such explosions are self-regulating. We also find a testable correlation between explosion energy and measures of explosion asymmetry, such as the ejecta energy and mass dipoles. While the correlations between two observables are roughly independent of the progenitor zero-age main-sequence (ZAMS) mass, the many correlations we derive with compactness cannot unambiguously be tied to a particular progenitor ZAMS mass. This relationship depends on the compactness/ZAMS mass mapping associated with the massive star progenitor models employed. Therefore, our derived correlations between compactness and observables may be more robust than with ZAMS mass but can nevertheless be used in the future once massive star modeling has converged.
A model for disk accretion by a rotating magnetic neutron star is proposed which includes a detailed description of matter flow in the transition region between the disk and the magnetosphere. It is shown that the disk plasma cannot be completely screened from the stellar magnetic field and that the resulting magnetic coupling between the star and the disk exerts a significant torque on the star. On the assumption that the distortion of the residual stellar field lines threading the disk is limited by reconnection, the total accretion torque on the star is calculated. The calculated torque gives period changes in agreement with those observed in the pulsating X-ray sources and provides a natural explanation of why a fast rotator like Her X-1 has a spin-up rate much below the conventional estimate for slow rotators. It is shown that for such fast rotators, fluctuations in the mass-accretion rate can produce fluctuations in the accretion torque about 100 times larger. For sufficiently fast rotators or, equivalently, for sufficiently low accretion rates, the star experiences a braking torque even while accretion continues and without any mass ejection from its vicinity.
A model for X-ray bursts from accreting neutron stars is developed by analogy with geomagnetic substorms. The essential steps in the substorm process are the nearly steady merging or reconnection of the magnetic field in the magnetosphere with the field in the stellar wind, the transport of some of the merged plasma into a magnetically controlled tail, and the explosive release of plasma from the tail into the magnetosphere. The strength of the magnetic field in the stellar wind required to drive a substorm is approximately 0.1 gauss. If the stellar wind is organized into large-scale magnetic sectors, as is the solar wind, topological dissipation will not occur, and the large-scale field will be available for merging at the magnetopause. Once the material is in the tail, the time scales for the Kruskal-Schwarzschild instability and the unidentified instability which drives terrestrial substorms may be comparable. Alternating periods of burst activity and quiescence could be caused by passage from one sector to another with opposite polarity, or be seasonal variations.
A neutron star in a close binary orbit around another neutron star (or stellar-mass black hole) spirals inward owing to gravitational radiation. We discuss the effects of tidal dissipation during this process. Tidal energy dissipated in the neutron star's core escapes mainly as neutrinos, but heating of the crust, and outward diffusion of photons, blows off the outer layers of the star. This photon-driven mass loss precedes the final coalescence. The presence of this eject material impedes the escape of gamma-rays created via neutrino interactions. If an e(+) - e(-) fireball, created in the late stages of coalescence, were loaded with (or surrounded by) material with the mean column density of the ejecta, it could not be an efficient source of gamma-rays. Models for cosmologically distant gamma-rays burst that involve neutron stars must therefore be anisotropic, so that the fireball expands preferentially in directions where the column density of previously blown-off material is far below the spherically averaged value which we have calculated. Some possible 'scenarios' along these lines are briefly discussed.