Engineering Papers⌕ Search

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.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Assessment of space environment induced microdamage in toughened composite materials

The effects of simulated space environments on the microdamage in a series of commercially available toughened matrix composite systems was determined. Low-earth orbit (LEO) exposures were simulated by thermal cycling; geosynchronous orbit (GEO) exposures were simulated by electron irradiation plus thermal cycling. Material response was characterized by assessing the induced microcracking and its influence on chemical and mechanical property changes. All materials, including several advanced, tough thermoplastics microcracked when exposed to the simulated LEO environment except a 177 C cured single phase toughened epoxy composite. The GEO simulated environment produced microdamage in all materials. The results suggest that increased matrix toughness may not be the overriding factor leading to improved durability in the space environment.

Sykes, George F.↗
NASA NTRS1987

Thomas precession, Wigner rotations and gauge transformations

The exact Lorentz kinematics of the Thomas precession is discussed in terms of Wigner's O(3)-like little group which describes rotations in the Lorentz frame in which the particle is at rest. A Lorentz-covariant form for the Thomas factor is derived. It is shown that this factor is a Lorentz-boosted rotation matrix, which becomes a gauge transformation in the infinite-momentum or zero-mass limit.

Han, D.↗
NASA NTRS1988

Modeling, model simplification and stability robustness with aeroelastic vehicles

The linearization and simplification of a nonlinear, literal model for flexible aircraft is highlighted. Areas of model fidelity that are critical if the model is to be used for control system synthesis are developed, and several simplification techniques that can deliver the necessary model fidelity are discussed. These techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle mass properties, vibrational characteristics, and rigid-body and aeroelastic stability derivatives, thus leading to the underlying causes for critical dynamic characteristics.

Schmidt, David K.↗
NASA NTRS1989

Flexible aircraft dynamic modeling for dynamic analysis and control synthesis

The linearization and simplification of a nonlinear, literal model for flexible aircraft is highlighted. Areas of model fidelity that are critical if the model is to be used for control system synthesis are developed and several simplification techniques that can deliver the necessary model fidelity are discussed. These techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle mass properties, vibrational characteristics, and rigid-body and aeroelastic stability derivatives, thus leading to the underlying causes for critical dynamic characteristics.

Schmidt, David K.↗
NASA NTRS1989

Efficient computer algebra algorithms for polynomial matrices in control design

The theory of polynomial matrices plays a key role in the design and analysis of multi-input multi-output control and communications systems using frequency domain methods. Examples include coprime factorizations of transfer functions, cannonical realizations from matrix fraction descriptions, and the transfer function design of feedback compensators. Typically, such problems abstract in a natural way to the need to solve systems of Diophantine equations or systems of linear equations over polynomials. These and other problems involving polynomial matrices can in turn be reduced to polynomial matrix triangularization procedures, a result which is not surprising given the importance of matrix triangularization techniques in numerical linear algebra. Matrices with entries from a field and Gaussian elimination play a fundamental role in understanding the triangularization process. In the case of polynomial matrices, matrices with entries from a ring for which Gaussian elimination is not defined and triangularization is accomplished by what is quite properly called Euclidean elimination. Unfortunately, the numerical stability and sensitivity issues which accompany floating point approaches to Euclidean elimination are not very well understood. New algorithms are presented which circumvent entirely such numerical issues through the use of exact, symbolic methods in computer algebra. The use of such error-free algorithms guarantees that the results are accurate to within the precision of the model data--the best that can be hoped for. Care must be taken in the design of such algorithms due to the phenomenon of intermediate expressions swell.

Baras, J. S.↗
NASA NTRS1991

Numerical and literal aeroelastic-vehicle-model reduction for feedback control synthesis

The simplification of a high-order, literal model for large flexible aircraft is discussed. Areas of model fidelity that are critical if the model is to be used for control law synthesis are presented. Several simplification techniques, some new and some widely available, that can deliver the necessary model fidelity are presented and applied to a model from the literature. The techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory, is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle vibrational characteristics and rigid-body and aeroelastic stability derivatives, thus providing insight into the underlying causes for critical dynamic characteristics.

Newman, Brett↗
NASA NTRS1995

Parallel CFD Algorithms for Aerodynamical Flow Solvers on Unstructured Meshes

