Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 469 records · Page 26

Solution of three-dimensional afterbody flow using reduced Navier-Stokes equations

The flow over afterbody geometries was investigated using the reduced Navier-Stokes (RNS) approximation. Both pressure velocity flux-split and composites velocity primitive variable formulations were considered. Pressure or pseudopotential relaxation procedures are combined with sparse matrix or coupled strongly implicit algorithms to form a three-dimensional solver for general non-orthogonal coordinates. Three-dimensional subsonic and transonic viscous/inviscid interacting flows were evaluated. Solutions with and without regions of recirculation were obtained.

Almahroos, H. M. H.↗

Asteroid orbital error analysis: Theory and application

We present a rigorous Bayesian theory for asteroid orbital error estimation in which the probability density of the orbital elements is derived from the noise statistics of the observations. For Gaussian noise in a linearized approximation the probability density is also Gaussian, and the errors of the orbital elements at a given epoch are fully described by the covariance matrix. The law of error propagation can then be applied to calculate past and future positional uncertainty ellipsoids (Cappellari et al. 1976, Yeomans et al. 1987, Whipple et al. 1991). To our knowledge, this is the first time a Bayesian approach has been formulated for orbital element estimation. In contrast to the classical Fisherian school of statistics, the Bayesian school allows a priori information to be formally present in the final estimation. However, Bayesian estimation does give the same results as Fisherian estimation when no priori information is assumed (Lehtinen 1988, and reference therein).

Muinonen, K.↗

Explicit and implicit compact high-resolution shock-capturing methods for multidimensional Euler equations 1: Formulation

Two classes of explicit compact high-resolution shock-capturing methods for the multidimensional compressible Euler equations for fluid dynamics are constructed. Some of these schemes can be fourth-order accurate away from discontinuities. For the semi-discrete case their shock-capturing properties are of the total variation diminishing (TVD), total variation bounded (TVB), total variation diminishing in the mean (TVDM), essentially nonoscillatory (ENO), or positive type of scheme for 1-D scalar hyperbolic conservation laws and are positive schemes in more than one dimension. These fourth-order schemes require the same grid stencil as their second-order non-compact cousins. One class does not require the standard matrix inversion or a special numerical boundary condition treatment associated with typical compact schemes. Due to the construction, these schemes can be viewed as approximations to genuinely multidimensional schemes in the sense that they might produce less distortion in spherical type shocks and are more accurate in vortex type flows than schemes based purely on one-dimensional extensions. However, one class has a more desirable high-resolution shock-capturing property and a smaller operation count in 3-D than the other class. The extension of these schemes to coupled nonlinear systems can be accomplished using the Roe approximate Riemann solver, the generalized Steger and Warming flux-vector splitting or the van Leer type flux-vector splitting. Modification to existing high-resolution second- or third-order non-compact shock-capturing computer codes is minimal. High-resolution shock-capturing properties can also be achieved via a variant of the second-order Lax-Friedrichs numerical flux without the use of Riemann solvers for coupled nonlinear systems with comparable operations count to their classical shock-capturing counterparts. The simplest extension to viscous flows can be achieved by using the standard fourth-order compact or non-compact formula for the viscous terms.

Yee, H. C.↗

Highly Reduced Forsterite and Enstatite from Stardust Track 61: Implications for Radial Transport of E Asteroid Material

The Stardust Mission returned a large fraction of high-temperature, crystalline material that was radially transported from the inner solar system to the Kuiper Belt [1,2]. The mineralogical diversity found in this single cometary collection points to an even greater number of source materials than most primitive chondrites. In particular, the type II olivine found in Wild 2 includes the three distinct Fe/Mn ratios found in the matrix and chondrules of carbonaceous chondrites (CCs) and unequilibrated ordinary chondrites (UOCs) [3]. We also find that low-Ca pyroxene is quite variable (approximately Fs3-29) and is usually indistinguishable from CC, UOC, and EH3 pyroxene as well. However, occasional olivine and pyroxene compositions are found in Wild 2 that are inconsistent with chondrites. The Stardust track 61 terminal particle (TP) is one such example and is the focus of this study. It s highly reduced forsterite and enstatite is consistent only with that in Aubrites, in which FeO is essentially absent from these phases (less than approximately 0.1 wt.% FeO) [4].

Frank, David R.↗

Linear and nonlinear dynamic analysis of redundant load path bearingless rotor systems

