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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Design of robust state feedback controllers via eigenvalue/eigenvector assignment

A new method of selecting a multivariable state feedback controller is presented. The resulting controller simultaneously realizes arbitrary closed-loop eigenvalues, approximates specified modal behavior and achieves low eigensystem sensitivity with respect to plant parameter variations. The method characterizes a vector space slightly larger than the null space of the closed-loop system matrix and projects a desired modal matrix onto this space. Sensitivity of eigenvalues and eigenvectors is then minimized local to the desired modal matrix using a gradient search technique. A tutorial example to illustrate the design procedure is given.

Mielke, R. R.↗
NASA NTRS1985

Lubrication of nonconformal contacts

Minimum film thickness results for piezoviscous-rigid regime of lubrication are developed for a compressible Newtonian fluid with Roelands viscosity. The results provide a basis for the analysis and design of a wide range of machine elements operating in the piezoviscous-rigid regime of lubrication. A new numerical method of calculating elastic deformation in contact stresses is developed using a biquadratic polynomial to approximate the pressure distribution on the whole domain analyzed. The deformation of every node is expressed as a linear combination of the nodal pressures whose coefficients can be combined into an influence coefficient matrix. This approach has the advantages of improved numerical accuracy, less computing time and smaller storage size required for influence matrix. The ideal elastohydrodynamic lubrication is extended to real bearing systems in order to gain an understanding of failure mechanisms in machine elements. The improved elastic deformation calculation is successfully incorporated into the EHL numerical scheme. Using this revised numerical technique and the flow factor model developed by Patir and Cheng (1978) the surface roughness effects on the elastohydrodynamic lubrication of point contact is considered. Conditions typical of an EHL contact in the piezoviscous-elastic regime entrained in pure rolling are investigated. Results are compared with the smooth surface solutions. Experiments are conducted to study the transient EHL effects in instrument ball bearings.

Jeng, Y. R.↗
NASA NTRS1987

Approximate polynomial preconditioning applied to biharmonic equations on vector supercomputers

Applying a finite difference approximation to a biharmonic equation results in a very ill-conditioned system of equations. This paper examines the conjugate gradient method used in conjunction with the generalized and approximate polynomial preconditionings for solving such linear systems. An approximate polynomial preconditioning is introduced, and is shown to be more efficient than the generalized polynomial preconditionings. This new technique provides a simple but effective preconditioning polynomial, which is based on another coefficient matrix rather than the original matrix operator as commonly used.

Wong, Yau Shu↗
NASA NTRS1997

MoSi2-Base Composite for Engine Applications

The intermetallic compound MoSi2 has long been known as a high temperature material that has excellent oxidation resistance and electrical/thermal conductivity. Also its low cost, high melting point (2023 C), relatively low density (6.2 g/cm versus 8 g/cm for current engine materials), and ease of machining make it an attractive structural material. However, the use of MoSi2 has been hindered because of the brittle nature of the material at low temperatures, inadequate creep resistance at high temperatures, accelerated oxidation (also known as 'pest' oxidation) at temperatures between approximately 400 and 500 C, and a coefficient of thermal expansion (CTE) that is relatively high in comparison to potential reinforcing fibers such as SiC. This CTE mismatch between the fiber and the matrix resulted in severe matrix cracking during thermal cycling.

Hebsur, Mohan G.↗
NASA NTRS1978

Thermal-vacuum facility with in-situ mechanical loading

The paper describes a thermal-vacuum space simulator used to assess property changes of fiber-reinforced polymer composite systems. The facility can achieve a vacuum of approximately .0000001 torr with temperatures ranging from -200 to +300 F. Some preliminary experimental results are presented for materials subjected to thermal loading up to 200 F. The tests conducted include the evaluation of matrix modulus and strength, coefficients of thermal expansion, and fracture toughness. Though the experimental program is at an early stage, the data appear to indicate that these parameters are influenced by hard vacuum.

Tennyson, R. C.↗
NASA NTRS1981

Determination of magnetic helicity in the solar wind and implications for cosmic ray propagation

Magnetic helicity (Hm) is the mean value of the correlation between a turbulent magnetic field and the magnetic vector potential. A technique is described for determining Hm and its 'reduced' spectrum from the two point magnetic correlation matrix. The application of the derived formalism to solar wind magnetic fluctuations is discussed, taking into account cases for which only single point measurements are available. The application procedure employs the usual 'frozen in approximation' approach. The considered method is applied to an analysis of several periods of Voyager 2 interplanetary magnetometer data near 2.8 AU. During these periods the correlation length, or energy containing length, was found to be approximately 3 x 10 to the 11th cm

Matthaeus, W. M.↗
NASA NTRS1986

