SEARCH · Search NASA

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Relative Navigation of Formation Flying Satellites

The Guidance, Navigation, and Control Center (GNCC) at Goddard Space Flight Center (GSFC) has successfully developed high-accuracy autonomous satellite navigation systems using the National Aeronautics and Space Administration's (NASA's) space and ground communications systems and the Global Positioning System (GPS). In addition, an autonomous navigation system that uses celestial object sensor measurements is currently under development and has been successfully tested using real Sun and Earth horizon measurements.The GNCC has developed advanced spacecraft systems that provide autonomous navigation and control of formation flyers in near-Earth, high-Earth, and libration point orbits. To support this effort, the GNCC is assessing the relative navigation accuracy achievable for proposed formations using GPS, intersatellite crosslink, ground-to-satellite Doppler, and celestial object sensor measurements. This paper evaluates the performance of these relative navigation approaches for three proposed missions with two or more vehicles maintaining relatively tight formations. High-fidelity simulations were performed to quantify the absolute and relative navigation accuracy as a function of navigation algorithm and measurement type. Realistically-simulated measurements were processed using the extended Kalman filter implemented in the GPS Enhanced Inboard Navigation System (GEONS) flight software developed by GSFC GNCC. Solutions obtained by simultaneously estimating all satellites in the formation were compared with the results obtained using a simpler approach based on differencing independently estimated state vectors.

Long, Anne↗

NASA NTRS2004

Linear State-Space Representation of the Dynamics of Relative Motion, Based on Restricted Three Body Dynamics

Precision Formation Flying is an enabling technology for a variety of proposed space-based observatories, including the Micro-Arcsecond X-ray Imaging Mission (MAXIM) , the associated MAXIM pathfinder mission, Stellar Imager (SI) and the Terrestrial Planet Finder (TPF). An essential element of the technology is the control algorithm, requiring a clear understanding of the dynamics of relative motion. This paper examines the dynamics of relative motion in the context of the Restricted Three Body Problem (RTBP). The natural dynamics of relative motion are presented in their full nonlinear form. Motivated by the desire to apply linear control methods, the dynamics equations are linearized and presented in state-space form. The stability properties are explored for regions in proximity to each of the libration points in the Earth/Moon - Sun rotating frame. The dynamics of relative motion are presented in both the inertial and rotating coordinate frames.

Luquette,Richard J.↗

NASA NTRS2000

Lidar Measurements of Relative Humidity and Ice Supersaturation in the Upper Troposphere

We compute upper tropospheric relative humidity profiles using water vapor profiles measured by an airborne DIAL and a ground-based Raman lidar. LASE water vapor and MTP temperature profiles acquired from the NASA DC-8 aircraft during the recent Pacific Exploratory Mission Tropics B (PEM Tropics B) field mission in the tropical Pacific and the SAGE-III Ozone Loss and Validation Experiment (SOLVE) in the Arctic as well as water vapor profiles derived from the ground-based DOE ARM Southern Great Plains (SGP) CART Raman lidar are used. Comparisons of the lidar water vapor measurements with available in situ measurements show reasonable agreement for water vapor mixing ratios above 0.05 g/kg. Relative humidity frequency distributions computed using LASE data indicate that ice supersaturation occurred about 5-11% of the time when temperatures were below -35 C. While a higher frequency of ice supersaturation was observed during SOLVE, higher peak values of relative humidity were observed during PEM Tropics B. The relative humidity fields associated with cirrus clouds are also examined.

Ferrare, Richard A.↗

NASA NTRS1997

Relative rates for plasma homo- and copolymerizations of olefins in a homologous series of fluorinated ethylenes

