Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical stars”

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Hydrodynamic effects in the atmosphere of variable stars

Numerical models of variable stars are established, using a nonlinear radiative transfer coupled hydrodynamics code. The variable Eddington method of radiative transfer is used. Comparisons are for models of W Virginis, beta Doradus, and eta Aquilae. From these models it appears that shocks are formed in the atmospheres of classical Cepheids as well as W Virginis stars. In classical Cepheids, with periods from 7 to 10 days, the bumps occurring in the light and velocity curves appear as the result of a compression wave that reflects from the star's center. At the head of the outward going compression wave, shocks form in the atmosphere. Comparisons between the hydrodynamic motions in W Virginis and classical Cepheids are made. The strong shocks in W Virginis do not penetrate into the interior as do the compression waves formed in classical Cepheids. The shocks formed in W Virginis stars cause emission lines, while in classical Cepheids the shocks are weaker.

Davis, C. G., Jr.↗
DOE OSTI2023

The Evolution of Population III and Extremely Metal-poor Binary Stars

Numerical simulations have now shown that Population III (Pop III) stars can form in binaries and small clusters and that these stars can be in close proximity to each other. If so, they could be subject to binary interactions such as mass exchange that could profoundly alter their evolution, ionizing UV and Lyman–Werner photon emission and explosion yields, with important consequences for early cosmological reionization and chemical enrichment. Here we investigate the evolution of Pop III and extremely metal-poor binary stars with the MESA code. We find that interactions ranging from stable mass transfer to common envelope evolution can occur in these binaries for a wide range of mass ratios and initial separations. Mass transfer can nearly double UV photon yields in some of these binaries with respect to their individual stars by extending the life of the companion star, which in turn can enhance early cosmological reionization but also suppress the formation of later generations of primordial stars. Binary interactions can also have large effects on the nucleosynthetic yields of the stars by promoting or removing them into or out of mass ranges for specific SN types. We provide fits to total photon yields for the binaries in our study for use in cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

A search for spectroscopic binaries among the runaway O type stars

Numerous radial velocity measurements of medium dispersion were made for the 10 brighter stars given in Stone's list of very probable O type runaways. All plates were measured with the KPNO PDS microdensitometer, and a new iterative reductional analysis was used to derive plate velocities, which are estimated to be 1.6 times more accurate internally than those found by using the traditional method. Of thse stars, psi Per, alpha Cam, HD 188209, and 26 Cep are identified as probable velocity variables, while 9 Sge, lambda Cep, and HD 218915 are classed as possible variables. If the source of this variability is Keplerian rather than atmospheric, which cannot be established unequivocally from the observations of this paper, psi Per could be a spectroscopic binary with a black hole companion, and at least 1.2 solar mass. The detection of runaway binary systems from radial velocity measurements is discussed.

Stone, R. C.↗
DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Spherically symmetric near-critical accretion onto neutron stars

Numerical and approximate analytic solutions for time-independent, spherically symmetric, radiation pressure-dominated accretion flows are presented. For flows with luminosities at infinity, L-infinity, sufficiently close to the Eddington limit L-crit, the flow velocity profile is qualitatively different from the modified free-fall profile v(r) = (1 - L-infinity/L-crit)exp 1/2 (2GM/r)exp 1/2. Advective contributions to the comoving radiation flux decelerate the flow within a criical radius, and, in this settling region, the velocity of the flow decreases linearly with decreasing radius.

Miller, Guy S.↗
NASA NTRS2003

Terrestrial Planet Formation in Binary Star Systems

Most stars reside in multiple star systems; however, virtually all models of planetary growth have assumed an isolated single star. Numerical simulations of the collapse of molecular cloud cores to form binary stars suggest that disks will form within such systems. Observations indirectly suggest disk material around one or both components within young binary star systems. If planets form at the right places within such circumstellar disks, they can remain in stable orbits within the binary star systems for eons. We are simulating the late stages of growth of terrestrial planets within binary star systems, using a new, ultrafast, symplectic integrator that we have developed for this purpose. We show that the late stages of terrestrial planet formation can indeed take place in a wide variety of binary systems and we have begun to delineate the range of parameter space for which this statement is true. Results of our initial simulations of planetary growth around each star in the alpha Centauri system and other 'wide' binary systems, as well as around both stars in very close binary systems, will be presented.

Lissauer, Jack J.↗
NASA NTRS2003

Terrestrial Planet Formation Around Close Binary Stars

