Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix structures”

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

Investigation of Kevlar fabric-based materials for use with inflatable structures

Design, manufacture and testing of laminated and coated composite materials incorporating a structural matrix of Kevlar are reported. The practicality of using Kevlar in aerostat materials is demonstrated, and data are provided on practical weaves, lamination and coating particulars, rigidity, strength, weight, elastic coefficients, abrasion resistance, crease effects, peel strength, blocking tendencies, helium permeability, and fabrication techniques. Properties of the Kevlar-based materials are compared with conventional Dacron-reinforced counterparts. A comprehensive test and qualification program is discussed, and considerable quantitative biaxial tensile and shear test data are provided.

Niccum, R. J.↗

Variational calculations of vibrational properties of ozone

A variational method is used to obtain vibrational-rotational properties for ozone from an experimental quartic force field. Band positions, average structures, matrix elements for calculating infrared intensities, and effective rotational constants are reported for (O-16)3 and its O-18 isotopic species. Also, the degree to which the vibrational energies and properties are converged is investigated as a function of the basis set parameters and basis set size, and of the method of obtaining the Hamiltonian matrix elements. A convenient procedure for assigning the vibrational states is developed for variational vibrational wavefunction expansions.

Carney, G. D.↗

A high order accurate finite element algorithm for high Reynolds number flow prediction

A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.

Baker, A. J.↗

Shear rupture of a directionally solidified eutectic gamma/gamma prime - alpha (Mo) alloy

Directionally solidified Mo alloys are evaluated to determine the shear rupture strength and to possibly improve it by microstructural and heat treatment variations. Bars of the alloy containing nominally 5.7% Al and 33.5% Mo by weight with balance Ni were directionally solidified at rates between 10 and 100 mm per hour in furnaces with thermal gradients at the liquid-solid interface of 250 or 100 C per cm. A limited number of longitudinal shear rupture tests were conducted at 760 C and 207 MPa in the as - solidified and in several heat treated conditions. It is shown that shear rupture failures are partly transgranular and that resistance to failure is prompted by good fiber alignment and a matrix structure consisting mainly of gamma prime. Well aligned as - solidified specimens sustained the shear stress for an average of 81 hours. A simulated coating heat treatment appeared to increase the transformation of gamma to gamma prime and raised the average shear life of aligned specimens to 111 hours. However, heat treatments at 1245 C and especially at 1190 C appeared to be detrimental by causing partial solutioning of the gamma prime, and reducing lives to 47 and 10 hours, respectively.

Harf, F. H.↗

On the accuracy and convergence of implicit numerical integration of finite element generated ordinary differential equations

A study of accuracy and convergence of linear functional finite element solution to linear parabolic and hyperbolic partial differential equations is presented. A variable-implicit integration procedure is employed for the resultant system of ordinary differential equations. Accuracy and convergence is compared for the consistent and two lumped assembly procedures for the identified initial-value matrix structure. Truncation error estimation is accomplished using Richardson extrapolation.

Baker, A. J.↗

Utility of a finite element solution algorithm for initial-value problems

The Galerkin criterion within a finite element Weighted Residuals formulation is employed to establish an implicit solution algorithm for an initial-value partial differential equation. Numerical solutions of a transient parabolic and a hyperbolic equation, obtained using linear, quadratic and two cubic finite element basis functions, are employed to quantize accuracy and confirm and refine theoretical convergence rate estimates. The linear basis algorithm for the hyperbolic equation displays excellent accuracy on a coarse computational grid and a high-order convergence rate with discretization refinement. Good accuracy and a strong convergence rate in surface flux are determined for a nonhomogeneous Neumann boundary constraint applied to a parabolic equation. The results amply demonstrate the impact of the nondiagonal finite element initial-value matrix structure on solution accuracy and/or convergence rate.

Baker, A. J.↗

Advanced technology lightweight fuel cell program

The potential of the alkaline electrolyte fuel cell as the power source in a multi hundred kilowatt orbital energy storage system was studied. The total system weight of an electrolysis cell energy storage system was determined. The tests demonstrated: (1) the performance stability of a platinum on carbon anode catalyst configuration after 5000 hours of testing has no loss in performance; (2) capability of the alkaline fuel cell to operate to a cyclical load profile; (3) suitability of a lightweight graphite electrolyte reservoir plate for use in the alkaline fuel cell; (4) long life potential of a hybrid polysulfone cell edge frame construction; and (5) long term stability of a fiber reinforced potassium titanate matrix structure. The power section tested operates with passive water removal eliminating the requirement for a dynamic hydrogen pump water separator thereby allowing a powerplant design with reduced weight, lower parasite power, and a potential for high reliability and extended endurance. It is concluded that two perovskites are unsuitable for use as a catalyst or as a catalyst support at the cathode of an alkaline fuel cell.

