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

Gravitational Waves: Elusive Cosmic Messengers

The final merger of two black holes is expected to be the strongest g ravitational wave source for ground-based interferometers such as LIG O, VIRGO, and GE0600, as well as the space-based interferometer LISA. Observing these sources with gravitational wave detectors requires t hat we know the radiation waveforms they emit. Since these mergers ta ke place in regions of extreme gravity, we need to solve Einstein's equations of general relativity on a computer in order to calculate t hese waveforms. For more than 30 years, scientists have tried to comp ute 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 comple te even a single orbit. Within the past few years, however, this situ ation has changed dramatically, with a series of remarkable breakthro ughs. This talk will focus on new simulations that are revealing the dynamics and waveforms of binary black hole mergers, and their applic ations in gravitational wave detection, data analysis, and astrophysi cs.

Centrella, Joan↗
NASA NTRS2000

High Order Discontinuous Gelerkin Methods for Convection Dominated Problems with Application to Aeroacoustics

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. On the analysis side, we have studied the efficient and stable discontinuous Galerkin framework for small second derivative terms, for example in Navier-Stokes equations, and also for related equations such as the Hamilton-Jacobi equations. This is a truly local discontinuous formulation where derivatives are considered as new variables. On the applied side, we have implemented and tested the efficiency of different approaches numerically. Related issues in high order ENO and WENO finite difference methods and spectral methods have also been investigated. Jointly with Hu, we have presented a discontinuous Galerkin finite element method for solving the nonlinear Hamilton-Jacobi equations. This method is based on the RungeKutta discontinuous Galerkin finite element method for solving conservation laws. The method has the flexibility of treating complicated geometry by using arbitrary triangulation, can achieve high order accuracy with a local, compact stencil, and are suited for efficient parallel implementation. One and two dimensional numerical examples are given to illustrate the capability of the method. Jointly with Hu, we have constructed third and fourth order WENO schemes on two dimensional unstructured meshes (triangles) in the finite volume formulation. The third order schemes are based on a combination of linear polynomials with nonlinear weights, and the fourth order schemes are based on combination of quadratic polynomials with nonlinear weights. We have addressed several difficult issues associated with high order WENO schemes on unstructured mesh, including the choice of linear and nonlinear weights, what to do with negative weights, etc. Numerical examples are shown to demonstrate the accuracies and robustness of the methods for shock calculations. Jointly with P. Montarnal, we have used a recently developed energy relaxation theory by Coquel and Perthame and high order weighted essentially non-oscillatory (WENO) schemes to simulate the Euler equations of real gas. The main idea is an energy decomposition under the form epsilon = epsilon(sub 1) + epsilon(sub 2), where epsilon(sub 1) is associated with a simpler pressure law (gamma)-law in this paper) and the nonlinear deviation epsilon(sub 2) is convected with the flow. A relaxation process is performed for each time step to ensure that the original pressure law is satisfied. The necessary characteristic decomposition for the high order WENO schemes is performed on the characteristic fields based on the epsilon(sub l) gamma-law. The algorithm only calls for the original pressure law once per grid point per time step, without the need to compute its derivatives or any Riemann solvers. Both one and two dimensional numerical examples are shown to illustrate the effectiveness of this approach.

Shu, Chi-Wang↗
NASA NTRS2006

Gauge Conditions for Moving Black Holes Without Excision

Recent demonstrations of unexcised, puncture black holes traversing freely across computational grids represent a significant advance in numerical relativity. Stable an$ accurate simulations of multiple orbits, and their radiated waves, result. This capability is critically undergirded by a careful choice of gauge. Here we present analytic considerations which suggest certain gauge choices, and numerically demonstrate their efficacy in evolving a single moving puncture.

van Meter, James↗
NASA NTRS1988

Flowfield visualization for SSME hot gas manifold

The objective of this research, as defined by NASA-Marshall Space Flight Center, was two-fold: (1) to numerically simulate viscous subsonic flow in a proposed elliptical two-duct version of the fuel side Hot Gas Manifold (HGM) for the Space Shuttle Main Engine (SSME), and (2) to provide analytical support for SSME related numerical computational experiments, being performed by the Computational Fluid Dynamics staff in the Aerophysics Division of the Structures and Dynamics Laboratory at NASA-MSFC. Numerical results of HGM were calculations to complement both water flow visualization experiments and air flow visualization experiments and air experiments in two-duct geometries performed at NASA-MSFC and Rocketdyne. In addition, code modification and improvement efforts were to strengthen the CFD capabilities of NASA-MSFC for producing reliable predictions of flow environments within the SSME.

