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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Numerical Relativity, Black Hole Mergers, and Gravitational Waves: Part II

This series of 3 lectures will present recent developments in numerical relativity, and their applications to simulating black hole mergers and computing the resulting gravitational waveforms. In this second lecture, we focus on simulations of black hole binary mergers. We hig hlight the instabilities that plagued the codes for many years, the r ecent breakthroughs that led to the first accurate simulations, and the current state of the art.

Centrella, Joan↗
NASA NTRS2010

Numerical Relativity Simulations for Black Hole Merger Astrophysics

Massive black hole mergers are perhaps the most energetic astronomical events, establishing their importance as gravitational wave sources for LISA, and also possibly leading to observable influences on their local environments. Advances in numerical relativity over the last five years have fueled the development of a rich physical understanding of general relativity's predictions for these events. Z will overview the understanding of these event emerging from numerical simulation studies. These simulations elucidate the pre-merger dynamics of the black hole binaries, the consequent gravitational waveform signatures ' and the resulting state, including its kick velocity, for the final black hole produced by the merger. Scenarios are now being considered for observing each of these aspects of the merger, involving both gravitational-wave and electromagnetic astronomy.

Baker, John G.↗
DOE OSTI2023

Numerical relativity estimates of the remnant recoil velocity in binary neutron star mergers

We present, for the first time, recoil velocity estimates for binary neutron star mergers using data from numerical relativity simulations. We find that binary neutron star merger remnants can have recoil velocity of the order of a few tens of kilometers per second and as high as 150 kilometers per second in our dataset. These recoils are attained due to equivalent contributions from the anisotropic gravitational wave emission as well as the asymmetric ejection of dynamical matter during the merger. Here, we provide fits for net recoil velocity as well as its ejecta component as a function of the amount of ejected matter, which may be useful when constraints on the ejected matter are obtained through electromagnetic observations. We also estimate the mass and spin of the remnants and find them to be in the range [2.34,3.38]M ⊙ and [0.63, 0.82], respectively, for our dataset.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2009

Merging Black Holes, Gravitational Waves, and Numerical Relativity

The final merger of two black holes will emit more energy than all the stars in the observable universe combined. This energy will come in the form of gravitational waves, which are a key prediction of Einstein's general relativity and a new tool for exploring the universe. Observing these mergers with gravitational wave detectors, such as the ground-based LIGO and the space-based LISA, requires knowledge of the radiation waveforms. Since these mergers take place in regions of extreme gravity, we need to solve Einstein's equations of general relativity on a computer. For more than 30 years, scientists have tried to compute black hole mergers using the methods of numerical relativity. The resulting computer codes were long plagued by instabilities, causing them to crash well before the black holes in the binary could complete even a single orbit. Within the past few years, however, this situation has changed dramatically, with a series of remarkable breakthroughs. This talk will focus on new simulations that are revealing the dynamics and w ae~fo rms of binary black hole mergers, and their applications in gravitational wave detection, testing general relativity, and astrophysics.

Centrella, Joan M.↗
NASA NTRS2007

Binary Black Holes, Gravitational Waves, and Numerical Relativity

Massive black hole (MBH) binaries are found at the centers of most galaxies. MBH mergers trace galaxy mergers and are strong sources of gravitational waves. Observing these sources with gravitational wave detectors requires that we know the radiation waveforms they emit. Since these mergers take place in regions of very strong gravitational fields, we need to solve Einstein's equations of general relativity on a computer in order to calculate these waveforms. For more than 30 years, scientists have tried to compute these waveforms using the methods of numerical relativity. The resulting computer codes have been plagued by instabilities. causing them to crash well before the black hole:, in the binary could complete even a single orbit. Recently this situation has changed dramatically, with a series of amazing breakthroughs. This presentation shows how a spacetime is constructed on a computer to build a simulation laboratory for binary black hole mergers. Focus is on the recent advances that that reveal these waveforms, and the potential for discoveries that arises when these sources are observed by LIGO and LISA.

Centrella, Joan↗
NASA NTRS2007

Binary Black Holes, Gravitational Waves, and Numerical Relativity

