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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A three-dimensional pressure flux-split RNS application to sub/supersonic flow in inlets and ducts

The reduced Navier-Stokes (RNS) formulation is combined with a pressure based flux split procedure for the computational analysis of three dimensional flow. Only lowest order diffusion, that is required to satisfy the no-slip boundary condition is retained in this approximation. The governing equations are appropriately differenced such that the physical boundary conditions are closely coupled and combined to form a closed discrete system. No numerical or characteristic type boundary conditions are required. A global pressure relaxation procedure, i.e. multi-sweep PNS, is considered. A sparse matrix direct solver is applied to the crossplane solution. At the outflow boundary only the pressure or pressure gradient is prescribed. The resulting formulation has been applied to the solution of three-dimensional subsonic/supersonic internal flows in inlets and ducts.

Khosla, P. K.↗
NASA NTRS1997

Comparison of Several Dissipation Algorithms for Central Difference Schemes

Several algorithms for introducing artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical results are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme.

Swanson, R. C.↗
NASA NTRS2001

Carbon Nanotube Anodes Being Evaluated for Lithium Ion Batteries

The NASA Glenn Research Center is evaluating the use of carbon nanotubes as anode materials for thin-film lithium-ion (Li) batteries. The motivation for this work lies in the fact that, in contrast to carbon black, directed structured nanotubes and nanofibers offer a superior intercalation media for Li-ion batteries. Carbon lamellas in carbon blacks are circumferentially oriented and block much of the particle interior, rendering much of the matrix useless as intercalation material. Nanofibers, on the other hand, can be grown so as to provide 100-percent accessibility of the entire carbon structure to intercalation. These tubes can be visualized as "rolled-up" sheets of carbon hexagons (see the following figure). One tube is approximately 1/10,000th the diameter of a human hair. In addition, the high accessibility of the structure confers a high mobility to ion-exchange processes, a fundamental for the batteries to respond dynamically because of intercalation.

Raffaelle, Ryne P.↗
NASA NTRS1979

Vibrational dependence of pressure induced spectral linewidths and line shifts - Application of the infinite order sudden scattering approximation

The infinite order sudden (IOS) approximation to molecular rotation is applied to simplify the theory of linewidths and shifts in vibration-rotation spectra. This approximation is expected to be most accurate for hard, short-range collisions and is therefore complementary to Anderson theory which is best for weak, glancing collisions. The IOS approximation predicts identical linewidths and shifts for P- and R-branch transitions with the same line number. It also predicts zero line shifts for pure rotational spectra. The dependence of linewidths and shifts on vibrational band is seen to be due mainly to variation in diagonal vibrational matrix elements of the intermolecular potential. Calculations are performed for the 0-0, 0-1, and 0-2 bands of CO perturbed by He, using a theoretical interaction potential with no semiempirical or adjustable parameters; results are in satisfactory accord with experimental data.

Green, S.↗
NASA NTRS1982

Ordered mixed-layer structures in the Mighei carbonaceous chondrite matrix

High resolution transmission electron microscopy of the Mighei carbonaceous chondrite matrix has revealed the presence of a new mixed layer structure material. This mixed-layer material consists of an ordered arrangement of serpentine-type (S) and brucite-type (B) layers in the sequence SBBSBB. Electron diffraction and imaging techniques show that the basal periodicity is approximately 17 A. Discrete crystals of SBB-type material are typically curved, of small size (less than 1 micron) and show structural variations similar to the serpentine group minerals. Mixed-layer material also occurs in association with planar serpentine. Characteristics of SBB-type material are not consistent with known terrestrial mixed-layer clay minerals. Evidence for formation by a condensation event or by subsequent alteration of pre-existing material is not yet apparent.

Mackinnon, I. D. R.↗
NASA NTRS1982

ACCESS computer program for the synthesis of large structural systems

The structural-synthesis computer code ACCESS-3, originally limited to the optimization of relatively simple problems involving truss, membrane, and shear panel elements, has been restructured to make it applicable to structures with a larger number of degrees-of-freedom, design variables, and constraints. The modifications made allow the use of six degrees-of-freedom per node and permit the out-of-core solution of large matrix equations encountered at the analysis stage. Also, finite elements which permit pure bending as well as combined bending and membrane (axial) behavior have been added to satisfy the need for a more comprehensive set of elements. Like its predecessor, the new version of the program combines the approximation concepts with primal or dual formulations to provide an efficient optimization tool.

Fleury, C.↗
NASA NTRS1993

Control of nonlinear systems represented in quasilinear form