The bearingless rotorcraft offers reduced weight, less complexity and superior flying qualities. Almost all the current industrial structural dynamic programs of conventional rotors which consist of single load path rotor blades employ the transfer matrix method to determine natural vibration characteristics because this method is ideally suited for one dimensional chain like structures. This method is extended to multiple load path rotor blades without resorting to an equivalent single load path approximation. Unlike the conventional blades, it isk necessary to introduce the axial-degree-of-freedom into the solution process to account for the differential axial displacements in the different load paths. With the present extension, the current rotor dynamic programs can be modified with relative ease to account for the multiple load paths without resorting to the equivalent single load path modeling. The results obtained by the transfer matrix method are validated by comparing with the finite element solutions. A differential stiffness matrix due to blade rotation is derived to facilitate the finite element solutions.

Murthy, V. R.↗

Albedo in the ATIC Experiment

ATIC(Advanced Thin Ionization Calorimeter) is a balloon borne experiment designed to measure the cosmic ray composition for elements from hydrogen to iron and their energy spectra from approx.50 GeV to near 100 TeV. It consists of a Si-matrix detector to determine the charge of a CR particle, a scintillator hodoscope for tracking, carbon interaction targets and a fully active BGO calorimeter. ATIC had its first 16-day flight from McMurdo, Antarctica from 28/12/2000 to 13/01/2000. The ATIC flight collected approximately 25 million events. To measure charge of primary particle in presence of radiation scattered back from the interaction and subsequent shower development in the calorimeter a charge detector must be a mosaic of small detector pads so that the pad containing the signal from the incident particle has no additional signal from albedo particles. Therefore the silicon matrix was built of 4480 individual silicon pads each 2 cm x 1.5 cm. The matrix consists of four planes of detectors and the active detector area, in these planes are partially overlapped to completely cover the aperture. The lateral and amplitude distributions of albedo signals in Si-matrix are analyzed for different primary nuclei and different energy deposits in BGO calorimeter. The greater part of albedo signals has Q near 1, where Q = square root of Amplitude(MIP). The albedo distribution exponentially decreases up to Q near 8. These high values are produced by slow protons and plans. There are also a small number of signals of Q > 8, mainly for heavy nucleus primaries. These signals are apparently generated by neutrons. The comparison of the experimental data and simulations with GEANT 3-21 code using QGSM generator for nucleus-nucleus interactions is presented.

Sokolskaya, N. V.↗

Alloy NASA-HR-1

NASA-HR-1 is a high-strength Fe-Ni-base superalloy that resists high-pressure hydrogen environment embrittlement (HEE), oxidation, and corrosion. Originally derived from JBK-75, NASA-HR-1 has exceptional HEE resistance that can be attributed to its gamma-matrix and eta-free (Ni3Ti) grain boundaries. The chemistry was formulated using a design approach capable of accounting for the simultaneous effects of several alloy additions. This approach included: (1) Systematically modifying gamma-matrix compositions based on JBK-75; (2) Increasing gamma (Ni3(Al,Ti)) volume fraction and adding gamma-matrix strengthening elements to obtain higher strength; and (3) Obtaining precipitate-free grain boundaries. The most outstanding attribute of NASA-HR-1 is its ability to resist HEE while showing much improved strength. NASA-HR-1 has approximately 25% higher yield strength than JXK-75 and exhibits tensile elongation of more than 20% with no ductility loss in a hydrogen environment at 5 ksi, an achievement unparalleled by any other commercially available alloy. Its Cr and Ni contents provide exceptional resistance to environments that promote oxidation and corrosion. Microstructural stability was maintained by improved solid solubility of the gamma-matrix, along with the addition of alloying elements to retard eta (Ni3Ti) precipitation. NASA-HR-1 represents a new system that greatly extends the compositional ranges of existing HEE-resistant Fe-Ni-base superalloys.

Chen, Po-Shou↗

Flight service evaluation of PRD-49/epoxy composite panels in wide-bodied commercial transport aircraft

Fairing panels were fabricated to evaluate the fabrication characteristics and flight service performance of PRD-49 (Kevlar-49) a composite reinforcing material and to compare it with the fiberglass which is currently in use. Panel configurations were selected to evaluate the PRD-49 with two resin matrix materials in sandwich and solid laminate construction. Left and right hand versions of these configurations were installed on L-1011's which will accumulate approximately 3000 flight hours per year per aircraft. The direct substitution of PRD-49 for fiberglass produced a twenty-six percent weight reduction on the panel configurations. Examination of these panels revealed that there was no visible difference between the PRD-49 and adjacent fiberglass panels.

Wooley, J. H.↗

Vibration-rotation intensities of SiO

Dipole moment matrix elements have been computed for a large number of transitions of astrophysical interest for the more abundant isotopes of SiO. The wave functions utilized were obtained from a direct solution of the Schroedinger equation with an accurate RKR potential. The dipole moment function, in the form of a Pade approximant, was chosen to reproduce the experimental measurements near equilibrium, to have the proper united and separated atom limits, and to have the correct long-range asymptotic dependence on internuclear separation. Because of the large number of transitions involved, and to facilitate applications, the squares of the dipole moment matrix elements were fitted by a least-squares procedure to polynomials in v and J. In addition, Einstein A coefficients are given for observed maser transitions and for selected vibrational bands. These latter are compared with previous calculations, and it is concluded that for the higher Delta v transitions, the present results represent a significant improvement.

