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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Numerical studies of 2-dimensional flows

A formulation of the lambda scheme for the analysis of two dimensional inviscid, compressible, unsteady transonic flows is presented. The scheme uses generalized Riemann variables to determine the appropriate two point, one sided finite difference approximation for each derivative in the unsteady Euler equations. These finite differences are applied at the predictor and corrector levels with shock updating at each level. The weaker oblique shocks are captured, but strong near normal shocks are fitted into the flow using the Rankine-Hugoniot relations. This code is demonstrated with a numerical example of a duct flow problem with developing normal and oblique shock waves. The technique is implemented in a code which has been made efficient by streamlining to a minimal number of operations and by eliminating branch statements. The scheme is shown to provide an accurate analysis of the flow, including formation, motions, and interactions of shocks; the results obtained on a relatively coarse mesh are comparable to those obtained by other methods on much finer meshes.

Moretti, G.↗
NASA NTRS1977

Transition and laminar instability

The linear stability theory was applied to the problem of boundary layer transition in incompressible flow. The theory was put into a form suitable for three-dimensional boundary layers; both the temporal and spatial theories were examined; and a generalized Gaster relation for three-dimensional boundary layers was derived. Numerical examples include the stability characteristics of Falkner-Skan boundary layers, the accuracy of the two-dimensional Gaster relation for these boundary layers, and the magnitude and direction of the group velocity for oblique waves in the Blasius boundary layer. Available experiments which bear on the validity of stability theory and its relation to transition are reviewed and the stability theory is applied to transition prediction. The amplitude method is described in which the wide band disturbance amplitude in the boundary layer is estimated from stability theory and an interaction relation for the initial amplitude density of the most unstable frequency.

Mack, L. M.↗
NASA NTRS1977

Fast but accurate techniques for calculating radiative terms in numerical atmospheric models and in remote sensing applications

Methods were developed for calculating radiative terms with relatively high accuracy but with sufficient speed, so that they can be used in numerical atmospheric models or in high volume processing of satellite measured radiances for remote sensing of atmospheric and surface parameters. Comparison with commonly used methods in both types of applications indicate improvements in calculating transmittances of factors between two and three, and in calculating radiances and cooling rates of factors between two and seven.

Arking, A.↗
NASA NTRS1990

An analysis of the vapor flow and the heat conduction through the liquid-wick and pipe wall in a heat pipe with single or multiple heat sources

A numerical analysis is presented for the overall performance of heat pipes with single or multiple heat sources. The analysis includes the heat conduction in the wall and liquid-wick regions as well as the compressibility effect of the vapor inside the heat pipe. The two-dimensional elliptic governing equations in conjunction with the thermodynamic equilibrium relation and appropriate boundary conditions are solved numerically. The solutions are in agreement with existing experimental data for the vapor and wall temperatures at both low and high operating temperatures.

Chen, Ming-Ming↗
NASA NTRS1992

The American Geophysical Union Chapman Conference on Tectonics and Topography

The Chapman Conference on Tectonics and Topography was held 31 Aug. - 4 Sep. 1992. The conference was designed to bring together disparate groups of earth scientists who increasingly found themselves working on similar problems but in relative isolation. Thus, process geomorphologists found themselves face-to-face with numerical modelers and field geomorphologists, hydrologists encountered geologists, and tectonophysicists found people with related data. The keynote speakers represented a wide variety of disciplines, all of which were relevant to the interdisciplinary theme of the conference. One of the most surprising issues that surfaced was the relative dearth of data that exists about erosion--process and rates. This was exacerbated by a reminder that erosion is critical to the evaluation of surface uplift.

Source record↗
NASA NTRS2012

Merging Black Holes

The final merger of two black holes is expected to be the strongest source of gravitational waves for both ground-based detectors such as LIGO and VIRGO, as well as future. space-based detectors. Since the merger takes place in the regime of strong dynamical gravity, computing the resulting gravitational waveforms requires solving the full Einstein equations of general relativity on a computer. For many years, numerical codes designed to simulate black hole mergers were plagued by a host of instabilities. However, recent breakthroughs have conquered these instabilities and opened up this field dramatically. This talk will focus on.the resulting 'gold rush' of new results that is revealing the dynamics and waveforms of binary black hole mergers, and their applications in gravitational wave detection, testing general relativity, and astrophysics

