Engineering PapersSearch

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 91 records · Page 5

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.

G-DYN Multibody Dynamics Engine

G-DYN is a multi-body dynamic simulation software engine that automatically assembles and integrates equations of motion for arbitrarily connected multibody dynamic systems. The algorithm behind G-DYN is based on a primal-dual formulation of the dynamics that captures the position and velocity vectors (primal variables) of each body and the interaction forces (dual variables) between bodies, which are particularly useful for control and estimation analysis and synthesis. It also takes full advantage of the spare matrix structure resulting from the system dynamics to numerically integrate the equations of motion efficiently. Furthermore, the dynamic model for each body can easily be replaced without re-deriving the overall equations of motion, and the assembly of the equations of motion is done automatically. G-DYN proved an essential software tool in the simulation of spacecraft systems used for small celestial body surface sampling, specifically in simulating touch-and-go (TAG) maneuvers of a robotic sampling system from a comet and asteroid. It is used extensively in validating mission concepts for small body sample return, such as Comet Odyssey and Galahad New Frontiers proposals.

Acikmese, Behcet

Low-Power Architectures for Large Radio Astronomy Correlators

The architecture of a cross-correlator for a synthesis radio telescope with N greater than 1000 antennas is studied with the objective of minimizing power consumption. It is found that the optimum architecture minimizes memory operations, and this implies preference for a matrix structure over a pipeline structure and avoiding the use of memory banks as accumulation registers when sharing multiply-accumulators among baselines. A straw-man design for N = 2000 and bandwidth of 1 GHz, based on ASICs fabricated in a 90 nm CMOS process, is presented. The cross-correlator proper (excluding per-antenna processing) is estimated to consume less than 35 kW.

signal processing

Mapping Past, Present, and Future Climatic Suitability for Invasive Aedes Aegypti and Aedes Albopictus in the United States: A Process-Based Modeling Approach Using CMIP5 Downscaled Climate Scenarios

The ongoing spread of the mosquitoes, Aedes aegypti and Aedes albopictus, in the continental United States leaves new areas at risk for local transmission of dengue, chikungunya, and Zika viruses. All three viruses have caused major disease outbreaks in the Americas with infected travelers returning regularly to the U.S. The expanding range of these mosquitoes raises questions about whether recent spread has been enabled by climate change or other anthropogenic influences. In this analysis, we used downscaled climate scenarios from the NASA Earth Exchange Global Daily Downscaled Projections (NEX GDDP) dataset to model Ae. aegypti and Ae. albopictus population growth rates across the United States. We used a stage-structured matrix population model to understand past and present climatic suitability for these vectors, and to project future suitability under CMIP5 climate change scenarios. Our results indicate that much of the southern U.S. is suitable for both Ae. aegypti and Ae. albopictus year-round. In addition, a large proportion of the U.S. is seasonally suitable for mosquito population growth, creating the potential for periodic incursions into new areas. Changes in climatic suitability in recent decades for Ae. aegypti and Ae. albopictus have occurred already in many regions of the U.S., and model projections of future climate suggest that climate change will continue to reshape the range of Ae. aegypti and Ae. albopictus in the U.S., and potentially the risk of the viruses they transmit.

process-based modeling

Observations on the Computation of Eigenvalue and Eigenvector Jacobians

Many scientific and engineering problems benefit from analytic expressions for eigenvalue and eigenvector derivatives with respect to the elements of the parent matrix. While there exists extensive literature on the calculation of these derivatives, which take the form of Jacobian matrices, there are a variety of deficiencies that have yet to be addressed — including the need for both left and right eigenvectors, limitations on the matrix structure, and issues with complex eigenvalues and eigenvectors. This work addresses these deficiencies by proposing a new analytic solution for the eigenvalue and eigenvector derivatives. The resulting analytic Jacobian matrices are numerically efficient to compute and are valid for the general complex case. It is further shown that this new general result collapses to previously known relations for the special cases of real symmetric matrices and real diagonal matrices. Finally, the new Jacobian expressions are validated using forward finite differencing and performance is compared with another technique.

Jacobian

Space Weathering of Genesis Mission Solar-wind Collectors with Inferences for Weathering on Airless Bodies

Samples from the Genesis Solar Wind Sample Return (NASA Discovery 5) are a unique opportunity to study the initial stages of space weathering, i.e., the physical and chemical effects of solar-wind irradiation. Arrays of collectors containing multiple materials were each exposed to a different solar-wind regime (fast, slow, bulk, or coronal mass ejection) at the L1 point for long durations (years). Materials exposed to the solar wind included metals, semiconductors, and insulators. Although the time of exposure was obviously short relative to samples having extraterrestrial origins, optical properties, surface chemical properties, and matrix structure have changed in many collectors due to exposure to solar wind. The thickness of amorphous zones, where present, appears to correspond with the depth of the peak of the solar-wind H distribution in each regime. Damage from high-energy particles was negligible because the collectors were 700 μm or less in thickness and shielded from the back by the spacecraft. Micrometeorite impacts and sputtering were also negligible because of the short exposure times. Our current results are preliminary: We hope future workers will extend this study both to support Genesis characterization efforts and to further understand space-weathering processes on a geologic timescale.