The final merger of two black holes is expected to be the strongest gravitational wave source for ground-based interferometers such as LIGO, VIRGO, and GE0600, as well as the space-based interferometer LISA. Observing these sources with gravitational wave detectors requires that we know the radiation waveforms they emit. Since these mergers take place in regions of extreme gravity, we need to solve Einstein's equations of general relativity on a computer in order to calculate these waveforms. For more than 30 years, scientists have tried to compute black hole mergers using the methods of numerical relativity. The resulting computer codes have been plagued by instabilities, causing them to crash well before the black holes in the binary could complete even a single orbit. Within the past few years, however, this situation has changed dramatically, with a series of remarkable breakthroughs. This talk will focus on new simulations that are revealing the dynamics and waveforms of binary black hole mergers, and their applications in gravitational wave detection, data analysis, and astrophysics.

Centrella, John↗
NASA NTRS2010

Black-Hole Binaries, Gravitational Waves, and Numerical Relativity

Understanding the predictions of general relativity for the dynamical interactions of two black holes has been a long-standing unsolved problem in theoretical physics. Black-hole mergers are monumental astrophysical events ' releasing tremendous amounts of energy in the form of gravitational radiation ' and are key sources for both ground- and spacebased gravitational wave detectors. The black-hole merger dynamics and the resulting gravitational waveforms can only he calculated through numerical simulations of Einstein's equations of general relativity. For many years, numerical relativists attempting to model these mergers encountered a host of problems, causing their codes to crash after just a fraction of a binary orbit cnuld be simulated. Recently ' however, a series of dramatic advances in numerical relativity has ' for the first time, allowed stable / robust black hole merger simulations. We chronicle this remarkable progress in the rapidly maturing field of numerical relativity, and the new understanding of black-hole binary dynamics that is emerging. We also discuss important applications of these fundamental physics results to astrophysics, to gravitationalwave astronomy, and in other areas.

Kelly, Bernard J.↗
DOE OSTI2022

Numerical relativity simulations of prompt collapse mergers: Threshold mass and phenomenological constraints on neutron star properties after GW170817

We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite-temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, binary tidal parameter at the threshold (Λth), maximum mass of nonrotating NSs, and radii of reference mass NSs. We combine the EOS-insensitive relations, phenomenological constraints on NS properties, and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and a 1.6 M ⊙ NS of 9.81 and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability (Λ) of a 1.4 M ⊙ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters—fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We discuss an approach to constrain the upper as well as lower limit of NS maximum mass using future binary NS detections and their identification as prompt or delayed collapse. With future observations, it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2011

Numerical Relativity for Space-Based Gravitational Wave Astronomy

In the next decade, gravitational wave instruments in space may provide high-precision measurements of gravitational-wave signals from strong sources, such as black holes. Currently variations on the original Laser Interferometer Space Antenna mission concepts are under study in the hope of reducing costs. Even the observations of a reduced instrument may place strong demands on numerical relativity capabilities. Possible advances in the coming years may fuel a new generation of codes ready to confront these challenges.

Baker, John G.↗
DOE OSTI2022

Numerical relativity simulations of the neutron star merger GW190425: microphysics and mass ratio effects

ABSTRACT GW190425 was the second gravitational wave (GW) signal compatible with a binary neutron star (BNS) merger detected by the Advanced LIGO and Advanced Virgo detectors. Since no electromagnetic counterpart was identified, whether the associated kilonova was too dim or the localization area too broad is still an open question. We simulate 28 BNS mergers with the chirp mass of GW190425 and mass ratio 1 ≤ q ≤ 1.67, using numerical-relativity simulations with finite-temperature, composition dependent equations of state (EOS) and neutrino radiation. The energy emitted in GWs is $\lesssim 0.083\mathrm{\, M_\odot }c^2$ with peak luminosity of 1.1–$2.4\times ~10^{58}/(1+q)^2\, {\rm {erg \, s^{-1}}}$. Dynamical ejecta and disc mass range between 5 × 10−6–10−3 and 10−5–$0.1 \mathrm{\, M_\odot }$, respectively. Asymmetric mergers, especially with stiff EOSs, unbind more matter and form heavier discs compared to equal mass binaries. The angular momentum of the disc is 8–$10\mathrm{\, M_\odot }~GM_{\rm {disc}}/c$ over three orders of magnitude in Mdisc. While the nucleosynthesis shows no peculiarity, the simulated kilonovae are relatively dim compared with GW170817. For distances compatible with GW190425, AB magnitudes are always dimmer than ∼20 mag for the B, r, and K bands, with brighter kilonovae associated to more asymmetric binaries and stiffer EOSs. We suggest that, even assuming a good coverage of GW190425’s sky location, the kilonova could hardly have been detected by present wide-field surveys and no firm constraints on the binary parameters or EOS can be argued from the lack of the detection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
DOE OSTI2023