Most stars reside in multiple star systems; however, virtually all models of planetary growth have assumed an isolated single star. Numerical simulations of the collapse of molecular cloud cores to form binary stars suggest that disks will form within such systems. Observations indirectly suggest disk material around one or both components within young binary star systems. If planets form at the right places within such circumstellar disks, they can remain in stable orbits within the binary star systems for eons. We are simulating the late stages of growth of terrestrial planets around close binary stars, using a new, ultrafast, symplectic integrator that we have developed for this purpose. The sum of the masses of the two stars is one solar mass, and the initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet growth within our Solar System and in the Alpha Centauri wide binary star system. Giant planets &are included in the simulations, as they are in most simulations of the late stages of terrestrial planet accumulation in our Solar System. When the stars travel on a circular orbit with semimajor axis of up to 0.1 AU about their mutual center of mass, the planetary embryos grow into a system of terrestrial planets that is statistically identical to those formed about single stars, but a larger semimajor axis and/or a significantly eccentric binary orbit can lead to significantly more dynamically hot terrestrial planet systems.

Lissauer, Jack J.↗
NASA NTRS1994

A high-resolution survey of interstellar Na I D1 lines

We present high-resolution (0.5 km/s) spectra, obtained with the McDonald Observatory 2.7 m coude echelle spectrograph, of interstellar Na I D1 absorption toward 38 bright stars. Numerous narrow, closely blended absorption components, showing resolved Na I hyperfine structure, are evident in these spectra; such narrow components appear in both low halo and quite local gas, as well as in gas toward more distant disk stars. We have used the method of profile fitting in an attempt to determine column densities, line widths, and velocities for the individual interstellar clouds contributing to the observed absorption lines. The resulting sample of 276 clouds is significantly larger, and likely more complete, than several previous samples of 'individual' interstellar clouds, and allows more precise investigation of various statistical properties. We find that the cloud column density (N) and line width parameter (b) are not correlated, for 0.3 km/s approximately less than b approximately less than 1.5 km/s and 10.0/sq cm approximately less than log (N(Na I)) approximately less than 11.6/sq cm. The median b is about 0.73 km/s, the median log N is about 11.09/sq cm, and the median separation between adjacent components is about 2.0 km/s. All these are overestimates of the true median values, however, due to our inability to completely resolve all the component structure present in some cases; even at a resolution of 0.5 km/s, we may have discerned only 60% of the full number of individual components actually present. The one-dimensional dispersion of component velocities, in the local standard of rest, is approximately 8.6 km/s; the distribution of velocities is broader and displaced to more negative velocities for the weaker components. If 80 K is a representative temperature for the interstellar clouds seen in Na I absorption, then at least 38% (and probably the majority) of the clouds have subsonic internal turbulent motions. The range in N(H I) observed at a given N(Na I) increases as N(Na I) decreases below about 10(exp 11)/sq cm, so that N(Na I) becomes a less reliable predictor of N(H I) at low column densities. These spectra will be quite useful in future detailed studies of these lines of sight with the Hubble Space Telescope (HST) Goddard High Resolution Spectrograph (GHRS) echelle -- to determine accurate velocities and to reveal the detailed interstellar component structures that cannot be discerned at the 3.5 km/s resolution available with the GHRS.

Welty, Daniel E.↗
DOE OSTI2021

Wide Dust Gaps in Protoplanetary Disks Induced by Eccentric Planets: A Mass-eccentricity Degeneracy

The tidal perturbation of embedded protoplanets on their natal disks has been widely attributed to be the cause of gap-ring structures in submillimeter images of protoplanetary disks around T Tauri stars. Numerical simulations of this process have been used to propose scaling of characteristic dust-gap width/gap-ring distance with respect to planet mass. Applying such scaling to analyze observed gap samples yields a continuous mass distribution for a rich population of hypothetical planets in the range of several Earth to Jupiter masses. In contrast, the conventional core-accretion scenario of planet formation predicts a bimodal mass function due to (1) the onset of runaway gas accretion above ~20 Earth masses and (2) suppression of accretion induced by gap opening. Here, we examine the dust disk response to the tidal perturbation of eccentric planets as a possible resolution of this paradox. Based on simulated gas and dust distributions, we show the gap-ring separation of Neptune-mass planets with small eccentricities might become comparable to that induced by Saturn-mass planets on circular orbits. Furthermore, this degeneracy may obliterate the discrepancy between the theoretical bimodal mass distribution and the observed continuous gap width distribution. Despite damping due to planet–disk interaction, modest eccentricity may be sustained either in the outer regions of relatively thick disks or through resonant excitation among multiple super Earths. Moreover, the ring-like dust distribution induced by planets with small eccentricities is axisymmetric even in low viscosity environments, consistent with the paucity of vortices in Atacama Large Millimeter/submillimeter Array images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Evolution of Inclined Binary Black Holes in the Disks of Active Galactic Nuclei

