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 55 records · Page 3

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.↗

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↗

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.↗

Binary Black Holes, Numerical Relativity, and Gravitational Waves

The final merger of two black holes releases a tremendous amount of energy, more than the combined light from all the stars in the visible universe. This energy is emitted in the form of gravitational waves, and observing these sources with gravitational wave detectors such as LISA requires that we know the pattern or fingerprint of the radiation emitted. Since black hole mergers take place in regions of extreme gravitational fields, we need to solve Einstein's equations of general relativity on a computer in order to calculate these wave patterns. For more than 30 years, scientists have tried to compute these wave patterns. However, their computer codes have been plagued by problems that caused them to crash. This situation has changed dramatically in the past 2 years, with a series of amazing breakthroughs. This talk will take you on this quest for these gravitational wave patterns, showing how a spacetime is constructed on a computer to build a simulation laboratory for binary black hole mergers. We will focus on the recent advances that are revealing these waveforms, and the dramatic new potential for discoveries that arises when these sources will be observed by LISA

Centrella, Joan↗

Results from Numerical General Relativity

For several years numerical simulations have been revealing the details of general relativity's predictions for the dynamical interactions of merging black holes. I will review what has been learned of the rich phenomenology of these mergers and the resulting gravitational wave signatures. These wave forms provide a potentially observable record of the powerful astronomical events, a central target of gravitational wave astronomy. Asymmetric radiation can produce a thrust on the system which may accelerate the single black hole resulting from the merger to high relative velocity.

Baker, John G.↗

Simulating Gravitational Radiation from Binary Black Holes Mergers as LISA Sources

A viewgraph presentation on the simulation of gravitational waves from Binary Massive Black Holes with LISA observations is shown. The topics include: 1) Massive Black Holes (MBHs); 2) MBH Binaries; 3) Gravitational Wavws from MBH Binaries; 4) Observing with LISA; 5) How LISA sees MBH binary mergers; 6) MBH binary inspirals to LISA; 7) Numerical Relativity Simulations; 8) Numerical Relativity Challenges; 9) Recent Successes; 10) Goddard Team; 11) Binary Black Hole Simulations at Goddard; 12) Goddard Recent Advances; 13) Baker, et al.:GSFC; 13) Starting Farther Out; 14) Comparing Initial Separation; 15) Now with AMR; and 16) Conclusion.

Baker, John↗

Statistical and systematic uncertainties in extracting the source properties of neutron star-black hole binaries with gravitational waves

Gravitational waves emitted by neutron star black hole mergers encode key properties of neutron stars - such as their size, maximum mass and spins - and black holes. However, the presence of matter and the high mass ratio makes generating long and accurate waveforms from these systems hard to do with numerical relativity, and not much is known about systematic uncertainties due to waveform modeling. In this work, we simulate gravitational waves from neutron star black hole mergers by hybridizing numerical relativity waveforms produced with the SpEC code with a recent numerical relativity surrogate NRHybSur3dq8Tidal. These signals are analyzed using a range of available waveform families, and statistical and systematic errors are reported. We find that at a network signal-to-noise ratio (SNR) of 30, statistical uncertainties are usually larger than systematic offsets, while at an SNR of 70 the two become comparable. The individual black hole and neutron star masses, as well as the mass ratios, are typically measured very precisely, though not always accurately at high SNR. At a SNR of 30 the neutron star tidal deformability can only be bound from above, while for louder sources it can be measured and constrained away from zero. All neutron stars in our simulations are non-spinning, but in no case we can constrain the neutron star spin to be smaller than ~0.4 (90% credible interval). Waveform families whose late inspiral has been tuned specifically for neutron star black hole signals typically yield the most accurate characterization of the source parameters. Their measurements are in tension with those obtained using waveform families tuned against binary neutron stars, even for mass ratios that could be relevant for both binary neutron stars and neutron star black holes mergers.

79 ASTRONOMY AND ASTROPHYSICS↗

Observing Mergers of Nonspinning Black Hole Binaries with LISA

