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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A new moment-based general-relativistic neutrino-radiation transport code: Methods and first applications to neutron star mergers

ABSTRACT We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Addition of tabulated equation of state and neutrino leakage support to illinoisgrmhd

Here we have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to illinoisgrmhd, an open-source GRMHD code for dynamical spacetimes in the einstein toolkit. These new features are provided by two new, nrpy+-based codes: nrpyeos, which performs highly efficient EOS table lookups and interpolations, and nrpyleakage, which implements a new, adaptive mesh refinement (AMR)-capable neutrino leakage scheme in the einstein toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of illinoisgrmhd, as well as nrpyeos and nrpyleakage, is pedagogically documented in jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the einstein toolkit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cocoon shock breakout emission from binary neutron star mergers

Shock breakout emission is among the first observable signals in a wide variety of astrophysical phenomena, including neutron star (NS) mergers, and it can be the dominant component in low-luminosity short γ-ray bursts, as exemplified by GRB 170817A. In this work, we investigate the cocoon shock breakout emission in NS mergers and how its signal depends on the outermost layers of the ejecta profile, which we derive from general relativistic radiation hydrodynamic simulations. We study the formation of the cocoon as a consequence of a relativistic jet propagating through the ejecta. To explore the influence of the outermost layers of the ejecta on the breakout emission, we explore cases where the ejecta has a sharp cutoff or an extended smooth tail. We find that the shock breakout emission is strongly influenced by the shape of the ejecta outer layers, with extended tails yielding results consistent with the observed properties of GRB 170817A, whereas sharp cutoffs overestimate the radiated energy. Using a Bayesian analysis, we estimate the best fit parameters of the central engine, considering both accreting black hole and magnetized neutron star scenarios. Our findings indicate a slight preference for the scenarios where the engine is a black hole. Furthermore, our work probes the nature of neutron star mergers and highlights the importance of the shape of the ejecta profile in modeling early electromagnetic counterparts to these mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants

We study out-of-thermodynamic-equilibrium effects in neutron-star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During mergers, the cores of the neutron stars remain cold (T ~ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ~2 to 3 ms matter and neutrinos reach equilibrium everywhere in the remnant massive neutron star. Furthermore, our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Neutrino emission from binary neutron star mergers: characterising light curves and mean energies

Abstract Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resulting r -process nucleosynthesis and kilonova emission. In this work, we analyse neutrino emission from a large sample of binary neutron star merger simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several $$10^{53}{\hbox {erg}~{\hbox {s}}^{-1}}$$ 10 53 erg s - 1 ) decreases as the reduced tidal deformability increases, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterisation of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Preliminary Thermal-Hydraulic Analysis of the UTA–2 Subcritical Linear Accelerator–Driven System

The National Nuclear Security Administration’s mission includes establishing a reliable supply of 99 Mo without highly enriched uranium. Oak Ridge National Laboratory (ORNL) supports this objective through collaborative research and development with industrial partners. Niowave Inc., a current partner, is currently designing a subcritical linear accelerator-driven system (ADS) and preparing for the US Nuclear Regulatory Commission’s licensing process. Niowave’s system has the potential to efficiently supply medical radioisotopes. The technology includes a superconducting electron accelerator and a pile of both natural and low-enriched uranium (LEU) targets. The process fissions uranium and many valuable isotopes can then be extracted from the targets. Niowave is currently iterating through conceptual and detailed design processes for several system sizes. This report discusses UTA–2, which is at the demonstration stage. UTA–2 will validate numerical modeling results with experimental measurements before progressing to the detailed design of UTA–3, the commercial-sized ADS. The thermal-hydraulic behavior of the UTA–2 core design was numerically investigated using STAR-CCM+, as described in this report. STAR-CCM+, a state-of-the-art computational fluid dynamics (CFD) software that was commercially developed by Siemens, has an extensive user base and a set of validation studies. It is also compliant with the American Society of Mechanical Engineers’ Nuclear Quality Assurance 1 standard. Three cases are investigated in this report. Case 1 quantified the temperature field within the UTA–2 assembly and water tank using only conduction as the method of thermal energy transport. This simplified approach was overly conservative and yielded wetted cladding temperatures above the coolant saturation temperature. In this case, the maximum temperature of the wetted cladding surface of the highest power rod exceeded the saturation temperature of water by 63.8°C. Because of the overly conservative approach taken in Case 1 and its negative subcooled margin, buoyancy-driven natural circulation flow physics were implemented in Case 2. Adding coolant motion significantly distributed the thermal energy of the system through convective heat transfer. This relatively small amount of convective heat transfer significantly reduced system temperatures and increased the subcooled margin from –63.8 to 61.3°C. This margin confirmed that no boiling was expected during normal operation of UTA–2 at 230 W. Case 3 had no additional physics models but considered an overpower event during which the power of each LEU and natural uranium rod was at its respective peak values. This resulted in a study with the total assembly power equal to 176% of the nominal power of 230 W considered in Cases 1 and 2. The resulting natural circulation flows were slightly enhanced. The subcooled margin decreased slightly to 46.3°C, which is still a significant margin to local boiling of the water within the tank. This margin confirmed that no boiling was expected during an abnormal operation of UTA–2 at 406 W.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Collapsar Black Holes Are Likely Born Slowly Spinning