The Advisory Group for Aerospace Research and Development (AGARD) has requested my participation in the lecture series entitled Parallel Computing in Computational Fluid Dynamics to be held at the von Karman Institute in Brussels, Belgium on May 15-19, 1995. In addition, a request has been made from the US Coordinator for AGARD at the Pentagon for NASA Ames to hold a repetition of the lecture series on October 16-20, 1995. I have been asked to be a local coordinator for the Ames event. All AGARD lecture series events have attendance limited to NATO allied countries. A brief of the lecture series is provided in the attached enclosure. Specifically, I have been asked to give two lectures of approximately 75 minutes each on the subject of parallel solution techniques for the fluid flow equations on unstructured meshes. The title of my lectures is "Parallel CFD Algorithms for Aerodynamical Flow Solvers on Unstructured Meshes" (Parts I-II). The contents of these lectures will be largely review in nature and will draw upon previously published work in this area. Topics of my lectures will include: (1) Mesh partitioning algorithms. Recursive techniques based on coordinate bisection, Cuthill-McKee level structures, and spectral bisection. (2) Newton's method for large scale CFD problems. Size and complexity estimates for Newton's method, modifications for insuring global convergence. (3) Techniques for constructing the Jacobian matrix. Analytic and numerical techniques for Jacobian matrix-vector products, constructing the transposed matrix, extensions to optimization and homotopy theories. (4) Iterative solution algorithms. Practical experience with GIVIRES and BICG-STAB matrix solvers. (5) Parallel matrix preconditioning. Incomplete Lower-Upper (ILU) factorization, domain-decomposed ILU, approximate Schur complement strategies.

Barth, Timothy J.↗
NASA NTRS2001

Maximally Informative Statistics for Localization and Mapping

This paper presents an algorithm for localization and mapping for a mobile robot using monocular vision and odometry as its means of sensing. The approach uses the Variable State Dimension filtering (VSDF) framework to combine aspects of Extended Kalman filtering and nonlinear batch optimization. This paper describes two primary improvements to the VSDF. The first is to use an interpolation scheme based on Gaussian quadrature to linearize measurements rather than relying on analytic Jacobians. The second is to replace the inverse covariance matrix in the VSDF with its Cholesky factor to improve the computational complexity. Results of applying the filter to the problem of localization and mapping with omnidirectional vision are presented.

Deans, Matthew C.↗
NASA NTRS2010

Mapping Successional Stages in a Wet Tropical Forest Using Landsat ETM+ and Forest Inventory Data

In this study, we test whether an existing classification technique based on the integration of Landsat ETM+ and forest inventory data enables detailed characterization of successional stages in a wet tropical forest site. The specific objectives were: (1) to map forest age classes across the La Selva Biological Station in Costa Rica; and (2) to quantify uncertainties in the proposed approach in relation to field data and existing vegetation maps. Although significant relationships between vegetation height entropy (a surrogate for forest age) and ETM+ data were detected, the classification scheme tested in this study was not suitable for characterizing spatial variation in age at La Selva, as evidenced by the error matrix and the low Kappa coefficient (12.9%). Factors affecting the performance of the classification at this particular study site include the smooth transition in vegetation structure between intermediate and advanced successional stages, and the low sensitivity of NDVI to variations in vertical structure at high biomass levels.

forest structure↗
NASA NTRS1993

Proposed framework for thermomechanical life modeling of metal matrix composites

