Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATRICES”

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 613 records · Page 34

Aqueous alteration in carbonaceous chondrites - Mass balance constraints on matrix mineralogy

Bulk chemical compositions of matrix material were determined in eleven Antarctic CM chondrites and five non-Antarctic CM and CI chondrites, using microprobe defocused-beam technique. The results, along with previously published data, are used to provide mass balance constraints on the relative proportions of intergrown and intermixed phyllosilicate phases in carbonaceous chondrite matrices. In terms of Fe, Si, and Mg, the CM matrix bulk compositions plot within a triangle defined by the compositions of PCP (a mixture of 25 percent tochilinite and 75 percent cronstedtite), Mg-rich serpentine, and Fe-rich serpentine. Results indicate differing amounts of PCP and serpentines in individual CM matrices. The discrepancies found between predicted and measured S and Ni values require additional sulfide phases. CI matrices were found to contain little, if any, PCP; they consist mostly of serpentine and montmorillonite.

Mcsween, Harry Y., Jr.↗

Bismaleimide-vinylpolystyrylpyridine graphite composites

Flammability, thermal, and selected mechanical properties of composites fabricated with epoxy and other thermally stable resin matrices are described. Properties which were measured included limiting-oxygen index, smoke evolution, thermal degradation products, total-heat release, heat-release rates, mass loss, flame spread, ignition resistance, thermogravimetric analysis and selected mechanical properties. This paper describes the properties of eight different graphite composite panels fabricated using four different resin matrices and two types of graphite reinforcement. The resin matrices included: XU71775/H795, a blend of vinylpolystyrylpyridine and bismaleimide; H795, a bismaleimide; Cycom 6162, a phenolic; and PSP 6022M, a polystyrylpyridine. The graphite fiber used was AS-4 in the form of either tape or fabric. The properties of these composites were compared with epoxy composites. It was determined that the blend of vinylpolystyrylpyridine and bismaleimide (XU71775/H795) with the graphite tape was the optimum design giving the lowest heat release rate.

Kourtides, D. A.↗

FREQ: A computational package for multivariable system loop-shaping procedures

Many approaches in the field of linear, multivariable time-invariant systems analysis and controller synthesis employ loop-sharing procedures wherein design parameters are chosen to shape frequency-response singular value plots of selected transfer matrices. A software package, FREQ, is documented for computing within on unified framework many of the most used multivariable transfer matrices for both continuous and discrete systems. The matrices are evaluated at user-selected frequency-response values, and singular values against frequency. Example computations are presented to demonstrate the use of the FREQ code.

Giesy, Daniel P.↗

SPARSKIT: A basic tool kit for sparse matrix computations

Presented here are the main features of a tool package for manipulating and working with sparse matrices. One of the goals of the package is to provide basic tools to facilitate the exchange of software and data between researchers in sparse matrix computations. The starting point is the Harwell/Boeing collection of matrices for which the authors provide a number of tools. Among other things, the package provides programs for converting data structures, printing simple statistics on a matrix, plotting a matrix profile, and performing linear algebra operations with sparse matrices.

Saad, Youcef↗

Dimensionally Stable Graphite-Fiber/Glass Composites

Method of making composites of glass matrices reinforced by graphite fibers provides for control of proportions, orientations, and distributions of fibers in matrices and for fused bonds between fibers and matrices. Enables fabrication of composites of high specific strength and dimensional stability. Method particularly suitable for making low-thermal-expansion platforms for optical instruments.

Harris, Robert↗

Time and wavelength domain algorithms for chemical analysis by laser radar