Collapsing stars constitute the main black hole (BH) formation channel, and are occasionally associated with the launch of relativistic jets that power γ-ray bursts (GRBs). Thus, collapsars offer an opportunity to infer the natal (before spin-up/down by accretion) BH spin directly from observations. We show that once the BH saturates with a large-scale magnetic flux, the jet power is dictated by the BH spin and mass accretion rate. Core-collapse simulations by Halevi et al. and GRB observations favor stellar density profiles that yield an accretion rate of $\dot m\approx 10^{-2} M_{\odot}s^{-1}$, weakly dependent on time. This leaves the spin as the main factor that governs the jet power. By comparing the jet power to characteristic GRB luminosities, we find that the majority of BHs associated with jets are likely born slowly spinning with a dimensionless spin of a ≃ 0.2, or a ≃ 0.5 for wobbling jets, with the main uncertainty originating in the unknown γ-ray radiative efficiency. This result could be applied to the entire core-collapse BH population, unless an anticorrelation between the stellar magnetic field and angular momentum is present. In a companion paper, Jacquemin-Ide et al., we show that regardless of the natal spin, the extraction of BH rotational energy leads to spin-down to a ≲ 0.2, consistent with gravitational-wave observations. We verify our results by performing the first 3D general-relativistic magnetohydrodynamic simulations of collapsar jets with characteristic GRB energies, powered by slowly spinning BHs. We find that jets of typical GRB power struggle to escape from the star, providing the first numerical indication that many jets fail to generate a GRB.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Viscous effects in rapidly rotating stars with application to white-dwarf models. III - Further numerical results

Improved viscous evolutionary sequences of differentially rotating, axisymmetric, nonmagnetic, zero-temperature white-dwarf models are constructed using the relativistically corrected degenerate electron viscosity. The results support the earlier conclusion that angular momentum transport due to viscosity does not lead to overall uniform rotation in many interesting cases. Qualitatively different behaviors are obtained, depending on how the total mass M and angular momentum J compare with the M and J values for which uniformly rotating models exist. Evolutions roughly determine the region in M and J for which models with a particular initial angular momentum distribution can reach carbon-ignition densities in 10 b.y. Such models may represent Type I supernova precursors.

Durisen, R. H.↗
NASA NTRS1979

Occultation of Epsilon Geminorum by Mars. III - Temperature structure of the Martian upper atmosphere

High-quality photoelectric observations of the April 8, 1976, occultation of Epsilon Gem by Mars are analyzed which were obtained at five different observatories. Temperature profiles are derived by numerical inversion of the star's light curves and found to be in mutual agreement. This is taken as evidence for the presence of strong atmospheric waves in the Martial upper atmosphere between number densities of 10 to the 13th and 10 to the 15th per cu cm. It is estimated that the vertical wavelength is about 20 km, the horizontal wavelength is greater than 600 km, and the peak-to-peak wave amplitude is approximately 40 K. The results are shown to be qualitatively consistent with Zurek's (1976) model of thermally driven tidal modes. No compelling evidence for turbulence effects is found in the occultation data.

French, R. G.↗
NASA NTRS1982

IUE spectroscopic investigation of interacting binary systems