It is well known that the rate of plasma polymerization, or deposition rate, of a given monomer depends on various plasma process parameters, e.g., monomer flow rate, pressure, power, frequency (DC, rf or microwave), location of the substrate in the reactor, reactor geometry or configuration, and temperature. In contrast, little work has been done to relate deposition rates to monomer structures for a homologous series of monomers where the rates are obtained under identical plasma process parameters. For the particular series of fluorinated ethylenes (C2HxF4-x; x = 0-4), deposition rates were reported for ethylene (ET), vinyl fluoride, vinylidene fluoride and tetrafluoroethylene (TFE), but for plasma polymerizations carried out under different discharge conditions, e.g., pressure, current density, and electrode temperature. Apparently, relative deposition rates were reported for only two members of that series (ET, x = 4, and TFE, x = 0) for which the plasma polymerizations were conducted under identical conditions. We now present relative deposition rates for both homopolymerizations and copolymerizations of the entire series of fluorinated ethylenes (x = 0-4). Our interest in such rates stems from prior work on the plasma copolymerization of ET and TFE in which it was found that the deposition rates for the plasma copolymers, when plotted versus mol % TFE in the ET/TFE feed stock, followed a concave-downward curve situated above the straight line joining the deposition rates for the plasma homopolymers. This type of plot (observed also for an argon-ET/TFE plasma copolymerization) indicated a positive interaction between ET and TFE such that each monomer apparently "sensitized" the plasma copolymerization of the other. Since the shape of that plot is not altered if mol % TFE is replaced by F/C, the fluorine-to-carbon ratio, this paper aims (1) to show how the relative deposition rates for plasma copolymers drawn from all pairs of monomers in the C2HxF4-x series, as well as the deposition rates for the individual plasma homopolymers, vary with F/C ratios of the monomers or monomer blends, and (2) to see if those rates give rise to a common plot.

NASA Center ARC↗

NASA NTRS1999

Circadian regulation of human sleep and age-related changes in its timing, consolidation and EEG characteristics

The light-entrainable circadian pacemaker located in the suprachiasmatic nucleus of the hypothalamus regulates the timing and consolidation of sleep by generating a paradoxical rhythm of sleep propensity; the circadian drive for wakefulness peaks at the end of the day spent awake, ie close to the onset of melatonin secretion at 21.00-22.00 h and the circadian drive for sleep crests shortly before habitual waking-up time. With advancing age, ie after early adulthood, sleep consolidation declines, and time of awakening and the rhythms of body temperature, plasma melatonin and cortisol shift to an earlier clock hour. The variability of the phase relationship between the sleep-wake cycle and circadian rhythms increases, and in old age sleep is more susceptible to internal arousing stimuli associated with circadian misalignment. The propensity to awaken from sleep advances relative to the body temperature nadir in older people, a change that is opposite to the phase delay of awakening relative to internal circadian rhythms associated with morningness in young people. Age-related changes do not appear to be associated with a shortening of the circadian period or a reduction of the circadian drive for wake maintenance. These changes may be related to changes in the sleep process itself, such as reductions in slow-wave sleep and sleep spindles as well as a reduced strength of the circadian signal promoting sleep in the early morning hours. Putative mediators and modulators of circadian sleep regulation are discussed.

Review, Tutorial↗

NASA NTRS2004

Linear State-Space Representation of the Dynamics of Relative Motion, Based on Restricted Three Body Dynamics

Precision Formation Flying is an enabling technology for a variety of proposed space- based observatories, including the Micro-Arcsecond X-ray Imaging Mission (MAXIM), Stellar Imager (SI) and the Terrestrial Planet Finder (TPF). An essential element of the technology is the control algorithm, requiring a clear understanding of the dynamics of relative motion. This paper examines the dynamics of relative motion in the context of the Restricted Three Body Problem (RTBP). The natural dynamics of relative motion are presented in their full nonlinear form. Motivated by the desire to apply linear control methods, the dynamics equations are linearized and presented in state-space form. The stability properties are explored for regions in proximity to each of the libration points in the Earth/Moon - Sun rotating frame. The dynamics of relative motion are presented in both the inertial and rotating coordinate frames.

Luquette, Richard J.↗

NASA NTRS1994

Gender- and age-related differences in heart rate dynamics: are women more complex than men?

