Search NASA⌕ Search

SEARCH · Search NASA

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

An approach to rescheduling activities based on determination of priority and disruptivity

A constraint-based scheduling system called SPIKE is being used to create long term schedules for the Hubble Space Telescope. Feedback for the spacecraft or from other ground support systems may invalidate some scheduling decisions and those activities concerned must be reconsidered. A function rescheduling priority is defined which for a given activity performs a heuristic analysis and produces a relative numerical value which is used to rank all such entities in the order that they should be rescheduled. A function disruptivity is also defined that is used to place a relative numeric value on how much a pre-existing schedule would be changed in order to reschedule an activity. Using these functions, two algorithms (a stochastic neural network approach and an exhaustive search approach) are proposed to find the best place to reschedule an activity. Prototypes were implemented and preliminary testing reveals that the exhaustive technique produces only marginally better results at much greater computational cost.

Sponsler, Jeffrey L.↗
DOE OSTI2024

Detectability of QCD phase transitions in binary neutron star mergers: Bayesian inference with the next generation gravitational wave detectors

We study the detectability of postmerger QCD phase transitions in neutron star binaries with next-generation gravitational-wave detectors Cosmic Explorer and Einstein Telescope. We perform numerical relativity simulations of neutron star mergers with equations of state that include a quark deconfinement phase transition through either a Gibbs or Maxwell construction. These are followed by Bayesian parameter estimation of the associated gravitational-wave signals using the nrpmw waveform model, with priors inferred from the analysis of the inspiral signal. We assess the ability of the model to measure the postmerger peak frequency $f$$^{peak}_{2}$ and identify aspects that should be improved in the model. We show that, even at postmerger signal to noise ratios as low as 10, the model can distinguish (at the 90% level) $f$$^{peak}_{2}$ between binaries with and without a phase transition in most cases. Phase-transition induced deviations in the $f$$^{peak}_{2}$ from the predictions of equation-of-state insensitive relations can also be detected if they exceed 1.6⁢σ. Our results suggest that next-generation gravitational wave detectors can measure phase transition effects in binary neutron star mergers. Furthermore, unless the phase transition is “strong,” disentangling it from other hadronic physics uncertainties will require significant theory improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2007

Computing Spacetimes: From Cosmology to Black Holes

Numerical relativity, the solution of the Einstein equations on a computer, is one of the most challenging and exciting areas of physics. Richard Matzner has played a key role in this subject from its birth, roughly 3 decades ago, to the present. This talk will present some of the highlights of Richard's work in numerical relativity.

Centrella, Joan↗
DOE OSTI2021

AT2017gfo: Bayesian inference and model selection of multicomponent kilonovae and constraints on the neutron star equation of state

The joint detection of the gravitational wave GW170817, of the short γ-ray burst GRB170817A and of the kilonova AT2017gfo, generated by the the binary neutron star (NS) merger observed on 2017 August 17, is a milestone in multimessenger astronomy and provides new constraints on the NS equation of state. We perform Bayesian inference and model selection on AT2017gfo using semi-analytical, multicomponents models that also account for non-spherical ejecta. Observational data favour anisotropic geometries to spherically symmetric profiles, with a log-Bayes’ factor of ~104, and favour multicomponent models against single-component ones. The best-fitting model is an anisotropic three-component composed of dynamical ejecta plus neutrino and viscous winds. Using the dynamical ejecta parameters inferred from the best-fitting model and numerical–relativity relations connecting the ejecta properties to the binary properties, we constrain the binary mass ratio to q < 1.54 and the reduced tidal parameter to $120\lt \tilde{\Lambda }\lt 1110$. Finally, we combine the predictions from AT2017gfo with those from GW170817, constraining the radius of a NS of 1.4 M ⊙ to 12.2 ± 0.5 km (1σ level). This prediction could be further strengthened by improving kilonova models with numerical-relativity information.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2010

The Final Merger of Black-Hole Binaries

Recent breakthroughs in the field of numerical relativity have led to dramatic progress in understanding the predictions of General Relativity for the dynamical interactions of two black holes in the regime of very strong gravitational fields. Such black-hole binaries are important astrophysical systems and are a key target of current and developing gravitational-wave detectors. The waveform signature of strong gravitational radiation emitted as the black holes fall together and merge provides a clear observable record of the process. After decades of slow progress / these mergers and the gravitational-wave signals they generate can now be routinely calculated using the methods of numerical relativity. We review recent advances in understanding the predicted physics of events and the consequent radiation, and discuss some of the impacts this new knowledge is having in various areas of astrophysics

Kelly, Bernard J.↗
DOE OSTI2026

AthenaK: A Performance-portable Version of the Athena++ Adaptive Mesh Refinement Framework