Recent advances in the field of numerical relativity now make it possible to calculate the final, most powerful merger phase of binary black hole coalescence. We present the application of nonspinning numerical relativity waveforms to the search for and precision measurement of black hole binary coalescences using LISA. In particular, we focus on the advances made in moving beyond the equal mass, nonspinning case into other regions of parameter space, focusing on the case of nonspinning holes with ever-increasing mass ratios. We analyze the available unequal mass merger waveforms from numerical relativity, and compare them to two models, both of which use an effective one body treatment of the inspiral, but which use fundamentally different approaches to the treatment of the merger-ringdown. We confirm the expected mass ratio scaling of the merger, and investigate the changes in waveform behavior and their observational impact with changing mass ratio. Finally, we investigate the potential contribution from the merger portion of the waveform to measurement uncertainties of the binary's parameters for the unequal mass case.

McWilliams S.↗

Canopy tree mortality and crown exposure data from the Amacayacu Forest Dynamics Plot, Northwestern Amazon

Data on the mortality of 984 canopy trees, their crown exposure to light (relative to total crown area), growth deviations (relative to conspecifics), tree size, and species’ wood density collected between 2013 and 2019 in 18 ha of the Amacayacu Forest Dynamics Plot, Northwestern Amazon. This dataset contains a single CSV data file. Variable definitions: 1. Species: [character] species identification 2. Family: [character] family of the species 3. tag: [character] unique consecutive for the tree 4. status: [character] status of the tree in the third census (ALIVE or DEAD) 5. wsg: [numeric]: species’ wood density (g cm-3) 6. growth_r1: [numeric]: annual growth rate of the tree between the first and second census (cm y-1) 7. growth_r2: [numeric]: annual growth rate of the tree between the second and third census (cm y-1) 8. gt1: [numeric] modulus transformed growth rate with a lambda of 0.4 for the annual growth rate of the tree between the first and second census (cm y-1) 9. gt2: [numeric]: modulus transformed growth rate with a lambda of 0.4 for the annual growth rate of the tree between the second and third census (cm y-1) 10. rgr1: [numeric] relative growth rate between the first and second census (cm y-1) 11. rgr2: [numeric] relative growth rate between the second and third census (cm y-1) 12. sa_gt: [numeric] species-adjusted modulus transformed growth rate 13. sa_rgr: [numeric] species-adjusted relative growth rate 14. gr_n: [numeric] number of individuals of the species used to calculate the mean and species modulus transformed growth rate 15. dbh_flight: [numeric] diameter at breast height (1.3 m) estimated at the time of the drone flight (cm) 16. eca: [numeric] exposed crown area calculated as the area of the crown polygon delineated in the orthomosaic (m2) 17. total_ca: [numeric] total crown area estimated from a crown area model (m2) 18. rcel: [numeric] relative crown exposure to light (m2) 19. time_flight_census3: [numeric] time in years from the drone flight date to the third census for that tree (yr)

54 ENVIRONMENTAL SCIENCES↗

GR-Athena++: Puncture Evolutions on Vertex-centered Oct-tree Adaptive Mesh Refinement

Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multimessenger astronomy call for highly performant and scalable codes on modern massively parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical spacetimes, GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong-scaling efficiencies above 95% for up to ~1.2 × 10 4 CPUs and excellent weak scaling is shown up to ~10 5 CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path toward numerical relativity at exascale.

79 ASTRONOMY AND ASTROPHYSICS↗

Reliable algorithm for modal decomposition

This paper describes a reliable, general algorithm for modal decomposition in real arithmetic and its use in analyzing and synthesizing control logic for linear dynamic systems. The numerical difficulties are described associated with computing the Jordan canonical form when the system has repeated, or nearly repeated, eigenvalues. A new algorithm is described that satisfactorily solves these numerical difficulties. The relation and extension to related numerical analysis research are discussed to clarify the reliability of the techniques. Finally, its implementation as a practical modal decomposition method for efficiently computing the matrix exponential, transfer functions, and frequency response is also described.