Effects of noise on modal parameters identified by the Eigensystem Realization Algorithm

The basic concept of the Eigensystem Realization Algorithm for modal parameter identification and model reduction is extended to minimize the distortion of the identified parameters caused by noise. The mathematical foundation for the properties of accuracy indicators, such as the singular values of the data matrix and modal amplitude coherence, is provided, based on knowledge of the noise characteristics. These indicators quantitatively discriminate noise from system information and are used to reduce the realized system model to a better approximation of the true model. Monte Carlo Simulations are included to support the analytical studies.

Juang, J.-N.↗
NASA NTRS1986

Effects of noise on ERA-identified modal parameters

The basic concept of the Eigensystem Realization Algorithm (ERA) for modal parameter identification and model reduction is extended to minimize the distortion of the identified parameters caused by noise. The mathematical foundation for the properties of accuracy indicators such as the singular values of the data matrix and modal amplitude coherence is provided based on knowledge of the noise characteristics. These indicators quantitatively discriminate noise from system information and are used to reduce the realized system model to a better approximation of the true model. Monte Carlo simulations are included to support the analytical studies.

Juang, J.-N.↗
NASA NTRS1989

Numerical solution of large Lyapunov equations

A few methods are proposed for solving large Lyapunov equations that arise in control problems. The common case where the right hand side is a small rank matrix is considered. For the single input case, i.e., when the equation considered is of the form AX + XA(sup T) + bb(sup T) = 0, where b is a column vector, the existence of approximate solutions of the form X = VGV(sup T) where V is N x m and G is m x m, with m small is established. The first class of methods proposed is based on the use of numerical quadrature formulas, such as Gauss-Laguerre formulas, applied to the controllability Grammian. The second is based on a projection process of Galerkin type. Numerical experiments are presented to test the effectiveness of these methods for large problems.

Saad, Youcef↗
NASA NTRS1993

On the cross-stream spectral method for the Orr-Sommerfeld equation

Cross-stream models are defined as solutions to the Orr-Sommerfeld equation which are propagating normal to the flow direction. These models are utilized as a basis for a Hilbert space to approximate the spectrum of the Orr-Sommerfeld equation with plane Poiseuille flow. The cross-stream basis leads to a standard eigenvalue problem for the frequencies of Poiseuille flow instability waves. The coefficient matrix in the eigenvalue problem is shown to be the sum of a real matrix and a negative-imaginary diagonal matrix which represents the frequencies of the cross-stream modes. The real coefficient matrix is shown to approach a Toeplitz matrix when the row and column indices are large. The Toeplitz matrix is diagonally dominant, and the diagonal elements vary inversely in magnitude with diagonal position. The Poiseuille flow eigenvalues are shown to lie within Gersgorin disks with radii bounded by the product of the average flow speed and the axial wavenumber. It is shown that the eigenvalues approach the Gersgorin disk centers when the mode index is large, so that the method may be used to compute spectra with an essentially unlimited number of elements. When the mode index is large, the real part of the eigenvalue is the product of the axial wavenumber and the average flow speed, and the imaginary part of the eigen value is identical to the corresponding cross-stream mode frequency. The cross-stream method is numerically well-conditioned in comparison to Chebyshev based methods, providing equivalent accuracy for small mode indices and superior accuracy for large indices.

Zorumski, William E.↗
NASA NTRS1995

Axtrell, a new CV3 chondrite find from Texas

We describe a previously unreported meteorite found in Axtell, Texas, in 1943. Based on the mineralogical composition and texture of its matrix and the sizes and abundance of chondrules, we classify it as a CV3 carbonaceous chondrite. The dominant opaque phase in the chondrules is magnetite, and that in refractory inclusions is Ni-rich metal (awaruite). Axtell, therefore, belongs to the oxidized subgroup of CV3 chondrites, although unlike Allende it escaped strong sulfidation. The meteorite bears a strong textural resemblance to Allende, and its chondrule population and matrix appear to be quite similar to those of Allende, but its refractory inclusions, thermoluminescence properties, and cosmogenic Co-60 abundances are not. Our data are consistent with a terrestrial age for Axtell of approximately 100 years and a metamorphic grade slightly lower than that of Allende.

Simon, S. B.↗
NASA NTRS2014

Assessment of the Accuracy of the Conventional Ray-Tracing Technique: Implications in Remote Sensing and Radiative Transfer Involving Ice Clouds.