We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic, and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers, we describe (1) a new solver for the Einstein equations based on the Z4c formalism, and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging Advanced RISC Machine processors. AthenaK shows excellent performance and weak scaling, achieving over 1 billion cell updates per second for hydrodynamics in three dimensions on a single NVIDIA Grace Hopper processor. It does this with a typical parallel efficiency of 80% on 65,536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Neutrino Flavor Transformation in Neutron Star Mergers

Here, we present the first numerical relativity simulations including neutrino flavor transformations that could result from flavor instabilities, quantum many-body effects, or potential beyond standard model physics in neutron star mergers. We find that neutrino flavor transformations impact the composition and structure of the remnant, potentially leaving an imprint on the postmerger gravitational-wave signal. They also have a significant impact on the composition and nucleosynthesis yields of the ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Conservative Numerical Schemes with Optimal Dispersive Wave Relations: Part II. Numerical Evaluations

A new energy and enstrophy conserving scheme (EEC) for the shallow water equations is proposed and evaluated using a suite of test cases over the global spherical or bounded domain. The evaluation is organized around a set of pre-defined properties: accuracy of individual operators, accuracy of the whole scheme, conservation of key quantities, control of the divergence variable, representation of the energy and enstrophy spectra, and simulation of nonlinear dynamics. The results confirm that the scheme is between the first and second order accurate, and conserves the total energy and potential enstrophy up to the time truncation errors. Here, the scheme is capable of producing more physically realistic energy and enstrophy spectra, indicating that it can help prevent the unphysical energy cascade towards the finest resolvable scales. With an optimal representation of the dispersive wave relations, the scheme is able to keep the flow close to being non-divergent, and maintain the geostrophically balanced structures with large-scale geophysical flows over long-term simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Oscillon preheating in full general relativity

Abstract Oscillons are dense nonperturbative objects that may be copiously produced in the preheating period after inflation. Oscillon preheatings are usually simulated with nonlinear matter interactions but in a rigid FLRW background, without taking into account the gravitational backreactions. We investigate the oscillon preheating scenario in full general relativity with a full numerical relativity scheme, and compare the general-relativistic simulations with the traditional ones done in an FLRW background. We find that in certain parameter regions the general-relativistic corrections are significant, producing more and denser oscillons, and can be strong enough to collapse the oscillons to black holes.

Astronomy & Astrophysics↗
DOE OSTI2023

Second release of the CoRe database of binary neutron star merger waveforms

Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1967

Ordering of implicants.

Correspondence between logical relations of Boolean function implicants and numerical relations between identifiers

BOOLEAN ALGEBRA↗
DOE OSTI2021

Signatures of deconfined quark phases in binary neutron star mergers

We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We employ a new finite-temperature composition-dependent equation of state (EOS) with a first order phase transition between hadrons and deconfined quarks to perform numerical relativity simulations of BNS mergers. Here, the softening of the EOS due to the phase transition causes the merger remnants to be more compact and to collapse to a black hole (BH) at earlier times. The phase transition is imprinted on the postmerger gravitational wave (GW) signal duration, amplitude, and peak frequency. However, this imprint is only detectable for binaries with sufficiently long-lived remnants. Moreover, the phase transition does not result in significant deviations from quasi-universal relations for the postmerger GW peak frequency. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, r-process nucleosynthetic yields and associated electromagnetic (EM) counterparts. While there are differences in the EM counterparts and nucleosynthesis yields between the purely hadronic models and the models with phase transitions, these can be primarily ascribed to the difference in remnant collapse time between the two. An exception is the non-thermal afterglow caused by the interaction of the fastest component of the dynamical ejecta and the interstellar medium, which is systematically boosted in the binaries with phase transition as a consequence of the more violent merger they experience.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reverse phase transitions in binary neutron-star systems with exotic-matter cores

Multimessenger observations of binary neutron-star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter—e.g., quark matter—the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open question if such phase transitions occur inside isolated neutron stars or during binary neutron-star mergers, or if they appear at even higher densities that are not realized in the cosmos. Here, in this article, we perform numerical relativity simulations of neutron-star mergers and investigate scenarios in which the onset density of such a phase transition is exceeded in at least one inspiraling binary component. Our simulations reveal that shortly before the merger, it is possible that such stars undergo a “reverse phase transition”—i.e., densities decrease and the quark core inside the star disappears, leaving a purely hadronic star at merger. After the merger, when densities increase once more, the phase transition occurs again and leads, in the cases considered in this work, to the rapid formation of a black hole. We compute the gravitational-wave signal and the mass ejection for our simulations of such scenarios and find clear signatures that are related to the postmerger phase transition—e.g., smaller ejecta masses due to the softening of the equation of state through the quark core formation. Unfortunately, we do not find measurable imprints of the reverse phase transition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dark magnetohydrodynamics: Black hole accretion in superradiant dark photon clouds