Tipping, R. H.↗

Numerical algorithms for steady and unsteady incompressible Navier-Stokes equations

The numerical analysis of the incompressible Navier-Stokes equations are becoming important tools in the understanding of some fluid flow problems which are encountered in research as well as in industry. With the advent of the supercomputers, more realistic problems can be studied with a wider choice of numerical algorithms. An alternative formulation is presented for viscous incompressible flows. The incompressible Navier-Stokes equations are cast in a velocity/vorticity formulation. This formulation consists of solving the Poisson equations for the velocity components and the vorticity transport equation. Two numerical algorithms for the steady two-dimensional laminar flows are presented. The first method is based on the actual partial differential equations. This uses a finite-difference approximation of the governing equations on a staggered grid. The second method uses a finite element discretization with the vorticity transport equation approximated using a Galerkin approximation and the Poisson equations are obtained using a least squares method. The equations are solved efficiently using Newton's method and a banded direct matrix solver (LINPACK). The method is extended to steady three-dimensional laminar flows and applied to a cubic driven cavity using finite difference schemes and a staggered grid arrangement on a Cartesian mesh. The equations are solved iteratively using a plane zebra relaxation scheme. Currently, a two-dimensional, unsteady algorithm is being developed using a generalized coordinate system. The equations are discretized using a finite-volume approach. This work will then be extended to three-dimensional flows.

Hafez, Mohammed↗

Computers boost structural technology

Derived from matrix methods of structural analysis and finite element methods developed over the last three decades, computational structures technology (CST) blends computer science, numerical analysis, and approximation theory into structural analysis and synthesis. Recent significant advances in CST include stochastic-based modeling, strategies for performing large-scale structural calculations on new computing systems, and the integration of CST with other disciplinary modules for multidisciplinary analysis and design. New methodologies have been developed at NASA for integrated fluid-thermal structural analysis and integrated aerodynamic-structure-control design. The need for multiple views of data for different modules also led to the development of a number of sophisticated data-base management systems. For CST to play a role in the future development of structures technology and in the multidisciplinary design of future flight vehicles, major advances and computational tools are needed in a number of key areas.

Noor, Ahmed K.↗

Discretized Newton-relaxation solution of high resolution flux-difference split schemes

A method is presented for the implicit solution of high resolution flux-difference split schemes for which it is impractical to obtain the Jacobian matrix analytically. The method is used frequently in numerical schemes for solving nonlinear systems of equations, but apparently has not been used to date for approximating the Jacobian of the Roe numerical flux vectors. The method consists of obtaining a discretized Jacobian by using simple finite-differences of the flux vectors, and then using this discretized Jacobian in a Newton-relaxation solution method. The method is applied to both compressible and incompressible flows. Numerical results are presented that demonstrate the method works for both inviscid and viscous flows.

Whitfield, David L.↗

Quiet-Time Spectra and Abundances of Energetic Particles During the 1996 Solar Minimum

We report the energy spectra and abundances of ions with atomic number, Z, in the interval Z is greater than or equal to 2 and Z is less than or equal to 36 and energies approximately 3-20 MeV/amu for solar and interplanetary quiet periods between 1994 November and 1998 April as measured by the large-geometry Low Energy Matrix Telescope (LEMT) telescope on the Wind spacecraft near Earth. The energy spectra show the presence of galactic (GCR) and "anomalous" cosmic ray (ACR) components, depending on the element. ACR components are reported for Mg and Si for the first time at 1 AU and the previous observation of S and Ar is confirmed. However, only GCR components are clearly apparent for the elements Ca, Ti, Cr, Fe, as well as for C. New limits are placed on a possible ACR contribution for other elements, including Kr.

Reames, Donald V.↗

Formation of Metal and Silicate Globules in Gujba: A New Bencubbin-like Meteorite Fall