Walker, Robert A.↗

Numerical experiments related to the summer 1980 U.S. heat wave

The NASA Goddard Laboratory for Atmospheres fourth-order GCM was used in a series of medium-range numerical forecast experiments in order to improve understanding of the severe summer 1980 heat wave over the U.S. The results show that the derived soil moisture anomalies in the summer of 1980 contributed positively to the model's simulation of the heat wave maintenance, and suggest that once a region of reduced soil moisture is established, it tends to persist and maintain warmer and drier conditions. The lower soil moisture values resulted in reduced evaporation, higher ground temperatures, increased sensible heat flux from ground to air, higher surface air temperature, lower sea-level pressure, and higher 500-mb height. The effects of North Pacific sea-surface-temperature anomalies were mostly opposite to those of the soil-moisture anomalies: enhanced northerly flow of cooler drier air, increased evaporation, lower ground and air temperature, higher sea level pressure, and lower 500 mb heights over the Great Plains.

Wolfson, N.↗

Anti-Ultralocality and Plateau Models of Inflation

Anti-ultralocality refers to the growth of spatial gradient terms relative to velocity terms in the coupled Einstein--scalar field equations. It is a characteristic feature of decelerated expansion before the onset of inflation. Previous numerical relativity studies have shown that anti-ultralocality prevents the onset of inflation in models with power-law inflaton potentials. In this paper, we show that models with plateau-shaped inflaton potentials, which are considered to be the simplest way to generate a tensor-to-scalar ratio below current observational upper limits, are especially vulnerable to anti-ultralocality effects. The reasons are the flatness of the plateau and the energy density gap of $\sim 10$ orders of magnitude between the Planck density and the plateau potential energy. To study the problem, we develop a protocol for assessing the viability of inflationary models in general, and we apply it to a plateau potential using a previously validated numerical relativity code. We find that, starting from generic initial conditions, the growth of gradient terms in the Einstein equations relative to non-gradient terms either prevents inflation from lasting for enough $e$-folds or triggers a phase of quantum runaway. We show that the fine-tuning of initial conditions necessary to avoid these issues becomes more severe as the energy scale of inflation is made smaller, disfavoring common approaches for reducing the tensor-to-scalar ratio.

FOS: Physical sciences↗

Accurate Waveforms for Non-spinning Binary Black Holes using the Effective-one-body Approach

Using numerical relativity as guidance and the natural flexibility of the effective-one-body (EOB) model, we extend the latter so that it can successfully match the numerical relativity waveforms of non-spinning binary black holes during the last stages of inspiral, merger and ringdown. Here, by successfully, we mean with phase differences < or approx. 8% of a gravitational-wave cycle accumulated until the end of the ringdown phase. We obtain this result by simply adding a 4 post-Newtonian order correction in the EOB radial potential and determining the (constant) coefficient by imposing high-matching performances with numerical waveforms of mass ratios m1/m2 = 1,2/3,1/2 and = 1/4, m1 and m2 being the individual black-hole masses. The final black-hole mass and spin predicted by the numerical simulations are used to determine the ringdown frequency and decay time of three quasi-normal-mode damped sinusoids that are attached to the EOB inspiral-(plunge) waveform at the light-ring. The accurate EOB waveforms may be employed for coherent searches of gravitational waves emitted by non-spinning coalescing binary black holes with ground-based laser-interferometer detectors.

Buonanno, Alessandra↗

Observing the Final Moments of Massive Black Hole Mergers with LISA.

Mergers of binary black hole systems are one of the strongest sources of gravitational radiation expected to be observed by LISA. Recent advances in modeling the final merger and ringdown of comparable-mass systems, particularly via numerical relativity simulations, are dramatically expanding our understanding of these systems and the radiation they generate. We summarize recent modeling results, highlighting the work of Goddard s numerical relativity group, and apply this emerging knowledge to the problem of observing the final moments of binary black hole mergers with LISA.

Baker, John↗