Methods to synthesize controllers for nonlinear systems are developed by exploiting the fact that under mild differentiability conditions, systems of the form: x-dot = f(x) + G(x)u can be represented in quasilinear form, viz: x-dot = A(x)x + B(x)u. Two classes of control methods are investigated. The first is zero-look-ahead control, where the control input depends only on the current values of A(x) and B(x). For this case the control input is computed by continuously solving a matrix Riccati equation as the system progresses along a trajectory. The second is controllers with look-ahead, where the control input depends on the future behavior of A(x) and B(x). These controllers use the similarity between quasilinear systems and linear time varying systems to find approximate solutions to optimal control type problems. The methods that are developed are not guaranteed to be globally stable. However in simulation studies they were found to be useful alternatives for synthesizing control laws for a general class of nonlinear systems.

Coetsee, Josef A.↗
NASA NTRS1994

Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

Slack, David C.↗
NASA NTRS2002

Efficient Reformulation of the Thermoelastic Higher-order Theory for Fgms

Functionally graded materials (FGMs) are characterized by spatially variable microstructures which are introduced to satisfy given performance requirements. The microstructural gradation gives rise to continuously or discretely changing material properties which complicate FGM analysis. Various techniques have been developed during the past several decades for analyzing traditional composites and many of these have been adapted for the analysis of FGMs. Most of the available techniques use the so-called uncoupled approach in order to analyze graded structures. These techniques ignore the effect of microstructural gradation by employing specific spatial material property variations that are either assumed or obtained by local homogenization. The higher-order theory for functionally graded materials (HOTFGM) is a coupled approach developed by Aboudi et al. (1999) which takes the effect of microstructural gradation into consideration and does not ignore the local-global interaction of the spatially variable inclusion phase(s). Despite its demonstrated utility, however, the original formulation of the higher-order theory is computationally intensive. Herein, an efficient reformulation of the original higher-order theory for two-dimensional elastic problems is developed and validated. The use of the local-global conductivity and local-global stiffness matrix approach is made in order to reduce the number of equations involved. In this approach, surface-averaged quantities are the primary variables which replace volume-averaged quantities employed in the original formulation. The reformulation decreases the size of the global conductivity and stiffness matrices by approximately sixty percent. Various thermal, mechanical, and combined thermomechanical problems are analyzed in order to validate the accuracy of the reformulated theory through comparison with analytical and finite-element solutions. The presented results illustrate the efficiency of the reformulation and its advantages in analyzing functionally graded materials.

Bansal, Yogesh↗
NASA NTRS1998

On Some Numerical Dissipation Schemes

Several schemes for introducing an artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar dissipation and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical solutions are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme. The coarse-grid accuracy for the original CUSP scheme is improved by modifying the limiter function used with the scheme, giving comparable accuracy to that obtained with the MATD scheme. The modifications reduce the background dissipation and provide control over the regions where the scheme can become first order.

Swanson, R. C.↗
NASA NTRS2018

Radiative Transfer in a Discrete Random Medium Adjacent to a Half-Space with a Rough Interface

For a macroscopically plane-parallel discrete random medium, the boundary conditions for the specific coherency dyadic at a rough interface are derived. The derivation is based on a modification of the Twersky approximation for a scattering system consisting of a group of particles and the rough surface, and reduces to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave propagating in a discrete random medium with non-scattering boundaries. In a matrix-form setting, the boundary conditions for the specific coherency dyadic imply the boundary conditions for specific intensity column vectors which in turn, yield the expressions for the reflection and transmission matrices. The derived expressions are shown to be identical to those obtained by applying a phenomenological approach based on a facet model to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave.

Adrian Doicu↗
NASA NTRS1979

High temperature dynamic modulus and damping of aluminum and titanium matrix composites

Dynamic modulus and damping capacity property data were measured from 20 to over 500 C for unidirectional B/Al (1100), B/Al (6061), B/SiC/Al (6061), Al2O3/Al, SiC/Ti-6Al-4V, and SiC/Ti composites. The measurements were made under vacuum by the forced vibration of composite bars at free-free flexural resonance near 2000 Hz and at amplitudes below 0.000001. Whereas little variation was observed in the dynamic moduli of specimens with approximately the same fiber content (50 percent), the damping of B/Al composites was found at all temperatures to be significantly greater than the damping of the Al2O3/Al and SiC/Ti composites. For those few situations where slight deviations from theory were observed, the dynamic data were examined for information concerning microstructural changes induced by composite fabrication and thermal treatment. The 270 C damping peak observed in B/Al (6061) composites after heat treatment above 460 C appears to be the result of a change in the 6061 aluminum alloy microstructure induced by interaction with the boron fibers. The growth characteristics of the damping peak suggest its possible value for monitoring fiber strength degration caused by excess thermal treatment during B/Al (6061) fabrication and use.