Centrella, Joan↗
NASA NTRS1976

A higher order panel method for general analysis and design applications in subsonic flow

A higher-order panel method is described for numerical solution of boundary-value problems relating to steady inviscid irrotational incompressible subsonic fluid flow in a domain. Both Neumann and Dirichlet boundary conditions are treated; two types of auxiliary conditions are used to remove the degrees of freedom that arise from specifying only the derivative of the perturbation velocity potential. Four general network types and two expansions of the induced potential kernel are employed in the numerical solution. Some results are presented which illustrate the modeling options and numerical characteristics of the method.

Johnson, F. T.↗
NASA NTRS2010

Merging Black Holes

The final merger of two black holes is expected to be the strongest source of gravitational waves for both ground-based detectors such as LIGO and VIRGO, as well as the space-based LISA. Since the merger takes place in the regime of strong dynamical gravity, computing the resulting gravitational waveforms requires solving the full Einstein equations of general relativity on a computer. For many years, numerical codes designed to simulate black hole mergers were plagued by a host of instabilities. However, recent breakthroughs have conquered these instabilities and opened up this field dramatically. This talk will focus on the resulting gold rush of new results that are revealing the dynamics and waveforms of binary black hole mergers, and their applications in gravitational wave detection, testing general relativity, and astrophysics.

Centrella, Joan↗
NASA NTRS2010

Merging Black Holes

The final merger of two black holes is expected to be the strongest source of gravitational waves for both ground-based detectors such as LIGO and VIRGO, as well as the space-based LISA. Since the merger takes place in the regime of strong dynamical gravity, computing the resulting gravitational waveforms requires solving the full Einstein equations of general relativity on a computer. For many years, numerical codes designed to simulate black hole mergers were plagued by a host of instabilities. However, recent breakthroughs have conquered these instabilities and opened up this field dramatically. This talk will focus on the resulting gold rush of new results that are revealing the dynamics and waveforms of binary black hole mergers, and their applications in gravitational wove detection, testing general relativity, and astrophysics.

Centrella, Joan↗
NASA NTRS1988

Hellas Basin ejecta: A prime stratigraphic marker and sampling target

The Hellas Basin and associated ejecta deposits dominate a large portion of the southern highlands of Mars. These deposits have been particularly well-documented at the scale 1:5,000,000 in the Sinus Sabaeus (MC-20), Iapygia (MC-21) and Mare Tyrrhenum (MC-22) quadrangles (1,2,3) where they have been mapped as ancient hilly and cratered material. In MC-22 the ejecta deposit is interpreted to thin from west to east, away from Hellas. In the eastern part of MC-22, the outlines of some large older craters, such as Hershel, can be seen through the thin ejecta cover and some of these buried craters even have enough relief to act as depositional traps for later mostly eolian sediments. Superposed on the ejecta is a population of numerous intermediate-sized to smaller craters which indicates the relatively old age of the ejecta deposit and causes the unit to have rather rough surface topography. Numerous small channels on steeply sloping crater rims and regionally sloping surface indicate erosion, transportation and redeposition of substantial amounts of ejecta materials, probably by aqueous fluvial processes. Numerous more recent volcanic constructs and flows also are superposed on the ejecta.

King, Elbert A.↗
DOE OSTI2022

Role of heterogeneous surface wettability on dynamic immiscible displacement, capillary pressure, and relative permeability in a CO 2 -water-rock system

Numerical examination of CO 2 displacement in brine filled rock pores with varying wettability. Results show that as the CO 2 contact angles on the rock surface increased (intermediate water wetting conditions) negative entry pressures for CO 2 displacement were observed, indicating that when heterogeneous surface wettability is accounted for, greater mobility of CO 2 into rock pores is expected. In addition, simulated relative permeability CO 2 curves with heterogeneous intermediate water wetting conditions exhibited a greater amount of variation, indicative of a complex interrelation of geometric capillary restrictions and surface wettability.

58 GEOSCIENCES↗
NASA NTRS2008