HPC-driven computational reproducibility in numerical relativity codes: a use case study with IllinoisGRMHD

Abstract Reproducibility of results is a cornerstone of the scientific method. Scientific computing encounters two challenges when aiming for this goal. Firstly, reproducibility should not depend on details of the runtime environment, such as the compiler version or computing environment, so results are verifiable by third-parties. Secondly, different versions of software code executed in the same runtime environment should produceconsistent numerical results for physical quantities. In this manuscript, we test the feasibility of reproducing scientific results obtained using theIllinoisGRMHDcode that is part of an open-source community software for simulation in relativistic astrophysics, theEinstein Toolkit. We verify that numerical results of simulating a single isolated neutron star withIllinoisGRMHDcan be reproduced, and compare them to results reported by the code authors in 2015. We use two different supercomputers: Expanse at SDSC, and Stampede2 at TACC. By compiling the source code archived along with the paper on both Expanse and Stampede2, we find thatIllinoisGRMHDreproduces results published in its announcement paper up to errors comparable to round-off level changes in initial data parameters. We also verify that a current version ofIllinoisGRMHDreproduces these results once we account for bug fixes which have occurred since the original publication.

Astronomy & Astrophysics↗
NASA NTRS2006

Consistency of Post-Newtonian Waveforms with Numerical Relativity

General relativity predicts the gravitational radiation signatures of mergers of compact binaries, such as coalescing binary black hole systems. Derivations of waveform predictions for such systems are required for optimal scientific analysis of observational gravitational wave data, and have so far been achieved primarily with the aid of the post-Newtonian (PN) approximation. The quality of this treatment is unclear, however, for the important late inspiral portion. We derive late-inspiral waveforms via a complementary approach, direct numerical simulation of Einstein's equations, which has recently matured sufficiently for such applications. We compare waveform phasing from simulations covering the last approximately 14 cycles of gravitational radiation from an equal-mass binary system of nonspinning black holes with the corresponding 3PN and 3.5PN orbital phasing. We find agreement consistent with internal error estimates based on either approach at the level of one radian over approximately 10 cycles. The result suggests that PN waveforms for this system are effective roughly until the system reaches its last stable orbit just prior to the final merger/

Baker, John G.↗
NASA NTRS2007

Consistency of Post-Newtonian Waveforms with Numerical Relativity

General relativity predicts the gravitational radiation signatures of mergers of compact binaries,such as coalescing binary black hole systems. Derivations of waveform predictions for such systems are required for optimal scientific analysis of observational gravitational wave data, and have so far been achieved primarily with the aid of the post-Newtonian (PN) approximation. The quaIity of this treatment is unclear, however, for the important late inspiral portion. We derive late-inspiral wave forms via a complementary approach, direct numerical simulation of Einstein's equations, which has recently matured sufficiently for such applications. We compare waveform phasing from simulations covering the last approximately 14 cycles of gravitational radiation from an equal-mass binary system of nonspinning black holes with corresponding 3PN and 3.5PN waveforms. We find phasing agreement consistent with internal error estimates based in either approach, at the level of one radian over approximately 10 cycles. The result suggests that PN waveforms for this system are effective roughly until the system reaches its last stable orbit just prior to the final merger.

Baker, John G.↗
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↗
NASA NTRS1986

Analysis of coaxial spray combustion flames and related numerical issues

An approach to the simulation of strongly coupled multiphase flows in combustion hardware is sketched and its unique requirements highlighted. An example of a successful application to a coaxial injector flame is presented. Furthermore, several numerical issues that tend to interact with the physics of the problem are discussed with special regard to their potential impact on the choices of numerical parameters by the analyst. These include the issues of stability, numerical diffusivity, stiffness, and boundary conditions. The theme of this paper focuses on the intriguing relationships among the grid, the solution algorithm, and the actual physical mechanisms themselves.

Liang, P. Y.↗