Laser-induced fluorescence (LIF) is a promising technique for laser radar applications. Laser radar using LIF has already been applied to algae blooms and oil slicks. Laser radar using LIF has great potential for remote chemical analysis because LIF spectra are extremely sensitive to chemical composition. However, most samples in the real world contain mixtures of fluorescing components, not merely individual components. Multicomponent analysis of laser radar returns from mixtures is often difficult because LIF spectra from solids and liquids are very broad and devoid of line structure. Therefore, algorithms for interpreting LIF spectra from laser radar returns must be able to analyze spectra that overlap in multicomponent systems. This paper analyzes the possibility of using factor analysis-rank annihilation (FARA) to analyze emission-time matrices (ETM) from laser radar returns instead of excitation-emission matrices (EEM). The authors here define ETM as matrices where the rows (or columns) are emission spectra at fixed times and the columns (or rows) are temporal profiles for fixed emission wavelengths. Laser radar usually uses pulsed lasers for ranging purposes, which are suitable for measuring temporal profiles. Laser radar targets are hard instead of diffuse; that is, a definite surface emits the fluorescence instead of an extended volume. A hard target would not broaden the temporal profiles as would a diffuse target. Both fluorescence lifetimes and emission spectra are sensitive to chemical composition. Therefore, temporal profiles can be used instead of excitation spectra in FARA analysis of laser radar returns. The resulting laser radar returns would be ETM instead of EEM.

Rosen, David L.↗

Krylov Subspace Methods for Complex Non-Hermitian Linear Systems

We consider Krylov subspace methods for the solution of large sparse linear systems Ax = b with complex non-Hermitian coefficient matrices. Such linear systems arise in important applications, such as inverse scattering, numerical solution of time-dependent Schrodinger equations, underwater acoustics, eddy current computations, numerical computations in quantum chromodynamics, and numerical conformal mapping. Typically, the resulting coefficient matrices A exhibit special structures, such as complex symmetry, or they are shifted Hermitian matrices. In this paper, we first describe a Krylov subspace approach with iterates defined by a quasi-minimal residual property, the QMR method, for solving general complex non-Hermitian linear systems. Then, we study special Krylov subspace methods designed for the two families of complex symmetric respectively shifted Hermitian linear systems. We also include some results concerning the obvious approach to general complex linear systems by solving equivalent real linear systems for the real and imaginary parts of x. Finally, numerical experiments for linear systems arising from the complex Helmholtz equation are reported.

Freund, Roland W.↗

A transfer matrix approach to vibration localization in mistuned blade assemblies

A study of mode localization in mistuned bladed disks is performed using transfer matrices. The transfer matrix approach yields the free response of a general, mono-coupled, perfectly cyclic assembly in closed form. A mistuned structure is represented by random transfer matrices, and the expansion of these matrices in terms of the small mistuning parameter leads to the definition of a measure of sensitivity to mistuning. An approximation of the localization factor, the spatially averaged rate of exponential attenuation per blade-disk sector, is obtained through perturbation techniques in the limits of high and low sensitivity. The methodology is applied to a common model of a bladed disk and the results verified by Monte Carlo simulations. The easily calculated sensitivity measure may prove to be a valuable design tool due to its system-independent quantification of mistuning effects such as mode localization.

Ottarson, Gisli↗

Applications of multiple-constraint matrix updates to the optimal control of large structures

Low-authority control or vibration suppression in large, flexible space structures can be formulated as a linear feedback control problem requiring computation of displacement and velocity feedback gain matrices. To ensure stability in the uncontrolled modes, these gain matrices must be symmetric and positive definite. In this paper, efficient computation of symmetric, positive-definite feedback gain matrices is accomplished through the use of multiple-constraint matrix update techniques originally developed for structural identification applications. Two systems were used to illustrate the application: a simple spring-mass system and a planar truss. From these demonstrations, use of this multiple-constraint technique is seen to provide a straightforward approach for computing the low-authority gains.

Smith, S. W.↗

Computation of vibration mode elastic-rigid and effective weight coefficients from finite-element computer program output

Post-processing algorithms are given to compute the vibratory elastic-rigid coupling matrices and the modal contributions to the rigid-body mass matrices and to the effective modal inertias and masses. Recomputation of the elastic-rigid coupling matrices for a change in origin is also described. A computational example is included. The algorithms can all be executed by using standard finite-element program eigenvalue analysis output with no changes to existing code or source programs.