A fundamental problem in remote sensing and radiative transfer simulations involving ice clouds is the ability to compute accurate optical properties for individual ice particles. While relatively simple and intuitively appealing, the conventional geometric-optics method (CGOM) is used frequently for the solution of light scattering by ice crystals. Due to the approximations in the ray-tracing technique, the CGOM accuracy is not well quantified. The result is that the uncertainties are introduced that can impact many applications. Improvements in the Invariant Imbedding T-matrix method (II-TM) and the Improved Geometric-Optics Method (IGOM) provide a mechanism to assess the aforementioned uncertainties. The results computed by the II-TMþIGOM are considered as a benchmark because the IITM solves Maxwell's equations from first principles and is applicable to particle size parameters ranging into the domain at which the IGOM has reasonable accuracy. To assess the uncertainties with the CGOM in remote sensing and radiative transfer simulations, two independent optical property datasets of hexagonal columns are developed for sensitivity studies by using the CGOM and the II-TMþIGOM, respectively. Ice cloud bulk optical properties obtained from the two datasets are compared and subsequently applied to retrieve the optical thickness and effective diameter from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Additionally, the bulk optical properties are tested in broadband radiative transfer (RT) simulations using the general circulation model (GCM) version of the Rapid Radiative Transfer Model (RRTMG) that is adopted in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM, version 5.1). For MODIS retrievals, the mean bias of uncertainties of applying the CGOM in shortwave bands (0.86 and 2.13 micrometers) can be up to 5% in the optical thickness and as high as 20% in the effective diameter, depending on cloud optical thickness and effective diameter. In the MODIS infrared window bands centered at 8.5, 11, and 12 micrometers biases in the optical thickness and effective diameter are up to 12% and 10%, respectively. The CGOM-based simulation errors in ice cloud radiative forcing calculations are on the order of 10Wm(exp 2).

remote sensing↗
NASA NTRS2018

Fuzzy Modeling and Parallel Distributed Compensation for Aircraft Flight Control from Simulated Flight Data

A method is described that combines fuzzy system identification techniques with Parallel Distributed Compensation (PDC) to develop nonlinear control methods for aircraft using minimal a priori knowledge, as part of NASA’s Learn-to-Fly initiative. A fuzzy model was generated with simulated flight data, and consisted of a weighted average of multiple linear time invariant state-space cells having parameters estimated using the equation-error approach and a least-squares estimator. A compensator was designed for each subsystem using Linear Matrix Inequalities (LMI) to guarantee closed-loop stability and performance requirements. This approach is demonstrated using simulated flight data to automatically develop a fuzzy model and design control laws for a simplified longitudinal approximation of the F-16 nonlinear flight dynamics simulation. Results include a comparison of flight data with the estimated fuzzy models and simulations that illustrate the feasibility and utility of the combined fuzzy modeling and control approach.

Weinstein, Rose↗
NASA NTRS2005

Effects of Control Mode and R-Ratio on the Fatigue Behavior of a Metal Matrix Composite

