Engineering PapersSearch

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

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

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.

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

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.

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

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

Ordering of implicants.

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

BOOLEAN ALGEBRA

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

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

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

Middle atmospheric electrodynamic modification by particle precipitation at the South Atlantic magnetic anomaly

Evidence for a localized middle atmospheric electrodynamic modification at low latitudes (southern Brazilian coast) of the South Atlantic Magnetic Anomaly (SAMA), in association with enhanced geomagnetic activity, are presented in a unified way combining recent observational efforts and related numerical studies. They involve a distortion effect in the fair weather electric field at balloon altitudes. This effect is attributed to a local intensification of energetic electron precipitation through a related middle atmospheric ionization enhancement and is elucidated by numeric simulation. From the electric field measurements and the numeric simulation, the intensification of precipitation is considered to occur in fairly narrow regions at the observed low L values (around L = 1.13) of the SAMA, with horizontal extensions of the order of a few hundred kilometers. A physical mechanism that could be responsible for this sort of intensification is suggested. Furthermore, a comparison of the phenomenon of middle atmospheric electrodynamic modification at the SAMA with a similar one at auroral latitudes, in response to enhanced solar and geomagnetic activity, is also given.

Gonzalez, W. D.

Binary Black Hole Late Inspiral: Simulations for Gravitational Wave Observations

Coalescing binary black hole mergers are expected to be the strongest gravitational wave sources for ground-based interferometers, such as the LIGO, VIRGO, and GEO600, as well as the spacebased interferometer LISA. Until recently it has been impossible to reliably derive the predictions of General Relativity for the final merger stage, which takes place in the strong-field regime. Recent progress in numerical relativity simulations is, however, revolutionizing our understanding of these systems. We examine here the specific case of merging equal-mass Schwarzschild black holes in detail, presenting new simulations in which the black holes start in the late inspiral stage on orbits with very low eccentricity and evolve for approximately 1200M through approximately 7 orbits before merging. We study the accuracy and consistency of our simulations and the resulting gravitational waveforms, which encompass approximately 14 cycles before merger, and highlight the importance of using frequency (rather than time) to set the physical reference when comparing models. Matching our results to PN calculations for the earlier parts of the inspiral provides a combined waveform with less than half a cycle of accumulated phase error through the entire coalescence. Using this waveform, we calculate signal-to-noise ratios (SNRs) for iLIGO, adLIGO, and LISA, highlighting the contributions from the late-inspiral and merger-ringdown parts of the waveform which can now be simulated numerically. Contour plots of SNR as a function of z and M show that adLIGO can achieve SNR 2 10 for some IMBBHs out to z approximately equals 1, and that LISA can see MBBHs in the range 3 x 10(exp 4) approximately < M/Mo approximately < 10(exp 7) at SNR > 100 out to the earliest epochs of structure formation at z > 15.

Baker, John G.

Black Hole Mergers and Gravitational Waves: Opening the New Frontier

The final merger of two black holes produces a powerful burst of gravitational waves, emitting more energy than all the stars in the observable universe combined. Since these mergers take place in the regime of strong dynamical gravity, computing the gravitational waveforms requires solving the full Einstein equations of general relativity on a computer. For more than 30 years, scientists tried to simulate these mergers using the methods of numerical relativity. The resulting computer codes were plagued by instabilities, causing them to crash well before the black holes in the binary could complete even a single orbit. In the past several years, this situation has changed dramatically, with a series of remarkable breakthroughs. This talk will highlight these breakthroughs and the resulting 'gold rush' of new results that is revealing the dynamics of binary black hole mergers, and their applications in gravitational wave detection, testing general relativity, and astrophysics.

Centrella, Joan

Solution of quadratic matrix equations for free vibration analysis of structures.

An efficient digital computer procedure and the related numerical algorithm are presented herein for the solution of quadratic matrix equations associated with free vibration analysis of structures. Such a procedure enables accurate and economical analysis of natural frequencies and associated modes of discretized structures. The numerically stable algorithm is based on the Sturm sequence method, which fully exploits the banded form of associated stiffness and mass matrices. The related computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be substantially more accurate and economical than other existing procedures of such analysis. Numerical examples are presented for two structures - a cantilever beam and a semicircular arch.

Gupta, K. K.

Gravitational Rocket from the Merging Massive Black Hole Binaries

Coalescing massive black hole binaries are expected to be among the most fascinating gravitational wave sources, observable by the NASA/ESA LISA detector. Not only will the merger events reveal the rich phenomenology of extremely strong and dynamical gravity deep inside the potential wells at the centers of galaxies (thus providing an excellent testing ground for general relativity), it will also make important contributions to the astrophysics of massive black hole evolutions. Typical black hole mergers involve asymmetric radiation of gravitational waves and lose linear momentum as well as energy and angular momentum. As a result, the merger remnant receives a kick from the GW emission: a gravitational rocket effect. High kick velocities (higher than the escape velocities of the host structure) would have a strong impact on our understanding of how massive black holes have evolved over time and, in particular, on the estimates of the merger rate for LISA. The main difficulties in calculations of the kick velocities has been in the last moments of the merger where the full theory of general relativity must be employed to accurately model the black hole dynamics. I describe a recent calculation of the kick velocities from numerical relativity simulations of the merging black hole binaries.

Choi, Dale

Modelling Gravitational Radiation from Binary Black Holes

The final merger and coalescence of binary black holes is a key source of strong gravitational waves for the LISA mission. Observing these systems will allow us to probe the formation of cosmic structure to high redshifts and test general relativity directly in the strong-field, dynamical regime. Recently, major breakthroughs have been made in modeling black hole mergers using numerical relativity. This talk will survey these exciting developments, focusing on the gravitational waveforms and the recoil kicks produced from non-equal mass mergers.

Centrella, Joan