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 NTRS1990

Vibration and wave propagation characteristics of multisegmented elastic beams

Closed form analytical solutions are derived for the vibration and wave propagation of multisegmented elastic beams. Each segment is modeled as a Timoshenko beam with possible inclusion of material viscosity, elastic foundation and axial forces. Solutions are obtained by using transfer matrix methods. According to these methods formal solutions are first constructed which relate the deflection, slope, moment and shear force of one end of the individual segment to those of the other. By satisfying appropriate continuity conditions at segment junctions, a global 4x4 matrix results which relates the deflection, slope, moment and shear force of one end of the beam to those of the other. If any boundary conditions are subsequently invoked on the ends of the beam one gets the appropriate characteristic equation for the natural frequencies. Furthermore, by invoking appropriate periodicity conditions the dispersion relation for a periodic system is obtained. A variety of numerical examples are included.

Nayfeh, Adnan H.↗
NASA NTRS1991

The first ANDES elements: 9-DOF plate bending triangles

New elements are derived to validate and assess the assumed natural deviatoric strain (ANDES) formulation. This is a brand new variant of the assumed natural strain (ANS) formulation of finite elements, which has recently attracted attention as an effective method for constructing high-performance elements for linear and nonlinear analysis. The ANDES formulation is based on an extended parametrized variational principle developed in recent publications. The key concept is that only the deviatoric part of the strains is assumed over the element whereas the mean strain part is discarded in favor of a constant stress assumption. Unlike conventional ANS elements, ANDES elements satisfy the individual element test (a stringent form of the patch test) a priori while retaining the favorable distortion-insensitivity properties of ANS elements. The first application of this formulation is the development of several Kirchhoff plate bending triangular elements with the standard nine degrees of freedom. Linear curvature variations are sampled along the three sides with the corners as gage reading points. These sample values are interpolated over the triangle using three schemes. Two schemes merge back to conventional ANS elements, one being identical to the Discrete Kirchhoff Triangle (DKT), whereas the third one produces two new ANDES elements. Numerical experiments indicate that one of the ANDES element is relatively insensitive to distortion compared to previously derived high-performance plate-bending elements, while retaining accuracy for nondistorted elements.

Militello, Carmelo↗
NASA NTRS1992

The distribution of anorthosite on the nearside of the moon

Remote sensing studies of lunar basin and crater deposits were conducted in order to determine the composition of surface units and to investigate the stratigraphy of the lunar crust. These studies have combined both visible and near-IR spectral observations with multispectral imaging in order to determine the lithology of relatively small portions (2-10 km) of the lunar surface. Numerous deposits of pure anorthosite (plagioclase greater than 90 percent) were identified and an interesting pattern has emerged. The purposes of this report are as follows: (1) to summarize the results of our previous studies of the distribution of anorthosite; (2)to present new findings concerning the modes of occurrence of lunar anorthosite; and (3) to assess the implication for the magma ocean hypothesis.

Hawke, B. R.↗
NASA NTRS1991

A sensitive method for measuring atmospheric concentrations of sulfur dioxide

A new method for measuring tropospheric sulfur dioxide concentrations is proposed which is based on the mist chamber sampling method. At the present stage of development, the detection limit of the method is approximately 20 parts per trillion for a 45-min sampling time, with lower concentrations detectable with lower precision. The overall reproducibility of the method (+/-95 percent confidence intervals) is estimated at +/-10 percent. The technique is relatively simple, inexpensive, and lightweight, making it ideally suited for numerous field applications in atmospheric chemistry and biogeochemical studies from both ground-based and airborne platforms.

Klemm, O.↗
NASA NTRS1992

Compression of color-mapped images

In a standard image coding scenario, pixel-to-pixel correlation nearly always exists in the data, especially if the image is a natural scene. This correlation is what allows predictive coding schemes (e.g., DPCM) to perform efficient compression. In a color-mapped image, the values stored in the pixel array are no longer directly related to the pixel intensity. Two color indices which are numerically adjacent (close) may point to two very different colors. The correlation still exists, but only via the colormap. This fact can be exploited by sorting the color map to reintroduce the structure. The sorting of colormaps is studied and it is shown how the resulting structure can be used in both lossless and lossy compression of images.

Hadenfeldt, A. C.↗
NASA NTRS1991

Fractal interfaces and product generation in the two dimensional mixing layer

The dependence of product generation on Peclet and Reynolds numbers in a numerically simulated, reacting, two dimensional, temporally growing mixing layer is related theoretically to the fractal dimension of the passive scalar interfaces. This relation is verified using product generation measurements and dimensions derived from a standard box counting technique. A transition from a low initial dimension to a higher one of approximately 5/3 is identified and shown to be associated to the kinematic distortion on the flow field during the first pairing interaction. It is suggested that the structures responsible for this transition are non-deterministic, non-random, inhomogeneous fractals. Only the large scales are involved. No further transitions, either in the spectra of the vorticity field or in the mixing behavior, are found for Reynolds numbers up to 90,000.