space weathering

Structural Analysis and Matrix Interpretive System /SAMIS/

SAMIS digital computer program simplifies automated structural analysis and eliminates reprogramming for problem changes. Program objectives are achieved by standardizing, by providing a modular program, and by programming for intermediate-size problems.

Source record

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE

Cooled Ceramic Matrix Composite Propulsion Structures Demonstrated

NASA's Next Generation Launch Technology (NGLT) Program has successfully demonstrated cooled ceramic matrix composite (CMC) technology in a scramjet engine test. This demonstration represented the world s largest cooled nonmetallic matrix composite panel fabricated for a scramjet engine and the first cooled nonmetallic composite to be tested in a scramjet facility. Lightweight, high-temperature, actively cooled structures have been identified as a key technology for enabling reliable and low-cost space access. Tradeoff studies have shown this to be the case for a variety of launch platforms, including rockets and hypersonic cruise vehicles. Actively cooled carbon and CMC structures may meet high-performance goals at significantly lower weight, while improving safety by operating with a higher margin between the design temperature and material upper-use temperature. Studies have shown that using actively cooled CMCs can reduce the weight of the cooled flow-path component from 4.5 to 1.6 lb/sq ft and the weight of the propulsion system s cooled surface area by more than 50 percent. This weight savings enables advanced concepts, increased payload, and increased range. The ability of the cooled CMC flow-path components to operate over 1000 F hotter than the state-of-the-art metallic concept adds system design flexibility to space-access vehicle concepts. Other potential system-level benefits include smaller fuel pumps, lower part count, lower cost, and increased operating margin.

Jaskowiak, Martha H.

HITCAN for actively cooled hot-composite thermostructural analysis

A computer code, high temperature composite analyzer (HITCAN), was developed to analyze/design hot metal matrix composite structures. HITCAN is a general purpose code for predicting the global structural and local stress-strain response of multilayered (arbitrarily oriented) metal matrix structures both at the constituent (fiber, matrix, and interphase) and the structural level, including the fabrication process effects. The thermomechanical properties of the constituents are considered to be nonlinearly dependent on several parameters, including temperature, stress, and stress rate. The computational procedure employs an incremental iterative nonlinear approach utilizing a multifactor-interaction material behavior model, i.e., the material properties are expressed in terms of a product of several factors that affect the properties. HITCAN structural analysis capabilities (static, load stepping - a multistep static analysis with material properties updated at each step, modal, and buckling) for cooled hot structures are demonstrated through a specific example problem.

Chamis, C. C.

HITCAN for actively cooled hot-composite thermostructural analysis

A computer code, high temperature composite analyzer (HITCAN), was developed to analyze/design hot metal matrix composite structures. HITCAN is a general purpose code for predicting the global structural and local stress-strain response of multilayered (arbitrarily oriented) metal matrix structures both at the constituent (fiber, matrix, and interphase) and the structural level, including the fabrication process effects. The thermomechanical properties of the constituents are considered to be nonlinearly dependent on several parameters, including temperature, stress, and stress rate. The computational procedure employs an incremental iterative nonlinear approach utilizing a multifactor-interaction material behavior model, i.e., the material properties are expressed in terms of a product of several factors that affect the properties. HITCAN structural analysis capabilities (static, load stepping - a multistep static analysis with material properties updated at each step, modal, and buckling) for cooled hot structures are demonstrated through a specific example problem.

Chamis, C. C.

HITCAN: High temperature composite analyzer

A computer code, HITCAN (High Temperature Composite Analyzer) was developed to analyze/design metal matrix composite structures. HITCAN is based on composite mechanics theories and computer codes developed at NASA LeRC over the last two decades. HITCAN is a general purpose code for predicting the global structural and local stress-strain response of multilayered (arbitrarily oriented) metal matrix structures both at the constituent (fiber, matrix, and interphase) and the structure level and including the fabrication process effects. The thermomechanical properties of the constituents are considered to be nonlinearly dependent on several parameters including temperature, stress, and stress rate. The computational procedure employs an incremental iterative nonlinear approach utilizing a multifactor-interaction material behavior model. HITCAN features and analysis capabilities (static, load stepping, modal, and buckling) are demonstrated through typical example problems.

Singhal, Surendra N.