Dicarlo, J. A.↗
NASA NTRS1972

Lumped parameter analysis of a stringer reinforced plate excited by band limited noise

The maximum root-mean-square response of a square clamped plate, subjected to a stationary random excitation, was determined both theoretically and empirically. For the tests, a 40 x 40 x 1/4 aluminum sheet was milled out to provide a thin waffle plate with 9 panels. The plate was acoustically excited, with the frequency range limited between 25 Hz and 500 Hz. The root-mean-square power of the random excitations, 149 decibels, was nearly constant for all natural frequencies of the plate system. Strain gage readings at 14 points were plotted as power spectral densities and root-mean-square displacements. A modal damping ratio matrix was determined from the response curves. Damping ratios were found to be a function of frequency, but not of location on the plate. For comparisons with the theoretical response, the plate was subdivided into 14 modes and anlayzed as a damped, lumped parameter system by use of an approximate normal-mode method. This method gave the best predictions of the power spectral densities for the lower frequencies. That prediction errors increased for the higher frequencies is attributed to improper assumptions for mass distribution and the insufficient number of lumped mass points.

Bilyeu, D. J.↗
NASA NTRS1978

An implicit algorithm for the conservative transonic full potential equation using an arbitrary mesh

A new, implicit approximate factorization (AF) algorithm designed to solve the conservative full-potential equation for the transonic flow past arbitrary airfoils has been developed. The new algorithm uses an upwind bias of the density coefficient to provide stability in supersonic regions. This allows the simple two- and three-banded matrix form of the AF scheme to be retained over the entire flow field, even in regions of supersonic flow. A numerical transformation is used to establish an arbitrary body-fitted finite-difference mesh. Airfoil pressure distributions have been computed and are in good agreement with independent results.

Holst, T. L.↗
NASA NTRS1989

Kaon-nucleus scattering

Two kinds of number density distributions of the nucleus, harmonic well and Woods-Saxon models, are used with the t-matrix that is taken from the scattering experiments to find a simple optical potential. The parameterized two body inputs, which are kaon-nucleon total cross sections, elastic slope parameters, and the ratio of the real to imaginary part of the forward elastic scattering amplitude, are shown. The eikonal approximation was chosen as the solution method to estimate the total and absorptive cross sections for the kaon-nucleus scattering.

Hong, Byungsik↗
NASA NTRS1994

Nonperturbative methods in HZE propagation

An analytical solution to the perturbative multiple collision series of a fragmenting HZE ion beam has limited usefulness since the first collision term has several hundred contributions, the second collision term has tens of thousands of contributions, and each successive collision term progresses to unwieldy computational proportions. Our previous work has revealed the multiple collision terms in the straight-ahead approximation to be simple products of a spatially dependent factor times a linear energy-dependent factor of limited domain and unit normalization. The properties of these forms allow the development of the nonperturbative summation of the series to all orders assuming energy-independent nuclear cross sections as matrix products of a scaled Green's function described herein. This nonperturbative Green's function with multiple scattering correction factors compares well with experiments using 670 MeV/u neon-20 ion beams in thick water targets.

Non-NASA Center↗
NASA NTRS1972

Investigation of the plastic fracture of high strength steels

An investigation of the plastic fracture process to improve tensile strength in high strength steels is presented. Two generic types of steels are considered: a quenched and tempered grade and a maraging grade, in order to compare two different matrix microstructures. Each type of steel was studied in commercial grade purity and in special melted high purity form, low in residual and impurity elements. The specific alloys dealt with include AISI 4340 and 18 Ni, 200 grade maraging steel, both heat treated to the same yield strength level of approximately 200 ksi.

Cox, T. B.↗
NASA NTRS1977

Apollo 17 KREEPy basalt - A rock type intermediate between mare and KREEP basalts

The Apollo 17 KREEPy basalt is a unique lunar volcanic rock, observed only as clasts in the light friable breccia matrix (72275) of Boulder 1, Station 2 at Taurus-Littrow. Its status as a volcanic rock is confirmed by the absence of any meteoritic contamination, a lack of cognate inclusions or xenocrystal material, and low Ni contents in metal grains. The basalt was extruded 4.01 + or - 0.04 b.y. ago, approximately contemporaneously with the high-alumina mare basalts at Fra Mauro; shortly afterwards it was disrupted, probably by the Serenitatis impact, and its fragments emplaced in the South Massif. The basalt, which is quartz-normative and aluminous, is chemically and mineralogically intermediate between the Apollo 15 KREEP basalts and the high-alumina mare basalts in most respects. It consists mainly of plagioclase and pigeonitic pyroxene in approximately equal amounts, and 10-30% of mesostatis.

Ryder, G.↗