Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

High Temperature Mechanical Characterization and Analysis of Al2O3 /Al2O3 Composition

Sixteen ply unidirectional zirconia coated single crystal Al2O3 fiber reinforced polycrystalline Al2O3 was tested in uniaxial tension at temperatures to 1400 C in air. Fiber volume fractions ranged from 26 to 31%. The matrix has primarily open porosity of approximately 40%. Theories for predicting the Young's modulus, first matrix cracking stress, and ultimate strength were applied and evaluated for suitability in predicting the mechanical behavior of Al2O3/Al2O3 composites. The composite exhibited pseudo tough behavior (increased area under the stress/strain curve relative to monolithic alumina) from 22 to 1400 C. The rule-of-mixtures provides a good estimate of the Young's modulus of the composite using the constituent properties from room temperature to approximately 1200 C for short term static tensile tests in air. The ACK theory provides the best approximation of the first matrix cracking stress while accounting for residual stresses at room temperature. Difficulties in determining the fiber/matrix interfacial shear stress at high temperatures prevented the accurate prediction of the first matrix cracking stress above room temperature. The theory of Cao and Thouless, based on Weibull statistics, gave the best prediction for the composite ultimate tensile strength.

Gyekenyesi, John Z.↗
NASA NTRS1989

Simultaneous quaternion estimation (QUEST) and bias determination

Tests of a new method for the simultaneous estimation of spacecraft attitude and sensor biases, based on a quaternion estimation algorithm minimizing Wahba's loss function are presented. The new method is compared with a conventional batch least-squares differential correction algorithm. The estimates are based on data from strapdown gyros and star trackers, simulated with varying levels of Gaussian noise for both inertially-fixed and Earth-pointing reference attitudes. Both algorithms solve for the spacecraft attitude and the gyro drift rate biases. They converge to the same estimates at the same rate for inertially-fixed attitude, but the new algorithm converges more slowly than the differential correction for Earth-pointing attitude. The slower convergence of the new method for non-zero attitude rates is believed to be due to the use of an inadequate approximation for a partial derivative matrix. The new method requires about twice the computational effort of the differential correction. Improving the approximation for the partial derivative matrix in the new method is expected to improve its convergence at the cost of increased computational effort.

Markley, F. Landis↗
NASA NTRS1967

Renormalized Born approximation.

Renormalizing approximate wave functions so that amplitude is correct by means of matrix, with applications to Born approximation

BORN APPROXIMATION↗
NASA NTRS1995

Metamorphism and aqueous alteration in low petrographic type ordinary chondrites

In order to investigate the relative importance of dry metamorphism and aqueous alteration in the history of chondruies, chondruies were hand-picked from the Semarkona (petrographic type 3.0), Bishunpur (3. 1), Chainpur (3.4), Dhajala (3.8) and Allegan (5) chondrites, and matrix samples were extracted from the first three ordinary chondrites. The thermoluminescence (TL) properties of all the samples were measured, and appropriate subsets of the samples were analyzed by electron-microprobe and radiochemical neutron activation and the water and H-isotopic composition determined. The TL data for chondrules from Semarkona and Bishunpur scatter widely showing no unambiguous trends, although group B1 chondrules tend to have lower sensitivities and lower peak temperatures compared with group A5 chondrules. It is argued that these data reflect the variety of processes accompanying chondrule formation. The chondrules show remarkably uniform contents of the highly labile elements, indicating mineralogical control on abundance and volatile loss from silicates and loss and recondensation of mobile chalcophiles and siderophiles in some cases. Very high D/H values (up to approx. 8000% SMOW) are observed in certain Semarkona chondrules, a confirmation of earlier work. With increasing petrographic type, mean TL sensitivities of the chondrules increase, the spread of values within an individual meteorite decreases, and peak temperatures and peak widths show trends indicating that the TL is mainly produced by feldspar and that dry, thermal metamorphism is the dominant secondary process experienced by the chondrules. The TL sensitivities of matrix samples also increase with petrographic type. Chainpur matrix samples show the same spread of peak temperatures and peak widths as Chainpur chondruies, indicating metamorphism-related changes in the feldspar are responsible for the TL of the matrix. The TL data for the Semarkona and Bishunpur matrix samples provide, at best, only weak evidence for aqueous alteration, but the matrix contains H with approximately terrestrial D/H values, even though it contains much water. Secondary processes (probably aqueous alteration) presumably lowered the D/H of the matrix and certain chondrules. While chondrule properties appear to be governed primarily by formation processes and subsequent metamorphism, the matrix of Semarkona has a more complex history involving aqueous alteration as a meteorite-wide process.