Jimenez, Javier↗
NASA NTRS1951

Analysis of thrust augmentation of turbojet engines by water injection at compressor inlet including charts for calculating compression processes with water injection

A psychrometric chart having total pressure (sum of partial pressures of air and water vapor) as a variable, a Mollier diagram for air saturated with water vapor, and charts showing the thermodynamic properties of various air-water vapor and exhaust gas-water vapor mixtures are presented as aids in calculating the thrust augmentation of a turbojet engine resulting from the injection of water at the compressor inlet. Curves are presented that show the theoretical performance of the augmentation method for various amounts of water injected and the effects of varying flight Mach number, altitude, ambient-air temperature, ambient relative humidity, compressor pressure ratio, and inlet-diffuser efficiency. Numerical examples, illustrating the use of the psychrometric chart and the Mollier diagram in calculating both compressor-inlet and compressor-outlet conditions when water is injected at the compressor inlet, are presented.

Wilcox, E Clinton↗
NASA NTRS1993

Relaxation schemes for Chebyshev spectral multigrid methods

Two relaxation schemes for Chebyshev spectral multigrid methods are presented for elliptic equations with Dirichlet boundary conditions. The first scheme is a pointwise-preconditioned Richardson relaxation scheme and the second is a line relaxation scheme. The line relaxation scheme provides an efficient and relatively simple approach for solving two-dimensional spectral equations. Numerical examples and comparisons with other methods are given.

Kang, Yimin↗
NASA NTRS1994

Energy deposition and K-shell ionization of supernovae

We present explicit evaluation of K-shell ionization caused by photoelectric absorption, by Compton scattering, and by electron collisions in supernovae envelopes. Each process derives from the radioactivity that dominates the bolometric luminosity of the model supernovae. We include the ionization by bremsstrahlung. We find that electron-collision ionization with decelerating Compton electrons is the dominant process of K-shell ionization for light elements (lighter than Si) whereas photoelectric K-shell absorption of Comptonized gammas is the dominant process of K-shell ionization for heavy elements (like Fe). The relative importance of the ionization processes also depends upon composition and upon time, as well as upon supernova type. For Type Ia all four ionization sources are of importance. We present fitting formulas for the deposition of radioactive power in each class of models studied. The energy deposition of all mechanisms decreases with time, including the rate of K-shell ionizations of all elements throughout the supernovae model evaluated herein. We compare the relative contributions to the energy deposition and confirm with greater numerical detail results already established.

The, Lih-Sin↗
NASA NTRS1995

Material Measurements Using Groundplane Apertures

A technique for material parameter determination using an aperture in a groundplane is studied. The material parameters are found by relating the measured reflected field in the aperture to a numerical model. Two apertures are studied which can have a variety of different material configurations covering the aperture. The aperture cross-sections studied are rectangular and coaxial. The material configurations involved combinations of single layer and dual layers with or without a resistive exterior resistive sheet. The resistivity of the resistive sheet can be specified to simulate a perfect electric conductor (PEC) backing (0 Ohms/square) to a free space backing (infinity Ohms/square). Numerical parameter studies and measurements were performed to assess the feasibility of the technique.

Komisarek, K.↗
NASA NTRS1998

The Use of a Satellite Climatological Data Set to Infer Large Scale Three Dimensional Flow Characteristics

Ever since the first satellite image loops from the 6.3 micron water vapor channel on the METEOSAT-1 in 1978, there have been numerous efforts (many to a great degree of success) to relate the water vapor radiance patterns to familiar atmospheric dynamic quantities. The realization of these efforts is becoming evident with the merging of satellite derived winds into predictive models (Velden et al., 1997; Swadley and Goerss, 1989). Another parameter that has been quantified from satellite water vapor channel measurements is upper tropospheric relative humidity (UTH) (e.g., Soden and Bretherton, 1996; Schmetz and Turpeinen, 1988). These humidity measurements, in turn, can be used to quantify upper tropospheric water vapor and its transport to more accurately diagnose climate changes (Lerner et al., 1998; Schmetz et al. 1995a) and quantify radiative processes in the upper troposphere. Also apparent in water vapor imagery animations are regions of subsiding and ascending air flow. Indeed, a component of the translated motions we observe are due to vertical velocities. The few attempts at exploiting this information have been met with a fair degree of success. Picon and Desbois (1990) statistically related Meteosat monthly mean water vapor radiances to six standard pressure levels of the European Centre for Medium Range Weather Forecast (ECMWF) model vertical velocities and found correlation coefficients of about 0.50 or less. This paper presents some preliminary results of viewing climatological satellite water vapor data in a different fashion. Specifically, we attempt to infer the three dimensional flow characteristics of the mid- to upper troposphere as portrayed by GOES VAS during the warm ENSO event (1987) and a subsequent cold period in 1998.