Martin, R. E.↗

Higher order methods for convection-diffusion problems

This paper applies C1 cubic Hermite polynomials embedded in an orthogonal collocation scheme to the spatial discretization of the unsteady nonlinear Burgers equation as a model of the equations of fluid mechanics. The temporal discretization is carried out by means of either a noniterative finite difference or an iterative finite difference procedure. Results of this method are compared with those of a second-order finite difference scheme and a splined-cubic Taylor's series scheme. Stability limits are derived and the matrix structure of the several schemes are compared.

Murphy, J. D.↗

Estimation of geopotential from satellite-to-satellite range rate data: Numerical results

A technique for high-resolution geopotential field estimation by recovering the harmonic coefficients from satellite-to-satellite range rate data is presented and tested against both a controlled analytical simulation of a one-day satellite mission (maximum degree and order 8) and then against a Cowell method simulation of a 32-day mission (maximum degree and order 180). Innovations include: (1) a new frequency-domain observation equation based on kinetic energy perturbations which avoids much of the complication of the usual Keplerian element perturbation approaches; (2) a new method for computing the normalized inclination functions which unlike previous methods is both efficient and numerically stable even for large harmonic degrees and orders; (3) the application of a mass storage FFT to the entire mission range rate history; (4) the exploitation of newly discovered symmetries in the block diagonal observation matrix which reduce each block to the product of (a) a real diagonal matrix factor, (b) a real trapezoidal factor with half the number of rows as before, and (c) a complex diagonal factor; (5) a block-by-block least-squares solution of the observation equation by means of a custom-designed Givens orthogonal rotation method which is both numerically stable and tailored to the trapezoidal matrix structure for fast execution.

Thobe, Glenn E.↗

Improved performance in NASTRAN (R)

Three areas of improvement in COSMIC/NASTRAN, 1989 release, were incorporated recently that make the analysis program run faster on large problems. Actual log files and actual timings on a few test samples that were run on IBM, CDC, VAX, and CRAY computers were compiled. The speed improvement is proportional to the problem size and number of continuation cards. Vectorizing certain operations in BANDIT, makes BANDIT run twice as fast in some large problems using structural elements with many node points. BANDIT is a built-in NASTRAN processor that optimizes the structural matrix bandwidth. The VAX matrix packing routine BLDPK was modified so that it is now packing a column of a matrix 3 to 9 times faster. The denser and bigger the matrix, the greater is the speed improvement. This improvement makes a host of routines and modules that involve matrix operation run significantly faster, and saves disc space for dense matrices. A UNIX version, converted from 1988 COSMIC/NASTRAN, was tested successfully on a Silicon Graphics computer using the UNIX V Operating System, with Berkeley 4.3 Extensions. The Utility Modules INPUTT5 and OUTPUT5 were expanded to handle table data, as well as matrices. Both INPUTT5 and OUTPUT5 are general input/output modules that read and write FORTRAN files with or without format. More user informative messages are echoed from PARAMR, PARAMD, and SCALAR modules to ensure proper data values and data types being handled. Two new Utility Modules, GINOFILE and DATABASE, were written for the 1989 release. Seven rigid elements are added to COSMIC/NASTRAN. They are: CRROD, CRBAR, CRTRPLT, CRBE1, CRBE2, CRBE3, and CRSPLINE.

Chan, Gordon C.↗

How well can gravity be recovered using Topex and GPS data?

When Topex is launched in mid-1992 it will carry a high quality GPS receiver which will operate in concert with a worldwide network of precision GPS ground receivers. The data from these receivers can be used to recover new information about the earth's gravity field at longer wavelengths. Software and algorithms have been developed which will allow this gravity field information to be recovered with much greater efficiency than with traditional techniques. The basis for these algorithms is the gravity bin formulation and related filtering techniques that exploit the repeat orbit of Topex and the sparse matrix structure of the problem. This new software has been used to evaluate the expected improvement in the gravity field using multiple ten-day arcs of GPS data from Topex.

Bertiger, Willy↗

Constraint treatment techniques and parallel algorithms for multibody dynamic analysis