The accretion disks that fuel active galactic nuclei (AGNs) may house numerous stars and compact objects, formed in situ or captured from nearby star clusters. Embedded neutron stars and black holes may form binaries and eventually merge, emitting gravitational waves detectable by LIGO/VIRGO. AGN disks are a particularly promising environment for the production of high-mass gravitational-wave events involving black holes in the pair-instability mass gap, and may facilitate electromagnetic counterparts to black hole binary mergers. However, many orders of magnitude separate the typical length scales of binary formation and those on which gravitational waves can drive binary inspirals, making binary mergers inside the disk uncertain. Previous hydrodynamical simulations of binaries have either been restricted to two dimensions entirely, or focused on binaries aligned with the midplane of the disk. Herein we present the first three-dimensional, high-resolution, local-shearing-box, inviscid hydrodynamical simulations of disk-embedded binaries over a range of orbital inclinations. We find that retrograde binaries can shrink up to 4 times as quickly as prograde binaries, and that all binaries not perfectly aligned (or anti-aligned) with the AGN disk are driven into alignment. An important consequence of this is that initially retrograde binaries will traverse the inclinations where von Zeipel–Lidov–Kozai oscillations can drive binary eccentricities to large values, potentially facilitating mergers. We also find that interactions with the AGN disk may excite eccentricities in retrograde binaries and cause the orbits of embedded binaries to precess.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Mapping the Outcomes of Stellar Evolution in the Disks of Active Galactic Nuclei

The disks of active galactic nuclei (AGNs) are expected to be populated by numerous stars, either formed in the outer regions of the disk via gravitational instability or captured from the nearby nuclear star cluster. Regardless of their formation mechanism, these stars experience altered evolutionary paths, mostly shaped by the accretion of dense disk material. In this study, through the comparison of different timescales, we chart the evolutionary outcomes of these AGN stars as a function of disk radius and across a range of supermassive black hole masses, spanning from 10 6 to 10 9 M ⊙ , for two popular AGN disk models. We find that in the outer regions of the disk, stars evolve similarly to those in the interstellar medium, but in the inner and denser regions, accretion quickly turns low-mass stars into massive stars, and their fate depends on just how quickly they accrete. If accretion occurs at a faster rate than nuclear burning, they can reach a quasi-steady "immortal" state. If stars accrete faster than they can thermally adjust, runaway accretion occurs, potentially preventing a quasi-steady state and altering the disk structure. During the AGN lifetime, in the regions of the disk that produce massive stars, supernovae (SNe) and gamma-ray bursts (GRBs) may occur within the disk over a wide range of optical depths and ambient densities. Subsequently, in the final phase of the AGN, as the disk becomes depleted, formerly immortal stars will be unable to replenish their fuel, leading to additional SNe and GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

Velocity dispersions in galaxies: 1: The SO galaxy NGC 7332

A Coude spectrum of the SO galaxy NGC 7332 with 0.9 A resolution from 4186 to 4364 A was obtained with the SEC vidicon television camera and the Hale telescope. Comparisons with spectra of G and K giant stars, numerically broadened for various Maxwellian velocity distributions, give a dispersion velocity in the line of sight of 160 + or - 20 km/sec with the best fit at G8III. The dispersion appears to be constant within + or - 35 km/sec out to 1.4 kpc (H = 100 km/sec/mpc). After correction for projection, the rotation curve has a slope of 0.16 km/sec/pc at the center and a velocity of 130 km/sec at 1.4 kpc where it is still increasing. For an estimated effective radius of 3.5 kpc enclosing half the light, the virial theorem gives a mass of 1.4 x 10 to the 11th power solar masses if the mass-to-light ratio is constant throughout the galaxy. The photographic luminosity is 8.3 x 10 to the 9th power solar luminosities so that the M/L ratio is 17.

Morton, D. C.↗
NASA NTRS1972

Velocity dispersions in galaxies. I - The E7 galaxy NGC 7332.