Lerner, Jeffrey A.↗
NASA NTRS1998

Arctic Sea Ice Motion from Wavelet Analysis of Satellite Data

Wavelet analysis of DMSP SSM/I (Special Sensor Microwave/Imager) 85 GHz and 37 GHz radiance data, SMMR (Scanning Multichannel Microwave Radiometer) 37 GHz, and NSCAT (NASA Scatterometer) 13.9 GHZ data can be used to obtain daily sea ice drift information for both the northern and southern polar regions. The derived maps of sea ice drift provide both improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars (SAR). Examples of derived ice-drift maps in the Arctic illustrate large-scale circulation reversals within a period of a couple weeks. Comparisons with ice displacements derived from buoys show good quantitative agreement. NASA Scatterometer (NSCAT) 13.9 GHZ data have been also used for wavelet analysis to derive sea-ice drift. First, the 40' incidence-angle, sigma-zero (surface roughness) daily map of whole Arctic region with 25 km of pixel size from satellite's 600 km swath has been constructed. Then, the similar wavelet transform procedure to SSM/I data can be applied. Various scales of wavelet transform and threshold have been tested. By overlaying , neighbor filtering, and block-averaging the results of multiscale wavelet transforms, the final sea ice drift vectors are much smooth and representative to the sea ice motion. This wavelet analysis procedure is robust and can make a major contribution to the understanding of ice motion over large areas at relatively high temporal resolutions. The results of wavelet analysis of SSM/I and NSCAT images and buoy data can be merged by some data fusion techniques and will help to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model.

Liu, Antony K.↗
NASA NTRS1999

S-Wave Dispersion Relations: Exact Left Hand E-Plane Discontinuity from the Born Series

We show, for a superposition of Yukawa potentials, that the left hand cut discontinuity in the complex E plane of the (S-wave) scattering amplitude is given exactly, in an interval depending on n, by the discontinuity of the Born series stopped at order n. This also establishes an inverse and unexpected correspondence of the Born series at positive high energies and negative low energies. We can thus construct a viable dispersion relation (DR) for the partial (S-) wave amplitude. The high numerical precision achievable by the DR is demonstrated for the exponential potential at zero scattering energy. We also briefly discuss the extension of our results to Field Theory.

Bessis, D.↗
NASA NTRS2000

Sea Ice Motion from Wavelet Analysis of Satellite Data

Wavelet analysis of NASA scatterometer (NSCAT) backscatter and Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave/Imager (SSM/I) radiance data can be used to obtain daily sea ice drift information for the Arctic region. This technique provides improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars. Comparisons with ice motion derived from ocean buoys give good quantitative agreement. Both comparison results from NSCAT and SSM/I are compatible, and the results from NSCAT can definitely complement that from SSM/I when there are cloud or surface effects. Then three sea-ice drift daily results from NSCAT, SSM/I, and buoy data can be merged as a composite map by some data fusion techniques. The ice flow streamlines are highly correlated with surface air pressure contours. Examples of derived ice-drift maps in December 1996 illustrate large-scale circulation reversals over a period of four days. A method for deriving divergence and shear at the large-scale has been developed and comparison between buoys and satellite results shows a good agreement. These calibrated/validated results indicate that NSCAT, SSM/I merged daily ice motion are suitably accurate to identify and closely locate sea ice processes, and to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model. For demonstration purpose, the ice velocities derived from satellite data are compared with the ice velocities derived from a coupled ice-ocean interaction model. The comparison reveals that the general circulation patterns of the two are quite similar but the ice velocity differences between the two are quite significant. In order to quantify the wind effects on ice motion, empirical orthogonal functions (EOF) are used in the principal component analysis for both ice motion and pressure field. Some preliminary results of sea-ice motion from QuikScat will also be presented.

Liu, Antony K.↗
NASA NTRS2001

Calibration of Sea Ice Motion from QuikSCAT with those from SSM/I and Buoy