Gujba is a coarse-grained meteorite fall composed of 41 vol% large kamacite globules, 20 vol% large light-colored silicate globules with cryptocrystalline, barred pyroxene and barred olivine textures, 39 vol% dark-colored, silicate-rich matrix, and rare refractory inclusions. Gujba resembles Bencubbin and Weatherford in texture, oxygen-isotopic composition and in having high bulk delta N-15 values (approximately +685%0). The He-3 cosmic-ray exposure age of Gujba (26 +/- 7 Ma) is essentially identical to that of Bencubbin, suggesting that they were both reduced to meter-size fragments in the same parent-body collision. The Gujba metal globules exhibit metal-troilite quench textures and vary in their abundances of troilite and volatile siderophile elements. We suggest that the metal globules formed as liquid droplets either via condensation in an impact-generated vapor plume or by evaporation of preexisting metal particles in a plume. The lower the abundance of volatile elements in the metal globules, the higher the globule quench temperature. We infer that the large silicate globules also formed from completely molten droplets; their low volatile-element abundances indicate that they also formed at high temperatures, probably by processes analogous to those that formed the metal globules. The coarse-grained Bencubbin-Weatherford-Gujba meteorites may represent a depositional component from the vapor cloud enriched in coarse and dense particles. A second class of Bencubbin-like meteorites (represented by Hammadah a1 Hamra 237 and QUE 94411) may be a finer fraction derived from the same vapor cloud

Rubin, Alan E.↗

Structure and dynamics of amorphous water ice

Further insight into the structure and dynamics of amorphous water ice, at low temperatures, was obtained by trapping in it Ar, Ne, H2, and D2. Ballistic water-vapor deposition results in the growth of smooth, approximately 1 x 0.2 micrometer2, ice needles. The amorphous ice seems to exist in at least two separate forms, at T < 85 K and at 85 < T < 136.8 K, and transform irreversibly from one form to the other through a series of temperature-dependent metastable states. The channels formed by the water hexagons in the ice are wide enough to allow the free penetration of H2 and D2 into the ice matrix even in the relatively compact cubic ice, resulting in H2-(D2-) to-ice ratios (by number) as high as 0.63. The larger Ar atoms can penetrate only into the wider channels of amorphous ice, and Ne is an intermediate case. Dynamic percolation behavior explains the emergence of Ar and Ne (but not H2 and D2) for the ice, upon warming, in small and big gas jets. The big jets, each containing approximately 5 x 10(10) atoms, break and propel the ice needles. Dynamic percolation also explains the collapse of the ice matrix under bombardment by Ar , at a pressure exceeding 2.6 dyn cm-2, and the burial of huge amounts of gas inside the collapsed matrix, up to an Ar-to-ice of 3.3 (by number). The experimental results could be relevant to comets, icy satellites, and icy grain mantles in dense interstellar clouds.

NASA Program Exobiology↗

Computer solutions of Wu's equations for compressible flow through turbomachines

Two computers programs, known as Matrix Through-Flow and Matrix Blade-To-Blade, for analyzing the meridional and blade-to-blade flow patterns are described. The numerical solutions are obtained by finite difference approximations to the governing Poisson-type differential equations for the stream function. Solutions for several turbomachines, giving flow patterns and velocity distributions, are included.

Smith, D. J. L.↗

Elastic-plastic deformation of a metal-matrix composite coupon with a center slot

A comprehensive experimental analysis of deformations of the surface of a metal-matrix specimen is reported. The specimen is a 6-ply 0 + or - 45 sub s boron-aluminum tensile coupon with a central slot. Moire interferometry is used for high-sensitivity whole-field measurements of in-plane displacements. Normal and shear strains are calculated from displacement gradients. Displacement fields are analyzed at various load levels from 15% to 95% of the failure load. Deformations of the boron fibers could be distinguished from those of the matrix. Highly localized plastic slip zones occur tangent to the ends of the slot. Shear strains and concurrent transverse compressive strains in the slip zones reach approximately 10% and 1%, respectively. Upon unloading, elastic recovery in surrounding regions causes a reverse plastic shear strain in the slip zone of about 4%. Longitudinal normal strains on the unslotted ligament peak at the slot boundary at about 1% strain. The strain concentration factor at the end of the slot decreases with load level and the advance of plasticity.

Post, D.↗

Response of thick, notched laminates subjected to tension-compression cyclic loads

The fatigue response of a (0/24/90/-45)s4 T300-5208 graphite-epoxy laminate with a drilled centerhole subjected to fully reversed tension-compression (R=-1, T-C) constant amplitude loading was investigated. Damage evaluation techniques such as stiffness monitoring, penetrant-enchanced X-ray radiography, C-scan, laminate deeply and residual strength were used to establish the mechanisms of damage development as well as the relations between this damage and the stiffness, strength and life of the laminate. Damage initiated at the hole as matrix cracking parallel to the fibers in all plies. Matrix cracks had a significant effect on delamination initiation and growth. Delaminations initiated near the surface in the densely cracked region at the hole and grew along major matrix cracks. Compressive properties degraded more rapidly than tensile properties. At the stress levels used, residual tensile strength increased early in the fatigue life and remained approximately constant to near the end of life, when failure was precipitated by excessive laminate instability during the compressive portion of the loading.

Bakis, C. E.↗