Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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

Correspondence of NNGP Kernel and the Matérn Kernel

Kernels representing limiting cases of neural network architectures have recently gained popularity. However, the application and performance of these new kernels compared to existing options, such as the Matérn kernel, is not well studied. We take a practical approach to explore the neural network Gaussian process (NNGP) kernel and its application to data in Gaussian process regression. We first demonstrate the necessity of normalization to produce valid NNGP kernels and explore related numerical challenges. We further demonstrate that the predictions from this model are quite inflexible, and therefore do not vary much over the valid hyperparameter sets. We then demonstrate a surprising result that the predictions given from the NNGP kernel correspond closely to those given by the Matérn kernel under specific circumstances, which suggests a deep similarity between overparameterized deep neural networks and the Matérn kernel. Finally, we demonstrate the performance of the NNGP kernel as compared to the Matérn kernel on three benchmark data cases, and we conclude that for its flexibility and practical performance, the Matérn kernel is preferred to the novel NNGP in practical applications.

97 MATHEMATICS AND COMPUTING↗

Physics model enhancements in MARLEY version 2

Large next-generation neutrino detectors, such as those designed for accelerator-based oscillation experiments, are expected to have impressive sensitivity to MeV astrophysical neutrinos generated by the Sun and by supernovae. There are also numerous related prospects to search for MeV-scale phenomena beyond the Standard Model. In contrast to detectors constructed with water and hydrocarbon targets, which will be primarily sensitive to electron antineutrinos via inverse beta decay, the upcoming Deep Underground Neutrino Experiment (DUNE) will provide unparalleled capabilities to measure low-energy electron neutrinos from cosmic sources. Interpreting such future measurements in DUNE will require a high-quality simulation of MeV interaction physics in the argon nucleus. The MARLEY (Model of Argon Reaction Low Energy Yields) neutrino event generator has been designed with the goal of delivering the needed simulation capabilities for DUNE and other argon-based neutrino detectors. In this poster, I will present recent improvements to MARLEY's suite of interaction models, which will be released for general use by the neutrino community in version 2 of the software.