Levy, R.↗

Identification of a second-order mechanical system model from a state space realization

Many system identification algorithms produce models in modern state variable form. When a mechanical system is identified one knows that there must be a model in second order form. In structural dynamics one is interested in either a modal model of the system, or in a model including a mass matrix, stiffness matrix and a damping matrix. In this paper, algorithms are developed that convert a modern state space realization into the above two representations of interest in mechanical systems. The algorithm can be used with any identification methods that produce a modern state variable representation such as the Eigensystem Realization Algorithm (ERA), its modified version by data correlation (ERA/DC), and combinations of these with Observer/Kalman Filter Identification (OKID). An algorithm is developed that allows one to identify the damping matrix in a model representation, so that one determines how near to modal the damping is in the system, and can therefore understand the degree to which the modes are coupled by the damping in the system. A second algorithm produced the mass, damping, and stiffness matrices, given the input and output matrices. It is assumed that the number of sensors (or the number of actuators) is greater than or equal to the number of modes in the system, and it is proved that this is a necessary condition to be able to uniquely identify these matrices. Experience with examples suggests that a singular value truncation involved in these algorithms can help one determine the true system order. This truncation has the benefit of knowledge of the correct form for a mechanical system model, which is not present in the modern control identification algorithm.

Tseng, Dong-Huei↗

Rapid calibration of seven-hole probes

This paper summarizes the major conclusions and some of the key supporting analyses resulting from the calibration and application of two small seven hole probes at NASA Lewis Research Center. These probes can produce reasonably accurate and rapid surveys of unknown steady flow fields which may include flow angles up to 70 degrees and Mach numbers up to 0.8. The probes were calibrated with both 'complete' and 'reduced' test matrices. Both types of test matrices produced similar results suggesting the the reduced matrices are adequate for most purposes. The average accuracy fo the calibration was about the same as that achieved in previous seven hole probe calibrations. At the higher Mach numbers, the calibration was sensitive to the diameter of the free jet in the calibration facility. Over a narrow angular range at the higher Mach numbers, the system had serious repeatability problems. This lack or repeatability apparently results from aliasing of high frequency (20 to 40 Hz) noise with the data acquisition system sampling frequency of 10 Hz. Analyses show that these noise frequencies are probably not related to airflow dynamics in the connecting tubing.

Ericksen, A. L.↗

Non-local transport in turbulent MHD convection

The nonlocal non-diffusive transport of passive scalars in turbulent magnetohydrodynamic (MHD) convection is investigated using transilient matrices. These matrices describe the probability that a tracer particle beginning at one position in a flow will be advected to another position after some time. A method for the calculation of these matrices from simulation data which involves following the trajectories of passive tracer particles and calculating their transport statistics, is presented. The method is applied to study the transport in several simulations of turbulent, rotating, three dimensional compressible, penetrative MDH convection. Transport coefficients and other diagnostics are used to quantify the transport, which is found to resemble advection more closely than diffusion. Some of the results are found to have direct relevance to other physical problems, such as the light element depletion in sun-type stars. The large kurtosis found for downward moving particles at the base of the convection zone implies several extreme events.

Miesch, Mark↗

Mullite fiber reinforced reaction bonded Si3N4 composites

Fracture toughnesses of brittle ceramic materials have been improved by introducing reinforcements and carefully tailored interface layers. Silicon carbide and Si3N4 have been emphasized as matrices of structural composites intended for high temperature service because they combine excellent mechanical, chemical, thermal and physical properties. Both matrices have been successfully toughened with SiC fibers, whiskers and particles for ceramic matrix composite (CMC) parts made by sintering, hot pressing or reaction forming processes. These SiC reinforced CMCs have exhibited significantly improved toughnesses at low and intermediate temperature levels, as well as retention of properties at high temperatures for selected exposures; however, they are vulnerable to attack from elevated temperature dry and wet oxidizing atmospheres after the matrix has cracked. Property degradation results from oxidation of interface layers and/or reinforcements. The problem is particularly acute for small diameter (-20 tim) polymer derived SiC fibers used for weavable toes. This research explored opportunities for reinforcing Si3N4 matrices with fibers having improved environmental stability; the findings should also be applicable to SiC matrix CMCs.