Roger, Robert P.↗
NASA NTRS2010

Post-Newtonian Initial Data with Waves: Progress in Evolution

"In Kelly \et al [Phys. Rev. D, 76:024008, (2007)], we presented new binary black-hole initial data adapted to puncture evolutions in numerical relativity. This data satisfies the constraint equations to 2.5 post-Newtonian order, and contains a transverse-traceless "wavy" metric contribution, violating the standard assumption of conformal flatness. We report on progress in evolving this data with a modem moving-puncture implementation of the BSSN equations in several numerical codes. We discuss the effect of the new metric terms on junk radiation and continuity of physical radiation extracted."

Kelly, Bernard J.↗
NASA NTRS1994

Numerical aperture limits on efficient ball lens coupling of laser diodes to single-mode fibers with defocus to balance spherical aberration

The potential capabilities and limitations of single ball lenses for coupling laser diode radiation to single-mode optical fibers have been analyzed; parameters important to optical communications were specifically considered. These parameters included coupling efficiency, effective numerical apertures, lens radius, lens refractive index, wavelength, magnification in imaging the laser diode on the fiber, and defocus to counterbalance spherical aberration of the lens. Limiting numerical apertures in object and image space were determined under the constraint that the lens perform to the Rayleigh criterion of 0.25-wavelength (Strehl ratio = 0.80). The spherical aberration-defocus balance to provide an optical path difference of 0.25 wavelength units was shown to define a constant coupling efficiency (i.e., 0.56). The relative numerical aperture capabilities of the ball lens were determined for a set of wavelengths and associated fiber-core diameters of particular interest for single-mode fiber-optic communication. The results support general continuing efforts in the optical fiber communications industry to improve coupling links within such systems with emphasis on manufacturing simplicity, system packaging flexibility, relaxation of assembly alignment tolerances, cost reduction of opto-electronic components and long term reliability and stability.

Wilson, R. Gale↗
NASA NTRS2009

On the Importance of Small Ice Crystals in Tropical Anvil Cirrus

In situ measurements of ice crystal concentrations and sizes made with aircraft instrumentation over the past two decades have often indicated the presence of numerous relatively small (< 50 m diameter) crystals in cirrus clouds. Further, these measurements frequently indicate that small crystals account for a large fraction of the extinction in cirrus clouds. The fact that the instruments used to make these measurements, such as the Forward Scattering Spectrometer Probe (FSSP) and the Cloud Aerosol Spectrometer (CAS), ingest ice crystals into the sample volume through inlets has led to suspicion that the indications of numerous small ]crystals could be artifacts of large ]crystal shattering on the instrument inlets. We present new aircraft measurements in anvil cirrus sampled during the Tropical Composition, Cloud, and Climate Coupling (TC4) campaign with the 2 ] Dimensional Stereo (2D ]S) probe, which detects particles as small as 10 m. The 2D ]S has detector "arms" instead of an inlet tube. Since the 2D ]S probe surfaces are much further from the sample volume than is the case for the instruments with inlets, it is expected that 2D ]S will be less susceptible to shattering artifacts. In addition, particle inter ]arrival times are used to identify and remove shattering artifacts that occur even with the 2D ]S probe. The number of shattering artifacts identified by the 2D ]S interarrival time analysis ranges from a negligible contribution to an order of magnitude or more enhancement in apparent ice concentration over the natural ice concentration, depending on the abundance of large crystals and the natural small ]crystal concentration. The 2D ]S measurements in tropical anvil cirrus suggest that natural small ]crystal concentrations are typically one to two orders of magnitude lower than those inferred from CAS. The strong correlation between the CAS/2D ]S ratio of small ]crystal concentrations and large ]crystal concentration suggests that the discrepancy is likely caused by shattering of large crystals on the CAS inlet. We argue that past measurements with CAS in cirrus with large crystals present may contain errors due to crystal shattering, and past conclusions derived from these measurements may need to be revisited. Further, we present correlations between CAS spurious concentration and 2D ]S large ]crystal mass from spatially uniform anvil cirrus sampling periods as an approximate guide for estimating quantitative impact of large ]crystal shattering on CAS concentrations in previous datasets. We use radiative transfer calculations to demonstrate that in the maritime anvil cirrus sampled during TC4, small crystals indicated by 2D ]S contribute relatively little cloud extinction, radiative forcing, or radiative heating in the anvils, regardless of anvil age or vertical location in the clouds. While 2D ]S ice concentrations in fresh anvil cirrus may often exceed 1 cm.3, and are observed to exceed 10 cm.3 in turrets, they are typically ~0.1 cm.3 and rarely exceed 1 cm.3 (<1.4% of the time) in aged anvil cirrus. We hypothesize that isolated occurrences of higher ice concentrations in aged anvil cirrus may be caused by ice nucleation driven by either small ]scale convection or gravity waves. It appears that the numerous small crystals detrained from convective updrafts do not persist in the anvil cirrus sampled during TC ]4.