Xie, T.↗
NASA NTRS1988

Approximation theory for LQG (Linear-Quadratic-Gaussian) optimal control of flexible structures

An approximation theory is presented for the LQG (Linear-Quadratic-Gaussian) optimal control problem for flexible structures whose distributed models have bounded input and output operators. The main purpose of the theory is to guide the design of finite dimensional compensators that approximate closely the optimal compensator. The optimal LQG problem separates into an optimal linear-quadratic regulator problem and an optimal state estimation problem. The solution of the former problem lies in the solution to an infinite dimensional Riccati operator equation. The approximation scheme approximates the infinite dimensional LQG problem with a sequence of finite dimensional LQG problems defined for a sequence of finite dimensional, usually finite element or modal, approximations of the distributed model of the structure. Two Riccati matrix equations determine the solution to each approximating problem. The finite dimensional equations for numerical approximation are developed, including formulas for converting matrix control and estimator gains to their functional representation to allow comparison of gains based on different orders of approximation. Convergence of the approximating control and estimator gains and of the corresponding finite dimensional compensators is studied. Also, convergence and stability of the closed-loop systems produced with the finite dimensional compensators are discussed. The convergence theory is based on the convergence of the solutions of the finite dimensional Riccati equations to the solutions of the infinite dimensional Riccati equations. A numerical example with a flexible beam, a rotating rigid body, and a lumped mass is given.

Gibson, J. S.↗
NASA NTRS1988

An invariant asymptotic formula for solutions of second-order linear ODE's

An invariant-matrix technique for the approximate solution of second-order ordinary differential equations (ODEs) of form y-double-prime = phi(x)y is developed analytically and demonstrated. A set of linear transformations for the companion matrix differential system is proposed; the diagonalization procedure employed in the final stage of the asymptotic decomposition is explained; and a scalar formulation of solutions for the ODEs is obtained. Several typical ODEs are analyzed, and it is shown that the Liouville-Green or WKB approximation is a special case of the present formula, which provides an approximation which is valid for the entire interval (0, infinity).

Gingold, H.↗
NASA NTRS2001

Dynamic modeling of gene expression data

We describe the time evolution of gene expression levels by using a time translational matrix to predict future expression levels of genes based on their expression levels at some initial time. We deduce the time translational matrix for previously published DNA microarray gene expression data sets by modeling them within a linear framework by using the characteristic modes obtained by singular value decomposition. The resulting time translation matrix provides a measure of the relationships among the modes and governs their time evolution. We show that a truncated matrix linking just a few modes is a good approximation of the full time translation matrix. This finding suggests that the number of essential connections among the genes is small.

NASA Discipline Plant Biology↗
NASA NTRS1985

Adaptive Filtering for Large Space Structures: A Closed-Form Solution

In a previous paper Schaechter proposes using an extended Kalman filter to estimate adaptively the (slowly varying) frequencies and damping ratios of a large space structure. The time varying gains for estimating the frequencies and damping ratios can be determined in closed form so it is not necessary to integrate the matrix Riccati equations. After certain approximations, the time varying adaptive gain can be written as the product of a constant matrix times a matrix derived from the components of the estimated state vector. This is an important savings of computer resources and allows the adaptive filter to be implemented with approximately the same effort as the nonadaptive filter. The success of this new approach for adaptive filtering was demonstrated using synthetic data from a two mode system.

Rauch, H. E.↗
NASA NTRS1996

Practical Aspects of Krylov Subspace Iterative Methods in CFD

Implementation issues associated with the application of Krylov subspace iterative methods, such as Newton-GMRES, are presented within the framework of practical computational fluid dynamic (CFD) applications. This paper categorizes, evaluates, and contrasts the major ingredients (function evaluations, matrix-vector products, and preconditioners) of Newton-GMRES Krylov subspace methods in terms of their effect on the local linear and global nonlinear convergence, memory requirements, and accuracy. The discussion focuses on Newton-GMRES in both a structured multi-zone incompressible Navier-Stokes solver and an unstructured mesh finite-volume Navier-Stokes solver. Approximate versus exact matrix-vector products, effective preconditioners, and other pertinent issues are addressed.

