SEARCH · Engineering Papers

Results for “Matrix factorization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS2001

The Efficiency of Various Computers and Optimizations in Performing Finite Element Computations

With the advent of computers with many processors, it becomes unclear how to best exploit this advantage. For example, matrices can be inverted by applying several processors to each vector operation, or one processor can be applied to each matrix. The former approach has diminishing returns beyond a handful of processors, but how many processors depends on the computer architecture. Applying one processor to each matrix is feasible with enough ram memory and scratch disk space, but the speed at which this is done is found to vary by a factor of three depending on how it is done. The cost of the computer must also be taken into account. A computer with many processors and fast interprocessor communication is much more expensive than the same computer and processors with slow interprocessor communication. Consequently, for problems that require several matrices to be inverted, the best speed per dollar for computers is found to be several small workstations that are networked together, such as in a Beowulf cluster. Since these machines typically have two processors per node, each matrix is most efficiently inverted with no more than two processors assigned to it.

Marcus, Martin H.↗

NASA NTRS2014

Injury Surveillance Among NASA Astronauts Using the Barell Injury Diagnosis Matrix

Astronauts perform physically demanding tasks and risk incurring musculoskeletal injuries during both groundbased training and missions. Increased injury rates throughout the history of the U.S. space program have been attributed to numerous factors, including an aging astronaut corps, increased Weightless Environment Training Facility (WETF) and Neutral Buoyancy Laboratory (NBL) training to construct the International Space Station, and improved clinical operations that promote injury prevention and reporting. With NASA program changes through the years (including retirement of the Shuttle program) and an improved training environment (including a new astronaut gym), there is no surveillance program to systematically track injury rates. A limited number of research projects have been conducted over the past 20 years to evaluate musculoskeletal injuries: (1) to evaluate orthopedic injuries from 1987 to 1995, (2) to describe upper extremity injuries, (3) to evaluate EVA spacesuit training related injuries, and (4) to evaluate in-flight musculoskeletal injuries. Nevertheless, there has been no consistently performed comprehensive assessment of musculoskeletal injuries among astronauts. The Barell Injury Diagnosis Matrix was introduced at the 2001 meeting of the International Collaborative Effort (ICE) on Injury Statistics. The Matrix proposes a standardized method of classifying body region by nature of injury. Diagnoses are coded using the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) coding system. The purpose of this study is to assess the usefulness and complexity of the Barell Injury Diagnosis Matrix to classify and track musculoskeletal injuries among NASA astronauts.

Murray, J. D.↗

NASA NTRS1972

Study of fuel cell thermal control systems for advanced missions.

This study evaluated many heat rejection and thermal control concepts which could be applied to fuel cells for long term (600 hours) orbital and lunar surface missions. The concepts considered several types of radiators which utilized pumped gas, liquid and two phase working fluids and incorporated solid conduction fins as well as heat pipe (vapor chamber) fins. The comparison of the concepts was based on weight, area and other factors such as standby power, ability to accommodate heat load variation, control complexity, and meteoroid survival capability. A design selection matrix was established and an optimum (primary) and an alternate (secondary) heat rejection concept was chosen. Heat rejection techniques utilizing self-controlled heat pipe radiators dominate the results.

Caputo, R. S.↗

NASA NTRS1991

Parallel computation in a three-dimensional elastic-plastic finite-element analysis

CRAY 2 and CRAY Y-MP evaluation runs are undertaken for a three-dimensional, elastoplastic FEM code implementation of the 'autotasking' parallel-processing technique, which is used in all code components except the matrix-equation solver. For a typical example problem in the case of a four-processor CRAY 2, a speedup factor of 2.1 was found achievable in a dedicated environment, and 1.7 in a multiuser environment; for a three-processor CRAY Y-MP, speedup in a multiuser environment was about 2.4. Only minimal effort was required to implement autotasking in this implementation code.

Shivakumar, K. N.↗