Jensen, E. J.↗
NASA NTRS1994

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗
NASA NTRS2016

On the Susceptibility of Cold Tropical Cirrus to Ice Nuclei Abundance

Numerical simulations of cirrus formation in the tropical tropopause layer (TTL) are used to evaluate the impact heterogeneous ice nuclei (IN) abundance on cold cloud microphysical properties. The model includes homogeneous and heterogeneous ice nucleation, deposition growthsublimation, and sedimentation. Reanalysis temperature and wind fields are used to force the simulations, with addition of high-frequency wave temperature variability based on in situ measurements. The model results are constrained by recent in situ observations of TTL cirrus and relative humidity, as well as satellite measurements of cirrus bulk properties and occurrence frequencies. Temperature variability driven by high-frequency waves has a dominant influence on TTL cirrus microphysical properties and occurrence frequencies, and inclusion of these waves is required to produce agreement between the simulated and observed cloud properties as well as the abundance of TTL cirrus. With a composition-independent supersaturation threshold for homogeneous freezing of aqueous aerosols, the model produces excessive ice concentrations compared with in situ observations. Inclusion of relatively numerous heterogeneous nuclei ( 100L) in the simulations improves the agreement with observed ice concentrations. However, when IN contribute significantly to TTL cirrus ice nucleation, the occurrence frequency of large supersaturations with respect to ice is less than indicated by in situ measurements. The simulated TTL cirrus extinction statistics agree with observations (within uncertainties) from both in situ and remote-sensing measurements. We find that the sensitivity of TTL cirrus extinction and ice water content statistics to heterogeneous ice nuclei abundance is relatively weak. Likewise, the simulated occurrence frequencies of TTL cirrus are insensitive to ice nuclei abundance, both in terms of cloud frequency height distribution and regional distribution throughout the tropics.

occurences↗
NASA NTRS2010

Observing Mergers of Non-Spinning Black-Hole Binaries

Advances in the field of numerical relativity now make it possible to calculate the final, most powerful merger phase of binary black-hole coalescence for generic binaries. The state of the art has advanced well beyond the equal-mass case into the unequal-mass and spinning regions of parameter space. We present a study of the nonspinning portion of parameter space, primarily using an analytic waveform model tuned to available numerical data, with an emphasis on observational implications. We investigate the impact of varied m8BS ratio on merger signal-to-noise ratios (SNR~) for several detectors, and compare our results with expectations from the test-mass limit. We note a striking similarity of the waveform phasing of the merger waveform across the available mass ratios. Motivated by this, we calculate the match between our equal-mass and 4:1 mass-ratio waveforms during the merger as a function of location on the source sky, using a new formalism for the match that accounts for higher harmonics. This is an indicator of the amount of degeneracy in mass ratio for mergers of moderate mass ratio systems.

McWilliams, Sean T.↗
NASA NTRS2001

Modeling and Analysis of Wrinkled Membranes: An Overview