OBJECTIVES. This study aimed to quantify the complex dynamics of beat-to-beat sinus rhythm heart rate fluctuations and to determine their differences as a function of gender and age. BACKGROUND. Recently, measures of heart rate variability and the nonlinear "complexity" of heart rate dynamics have been used as indicators of cardiovascular health. Because women have lower cardiovascular risk and greater longevity than men, we postulated that there are important gender-related differences in beat-to-beat heart rate dynamics. METHODS. We analyzed heart rate dynamics during 8-min segments of continuous electrocardiographic recording in healthy young (20 to 39 years old), middle-aged (40 to 64 years old) and elderly (65 to 90 years old) men (n = 40) and women (n = 27) while they performed spontaneous and metronomic (15 breaths/min) breathing. Relatively high (0.15 to 0.40 Hz) and low (0.01 to 0.15 Hz) frequency components of heart rate variability were computed using spectral analysis. The overall "complexity" of each heart rate time series was quantified by its approximate entropy, a measure of regularity derived from nonlinear dynamics ("chaos" theory). RESULTS. Mean heart rate did not differ between the age groups or genders. High frequency heart rate power and the high/low frequency power ratio decreased with age in both men and women (p < 0.05). The high/low frequency power ratio during spontaneous and metronomic breathing was greater in women than men (p < 0.05). Heart rate approximate entropy decreased with age and was higher in women than men (p < 0.05). CONCLUSIONS. High frequency heart rate spectral power (associated with parasympathetic activity) and the overall complexity of heart rate dynamics are higher in women than men. These complementary findings indicate the need to account for gender-as well as age-related differences in heart rate dynamics. Whether these gender differences are related to lower cardiovascular disease risk and greater longevity in women requires further study.

NASA Discipline Cardiopulmonary↗

NASA NTRS2004

CDGPS-Based Relative Navigation for Multiple Spacecraft

This thesis investigates the use of Carrier-phase Differential GPS (CDGPS) in relative navigation filters for formation flying spacecraft. This work analyzes the relationship between the Extended Kalman Filter (EKF) design parameters and the resulting estimation accuracies, and in particular, the effect of the process and measurement noises on the semimajor axis error. This analysis clearly demonstrates that CDGPS-based relative navigation Kalman filters yield good estimation performance without satisfying the strong correlation property that previous work had associated with "good" navigation filters. Several examples are presented to show that the Kalman filter can be forced to create solutions with stronger correlations, but these always result in larger semimajor axis errors. These linear and nonlinear simulations also demonstrated the crucial role of the process noise in determining the semimajor axis knowledge. More sophisticated nonlinear models were included to reduce the propagation error in the estimator, but for long time steps and large separations, the EKF, which only uses a linearized covariance propagation, yielded very poor performance. In contrast, the CDGPS-based Unscented Kalman relative navigation Filter (UKF) handled the dynamic and measurement nonlinearities much better and yielded far superior performance than the EKF. The UKF produced good estimates for scenarios with long baselines and time steps for which the EKF would diverge rapidly. A hardware-in-the-loop testbed that is compatible with the Spirent Simulator at NASA GSFC was developed to provide a very flexible and robust capability for demonstrating CDGPS technologies in closed-loop. This extended previous work to implement the decentralized relative navigation algorithms in real time.

Mitchell, Megan Leigh↗

NASA NTRS2005

Controlling Reflections from Mesh Refinement Interfaces in Numerical Relativity

A leading approach to improving the accuracy on numerical relativity simulations of black hole systems is through fixed or adaptive mesh refinement techniques. We describe a generic numerical error which manifests as slowly converging, artificial reflections from refinement boundaries in a broad class of mesh-refinement implementations, potentially limiting the effectiveness of mesh- refinement techniques for some numerical relativity applications. We elucidate this numerical effect by presenting a model problem which exhibits the phenomenon, but which is simple enough that its numerical error can be understood analytically. Our analysis shows that the effect is caused by variations in finite differencing error generated across low and high resolution regions, and that its slow convergence is caused by the presence of dramatic speed differences among propagation modes typical of 3+1 relativity. Lastly, we resolve the problem, presenting a class of finite-differencing stencil modifications which eliminate this pathology in both our model problem and in numerical relativity examples.

Baker, John G.↗

NASA NTRS2005

Coping with Trial-to-Trial Variability of Event Related Signals: A Bayesian Inference Approach