A coude spectrum of the E7 galaxy NGC 7332 with 0.9 A-resolution from 4186 to 4364 A was obtained with the Princeton SEC vidicon television camera and the Hale telescope. Comparisons with spectra of G and K giant stars, numerically broadened for various Maxwellian velocity distributions, give a dispersion velocity in the line of sight of 160 (plus or minus 20) km/sec with the best fit at G8 III. The dispersion appears to be constant within plus or minus 35 km/sec out to 1.4 kpc. After correction for projection, the rotation curve has a slope of 0.18 km/sec per pc at the center and a velocity of 130 km/sec at 1.4 kpc where it is still increasing. For an estimated effective radius of 3.5 kpc enclosing half the light, the virial theorem gives a mass of 140 billion solar masses if the mass-to-light ratio is constant throughout the galaxy.

Morton, D. C.↗
NASA NTRS2006

The HST/STIS Next Generation Spectral Library

During Cycles 10, 12, and 13, we obtained STIS G230LB, G430L, and G750L spectra of 378 bright stars covering a wide range in abundance, effective temperature, and luminosity. This HST/STIS Next Generation Spectral Library was scheduled to reach its goal of 600 targets by the end of Cycle 13 when STIS came to an untimely end. Even at 2/3 complete, the library significantly improves the sampling of stellar atmosphere parameter space compared to most other spectral libraries by including the near-UV and significant numbers of metal poor and super-solar abundance stars. Numerous calibration challenges have been encountered, some expected, some not; these arise from the use of the E1 aperture location, non-standard wavelength calibration, and, most significantly, the serious contamination of the near-UV spectra by red light. Maximizing the utility of the library depends directly on overcoming or at least minimizing these problems, especially correcting the UV spectra.

Gregg, M. D.↗
NASA NTRS1987

Numerical models of star formation in X-ray cluster cooling flows

A grid of numerical models is used here to study how accretion flows are affected by star formation. The observables associated with these models are investigated, presenting their X-ray surface brightness profiles. An X-ray diagnostic is searched for that, in the absence of spatially resolved X-ray spectral data, can be an unambiguous indicator for ongoing star formation in cooling flows. The spatial distributions of newly formed stars are presented and compared to profiles of the accreting galaxy in each model. Whether the accretion population can form the luminous parts of the accreting galaxy or whether it forms a dark component is considered.

White, Raymond E., III↗
DOE OSTI2023

Study of the agreement between binary neutron star ejecta models derived from numerical relativity simulations

Neutron star mergers have recently become a tool to study extreme gravity, nucleosynthesis, and the chemical composition of the Universe. To date, there has been one joint gravitational and electromagnetic observation of a binary neutron star merger, GW170817, as well as a solely gravitational observation, GW190425. In order to accurately identify and interpret electromagnetic signals of neutron star mergers, better models of the matter outflows generated by these mergers are required. We compare a series of ejecta models to see where they provide strong constraints on the amount of ejected mass expected from a system, and where systematic uncertainties in current models prevent us from reliably extracting information from observed events. We also examine 2396 neutron star equations of state compatible with GW170817 to see whether a given ejecta mass could be reasonably produced with a neutron star of said equation of state, and whether different ejecta models provide consistent predictions. We find that the difference between models is often comparable to or larger than the error generally assumed for these models, implying better constraints on the models are needed. In conclusion, we also note that the extrapolation of outflow models outside of their calibration window, while commonly needed to analyze gravitational wave events, is extremely unreliable and occasionally leads to completely unphysical results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Error quantification and comparison of binary neutron star gravitational waveforms from numerical relativity codes

Future gravitational wave detections of merging binary neutron star systems have the possibility to tightly constrain the equation of state of dense nuclear matter. In order to extract such constraints, gravitational waveform models need to be calibrated to accurate numerical relativity simulations of the late inspiral and merger. In this work, we take an essential step toward classifying the error and potential systematics in current generation numerical relativity simulations of merging binary neutron stars. To this end, we perform a direct comparison of two codes (fil, spec), which differ in many aspects, including the numerical methods and discretizations used and equations solved. We find that despite these different approaches, the codes are—within current numerical resolution bounds—fully consistent, and broadly comparable in cost for a given accuracy level. Our results indicate that the error in the waveforms is primarily dominated by the hydrodynamic evolution, consistent with earlier findings in the literature. Finally, we also discuss current limitations and cost estimates for numerical relativity simulations to reach the accuracies required in the era of next-generation gravitational detectors.

gravitational waves↗