NASA NTRS1991

A comparison of experimental observations and numerical predictions for the initiation of fiber microbuckling in notched composite laminates

The results of the present experimental study of fiber microbuckling indicates that free surfaces, ply orientations, and resin-rich regions affect the strain level for the initiation of in-plane fiber microbuckling. Nonlinear matrix constitutive behavior, in conjunction with initial fiber waviness, cause significant reductions of fiber microbuckling initiation strain level; these two factors generate nonlinear global stress-strain responses in the models. The fiber microbuckling initiation strains correlate well with the applied strain associated with asymptotic local shear strains.

Guynn, E. G.↗

NASA NTRS1992

The 3D modeling of high numerical aperture imaging in thin films

A modelling technique is described which is used to explore three dimensional (3D) image irradiance distributions formed by high numerical aperture (NA is greater than 0.5) lenses in homogeneous, linear films. This work uses a 3D modelling approach that is based on a plane-wave decomposition in the exit pupil. Each plane wave component is weighted by factors due to polarization, aberration, and input amplitude and phase terms. This is combined with a modified thin-film matrix technique to derive the total field amplitude at each point in a film by a coherent vector sum over all plane waves. Then the total irradiance is calculated. The model is used to show how asymmetries present in the polarized image change with the influence of a thin film through varying degrees of focus.

Flagello, D. G.↗

NASA NTRS1995

Krylov methods for compressible flows

We investigate the application of Krylov methods to compressible flows, and the effect of implicit boundary conditions on the implicit solution of nonlinear problems. Two defect-correction procedures, namely, approximate factorization (AF) for structured grids and ILU/GMRES for general grids, are considered. Also considered here are Newton-Krylov matrix-free methods that we combined with the use of mixed discretization schemes in the implicitly defined Jacobian and its preconditioner. Numerical experiments that show the performance of our approaches are then presented.

Tidriri, M. D.↗

NASA NTRS1993

A spectral algorithm for envelope reduction of sparse matrices

A new algorithm for reducing the envelope of a sparse matrix is presented. This algorithm is based on the computation of eigenvectors of the Laplacian matrix associated with the graph of the sparse matrix. A reordering of the sparse matrix is determined based on the numerical values of the entries of an eigenvector of the Laplacian matrix. Numerical results show that the new reordering algorithm can in some cases reduce the envelope by more than a factor of two over the current standard algorithms such as Gibbs-Poole-Stockmeyer (GPS) or SPARSPAK's reverse Cuthill-McKee (RCM).

Barnard, Stephen T.↗

NASA NTRS1998

Numerical Procedures for Inlet/Diffuser/Nozzle Flows

Two primitive variable, pressure based, flux-split, RNS/NS solution procedures for viscous flows are presented. Both methods are uniformly valid across the full Mach number range, Le., from the incompressible limit to high supersonic speeds. The first method is an 'optimized' version of a previously developed global pressure relaxation RNS procedure. Considerable reduction in the number of relatively expensive matrix inversion, and thereby in the computational time, has been achieved with this procedure. CPU times are reduced by a factor of 15 for predominantly elliptic flows (incompressible and low subsonic). The second method is a time-marching, 'linearized' convection RNS/NS procedure. The key to the efficiency of this procedure is the reduction to a single LU inversion at the inflow cross-plane. The remainder of the algorithm simply requires back-substitution with this LU and the corresponding residual vector at any cross-plane location. This method is not time-consistent, but has a convective-type CFL stability limitation. Both formulations are robust and provide accurate solutions for a variety of internal viscous flows to be provided herein.

Rubin, Stanley G.↗

NASA NTRS2002

The Frequency Detuning Correction and the Asymmetry of Line Shapes: The Far Wings of H2O-H2O