Black holes threaded by massive vector fields can be subject to a superradiant instability, growing a cloud of massive vector particles around it. In this work, we consider what happens if such a dark matter candidate field mimicking a dark photon interacts with an accretion flow onto the black hole. By including a kinetic mixing term with the standard model photon, we extend the commonly used equations of general-relativistic magnetohydrodynamics to a dark photon constituent. The coupling to the dark photon then appears as an effective dynamo term together with a dark Lorentz force acting on the accreting matter. We numerically study the interactions between the superradiant dark photon cloud and the inner accretion flow by solving the coupled system in full numerical relativity. By parameterically varying the mixing parameter between the dark and standard model sector, we provide a first investigation of how the accretion flow could be modified. In conclusion, depending on the coupling strength, our solutions exhibit increased wind launching, as well as oscillation modes in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Impact of neutrino flavor conversions on neutron star merger dynamics, ejecta, nucleosynthesis, and multimessenger signals

We present numerical relativity simulations of binary neutron star mergers incorporating neutrino flavor transformations triggered by fast flavor instability, quantum many-body effects, or potential beyond standard model physics. In both long-lived and short-lived remnant scenarios, neutrino flavor conversions modify species-dependent neutrino luminosities and mean energies, and drive the matter towards more neutron rich conditions. They produce up to more neutron rich ejecta and significantly boost the r-process yields, especially in low-density, near-equatorial outflows. We identify regions unstable to fast flavor instabilities and find that these instabilities persist despite flavor conversions. We further test the sensitivity to the equilibration timescale of the flavor conversions, finding that slower flavor conversions can interact with thermodynamic equilibration, and increase the neutron richness of the ejecta. Flavor conversions may also contribute to stronger gravitational wave and neutrino emissions, pointing to a correlation between neutrino transport and merger dynamics. Furthermore, these results highlight the potential impact of flavor conversions while motivating future work to improve on theoretical understanding of flavor instabilities in global simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1937

Tests of related forward-camber airfoils in the variable-density wind tunnel

A recent investigation of numerous related airfoils indicated that positions of camber forward of the usual location resulted in an increase of the maximum lift. As an extension of this investigation, a series of forward-camber airfoils has been developed, the members of which show airfoil characteristics superior to those of the airfoils previously investigated. The primary object of this report is to present fully corrected results for airfoils in the useful range of shapes. With the data thus made available, an airplane designer may intelligently choose the best possible airfoil-section shape for a given application and may predict to a reasonable degree the aerodynamic characteristics to be expected in flight from the section shape chosen.

Jacobs, Eastman N↗
DOE OSTI2024

Do black holes remember what they are made of?

We study the ringdown signal of black holes formed in prompt-collapse binary neutron star mergers. We analyze data from 47 numerical relativity simulations. We show that the ($l$ = 2, $m$ = 2) and ($l$ = 2, $m$ = 1) multipoles of the gravitational wave signal are well fitted by decaying damped exponentials, as predicted by black-hole perturbation theory. We show that the ratio of the amplitude in the two modes depends on the progenitor binary mass ratio q and reduced tidal parameter $\overline{Λ}$. Unfortunately, the numerical uncertainty in our data is too large to fully quantify this dependency. If confirmed, these results will enable novel tests of general relativity in the presence of matter with next-generation gravitational-wave observatories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

SERS and the scientific career of Richard P. Van Duyne (1945–2019)

This article describes the scientific career of our dear friend Richard P. Van Duyne. Rick was the discoverer of surface-enhanced Raman spectroscopy (SERS), and over the period from the mid-1970s until his death in 2019, he served as both the intellectual leader of the field, and the technical genius who pushed the field forward through his numerous inventions and discoveries. Rick's work went through many phases, and even involved a significant period of inactivity that Rick was able to break out of, but especially after 2000, there was a period of incredible creativity that ultimately led to a detailed understanding of what the SERS enhancement factor is all about and also led to fabrication techniques for making optimum SERS substrates such that applications of SERS to a wide variety of real-world problems in sensing became possible. In addition, Rick's SERS work spawned numerous related technologies, including index of refraction sensing based on localized surface plasmon resonance (LSPR) measurements, and the extension of SERS to its tip-enhanced cousin TERS, to electrochemical versions of SERS and TERS, and to ultrafast nonlinear techniques. The practical applications he pursued included diagnostics for diseases such as cancer, hand-held glucose sensing, and SERS studies related to art conservation. Furthermore, Rick's legacy will probably be best known through the numerous students, postdocs, and collaborators that passed through his lab and who are now active researchers in areas related to SERS at universities and industrial labs around the world.

47 OTHER INSTRUMENTATION↗