Saleh, T.↗

Effect of Composite Fabrication on the Strength of Single Crystal Al2O3 Fibers in Two Fe-Base Alloy Composites

Continuous single-crystal Al2O3 fibers have been incorporated into a variety of metal and intermetallic matrices and the results have consistently indicated that the fiber strength had been reduced by 32 to 50% during processing. Two iron-based alloys, FeNiCoCrAl and FeAlVCMn, were chosen as matrices for Al2O3 fiber reinforced metal matrix composites (MMC) with the goal of maintaining Al2O3 fiber strength after composite processing. The feasibility of Al2O3/FeNiCoCrAl and Al2O3/FeAlVCMn composite systems for high temperature applications were assessed in terms of fiber-matrix chemical compatibility, interfacial bond strength, and composite tensile properties. The strength of etched-out fibers was significantly improved by choosing matrices containing less reactive elements. The ultimate tensile strength (UTS) values of the composites could generally be predicted with existing models using the strength of etched-out fibers. However, the UTS of the composites were less than desired due to a low fiber Weibull modulus. Acoustic emission analysis during tensile testing was a useful tool for determining the efficiency of the fibers in the composite and for determining the failure mechanism of the composites.

Draper, Susan L.↗

Optimal Frequency-Domain System Realization with Weighting

Several approaches are presented to identify an experimental system model directly from frequency response data. The formulation uses a matrix-fraction description as the model structure. Frequency weighting such as exponential weighting is introduced to solve a weighted least-squares problem to obtain the coefficient matrices for the matrix-fraction description. A multi-variable state-space model can then be formed using the coefficient matrices of the matrix-fraction description. Three different approaches are introduced to fine-tune the model using nonlinear programming methods to minimize the desired cost function. The first method uses an eigenvalue assignment technique to reassign a subset of system poles to improve the identified model. The second method deals with the model in the real Schur or modal form, reassigns a subset of system poles, and adjusts the columns (rows) of the input (output) influence matrix using a nonlinear optimizer. The third method also optimizes a subset of poles, but the input and output influence matrices are refined at every optimization step through least-squares procedures.

Juang, Jer-Nan↗

Optimal System Realization in Frequency Domain

Several approaches are presented to identify an experimental system model directly from frequency response data. The formulation begins with a matrix-fraction description as the model structure. Frequency weighting such as exponential weighting is introduced to solve a weighted least-squares problem to obtain the coefficient matrices for the matrix-fraction description. A multi-variable state-space model can then be formed using the coefficient matrices of the matrix-fraction description. An approach is introduced to fine-tune the model using non-linear programming methods to minimize the desired cost function. The method deals with the model in the real Schur or modal form and reassigns a subset of system poles using a nonlinear optimizer. At every optimization step, the input and output influence matrices are refined through least-squares procedures. The proposed approaches are used to identify an analytical model for a NASA testbed from experimental data.

Juang, Jer-Nan↗

User's Manual for PCSMS (Parallel Complex Sparse Matrix Solver)

PCSMS (Parallel Complex Sparse Matrix Solver) is a computer code written to make use of the existing real sparse direct solvers to solve complex, sparse matrix linear equations. PCSMS converts complex matrices into real matrices and use real, sparse direct matrix solvers to factor and solve the real matrices. The solution vector is reconverted to complex numbers. Though, this utility is written for Silicon Graphics (SGI) real sparse matrix solution routines, it is general in nature and can be easily modified to work with any real sparse matrix solver. The User's Manual is written to make the user acquainted with the installation and operation of the code. Driver routines are given to aid the users to integrate PCSMS routines in their own codes.

Reddy, C. J.↗