The IUE observations of binaries consisting of at least one early type star, to the Algol systems, and to W Serpentis systems are discussed. The OB close binaries almost invaribly show a stellar wind. The difference between an unevolved early type system and a slightly evolved system is difficult to discern observationally and probably not meaningful in terms of Roche lobes. the zeta Aur stars show complex gas motions and shocks. Many of the Algol systems show peculiarities in their ultraviolet continua and the resonance lines, particularly Si IV, are of abnormal and variable strength. Emission lines were detected in several systems during the total eclipse. The W Ser stars are characterized by numerous strong emission lines in their ultraviolet spectra. It appears that these systems may be related to beta Lyrae as they seem to be in a phase of rapid mass transfer and their spectra show similarities to that of beta Lyrae.

Mccluskey, G. E., Jr.↗
NASA NTRS1981

Thermonuclear processes on accreting neutron stars

Theoretical models for X-ray burst sources that invoke thermonuclear flashes on the surface layers of an accreting neutron star are discussed. The historical development of X-ray burst observation is summarized, and a physical picture of a neutron star undergoing accretion is drawn. Detailed numerical computations of the evolution of the surface layers of such a star are reviewed. The need for general relativistic corrections to the model is pointed out. Finally, comparisons are made with observations of X-ray bursts, the rapid burster, fast X-ray transients, X-ray pulsars, and gamma-ray burst sources.

Joss, P. C.↗
NASA NTRS1991

Cooling flows and the formation of massive halos in cD galaxies

X-ray brightness profiles of galaxy clusters frequently indicate the presence of a hot gaseous envelope, and theoretical models predict that numerous low-mass stars will be formed as these gaseous halos cool and contract toward the center of the cluster potential. In this Letter, the optical surface brightness profiles of the halos seen in six cooling flow galaxies are modeled, and it is concluded that the optical profiles and mass-to-light ratios are consistent with the cooling flow accretion hypothesis.

Prestwich, Andrea H.↗
NASA NTRS1994

Dense stellar cores in merger remnants

We use numerical models which include star formation to analyze the mass profiles of remnants formed by mergers of disk galaxies. During a merger, dissipation in gas and ensuing star formation leave behind a dense stellar core in the remnant. Rather than joining smoothly onto a de Vaucouleurs profile, the starburst population leads to a sharp break in the surface density profile at a few percent of the effective radius. While our results are preliminary, the lack of such signatures in most elliptical galaxies suggests that mergers of gas-rich disk galaxies may not have contributed greatly to the populations of present-day ellipticals.

Mihos, J. Christopher↗
NASA NTRS1998

Stellar Occultation Studies of the Solar System

Earth-based observations of stellar occultations provide extremely high spatial resolution for bodies in the outer solar system, about 10,000 times better than that of traditional imaging observations. Stellar occultation data can be used to establish the structure of atmospheres and rings of solar system bodies at high spatial resolution. Airborne occultation observations are particularly effective, since the controlled mobility of the observing platform allows the observer to fly within the optimum part of the occultation shadow for most events that are visible from Earth. Airborne observations are carried out above any clouds and are nearly free of scintillation noise from the Earth's atmosphere. KAO occultation observations resulted in the first detection of gravity waves in the Martian atmosphere, discovery of the Uranian rings, the first detection of Pluto's atmosphere, the first Earth-based investigations of Triton's atmosphere, and the discovery of narrow jets from Chiron's nucleus. The first SOFIA occultation opportunity will be an investigation of Pluto's atmospheric structure in November, 2002, and will resolve a problem that has lingered since the KAO discovery observation fourteen years earlier. We plan to continue our successful airborne occultation program with the greatly enhanced capability provided by SOFIA. We propose here to replace our KAO occultation photometer with one having twice the throughput, half the noise, a somewhat wider wavelength range, four times the field of view, and ten times the frame rate to optimize its performance and to capitalize on the larger collecting area offered by SOFIA. It will also allow for simultaneous visible and IR occultation observations, greatly enriching the results that we can obtain from occultations. We call this new imaging occultation photometer HOPI (High-speed Occultation Photometer and Imager). HOPI will provide a signal-to-noise ratio two to four times that of our present photometer for a given event, will permit useful observations of the more numerous events involving fainter stars, and will allow higher time resolution observations to be made. HOPI's characteristics also make it an ideal instrument for initial and ongoing evaluations of the SOFIA telescope's performance.

Elliot, James L.↗