In electro-neurophysiology, single-trial brain responses to a sensory stimulus or a motor act are commonly assumed to result from the linear superposition of a stereotypic event-related signal (e.g. the event-related potential or ERP) that is invariant across trials and some ongoing brain activity often referred to as noise. To extract the signal, one performs an ensemble average of the brain responses over many identical trials to attenuate the noise. To date, h s simple signal-plus-noise (SPN) model has been the dominant approach in cognitive neuroscience. Mounting empirical evidence has shown that the assumptions underlying this model may be overly simplistic. More realistic models have been proposed that account for the trial-to-trial variability of the event-related signal as well as the possibility of multiple differentially varying components within a given ERP waveform. The variable-signal-plus-noise (VSPN) model, which has been demonstrated to provide the foundation for separation and characterization of multiple differentially varying components, has the potential to provide a rich source of information for questions related to neural functions that complement the SPN model. Thus, being able to estimate the amplitude and latency of each ERP component on a trial-by-trial basis provides a critical link between the perceived benefits of the VSPN model and its many concrete applications. In this paper we describe a Bayesian approach to deal with this issue and the resulting strategy is referred to as the differentially Variable Component Analysis (dVCA). We compare the performance of dVCA on simulated data with Independent Component Analysis (ICA) and analyze neurobiological recordings from monkeys performing cognitive tasks.

Ding, Mingzhou↗

NASA NTRS2006

Are Ring Current Ions Lost in Electromagnetic Ion Cyclotron Wave Dispersion Relation?

Electromagnetic ion cyclotron (EMIC) waves are widely observed in the inner and outer magnetosphere, at geostationary orbit, at high latitudes along the plasmapause, and at the ionospheric altitudes. Interaction of the Ring Current (RC) ions and EMIC waves causes ion scattering into the loss cone and leads to decay of the RC, especially during the main phase of storms when the RC decay times of about one hour or less are observed. The oblique EMIC waves damp due to Landau resonance with the thermal plasmaspheric electrons, and subsequent transport of the dissipating wave energy into the ionosphere below causes an ionosphere temperature enhancement. Induced scattering of these waves by the plasmaspheric thermal ions leads to ion temperature enhancement, and forms a so-called hot zone near the plasmapause where the temperature of core plasma ions can reach tens of thousands of degrees. Relativistic electrons in the outer radiation belt also interact well with the EMIC waves, and during the main and/or recovery phases of the storms these electrons can easily be scattered into the loss cone over a time scale from several hours to a day. The plasma density distribution in the magnetosphere and the ion content play a critical role in EMIC wave generation and propagation, but the wave dispersion relation in the known RC-EMIC wave interaction models is assumed to be determined by the thermal plasma distribution only. In these models, the modification of the EMIC wave dispersion relation caused by the RC ions is not taken into account, and the RC ions are only treated as a source of free energy in order to generate EMIC waves. At the same time, the RC ions can dominate the thermal magnetospheric content in the night MLT sector at great L shells during the main and/or recovery storm phase. In this study, using our self-consistent RC-EMIC wave model [Khazanov et al., 2006], we simulate the May 1998 storm in order to quantify the global EMIC wave redistribution caused by taking into account the RC ions in the EMIC wave dispersion relation. The dramatic wave pattern redistribution is observed in the postdusk-predawn MLT sector (night sector) for L greater than 5. We found the intense EMIC waves (about a few nT) there during the main and early recovery phases of the storm. The observed wave generation in this sector is caused by taking into account the EMIC wave dispersion change due to the RC ions. There are no waves at these locations in our model if the RC ions are taken into account in the wave growth rate only, and the wave dispersion relation is only governed by the thermal plasmaspheric model.

Khazanov, G. V.↗

NASA NTRS2006

Effect of Ring Current Ions on Electromagnetic Ion Cyclotron Wave Dispersion Relation

Gamayunov, K. V.↗

NASA NTRS2007

Single-Frequency GPS Relative Navigation in a High Ionosphere Orbital Environment