Thin-film membranes are basic elements of a variety of space inflatable/deployable structures. Wrinkling degrades the performance and reliability of these membrane structures, and hence has been a topic of continued interest. Wrinkling analysis of membranes for general geometry and arbitrary boundary conditions is quite challenging. The objective of this presentation is two-fold. Firstly, the existing models of wrinkled membranes and related numerical solution methods are reviewed. The important issues to be discussed are the capability of a membrane model to characterize taut, wrinkled and slack states of membranes in a consistent and physically reasonable manner; the ability of a wrinkling analysis method to predict the formation and growth of wrinkled regions, and to determine out-of-plane deformation and wrinkled waves; the convergence of a numerical solution method for wrinkling analysis; and the compatibility of a wrinkling analysis with general-purpose finite element codes. According to this review, several opening issues in modeling and analysis of wrinkled membranes that are to be addressed in future research are summarized, The second objective of this presentation is to discuss a newly developed membrane model of two viable parameters (2-VP model) and associated parametric finite element method (PFEM) for wrinkling analysis are introduced. The innovations and advantages of the proposed membrane model and PFEM-based wrinkling analysis are: (1) Via a unified stress-strain relation; the 2-VP model treat the taut, wrinkled, and slack states of membranes consistently; (2) The PFEM-based wrinkling analysis has guaranteed convergence; (3) The 2-VP model along with PFEM is capable of predicting membrane out-of-plane deformations; and (4) The PFEM can be integrated into any existing finite element code. Preliminary numerical examples are also included in this presentation to demonstrate the 2-VP model and PFEM-based wrinkling analysis approach.

Yang, B.↗
NASA NTRS1997

The Stability of Radiatively Cooling Jets I. Linear Analysis

The results of a spatial stability analysis of a two-dimensional slab jet, in which optically thin radiative cooling is dynamically important, are presented. We study both magnetized and unmagnetized jets at external Mach numbers of 5 and 20. We model the cooling rate by using two different cooling curves: one appropriate to interstellar gas, and the other to photoionized gas of reduced metallicity. Thus, our results will be applicable to both protostellar (Herbig-Haro) jets and optical jets from active galactic nuclei. We present analytical solutions to the dispersion relations in useful limits and solve the dispersion relations numerically over a broad range of perturbation frequencies. We find that the growth rates and wavelengths of the unstable Kelvin-Helmholtz (K-H) modes are significantly different from the adiabatic limit, and that the form of the cooling function strongly affects the results. In particular, if the cooling curve is a steep function of temperature in the neighborhood of the equilibrium state, then the growth of K-H modes is reduced relative to the adiabatic jet. On the other hand, if the cooling curve is a shallow function of temperature, then the growth of K-H modes can be enhanced relative to the adiabatic jet by the increase in cooling relative to heating in overdense regions. Inclusion of a dynamically important magnetic field does not strongly modify the important differences between an adiabatic jet and a cooling jet, provided the jet is highly supermagnetosonic and not magnetic pressure-dominated. In the latter case, the unstable modes behave more like the transmagnetosonic magnetic pressure-dominated adiabatic limit. We also plot fluid displacement surfaces associated with the various waves in a cooling jet in order to predict the structures that might arise in the nonlinear regime. This analysis predicts that low-frequency surface waves and the lowest order body modes will be the most effective at producing observable features in the jet.

Hardee, Philip E.↗
NASA NTRS2006

Binary Black Hole Mergers and Recoil Kicks

Recent developments in numerical relativity have made it possible to follow reliably the coalescence of two black holes from near the innermost stable circular orbit to final ringdown. This opens up a wide variety of exciting astrophysical applications of these simulations. Chief among these is the net kick received when two unequal mass or spinning black holes merge. The magnitude of this kick has bearing on the production and growth of supermassive black holes during the epoch of structure formation, and on the retention of black holes in stellar clusters. Here we report the first accurate numerical calculation of this kick, for two nonspinning black holes in a 1.5:1 mass ratio, which is expected based on analytic considerations to give a significant fraction of the maximum possible recoil. We have performed multiple runs with different initial separations, orbital angular momenta, resolutions, extraction radii, and gauges. The full range of our kick speeds is 86-116 kilometers per second, and the most reliable runs give kicks between 86 and 97 kilometers per second. This is intermediate between the estimates from two recent post-Newtonian analyses and suggests that at redshifts z greater than 10, halos with masses less than 10(exp 9) M(sub SUN) will have difficulty retaining coalesced black holes after major mergers.

Centrella, Joan↗
NASA NTRS2011