The framework of a mechanics of materials model is proposed for thermomechanical fatigue (TMF) life prediction of unidirectional, continuous-fiber metal matrix composites (MMC's). Axially loaded MMC test samples are analyzed as structural components whose fatigue lives are governed by local stress-strain conditions resulting from combined interactions of the matrix, interfacial layer, and fiber constituents. The metallic matrix is identified as the vehicle for tracking fatigue crack initiation and propagation. The proposed framework has three major elements. First, TMF flow and failure characteristics of in situ matrix material are approximated from tests of unreinforced matrix material, and matrix TMF life prediction equations are numerically calibrated. The macrocrack initiation fatigue life of the matrix material is divided into microcrack initiation and microcrack propagation phases. Second, the influencing factors created by the presence of fibers and interfaces are analyzed, characterized, and documented in equation form. Some of the influences act on the microcrack initiation portion of the matrix fatigue life, others on the microcrack propagation life, while some affect both. Influencing factors include coefficient of thermal expansion mismatch strains, residual (mean) stresses, multiaxial stress states, off-axis fibers, internal stress concentrations, multiple initiation sites, nonuniform fiber spacing, fiber debonding, interfacial layers and cracking, fractured fibers, fiber deflections of crack fronts, fiber bridging of matrix cracks, and internal oxidation along internal interfaces. Equations exist for some, but not all, of the currently identified influencing factors. The third element is the inclusion of overriding influences such as maximum tensile strain limits of brittle fibers that could cause local fractures and ensuing catastrophic failure of surrounding matrix material. Some experimental data exist for assessing the plausibility of the proposed framework.

Halford, Gary R.↗
NASA NTRS1996

Some Experiences with Nonoverlapping Schur Complement Parallel Preconditioning for CFD Calculations

In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Barth, Timothy J.↗
NASA NTRS1997

Effect of Cyclic Thermal Loads on Fatigue Reliability in Polymer Matrix Composites

Technological solutions that will ensure the economic viability and environmental compatibility of a future High Speed Civil Transport plane are currently being sought. Lighter structural materials for both airframe primary structures and engine structure components are being investigated. We believe that such objectives can be achieved through the use of high-temperature composites as well as other conventional, lighter weight alloys. One of the prime issues for these structural components is assured long-term behavior with a specified reliability. An investigation was conducted to describe a computational simulation methodology for predicting fatigue life, reliability, and probabilistic long-term behavior of polymer matrix composites. A unified time-, stress-, and load-dependent Multi- Factor Interaction Equation (MFIE) model developed at the NASA Lewis Research Center was used to simulate the long-term behavior of polymer matrix composites.

Chamis, Christos C.↗
NASA NTRS1990

Identification of thermodynamically stable ceramic reinforcement materials for iron aluminide matrices

Aluminide-base intermetallic matrix composites are currently being considered as potential high-temperature materials. One of the key factors in the selection of a reinforcement material is its chemical stability in the matrix. In this study, chemical interactions between iron aluminides and several potential reinforcement materials, which include carbides, oxides, borides, and nitrides, are analyzed from thermodynamic considerations. Several chemically compatible reinforcement materials are identified for the iron aluminides with Al concentrations ranging from 40 to 50 at. pct.

Misra, Ajay K.↗
NASA NTRS1993

Interface considerations in Al2O3/NiAl composite

The fiber-matrix interface requirements in an Al2O3/NiAl composite were examined from theoretical considerations. Several factors that influence the interface bonding requirements were analyzed. These include: (1) residual stresses due to fiber-matrix coefficient of thermal expansion (CTE) mismatch; (2) matrix cracking stress at room temperature; (3) fracture toughness at room temperature; (4) load transfer from the matrix to the fiber and ultimate tensile strength at the use temperature; and (5) creep resistance at high temperature. A relatively weak fiber-matrix bond, with an interfacial shear strength of approximately 15-20 MPa, might be sufficient for attaining the desired mechanical properties in the fiber direction at the use temperature. A weak fiber-matrix bond is also beneficial for increasing the fracture toughness of the composite at room temperature. In contrast, a strong fiber-matrix bond is required to withstand some of the residual stresses resulting from the fiber-matrix CTE mismatch, which are not likely to be reduced significantly by interface coatings. A relatively strong bond is also beneficial in increasing the matrix cracking stress at room temperature. Various interface coating options to accommodate the conflicting bonding requirements were reviewed. One viable coating option is to incorporate a thick, ductile interface layer well bonded to both the fiber and the matrix.

Misra, Ajay K.↗
NASA NTRS2000

Growth and differentiation of human lens epithelial cells in vitro on matrix

PURPOSE: To characterize the growth and maturation of nonimmortalized human lens epithelial (HLE) cells grown in vitro. METHODS: HLE cells, established from 18-week prenatal lenses, were maintained on bovine corneal endothelial (BCE) extracellular matrix (ECM) in medium supplemented with basic fibroblast growth factor (FGF-2). The identity, growth, and differentiation of the cultures were characterized by karyotyping, cell morphology, and growth kinetics studies, reverse transcription-polymerase chain reaction (RT-PCR), immunofluorescence, and Western blot analysis. RESULTS: HLE cells had a male, human diploid (2N = 46) karyotype. The population-doubling time of exponentially growing cells was 24 hours. After 15 days in culture, cell morphology changed, and lentoid formation was evident. Reverse transcription-polymerase chain reaction (RT-PCR) indicated expression of alphaA- and betaB2-crystallin, fibroblast growth factor receptor 1 (FGFR1), and major intrinsic protein (MIP26) in exponential growth. Western analyses of protein extracts show positive expression of three immunologically distinct classes of crystallin proteins (alphaA-, alphaB-, and betaB2-crystallin) with time in culture. By Western blot analysis, expression of p57(KIP2), a known marker of terminally differentiated fiber cells, was detectable in exponential cultures, and levels increased after confluence. MIP26 and gamma-crystallin protein expression was detected in confluent cultures, by using immunofluorescence, but not in exponentially growing cells. CONCLUSIONS: HLE cells can be maintained for up to 4 months on ECM derived from BCE cells in medium containing FGF-2. With time in culture, the cells demonstrate morphologic characteristics of, and express protein markers for, lens fiber cell differentiation. This in vitro model will be useful for investigations of radiation-induced cataractogenesis and other studies of lens toxicity.

Non-NASA Center↗
NASA NTRS1974

Semiclassical S-matrix theory of vibrationally inelastic collisions between two diatomic molecules

We derive a semiclassical S matrix for vibrationally inelastic collisions between two diatomic molecules, assuming a collinear geometry. Our theory incorporates a quantum mechanical superposition principle with classical dynamics and, as such, is an extension of the atom-diatomic molecule theory of Miller. The several approximations to the S matrix differ in the complexity with which the interference between various classical trajectories is treated. We report numerical calculations for H2-D2 and D2-D2 collisions based on two different interaction potentials. The cruder approximations yield transition probabilities which agree with exact quantum mechanical results to within a factor of 2. More sophisticated approximations to the S matrix yield excellent quantitative agreement with the quantum calculations.

Cohen, S. C.↗
NASA NTRS2016

Hydrogen Isotopic Composition of Apatite in Northwest Africa 7034: A Record of the "Intermediate" H-Isotopic Reservoir in the Martian Crust?

Northwest Africa (NWA) 7034 and its pairings comprise a regolith breccia with a basaltic bulk composition [1] that yields a better match than any other martian meteorite to visible-infrared reflectance spectra of the martian surface measured from orbit [2]. The composition of the fine-grained matrix within NWA 7034 bears a striking resemblance to the major element composition estimated for the martian crust, with several exceptions. The NWA 7034 matrix is depleted in Fe, Ti, and Cr and enriched in Al, Na, and P [3]. The differences in Al and Fe are the most substantial, but the Fe content of NWA 7034 matrix falls within the range reported for the southern highlands crust [6]. It was previously suggested by [4] that NWA 7034 was sourced from the southern highlands based on the ancient 4.4 Ga ages recorded in NWA 7034/7533 zircons [4, 5]. In addition, the NWA 7034 matrix material is enriched in incompatible trace elements by a factor of 1.2-1.5 [7] relative to estimates of the bulk martian crust. The La/Yb ratio of the bulk martian crust is estimated to be approximately 3 [7], and the La/Yb of the NWA 7034 matrix materials ranges from approximately 3.9 to 4.4 [3, 8], indicating a higher degree of LREE enrichment in the NWA 7034 matrix materials. This elevated La/Yb ratio and enrichment in incompatible lithophile trace elements is consistent with NWA 7034 representing a more geochemically enriched crustal terrain than is represented by the bulk martian crust, which would be expected if NWA 7034 represents the bulk crust from the southern highlands. Given the similarities between NWA 7034 and the martian crust, NWA 7034 may represent an important sample for constraining the composition of the martian crust, particularly the ancient highlands. In the present study, we seek to constrain the H isotopic composition of the martian crust using Cl-rich apatite in NWA 7034. Usui et al., [9] recently proposed that a H isotopic reservoir exists within the martian crust that has a H-isotopic composition that is intermediate (δD of 1000-2000per mille) between an isotopically light mantle (Delta D is less than 275per mille [10]) and an isotopically heavy atmosphere (δD of 2500-6100per mille [11, 12]). Apatites in NWA 7034 occur in a number of lithologic domains, however apatites across all lithologic domains have been affected by a Pb-loss event at about 1.5 Ga before present [5], so they are unlikely to have retained a primary composition and are more likely to have equilibrated with fluids within the martian crust that may or may not have exchanged with the martian atmosphere. Equilibration of apatite with crustal fluids is further supported by the chlorine-rich compositions exhibited by apatites in NWA 7034 in comparison to apatites from other martian meteorites (Figure 1; [13]). Cl is more hydrophilic than F, which promotes formation of Cl-rich apatite compositions in fluid-rich systems [e.g., 14, 15-17].

McCubbin, F. M.↗
NASA NTRS1987

A direct algorithm for solution of incompressible three-dimensional unsteady Navier-Stokes equations

A direct, implicit, numerical solution algorithm, with second-order accuracy in space and time, is constructed for the three-dimensional unsteady incompressible Navier-Stokes equations formulated in terms of velocity and vorticity, using generalized orthogonal coordinates to achieve the accurate solution of complex viscous flow configurations. A numerically stable, efficient, direct inversion procedure is developed for the computationally intensive divergence-curl elliptic velocity problem. This overdetermined partial differential operator is first formulated as a uniquely determined, nonsingular matrix-vector problem; this aspect of the procedure is a unique feature of the present analysis. The three-dimensional vorticity-transport equation is solved by a modified factorization technique which completely eliminates the need for any block-matrix inversions and only scalar tridiagonal matrices need to be inverted. The method is applied to the test problem of the three-dimensional flow within a shear-driven cubical box. Coherent streamwise vortex structures are observed within the steady-state flow at Re = 100.

Osswald, G. A.↗