Composite Because of their high specific stiffness and strength at elevated temperatures, continuously reinforced metal matrix composites (MMC's) are under consideration for a future generation of aeropropulsion systems. Since components in aeropropulsion systems experience substantial cyclic thermal and mechanical loads, the fatigue behavior of MMC's is of great interest. Almost without exception, previous investigations of the fatigue behavior of MMC's have been conducted in a tension-tension, load-controlled mode. This has been due to the fact that available material is typically less than 2.5-mm thick and, therefore, unable to withstand high compressive loads without buckling. Since one possible use of MMC's is in aircraft skins, this type of testing mode may be appropriate. However, unlike aircraft skins, most engine components are thick. In addition, the transient thermal gradients experienced in an aircraft engine will impose tension-compression loading on engine components, requiring designers to understand how the MMC will behave under fully reversed loading conditions. The increased thickness of the MMC may also affect the fatigue life. Traditionally, low-cycle fatigue (LCF) tests on MMC's have been performed in load control. For monolithic alloys, low-cycle fatigue tests are more typically performed in strain control. Two reasons justify this choice: (1) the critical volume from which cracks initiate and grow is generally small and elastically constrained by the larger surrounding volume of material, and (2) load-controlled, low-cycle fatigue tests of monolithics invariably lead to unconstrained ratcheting and localized necking--an undesired material response because the failure mechanism is far more severe than, and unrelated to, the fatigue mechanism being studied. It is unknown if this is the proper approach to composite testing. However, there is a lack of strain-controlled data on which to base any decisions. Consequently, this study addresses the isothermal, LCF behavior of a [0]_32 MMC tested under strain- and load-controlled conditions for both zero-tension and tension-compression loading conditions. These tests were run at 427 C on thick specimens of SiC-reinforced Ti-15-3. For the fully-reversed tests, no difference was observed in the lives between the load- and strain-controlled tests. However, for the zero-tension tests, the strain-controlled tests had longer lives by a factor of 3 in comparison to the load-controlled tests. This was due to the fact that under strain-control the specimens cyclically softened, reducing the cracking potential. In contrast, the load-controlled tests ratcheted toward larger tensile strains leading to an eventual overload of the fibers. Fatigue tests revealed that specimens tested under fully-reversed conditions had lives approximately an order of magnitude longer than for those specimens tested under zero tension. When examined on a strain-range basis, the fully reversed specimens had similar, but still shorter lives than those of the unreinforced matrix material. However, the composite had a strain limitation at short lives because of the limited strain capacity of the brittle ceramic fiber. The composite also suffered at very high lives because of the lack of an apparent fatigue limit in comparison to the unreinforced matrix. The value of adding fibers to the matrix is apparent when the fatigue lives are plotted as a function of stress range. Here, the composite is far superior to the unreinforced matrix because of the additional load-carrying capacity of the fibers.

Source record↗
NASA NTRS1993

A test of the applicability of independent scattering to high albedo planetary regoliths

We show that 'independent scattering' is a useful approximation for high albedo particles whose size and packing density are similar to typical particles in the lunar regolith. Laboratory measurements of the intensity and linear polarization of light scattered from a laboratory sample of glass spheres of know size and composition are compared to radiative transfer calculations of the same observable quantities. Mie scattering is integrated over the size distribution of the particles to determine the mean phase and polarization functions, or phase matrix, of the particle. Assuming that the particles scatter independently, the 'doubling method' is used to rigorously calculate multiple scattering in an optically thick layer of these anisotropically scattering particles. All of the major features 'predicted' in the calculated intensity (double peaks at small phase angles) and polarization (negative branch at small phase angles, large positive peak near 20 degrees phase, and small polarization elsewhere) are observed in the laboratory measurements, with good quantitative agreement indicated at phase angles less than 90 degrees. Even though the particles are supported by physical contact with each other in the sample, as are the particles in planetary regoliths, the independent scattering calculation yields a good approximation to both the intensity and polarization. The physical parameters input to the calculation are only the size distribution of the particles and their complex index of refraction (composition). Significant advantages of this approach are that the phase matrix is calculated from basic physical principles and both the intensity and polarization are determined simultaneously. This model may have broad applications to the interpretation of photometry, spectroscopy, and polarimetry of the icy regoliths of high albedo satellites. The intent of this effort is to perform a controlled experiment that tests the utility of the independent scattering assumption

Goguen, J. D.↗
NASA NTRS1972

On a stochastic control method for weakly coupled linear systems

The stochastic control of two weakly coupled linear systems with different controllers is considered. Each controller only makes measurements about his own system; no information about the other system is assumed to be available. Based on the noisy measurements, the controllers are to generate independently suitable control policies which minimize a quadratic cost functional. To account for the effects of weak coupling directly, an approximate model, which involves replacing the influence of one system on the other by a white noise process is proposed. Simple suboptimal control problem for calculating the covariances of these noises is solved using the matrix minimum principle. The overall system performance based on this scheme is analyzed as a function of the degree of intersystem coupling.

Kwong, R. H.↗
NASA NTRS1972

The k-space formulation of the n-dimensional scattering problem

The n-dimensional scattering problem is solved by means of a k-space formulation of the field equations, thereby replacing the conventional integral equation formulation by a set of two algebraic equations in two unknowns in two spaces (the constitutive equation being an algebraic equation in x-space). These equations are solved by an iterative method with the aid of the fast Fourier transform (FFT) algorithm connecting the two spaces, requiring very simple initial approximations. Since algebraic and FFT equations are used, the number of arithmetic multiple-add operations and storage allocations required for a numerical solution are reduced from the order of N sq (for solving the matrix equations resulting from the conventional integral equations) to the order of N(log base 2 of N) and N, respectively (where N is the number of data points required for the specification of the problem). The advantage gained in speed and storage is thus of the order of N/log base 2 of N and N, respectively. This method is thus considerably more efficient than the conventional matrix method, and permits exact numerical solutions for much larger problems. Arguments are presented toward the view that the field equations are more fundamental in k-space. The details and some numerical results of the application of this method to the three-dimensional electromagnetic scattering problems are presented as an example.

Bojarski, N. N.↗
NASA NTRS1973

Nonequilibrium chemistry boundary layer integral matrix procedure

The development of an analytic procedure for the calculation of nonequilibrium boundary layer flows over surfaces of arbitrary catalycities is described. An existing equilibrium boundary layer integral matrix code was extended to include nonequilibrium chemistry while retaining all of the general boundary condition features built into the original code. For particular application to the pitch-plane of shuttle type vehicles, an approximate procedure was developed to estimate the nonequilibrium and nonisentropic state at the edge of the boundary layer.

Tong, H.↗