The Global Positioning System (GPS) provides a convenient source for space vehicle relative navigation measurements, especially for low Earth orbit formation flying and autonomous rendezvous mission concepts. For single-frequency GPS receivers, ionospheric path delay can be a significant error source if not properly mitigated. In particular, ionospheric effects are known to cause significant radial position error bias and add dramatically to relative state estimation error if the onboard navigation software does not force the use of measurements from common or shared GPS space vehicles. Results from GPS navigation simulations are presented for a pair of space vehicles flying in formation and using GPS pseudorange measurements to perform absolute and relative orbit determination. With careful measurement selection techniques relative state estimation accuracy to less than 20 cm with standard GPS pseudorange processing and less than 10 cm with single-differenced pseudorange processing is shown.

Conrad, Patrick R.↗

NASA NTRS2007

Impact of Ring Current Ions on Electromagnetic Ion Cyclotron Wave Dispersion Relation

Effect of the ring current ions in the real part of electromagnetic ion Cyclotron wave dispersion relation is studied on global scale. Recent Cluster observations by Engebretson et al. showed that although the temperature anisotropy of is energetic (> 10 keV) ring current protons was high during the entire 22 November 2003 perigee pass, electromagnetic ion cyclotron waves were observed only in conjunction with intensification of the ion fluxes below 1 keV by over an order of magnitude. To study the effect of the ring current ions on the wave dispersive properties and the corresponding global wave redistribution, we use a self-consistent model of interacting ring current and electromagnetic ion cyclotron waves, and simulate the May 1998 storm. The main findings of our simulation can be summarized as follows: First, the plasma density enhancement in the night MLT sector during the main and recovery storm phases is mostly caused by injection of suprathermal plasma sheet H + (approximately < 1 keV), which dominate the thermal plasma density. Second, during the recovery storm phases, the ring current modification of the wave dispersion relation leads to a qualitative change of the wave patterns in the postmidnight-dawn sector for L > 4.75. This "new" wave activity is well organized by outward edges of dense suprathermal ring current spots, and the waves are not observed if the ring current ions are not included in the real part of dispersion relation. Third, the most intense wave-induced ring current precipitation is located in the night MLT sector and caused by modification of the wave dispersion relation. The strongest precipitating fluxes of about 8 X 10(exp 6)/ (cm(exp 2) - s X st) are found near L=5.75, MLT=2 during the early recovery phase on 4 May. Finally, the nightside precipitation is more intense than the dayside fluxes, even if there are less intense waves, because the convection field moves ring current ions into the loss cone on the nightside, but drives them out of the loss cone on the dayside. So convection and wave scattering reinforce each other in the nightside, but interfere in the dayside sector.

Khazanov, G. V.↗

NASA NTRS2006

The Solar Wind and Geomagnetic Activity as a Function of Time Relative to Corotating Interaction Regions

Corotating interaction regions during the declining phase of the solar cycle are the cause of recurrent geomagnetic storms and are responsible for the generation of high fluxes of relativistic electrons. These regions are produced by the collision of a high-speed stream of solar wind with a slow-speed stream. The interface between the two streams is easily identified with plasma and field data from a solar wind monitor upstream of the Earth. The properties of the solar wind and interplanetary magnetic field are systematic functions of time relative to the stream interface. Consequently the coupling of the solar wind to the Earth's magnetosphere produces a predictable sequence of events. Because the streams persist for many solar rotations it should be possible to use terrestrial observations of past magnetic activity to predict future activity. Also the high-speed streams are produced by large unipolar magnetic regions on the Sun so that empirical models can be used to predict the velocity profile of a stream expected at the Earth. In either case knowledge of the statistical properties of the solar wind and geomagnetic activity as a function of time relative to a stream interface provides the basis for medium term forecasting of geomagnetic activity. In this report we use lists of stream interfaces identified in solar wind data during the years 1995 and 2004 to develop probability distribution functions for a variety of different variables as a function of time relative to the interface. The results are presented as temporal profiles of the quartiles of the cumulative probability distributions of these variables. We demonstrate that the storms produced by these interaction regions are generally very weak. Despite this the fluxes of relativistic electrons produced during those storms are the highest seen in the solar cycle. We attribute this to the specific sequence of events produced by the organization of the solar wind relative to the stream interfaces. We also show that there are large quantitative differences in various parameters between the two cycles.

McPherron, Robert L.↗

NASA NTRS2008