Efficient Simulation Budget Allocation for Selecting an Optimal Subset

We consider a class of the subset selection problem in ranking and selection. The objective is to identify the top m out of k designs based on simulated output. Traditional procedures are conservative and inefficient. Using the optimal computing budget allocation framework, we formulate the problem as that of maximizing the probability of correc tly selecting all of the top-m designs subject to a constraint on the total number of samples available. For an approximation of this corre ct selection probability, we derive an asymptotically optimal allocat ion and propose an easy-to-implement heuristic sequential allocation procedure. Numerical experiments indicate that the resulting allocatio ns are superior to other methods in the literature that we tested, and the relative efficiency increases for larger problems. In addition, preliminary numerical results indicate that the proposed new procedur e has the potential to enhance computational efficiency for simulation optimization.

Chen, Chun-Hung↗
NASA NTRS2014

Evaluating the Impact of Aerosols on Numerical Weather Predictions

This presentation focused on proposed Aerosol related case studies, and details were given about participating centers and modeling systems, as well as webpage and analysis tools. Preliminary results were given with regards to a case study on the dust storm over Egypt and another case study on the extreme pollution event in China.

Aerosol↗
NASA NTRS1993

A detailed model for simulation of catchment scale subsurface hydrologic processes

A catchment scale numerical model is developed based on the three-dimensional transient Richards equation describing fluid flow in variably saturated porous media. The model is designed to take advantage of digital elevation data bases and of information extracted from these data bases by topographic analysis. The practical application of the model is demonstrated in simulations of a small subcatchment of the Konza Prairie reserve near Manhattan, Kansas. In a preliminary investigation of computational issues related to model resolution, we obtain satisfactory numerical results using large aspect ratios, suggesting that horizontal grid dimensions may not be unreasonably constrained by the typically much smaller vertical length scale of a catchment and by vertical discretization requirements. Additional tests are needed to examine the effects of numerical constraints and parameter heterogeneity in determining acceptable grid aspect ratios. In other simulations we attempt to match the observed streamflow response of the catchment, and we point out the small contribution of the streamflow component to the overall water balance of the catchment.

Paniconi, Claudio↗
DOE OSTI2023

Temporal Subsampling Diminishes Small Spatial Scales in Recurrent Neural Network Emulators of Geophysical Turbulence

The immense computational cost of traditional numerical weather and climate models has sparked the development of machine learning (ML) based emulators. Because ML methods benefit from long records of training data, it is common to use data sets that are temporally subsampled relative to the time steps required for the numerical integration of differential equations. Here, we investigate how this often overlooked processing step affects the quality of an emulator's predictions. We implement two ML architectures from a class of methods called reservoir computing: (a) a form of Nonlinear Vector Autoregression (NVAR), and (b) an Echo State Network (ESN). Despite their simplicity, it is well documented that these architectures excel at predicting low dimensional chaotic dynamics. We are therefore motivated to test these architectures in an idealized setting of predicting high dimensional geophysical turbulence as represented by Surface Quasi-Geostrophic dynamics. In all cases, subsampling the training data consistently leads to an increased bias at small spatial scales that resembles numerical diffusion. Interestingly, the NVAR architecture becomes unstable when the temporal resolution is increased, indicating that the polynomial based interactions are insufficient at capturing the detailed nonlinearities of the turbulent flow. The ESN architecture is found to be more robust, suggesting a benefit to the more expensive but more general structure. Spectral errors are reduced by including a penalty on the kinetic energy density spectrum during training, although the subsampling related errors persist. Future work is warranted to understand how the temporal resolution of training data affects other ML architectures.

58 GEOSCIENCES↗
NASA NTRS1990

Analysis of spectral projectors in one-dimensional domains

A class of projection operators with values in a subspace of polynomials is analyzed. These projection operators are related to the Hilbert spaces involved in the numerical analysis of spectral methods. They are, in the first part of the paper, the standard Sobolev spaces and, in the second part, some weighted Sobolev spaces, the weight of which is related to the orthoronality relation satisfied by the Chebyshev polynomials. These results are used to study the approximation of a model fourth-order problem.

Maday, Y.↗