Pulliam, Thomas H.↗
NASA NTRS1978

Conservative implicit schemes for the full potential equation applied to transonic flows

Implicit approximate factorization techniques (AF) were investigated for the solution of matrix equations resulting from finite difference approximations to the full potential equation in conservation form. For transonic flows, an artificial viscosity, required to maintain stability in supersonic regions, was introduced by an upwind bias of the density. Two implicit AF procedures are presented and their convergence performance is compared with that of the standard transonic solution procedure, successive line overrelaxation (SLOR). Subcritical and supercritical test cases are considered. The results indicate that the AF schemes are substantially faster than SLOR.

Holst, T. L.↗
NASA NTRS2005

Kurtosis Approach Nonlinear Blind Source Separation

In this paper, we introduce a new algorithm for blind source signal separation for post-nonlinear mixtures. The mixtures are assumed to be linearly mixed from unknown sources first and then distorted by memoryless nonlinear functions. The nonlinear functions are assumed to be smooth and can be approximated by polynomials. Both the coefficients of the unknown mixing matrix and the coefficients of the approximated polynomials are estimated by the gradient descent method conditional on the higher order statistical requirements. The results of simulation experiments presented in this paper demonstrate the validity and usefulness of our approach for nonlinear blind source signal separation Keywords: Independent Component Analysis, Kurtosis, Higher order statistics.

Duong, Vu A.↗
NASA NTRS2005

Kurtosis Approach for Nonlinear Blind Source Separation

In this paper, we introduce a new algorithm for blind source signal separation for post-nonlinear mixtures. The mixtures are assumed to be linearly mixed from unknown sources first and then distorted by memoryless nonlinear functions. The nonlinear functions are assumed to be smooth and can be approximated by polynomials. Both the coefficients of the unknown mixing matrix and the coefficients of the approximated polynomials are estimated by the gradient descent method conditional on the higher order statistical requirements. The results of simulation experiments presented in this paper demonstrate the validity and usefulness of our approach for nonlinear blind source signal separation.

kurtosis↗
NASA NTRS2015

Theoretical Studies of Spectroscopic Line Mixing in Remote Sensing Applications

The phenomenon of collisional transfer of intensity due to line mixing has an increasing importance for atmospheric monitoring. From a theoretical point of view, all relevant information about the collisional processes is contained in the relaxation matrix where the diagonal elements give half-widths and shifts, and the off-diagonal elements correspond to line interferences. For simple systems such as those consisting of diatom-atom or diatom-diatom, accurate fully quantum calculations based on interaction potentials are feasible. However, fully quantum calculations become unrealistic for more complex systems. On the other hand, the semi-classical Robert-Bonamy (RB) formalism, which has been widely used to calculate half-widths and shifts for decades, fails in calculating the off-diagonal matrix elements. As a result, in order to simulate atmospheric spectra where the effects from line mixing are important, semi-empirical fitting or scaling laws such as the ECS (Energy-Corrected Sudden) and IOS (Infinite-Order Sudden) models are commonly used. Recently, while scrutinizing the development of the RB formalism, we have found that these authors applied the isolated line approximation in their evaluating matrix elements of the Liouville scattering operator given in exponential form. Since the criterion of this assumption is so stringent, it is not valid for many systems of interest in atmospheric applications. Furthermore, it is this assumption that blocks the possibility to calculate the whole relaxation matrix at all. By eliminating this unjustified application, and accurately evaluating matrix elements of the exponential operators, we have developed a more capable formalism. With this new formalism, we are now able not only to reduce uncertainties for calculated half-widths and shifts, but also to remove a once insurmountable obstacle to calculate the whole relaxation matrix. This implies that we can address the line mixing with the semi-classical theory based on interaction potentials between molecular absorber and molecular perturber. We have applied this formalism to address the line mixing for Raman and infrared spectra of molecules such as N2, C2H2, CO2, NH3, and H2O. By carrying out rigorous calculations, our calculated relaxation matrices are in good agreement with both experimental data and results derived from the ECS model.

Formalism↗
NASA NTRS1998