Comment on "An Analysis of VLF Electric Field Spectra Measured in Titan's Atmosphere by The Huygens Probe" By J. A. Morente et al.

Morente et al. have recently revisited the VLF electric field measurements made with the Permittivity, Wave and Altimetry (PWA) instrument during the descent of the Huygens Probe through the atmosphere of Titan. They assert that they have identified several harmonics of the transverse resonance mode of the surface?]ionosphere cavity, which would prove the existence of an electrical activity in the atmosphere of the largest satellite of Saturn. We refute this finding on the basis that it results from an artifact due to an improper analysis of the data set. [2] The investigators of the Permittivity, Wave and Altimetry (PWA) experiment on the Huygens Probe have reported the extremely low frequency (ELF) and very low frequency (VLF) electric signals recorded during the descent through the atmosphere of Titan. The PWA data are archived in the Planetary Science Archive (PSA) of ESA, and an extensive description of the instrument is at the disposal of the scientific community. Morente and his coworkers have revisited this data set and reported the results of their investigations in two papers. In a first paper, they claim that they have detected in the ELF range (0.100 Hz) several harmonics of a global resonance allegedly generated by lightning activity in the spherical cavity guide formed by the surface of Titan and the inner boundary of the ionosphere, a phenomenon similar to the Schumann resonance observed at EartH In the second paper dedicated to the VLF electric signal recorded by PWA, in the range 0.10 kHz, they argue that they can also bring out the transverse resonance and its harmonics, a more local phenomenon that develops around the excitation source and whose frequency is controlled by the separation between Titan?fs surface and the inner ionospheric boundary. [3] The PWA investigators have analyzed the narrowband ELF signal at about 36 Hz effectively observed during the entire descent. They have not endorsed, however, the alternative approach of Morente et al. that discloses additional ELF narrow-band signals and numerous related harmonics. Hamelin et al. and Beghin et al. have questioned the work of Morente and his coworkers and demonstrated unambiguously that their findings are mere artifacts entirely due to a faulty procedure. Morente et al. refuted this critical analysis of their ELF paper in a rebuttal. The present comment draws a parallelism between the ELF and VLF papers and demonstrates that the VLF signal carries no evidence whatever of any transverse resonance. The analyses of the ELF and VLF signals by Morente and his coworkers are flawed, and their conclusions are shown here to be invalid. A more comprehensive and thorough assessment of the numerical approaches proposed by Morente et al. is given by Berthelin et al..

Grard, Rejean↗
NASA NTRS1974

Overview of NASA aeronautical propulsion research and technology program

The program discussed is aimed at improving performance within an extended operating range, reducing weight, increasing service life, achieving greater cost effectiveness, reducing noise and exhaust pollution, and improving means of energy conservation. The program places emphasis on basic research in numerous related technical disciplines, system analysis, component technology, full scale propulsion system studies, and technology demonstrations. Much attention in the discussion is given to noise and pollution minimization research.

Johnson, H. W.↗
NASA NTRS1976

Prediction of the diffuse far infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far infrared studies of the galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far infrared flux distribution from the galactic plane as a function of galactic longitude alternate theta for 4 deg or = alternate theta or = 90 and the far infrared emissivity as a function of galactocentric distance is predicted. Future measurements of the galactic far infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the galaxy, particularly the inner region.

Fazio, G. G.↗
NASA NTRS1976

Prediction of the diffuse far-infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far-infrared studies of the Galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far-infrared flux distribution from the galactic plane is predicted as a function of galactic longitude for longitudes between 4 and 90 deg; the far-infrared emissivity is predicted as a function of galactocentric distance. Future measurements of the galactic far-infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the Galaxy, particularly the inner region.

Fazio, G. G.↗
NASA NTRS1976

Prediction of the diffuse far-infrared flux from the galactic plane

A basic model and simple numerical relations useful for future far-infrared studies of the Galaxy are presented. Making use of recent CO and other galactic surveys, the diffuse far-infrared flux distribution from the galactic planes is predicted as a function of galactic longitude for the region between 4 and 90 deg longitude; the far-infrared emissivity is predicted as a function of galactocentric distance. Future measurements of the galactic far-infrared flux would yield valuable information on the physical properties and distribution of dust and molecular clouds in the Galaxy, particularly the inner region.

Fazio, G. G.↗