QuikSCAT backscatter and DMSP SSM/I radiance data are used to derive sea ice motion for both the Arctic and Antarctic region using wavelet analysis method. This technique provides improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing satellite data from Synthetic Aperture Radar (SAR). Sea ice motion of the Arctic for the period from October 1999 to March 2000 derived from QuikSCAT and SSM/I data agrees well with that derived from ocean buoys quantitatively. Thus the ice tracking results from QuikSCAT and SSM/I are complement to each other, Then, three sea-ice drift daily results from QuikSCAT, SSM/I, and buoy data can be merged to generate composite maps with more complete coverage of sea ice motion than those from single data source. A series of composite sea ice motion maps for December 1999 show that the major circulation patterns of sea ice motion are changing and shifting significantly within every four days and they are dominated by wind forcing. Sea-ice drift in the summer can not be derived from NSCAT and SSM/I data. In later summer of 1999 (in September), however, QuikSCAT data can provide good sea ice motion information in the Arctic. QuiksCAT can also provide at least partial sea ice motion information until June 15 in early summer 1999. For the Antarctic, case study shows that sea ice motion derived from QuikSCAT data is predominantly forced by and is consistent with wind field derived from QuikSCAT around the polar region. These calibrated/validated results indicate that QuikSCAT, SSM/I, and buoy merged daily ice motion are suitably accurate to identify and closely locate sea ice processes, and to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model.

Liu, Antony K.↗
NASA NTRS2007

Core Physics and Kinetics Calculations for the Fissioning Plasma Core Reactor

Highly efficient, compact nuclear reactors would provide high specific impulse spacecraft propulsion. This analysis and numerical simulation effort has focused on the technical feasibility issues related to the nuclear design characteristics of a novel reactor design. The Fissioning Plasma Core Reactor (FPCR) is a shockwave-driven gaseous-core nuclear reactor, which uses Magneto Hydrodynamic effects to generate electric power to be used for propulsion. The nuclear design of the system depends on two major calculations: core physics calculations and kinetics calculations. Presently, core physics calculations have concentrated on the use of the MCNP4C code. However, initial results from other codes such as COMBINE/VENTURE and SCALE4a. are also shown. Several significant modifications were made to the ISR-developed QCALC1 kinetics analysis code. These modifications include testing the state of the core materials, an improvement to the calculation of the material properties of the core, the addition of an adiabatic core temperature model and improvement of the first order reactivity correction model. The accuracy of these modifications has been verified, and the accuracy of the point-core kinetics model used by the QCALC1 code has also been validated. Previously calculated kinetics results for the FPCR were described in the ISR report, "QCALC1: A code for FPCR Kinetics Model Feasibility Analysis" dated June 1, 2002.

Butler, C.↗
NASA NTRSDate not supplied

Benefit and Critical Factors for the Performance of the Boundary Layer Ingesting Propulsion

The benefit of the boundary layer ingestion (BLI) is described in the perspective of the propulsion and engine development. A power saving map of the BLI engines is derived based on the correlation of the shape factor of the ingested boundary layer and the propulsive efficiency. The ratio of the mass flow rate between BLI and non-BLI propulsors is introduced to quantify the power saving of the BLI engine relative to a corresponding clean inlet flow engine that generates the same amount of net thrust. The wake recovery factor from the jet flow is applied to find the optimum sizing of the engine for the given design requirement. The effects of the fan pressure ratio of the propulsor toward the power saving are also investigated to seek out the feasible range of the BLI propulsor design. The derived correlation is validated with CFD analyses. Simulation models using variously sized engines relative to an influencing body are developed so that a numerical experiment is carried out with the various shape factors. The effect of the propulsion airframe integration toward the aerodynamic forces on the fuselage and nacelle is investigated via parametric study of the engine sizing and fan pressure ratio for the wake ingestion models. The propulsor efficiency is quantified toward the saving in the actual shaft power and correlated with the efficiency penalty of the BLI propulsor. The correlation is also validated via comparison of the turbo-machinery CFD results for BLI2DTF and conventional transonic fan stages.

Boundary Layer Ingestion↗
NASA NTRS1989

Higher order dispersion in the propagation of a gravity wave packet

To the first order of approximation, the complex amplitude of a wave packet in an anisotropic and dispersive medium is convected with the group of velocity. However, a gravity wave is a vector wave. Its wave packet must be formed by superposition of various wave numbers with corresponding frequencies, as is the case for scalar waves, and additionally by superposing many eigenmodes which also depend on the wave number. To represent the vector wave packet self-consistently, it is found that a gradient term must be included in the expansion. For a Guassian wave packet, this gradient term is shown to have important implications on the velocity vector as represented by its hodograph. Numerical results show that the hodograph is influenced by the location of the relative position of interest from the center of a Gaussian pulse. Higher order expansion shows that an initial Gaussian wave packet will retain its Gaussian shape as it propagates, but the pulse will spread in all directions with its major axis undergoing a rotation. Numerical results indicate that these higher order dispersive effects may be marginally observable in the atmosphere.

Yeh, K. C.↗