A far-wing line shape theory which satisfies the detailed balance principle is applied to the H2O-H2O system. Within this formalism, two line shapes are introduced, corresponding to band-averages over the positive and negative resonance lines, respectively. Using the coordinate representation, the two line shapes can be obtained by evaluating 11-dimensional integrations whose integrands are a product of two factors. One depends on the interaction between the two molecules and is easy to evaluate. The other contains the density matrix of the system and is expressed as a product of two 3-dimensional distributions associated with the density matrices of the absorber and the perturber molecule, respectively. If most of the populated states are included in the averaging process, to obtain these distributions requires extensive computer CPU time, but only have to be computed once for a given temperature. The 11-dimensional integrations are evaluated using the Monte Carlo method, and in order to reduce the variance, the integration variables are chosen such that the sensitivity of the integrands on them is clearly distinguished.

Ma, Q.↗

NASA NTRS2006

Simulated Data for High Temperature Composite Design

The paper describes an effective formal method that can be used to simulate design properties for composites that is inclusive of all the effects that influence those properties. This effective simulation method is integrated computer codes that include composite micromechanics, composite macromechanics, laminate theory, structural analysis, and multi-factor interaction model. Demonstration of the method includes sample examples for static, thermal, and fracture reliability for a unidirectional metal matrix composite as well as rupture strength and fatigue strength for a high temperature super alloy. Typical results obtained for a unidirectional composite show that the thermal properties are more sensitive to internal local damage, the longitudinal properties degrade slowly with temperature, the transverse and shear properties degrade rapidly with temperature as do rupture strength and fatigue strength for super alloys.

Chamis, Christos C.↗

NASA NTRS1995

Collagen and Stretch Modulate Autocrine Secretion of Insulin-like Growth Factor-1 and Insulin-like Growth Factor Binding Proteins from Differentiated Skeletal Muscle Cells

Stretch-induced skeletal muscle growth may involve increased autocrine secretion of insulin-like growth factor-1 (IGF-1) since IGF-1 is a potent growth factor for skeletal muscle hypertrophy, and stretch elevates IGF-1 mRNA levels in vivo. In tissue cultures of differentiated avian pectoralis skeletal muscle cells, nanomolar concentrations of exogenous IGF-1 stimulated growth in mechanically stretched but not static cultures. These cultures released up to 100 pg of endogenously produced IGF-1/micro-g of protein/day, as well as three major IGF binding proteins of 31, 36, and 43 kilodaltons (kDa). IGF-1 was secreted from both myofibers and fibroblasts coexisting in the muscle cultures. Repetitive stretch/relaxation of the differentiated skeletal muscle cells stimulated the acute release of IGF-1 during the first 4 h after initiating mechanical activity, but caused no increase in the long-term secretion over 24-72 h of IGF-1, or its binding proteins. Varying the intensity and frequency of stretch had no effect on the long-term efflux of IGF-1. In contrast to stretch, embedding the differentiated muscle cells in a three-dimensional collagen (Type I) matrix resulted in a 2-5-fold increase in long-term IGF-1 efflux over 24-72 h. Collagen also caused a 2-5-fold increase in the release of the IGF binding proteins. Thus, both the extracellular matrix protein type I collagen and stretch stimulate the autocrine secretion of IGF-1, but with different time kinetics. This endogenously produced growth factor may be important for the growth response of skeletal myofibers to both types of external stimuli.

Perrone, Carmen E.↗

NASA NTRS1992

Damage development in titanium metal matrix composites subjected to cyclic loading

Several layups of SCS-6/Ti-15-3 composites were investigated. Fatigue tests were conducted and analyzed for both notched and unnotched specimens at room temperature and elevated temperatures. Thermo-mechanical fatigue results were analyzed. Test results indicated that the stress in the 0 degree fibers is the controlling factor in fatigue life. The static and fatigue strength of these materials is shown to be strongly dependent on the level of residual stresses and the fiber/matrix interfacial strength. Fatigue tests of notched specimens showed that cracks can initiate and grow many fiber spacings in the matrix materials without breaking fibers. Fiber bridging models were applied to characterize the crack growth behavior. The matrix cracks are shown to significantly reduce the residual strength of notched composites. The notch strength of these composites was accurately predicted using a micromechanics based methodology.