Computational procedures for kinematic and dynamic analysis of three-dimensional multibody dynamic (MBD) systems are developed from the differential-algebraic equations (DAE's) viewpoint. Constraint violations during the time integration process are minimized and penalty constraint stabilization techniques and partitioning schemes are developed. The governing equations of motion, a two-stage staggered explicit-implicit numerical algorithm, are treated which takes advantage of a partitioned solution procedure. A robust and parallelizable integration algorithm is developed. This algorithm uses a two-stage staggered central difference algorithm to integrate the translational coordinates and the angular velocities. The angular orientations of bodies in MBD systems are then obtained by using an implicit algorithm via the kinematic relationship between Euler parameters and angular velocities. It is shown that the combination of the present solution procedures yields a computationally more accurate solution. To speed up the computational procedures, parallel implementation of the present constraint treatment techniques, the two-stage staggered explicit-implicit numerical algorithm was efficiently carried out. The DAE's and the constraint treatment techniques were transformed into arrowhead matrices to which Schur complement form was derived. By fully exploiting the sparse matrix structural analysis techniques, a parallel preconditioned conjugate gradient numerical algorithm is used to solve the systems equations written in Schur complement form. A software testbed was designed and implemented in both sequential and parallel computers. This testbed was used to demonstrate the robustness and efficiency of the constraint treatment techniques, the accuracy of the two-stage staggered explicit-implicit numerical algorithm, and the speed up of the Schur-complement-based parallel preconditioned conjugate gradient algorithm on a parallel computer.

Chiou, Jin-Chern↗

Highly parallel sparse Cholesky factorization

Several fine grained parallel algorithms were developed and compared to compute the Cholesky factorization of a sparse matrix. The experimental implementations are on the Connection Machine, a distributed memory SIMD machine whose programming model conceptually supplies one processor per data element. In contrast to special purpose algorithms in which the matrix structure conforms to the connection structure of the machine, the focus is on matrices with arbitrary sparsity structure. The most promising algorithm is one whose inner loop performs several dense factorizations simultaneously on a 2-D grid of processors. Virtually any massively parallel dense factorization algorithm can be used as the key subroutine. The sparse code attains execution rates comparable to those of the dense subroutine. Although at present architectural limitations prevent the dense factorization from realizing its potential efficiency, it is concluded that a regular data parallel architecture can be used efficiently to solve arbitrarily structured sparse problems. A performance model is also presented and it is used to analyze the algorithms.

Gilbert, John R.↗

A Capability to Generate Physics-based Mass Estimating Relationships for Conceptual Space Vehicle Design

This paper is written in support of the on-going research into conceptual space vehicle design conducted at the Space Systems Design Laboratory (SSDL) at the Georgia Institute of Technology. Research at the SSDL follows a sequence of a number of the traditional aerospace disciplines. The sequence of disciplines and interrelationship among them is shown in the Design Structure Matrix (DSM). The discipline of Weights and Sizing occupies a central location in the design of a new space vehicle. Weights and Sizing interact, either in a feed forward or feed back manner, with every other discipline in the DSM. Because of this principle location, accuracy in Weights and Sizing is integral to producing an accurate model of a space vehicle concept. Instead of using conceptual level techniques, a simplified Finite Element Analysis (FEA) technique is described as applied to the problem of the Liquid Oxygen (LOX) tank bending loads applied to the forward Liquid Hydrogen (LH2) tank of the Georgia Tech Air Breathing Launch Vehicle (ABLV).

Olds, John R.↗

Mechanotransduction and the functional response of bone to mechanical strain

Mechanotransduction plays a crucial role in the physiology of many tissues including bone. Mechanical loading can inhibit bone resorption and increase bone formation in vivo. In bone, the process of mechanotransduction can be divided into four distinct steps: (1) mechanocoupling, (2) biochemical coupling, (3) transmission of signal, and (4) effector cell response. In mechanocoupling, mechanical loads in vivo cause deformations in bone that stretch bone cells within and lining the bone matrix and create fluid movement within the canaliculae of bone. Dynamic loading, which is associated with extracellular fluid flow and the creation of streaming potentials within bone, is most effective for stimulating new bone formation in vivo. Bone cells in vitro are stimulated to produce second messengers when exposed to fluid flow or mechanical stretch. In biochemical coupling, the possible mechanisms for the coupling of cell-level mechanical signals into intracellular biochemical signals include force transduction through the integrin-cytoskeleton-nuclear matrix structure, stretch-activated cation channels within the cell membrane, G protein-dependent pathways, and linkage between the cytoskeleton and the phospholipase C or phospholipase A pathways. The tight interaction of each of these pathways would suggest that the entire cell is a mechanosensor and there are many different pathways available for the transduction of a mechanical signal. In the transmission of signal, osteoblasts, osteocytes, and bone lining cells may act as sensors of mechanical signals and may communicate the signal through cell processes connected by gap junctions. These cells also produce paracrine factors that may signal osteoprogenitors to differentiate into osteoblasts and attach to the bone surface. Insulin-like growth factors and prostaglandins are possible candidates for intermediaries in signal transduction. In the effector cell response, the effects of mechanical loading are dependent upon the magnitude, duration, and rate of the applied load. Longer duration, lower amplitude loading has the same effect on bone formation as loads with short duration and high amplitude. Loading must be cyclic to stimulate new bone formation. Aging greatly reduces the osteogenic effects of mechanical loading in vivo. Also, some hormones may interact with local mechanical signals to change the sensitivity of the sensor or effector cells to mechanical load.

Non-NASA Center↗

A fiber matrix model for fluid flow and streaming potentials in the canaliculi of an osteon

A theoretical model is developed to predict the fluid shear stress and streaming potential at the surface of osteocytic processes in the lacunar-canalicular porosity of an osteon when the osteon is subject to mechanical loads that are parallel or perpendicular to its axis. The theory developed in Weinbaum et al. (31) for the flow through a proteoglycan matrix in a canaliculus is employed in a poroelastic model for the osteon. Our formulation is a generalization of that of Petrov et al. (17). Our model predicts that, in order to satisfy the measured frequency dependence of the phase and magnitude of the SGP in macroscopic bone samples, the fiber spacing in the fluid annulus must lie in the narrow range 6-7 nm typical of the spacing of GAG sidechains along a protein monomer. The model predictions for the local SGP profiles in the osteon agree with the experimental observations of Starkebaum et al. (24). The theory predicts that the pore pressure relaxation time, tau d, for a 150-300 microns diameter osteon with the foregoing matrix structure is approximately 0.03-0.13 sec, and that the amplitude of the mean fluid shear stress on the membrane of the osteocytic process at the mean areal radius of the osteon has a maximum at 28 Hz if tau d = 0.06 sec. This maximum, which is independent of the magnitude of the loading, could be important in vivo since the recent experiments of Turner et al. (28) and McLeod et al. (15) have a peak in the strain frequency spectrum between 20 and 30 Hz that also appears to be independent of the type (magnitude) of loading. Numerical predictions for the amplitude of the average fluid shear stress on the osteocytic membrane at the mean areal radius of the osteon show that the fluid shear stress associated with the low amplitude 20-30 Hz spectral strain component is at least as large as the average fluid shear stress associated with the high amplitude 1 Hz stride component, although the latter loading is an order of magnitude larger, and has a magnitude that lies within the middle of the range, 6-30 dynes/cm2, where fluid shear stresses in tissue culture studies with osteoblast monolayers have elicited an intracellular Ca++ response (31). The implications of these results for intracellular electrical communication are discussed.

Non-NASA Center↗

Silent Aircraft Initiative Concept Risk Assessment

A risk assessment of the Silent Aircraft Initiative's SAX-40 concept design for extremely low noise has been performed. A NASA team developed a list of 27 risk items, and evaluated the level of risk for each item in terms of the likelihood that the risk would occur and the consequences of the occurrence. The following risk items were identified as high risk, meaning that the combination of likelihood and consequence put them into the top one-fourth of the risk matrix: structures and weight prediction; boundary-layer ingestion (BLI) and inlet design; variable-area exhaust and thrust vectoring; displaced-threshold and continuous descent approach (CDA) operational concepts; cost; human factors; and overall noise performance. Several advanced-technology baseline concepts were created to serve as a basis for comparison to the SAX-40 concept. These comparisons indicate that the SAX-40 would have significantly greater research, development, test, and engineering (RDT&E) and production costs than a conventional aircraft with similar technology levels. Therefore, the cost of obtaining the extremely low noise capability that has been estimated for the SAX-40 is significant. The SAX-40 concept design proved successful in focusing attention toward low noise technologies and in raising public awareness of the issue.

Nickol, Craig L.↗

Low Cost Nuclear Thermal Rocket Cermet Fuel Element Environment Testing

Deep space missions with large payloads require high specific impulse and relatively high thrust to achieve mission goals in reasonable time frames.1,2 Conventional storable propellants produce average specific impulse. Nuclear thermal rockets capable of producing high specific impulse are proposed. Nuclear thermal rockets employ heat produced by fission reaction to heat and therefore accelerate hydrogen, which is then forced through a rocket nozzle providing thrust. Fuel element temperatures are very high (up to 3000 K), and hydrogen is highly reactive with most materials at high temperatures. Data covering the effects of high-temperature hydrogen exposure on fuel elements are limited.3 The primary concern is the mechanical failure of fuel elements that employ high-melting-point metals, ceramics, or a combination (cermet) as a structural matrix into which the nuclear fuel is distributed. The purpose of the testing is to obtain data to assess the properties of the non-nuclear support materials, as-fabricated, and determine their ability to survive and maintain thermal performance in a prototypical NTR reactor environment of exposure to hydrogen at very high temperatures. The fission process of the planned fissile material and the resulting heating performance is well known and does not therefore require that active fissile material be integrated in this testing. A small-scale test bed designed to heat fuel element samples via non-contact radio frequency heating and expose samples to hydrogen is being developed to assist in optimal material and manufacturing process selection without employing fissile material. This paper details the test bed design and results of testing conducted to date.

Bradley, D. E.↗