High Temperature Mechanical Characterization of Ceramic Matrix Composites

A high temperature mechanical characterization laboratory has been assembled at NASA Lewis Research Center. One contribution of this work is to test ceramic matrix composite specimens in tension in environmental extremes. Two high temperature tensile testing systems were assembled. The systems were assembled based on the performance and experience of other laboratories and meeting projected service conditions for the materials in question. The systems use frames with an electric actuator and a center screw. A PC based data acquisition and analysis system is used to collect and analyze the data. Mechanical extensometers are used to measure specimen strain. Thermocouples, placed near the specimen, are used to measure the specimen gage section temperature. The system for testing in air has a resistance element furnace with molybdenum disilicide elements and pneumatic grips with water cooling attached to hydraulic alignment devices. The system for testing in an inert gas has a graphite resistance element furnace in a chamber with rigidly mounted, water cooled, hydraulically actuated grips. Unidirectional SiC fiber reinforced reaction bonded Si3N4 and triaxially woven, two dimensional, SiC fiber reinforced enhanced SiC composites were tested in unidirectional tension. Theories for predicting the Young's modulus, modulus near the ultimate strength, first matrix cracking stress, and ultimate strength were applied and evaluated for suitability in predicting the mechanical behavior of SiC/RBSN and enhanced SiC/SiC composites. The SiC/RBSN composite exhibited pseudo tough behavior (increased area under the stress/strain curve) from 22 C to 1500 C. The rule of mixtures provides a good estimate of the Young's modulus of the SiC/RBSN composite using the constituent properties from room temperature to 1440 C for short term static tensile tests in air or nitrogen. The rule of mixtures significantly overestimates the secondary modulus near the ultimate strength. The ACK theory provides the best approximation of the first matrix cracking stress when residual stresses are ignored. The theory of Cao and Thouless, based on Weibull statistics, gave the best prediction for the composite ultimate strength. The enhanced SiC/SiC composite exhibited nonlinear stress/strain behavior from 24 C to 1370 C in air with increased ultimate strain when compared to monolithic SiC. The theory of Yang and Chou with the assumption of a frictional fiber/matrix interface provided the best estimate of the Young's modulus. The theory of Cao and Thouless gave the best estimate for the ultimate strength.

Gyekenyesi, John Z.↗
NASA NTRS2008

Modeling State-Space Aeroelastic Systems Using a Simple Matrix Polynomial Approach for the Unsteady Aerodynamics

A simple matrix polynomial approach is introduced for approximating unsteady aerodynamics in the s-plane and ultimately, after combining matrix polynomial coefficients with matrices defining the structure, a matrix polynomial of the flutter equations of motion (EOM) is formed. A technique of recasting the matrix-polynomial form of the flutter EOM into a first order form is also presented that can be used to determine the eigenvalues near the origin and everywhere on the complex plane. An aeroservoelastic (ASE) EOM have been generalized to include the gust terms on the right-hand side. The reasons for developing the new matrix polynomial approach are also presented, which are the following: first, the "workhorse" methods such as the NASTRAN flutter analysis lack the capability to consistently find roots near the origin, along the real axis or accurately find roots farther away from the imaginary axis of the complex plane; and, second, the existing s-plane methods, such as the Roger s s-plane approximation method as implemented in ISAC, do not always give suitable fits of some tabular data of the unsteady aerodynamics. A method available in MATLAB is introduced that will accurately fit generalized aerodynamic force (GAF) coefficients in a tabular data form into the coefficients of a matrix polynomial form. The root-locus results from the NASTRAN pknl flutter analysis, the ISAC-Roger's s-plane method and the present matrix polynomial method are presented and compared for accuracy and for the number and locations of roots.

Pototzky, Anthony S.↗
NASA NTRS2010

Multi-functional layered structure having structural and radiation shielding attributes

A cosmic and solar radiation shielding structure that also has structural attributes is comprised of three layers. The first layer is 30-42 percent by volume of ultra-high molecular weight (UHMW) polyethylene fibers, 18-30 percent by volume of graphite fibers, and a remaining percent by volume of an epoxy resin matrix. The second layer is approximately 68 percent by volume of UHMW polyethylene fibers and a remaining percent by volume of a polyethylene matrix. The third layer is a ceramic material.

Kaul, Raj K.↗