Johnson, W. S.↗

NASA NTRS1993

Damage development in titanium metal-matrix composites subjected to cyclic loading

Johnson, W. S.↗

NASA NTRS1990

A New Approximate Fracture Mechanics Analysis Methodology for Composites with a Crack or Hole

A new approximate theory which links the inherent flaw concept with the theory of crack tip stress singularities at a bi-material interface was developed. Three assumptions were made: (1) the existence of inherent flaw (i.e., damage zone) at the tip of the crack, (2) a fracture of the filamentary composites initiates at a crack lying in the matrix material at the interface of the matrix/filament, and (3) the laminate fails whenever the principal load-carrying laminae fails. This third assumption implies that for a laminate consisting of 0 degree plies, cracks into matrix perpendicular to the 0 degree filaments are the triggering mechanism for the final failure. Based on this theory, a parameter bar K sub Q which is similar to the stress intensity factor for isotropic materials but with a different dimension was defined. Utilizing existing test data, it was found that bar K sub Q can be treated as a material constant. Based on this finding a fracture mechanics analysis methodology was developed. The analytical results are correlated well with test results. This new approximate theory can apply to both brittle and metal matrix composite laminates with crack or hole.

Tsai, H. C.↗

NASA NTRS1976

Monotectic and syntectic alloys experiment MA-044

The flight test furnace run of the monotectic and syntectic alloys experiment was successfully performed during the Apollo Soyuz Test Project mission. All three experiment cartridges were returned intact. A second ground based test was conducted after the mission to provide one g samples with an equivalent thermal history. The examinations of all flight and ground based cartridges and ampoules have been completed. Encouraging results have been obtained in the preliminary analysis of one flight aluminum antimonide sample as compared to ground based and prototype test counterparts. In terms of reduction of the unwanted second phase, the low g sample shows a homogeneity improvement by a factor of 5 to 20. Preliminary examination of one lead zinc flight sample shows a dispersion of lead in the zinc matrix. However, complete homogenization of the material did not occur. A thermal analysis of flight data shows that the average soak temperature for the head zinc ampoules was 20 C lower than planned but 35 C above the published consolute temperature..

Ang, C. Y.↗

NASA NTRS1976

A numerical comparison of discrete Kalman filtering algorithms - An orbit determination case study

An improved Kalman filter algorithm based on a modified Givens matrix triangularization technique is proposed for solving a nonstationary discrete-time linear filtering problem. The proposed U-D covariance factorization filter uses orthogonal transformation technique; measurement and time updating of the U-D factors involve separate application of Gentleman's fast square-root-free Givens rotations. Numerical stability and accuracy of the algorithm are compared with those of the conventional and stabilized Kalman filters and the Potter-Schmidt square-root filter, by applying these techniques to a realistic planetary navigation problem (orbit determination for the Saturn approach phase of the Mariner Jupiter-Saturn Mission, 1977). The new algorithm is shown to combine the numerical precision of square root filtering with the efficiency of the original Kalman algorithm.

Thornton, C. L.↗

NASA NTRS1982

Refraction at a curved dielectric interface - Geometrical optics solution

The transmission of a spherical or plane wave through an arbitrarily curved dielectric interface is solved by the geometrical optics theory. The transmitted field is proportional to the product of the conventional Fresnel's transmission coefficient and a divergence factor (DF), which describes the cross-sectional variation (convergence or divergence) of a ray pencil as the latter propagates in the transmitted region. The factor DF depends on the incident wavefront, the curvatures of the interface, and the relative indices of the two media. Explicit matrix formulas for calculating DF are given, and its physical significance is illustrated via examples.

Lee, S.-W.↗