Gardiner, Steven James [Fermilab] (ORCID:000000028↗

Physics model enhancements in MARLEY version 2

Large next-generation neutrino detectors, such as those designed for accelerator-based oscillation experiments, are expected to have impressive sensitivity to MeV astrophysical neutrinos generated by the Sun and by supernovae. There are also numerous related prospects to search for MeV-scale phenomena beyond the Standard Model. In contrast to detectors constructed with water and hydrocarbon targets, which will be primarily sensitive to electron antineutrinos via inverse beta decay, the upcoming Deep Underground Neutrino Experiment (DUNE) will provide unparalleled capabilities to measure low-energy electron neutrinos from cosmic sources. Interpreting such future measurements in DUNE will require a high-quality simulation of MeV interaction physics in the argon nucleus. The MARLEY (Model of Argon Reaction Low Energy Yields) neutrino event generator has been designed with the goal of delivering the needed simulation capabilities for DUNE and other argon-based neutrino detectors. In this talk, I will present recent improvements to MARLEY's physics models, which will be released for general use by the neutrino community in version 2 of the software.

Gardiner, Steven James [Fermilab] (ORCID:000000028↗

Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗

Electromagnetic Signatures from the Tidal Tail of a Black Hole—Neutron Star Merger

Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle because of two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of ~3 with orientation while the peak in the optical bands varies by ~3 magnitudes. The spectrum is blueshifted at viewing angles along the bulk motion, which increases the V-band peak magnitude to ~-14 despite the lanthanide-rich composition.

79 ASTRONOMY AND ASTROPHYSICS↗

Production of Very Light Elements and Strontium in the Early Ejecta of Neutron Star Mergers

Abstract We study the production of very light elements ( Z < 20) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition-dependent nuclear equations of state, and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, ∼0.5–2 × 10 −6 M ⊙ of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10 −6 M ⊙ of helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavy r -process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including nonlocal thermodynamic equilibrium effects. The latter could be crucial to observe helium lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long-lived remnant) the one required to explain early spectral features in the kilonova of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗

Tight Bound on the Neutron Star Radius with Quasiperiodic Oscillations in Short Gamma-Ray Bursts

Abstract Quasiperiodic oscillations (QPOs) have been recently discovered in the short gamma-ray bursts (GRBs) 910711 and 931101B. Their frequencies are consistent with those of the quasiradial and quadrupolar oscillations of binary neutron star (BNS) merger remnants, as obtained in numerical relativity simulations. These simulations reveal quasi-universal relations between the remnant oscillation frequencies and the tidal coupling constant of the binaries. Under the assumption that the observed QPOs are due to these postmerger oscillations, we use the frequency–tide relations in a Bayesian framework to infer the source redshift, as well as the chirp mass and the binary tidal deformability of the BNS progenitors for GRBs 910711 and 931101B. We further use this inference to estimate bounds on the mass–radius relation for neutron stars. By combining the estimates from the two GRBs, we find a 68% credible range R 1.4 = 12.4 8 − 0.40 + 0.41 km for the radius of a neutron star with mass M = 1.4 M ⊙ , which is one of the tightest bounds to date.

Guedes, Victor↗

Evidence for X-Ray Emission in Excess to the Jet-afterglow Decay 3.5 yr after the Binary Neutron Star Merger GW 170817: A New Emission Component

For the first ~3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n < 0.01 cm -3 . We report on observational evidence for an excess of X-ray emission at δt > 900 days after the merger. With L x ≈ 5 × 10 38 erg s -1 at 1234 days, the recently detected X-ray emission represents a ≥ 3.2σ (Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance ≥ 3.1σ, depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥ 3.7σ. A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented

79 ASTRONOMY AND ASTROPHYSICS↗

Thermal Effects in Binary Neutron Star Mergers

Abstract We study the impact of finite-temperature effects in numerical-relativity simulations of binary neutron star mergers with microphysical equations of state and neutrino transport in which we vary the effective nucleon masses in a controlled way. We find that, as the specific heat is increased, the merger remnants become colder and more compact due to the reduced thermal pressure support. Using a full Bayesian analysis, we demonstrate that this effect will be measurable in the postmerger gravitational wave signal with next-generation observatories at signal-to-noise ratios of 15.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hours-long Near-UV/Optical Emission from Mildly Relativistic Outflows in Black Hole–Neutron Star Mergers

The ongoing LIGO–Virgo–KAGRA observing run O4 provides an opportunity to discover new multimessenger events, including binary neutron star (BNS) mergers such as GW170817 and the highly anticipated first detection of a multimessenger black hole–neutron star (BH–NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH–NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH–NS mergers. Adopting optimal binary properties, a mass ratio of q = 2, and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary's evolution from premerger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of ~–15 a few hours after the merger, longer than previous estimates, which did not consider the first principles–based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of ≳1 Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within ~200 Mpc. The BH–NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals.

79 ASTRONOMY AND ASTROPHYSICS↗

What Can We Learn about the Unstable Equation-of-state Branch from Neutron Star Mergers?

Abstract The equation of state (EOS) of dense strongly interacting matter can be probed by astrophysical observations of neutron stars (NS), such as X-ray detections of pulsars or the measurement of the tidal deformability of NSs during the inspiral stage of NS mergers. These observations constrain the EOS at most up to the density of the maximum-mass configuration, n TOV , which is the highest density that can be explored by stable NSs for a given EOS. However, under the right circumstances, binary neutron star (BNS) mergers can create a postmerger remnant that explores densities above n TOV . In this work, we explore whether the EOS above n TOV can be measured from gravitational-wave or electromagnetic observations of the postmerger remnant. We perform a total of 25 numerical-relativity simulations of BNS mergers for a range of EOSs and find no case in which different descriptions of the matter above n TOV have a detectable impact on postmerger observables. Hence, we conclude that the EOS above n TOV can likely not be probed through BNS merger observations for the current and next generation of detectors.

79 ASTRONOMY AND ASTROPHYSICS↗

Element Formation in Radiation-hydrodynamics Simulations of Kilonovae

Abstract Understanding the details of r -process nucleosynthesis in binary neutron star merger (BNSM) ejecta is key to interpreting kilonova observations and identifying the role of BNSMs in the origin of heavy elements. We present a self-consistent, two-dimensional, ray-by-ray radiation-hydrodynamic evolution of BNSM ejecta with an online nuclear network (NN) up to a timescale of days. For the first time, an initial numerical relativity ejecta profile composed of the dynamical component and spiral-wave and disk winds is evolved including detailed r -process reactions and nuclear heating effects. A simple model for the jet energy deposition is also included. Our simulation highlights that the common approach of relating in postprocessing the final nucleosynthesis yields to the initial thermodynamic profile of the ejecta can lead to inaccurate predictions. Moreover, we find that neglecting the details of the radiation-hydrodynamic evolution of the ejecta in nuclear calculations can introduce deviations of up to 1 order of magnitude in the final abundances of several elements, including very light and second r -process peak elements. The presence of a jet affects element production only in the innermost part of the polar ejecta, and it does not alter the global nucleosynthesis results. Overall, our analysis shows that employing an online NN improves the reliability of nucleosynthesis and kilonova light-curve predictions.

Magistrelli, Fabio (ORCID:0009000509767851)↗

Impact of Magnetic-field-driven Anisotropies on the Equation of State Probed in Neutron Star Mergers

Binary neutron star mergers can produce extreme magnetic fields, some of which can lead to strong magnetar-like remnants. While strong magnetic fields have been shown to affect the dynamics of outflows and angular momentum transport in the remnant, they can also crucially alter the properties of nuclear matter probed in the merger. In this work, we provide a first assessment of the latter, determining the strength of the pressure anisotropy caused by Landau-level quantization and the anomalous magnetic moment. To this end, we perform the first numerical relativity simulation with a magnetic polarization tensor and a magnetic-field-dependent equation of state using a new algorithm we present here, which also incorporates a mean-field dynamo model to control the magnetic field strength present in the merger remnant. Our results show that—in the most optimistic case—corrections to the anisotropy can be in excess of 10% and are potentially largest in the outer layers of the remnant. This work paves the way for a systematic investigation of these effects.

General relativity↗

Overview of NASA aeronautical propulsion research and technology program

The program discussed is aimed at improving performance within an extended operating range, reducing weight, increasing service life, achieving greater cost effectiveness, reducing noise and exhaust pollution, and improving means of energy conservation. The program places emphasis on basic research in numerous related technical disciplines, system analysis, component technology, full scale propulsion system studies, and technology demonstrations. Much attention in the discussion is given to noise and pollution minimization research.

Johnson, H. W.↗

Prediction of the diffuse far infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far infrared studies of the galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far infrared flux distribution from the galactic plane as a function of galactic longitude alternate theta for 4 deg or = alternate theta or = 90 and the far infrared emissivity as a function of galactocentric distance is predicted. Future measurements of the galactic far infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the galaxy, particularly the inner region.

Fazio, G. G.↗

Prediction of the diffuse far-infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far-infrared studies of the Galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far-infrared flux distribution from the galactic plane is predicted as a function of galactic longitude for longitudes between 4 and 90 deg; the far-infrared emissivity is predicted as a function of galactocentric distance. Future measurements of the galactic far-infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the Galaxy, particularly the inner region.

Fazio, G. G.↗

Prediction of the diffuse far-infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far-infrared studies of the Galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far-infrared flux distribution from the galactic planes is predicted as a function of galactic longitude for the region between 4 and 90 deg longitude; the far-infrared emissivity is predicted as a function of galactocentric distance. Future measurements of the galactic far-infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the Galaxy, particularly the inner region.

Fazio, G. G.↗

Terrestrial photovoltaic measurements, 2

The following major topics are discussed; (1) Terrestrial solar irradiance; (2) Solar simulation and reference cell calibration; and (3) Cell and array measurement procedures. Numerous related subtopics are also discussed within each major topic area.

Source record↗