Relating a Spectral Index from MODIS and Tower-based Measurements to Ecosystem Light Use Efficiency for a Fluxnet-Canada Coniferous Forest

As part of the North American Carbon Program effort to quantify the terrestrial carbon budget of North America, we have been examining the possibility of retrieving ecosystem light use efficiency (LUE, the carbon sequestered per unit photosynthetically active radiation) directly from satellite observations. Our novel approach has been to compare LUE derived from tower fluxes with LUE estimated using spectral indices computed from MODIS satellite observations over forests in the Fluxnet-Canada Research Network, using the MODIS narrow ocean bands acquired over land. We matched carbon flux data collected around the time of the MODIS mid-day overpass for over one hundred relatively clear days in five years (2001-2006) from a mature Douglas fir forest in British Columbia. We also examined hyperspectral reflectance data collected diurnally from the tower in conjunction with the eddy correlation fluxes and meteorological measurements made throughout the 2006 growing season at this site. The tower-based flux data provided an opportunity to examine diurnal and seasonal LUE processes and their relationship to spectral indices at the scale of the forest stand. We evaluated LUE in conjunction with the Photochemical Reflectance Index (PRI), a normalized difference spectral index that uses 531 nm and a reference band to capture responses to high light induced stress afforded by the xanthophyll cycle. Canopy structure information, retrieved from airborne laser scanning radar (LiDAR) observations, was used to partition the forest canopy into sunlit and shaded fractions throughout the day, on numerous days during 2006. At each observation period throughout a day, the PRI was examined for the sunlit, shaded, and intermediate canopy segments defined by their instantaneous position relative to the solar principal plane (SPP). The sunlit sector was associated with the illumination "hotspot" (the reflectance backscatter maximum), the shaded sector with the "cold or dark spot" (the reflectance forward scatter minimum), while the intermediate, mixed sunlit/shade sector was located in the cross-plane to the SPP. The PRI indices clearly captured the differences in leaf groups, with sunlit foliage exhibiting the lowest values on sunny days throughout the 2006 season. When tower-based canopy-level LUE was recalculated to estimate foliage-based values (LUE(sub foilage) for the three foliage groups under their incident light environments, a strong linear relationship for PRI:LUE(sub foilage) was demonstrated (0.6 less than or equal to r(sup 2) less than or equal to 0.8, n=822, P<0.0001). The MODIS data represent relatively large areas when acquired at nadir (approx.1 sq km) or at variable off-nadir view angles (greater than or equal to 1 sq km) looking forward or aft. Nevertheless, a similar relationship between MODIS PRI and tower-based LUE was obtained from satellite observations (r(sup 2) = 0.76, n=105, P= 0.026) when the azimuth offsets from the SPP for off-nadir observations were considered. At this relatively high latitude of 50 degrees, the MODIS directional observations were offset from the SPP by approximately 50 degrees, but still represented backscatter or forward scatter sectors of the bidirectional reflectance distribution function (BRDF). The backscatter observations sampled the sunlit forest and provided lower PRI values, in general, than the forward scatter observations from the shaded forest. Since the hotspot and darkspot were not typically directly observed, the dynamic range for MODIS PRI was less than that observed in the SPP at the canopy level; therefore, MODIS PRI values were more similar to those observed in sifu in the BRDF cross-plane. While not ideal in terms of spatial resolution or optimal viewing configuration, the MODIS observations nevertheless provide a means to monitor forest under stress using narrow spectral band indices and off-nadir observations. This research has stimulated several spin-off studies for remote sensinf LUE, and demonstrates the importance of the connection between ecosystem structure and physiological function.

Middleton, Elizabeth M.↗

NASA NTRS2007

Binary Black Holes, Gravitational Waves, and Numerical Relativity

This viewgraph presentation reviews the massive black hole (MBH) binaries that are found at the center of most galaxies, "astronomical messenger", gravitational waves (GW), and the use of numerical relativity understand the features of these phenomena. The final merger of two black holes releases a tremendous amount of energy and is one of the brightest sources in the gravitational wave sky. 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.. This talk will take you on this quest for the holy grail of numerical relativity, 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 LIGO and LISA.

Centrella, Joan↗

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