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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

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↗
NASA NTRS2011

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

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↗
NASA NTRS2008

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

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↗
NASA NTRS1990

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

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

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↗
NASA NTRS2007

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↗
NASA NTRS2006

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