Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix decomposition”

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 199 records · Page 11

In situ deposits of copper and copper oxide containing condensation polyimide films

Novel copper-polyimide composites have been synthesized via simultaneous thermal decomposition of solid solutions of bis (trifluoroacetylacetonato) copper (II) and thermal cyclodehydration of polyimide acid. In contrast to conventional filled polymer composites which are prepared by dispersion of particles or fibers in a polymer matrix this study has yielded in general uniform Cu or CuO dispersions of very small particle size that reside near the film surface that was exposed to the atmosphere during curing. The nature of the copper deposit, the thickness of the copper deposit, and the polyimide overlayer which bonds the copper to the polymer substrate depend on the curing atmosphere used. A variety of analytical surface methods along with thermogravimetric analysis and variable temperature (surface and volume) electrical resistivity measurements have been used to characterize these thin, flexible copper doped polyimide films.

Porta, G. M.↗

Quaternion Averaging

Many applications require an algorithm that averages quaternions in an optimal manner. For example, when combining the quaternion outputs of multiple star trackers having this output capability, it is desirable to properly average the quaternions without recomputing the attitude from the the raw star tracker data. Other applications requiring some sort of optimal quaternion averaging include particle filtering and multiple-model adaptive estimation, where weighted quaternions are used to determine the quaternion estimate. For spacecraft attitude estimation applications, derives an optimal averaging scheme to compute the average of a set of weighted attitude matrices using the singular value decomposition method. Focusing on a 4-dimensional quaternion Gaussian distribution on the unit hypersphere, provides an approach to computing the average quaternion by minimizing a quaternion cost function that is equivalent to the attitude matrix cost function Motivated by and extending its results, this Note derives an algorithm that deterniines an optimal average quaternion from a set of scalar- or matrix-weighted quaternions. Rirthermore, a sufficient condition for the uniqueness of the average quaternion, and the equivalence of the mininiization problem, stated herein, to maximum likelihood estimation, are shown.

Markley, F. Landis↗

Optical matrix-vector processing for computational fluid dynamics

An optical processor to solve partial differential equations for computational fluid dynamics applications is considered. This application is new and original for optical processors. The algorithms that are used are optical realizations of the Newton-Raphson method for nonlinear equations and a new optical LU direct decomposition and Gauss-Seidel iterative solution to the resultant linear algebraic equations. These algorithms are used to solve Burger's equation (a specific form of the momentum equation in fluid dynamics). The nonlinear equations provide 1-D velocity data at each time step. Simulation results of optical processing with these algorithms on computational fluid dynamics data is included.

Perlee, Caroline J.↗

The High-Resolution Wave-Propagation Method Applied to Meso- and Micro-Scale Flows

The high-resolution wave-propagation method for computing the nonhydrostatic atmospheric flows on meso- and micro-scales is described. The design and implementation of the Riemann solver used for computing the Godunov fluxes is discussed in detail. The method uses a flux-based wave decomposition in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The two advantages of the technique are: 1) the need for an explicit definition of the Roe matrix is eliminated and, 2) the inclusion of source term due to gravity does not result in discretization errors. The resulting flow solver is conservative and able to resolve regions of large gradients without introducing dispersion errors. The methodology is validated against exact analytical solutions and benchmark cases for non-hydrostatic atmospheric flows.

Ahmad, Nashat N.↗

Parallel pseudospectral domain decomposition techniques

The influence of interface boundary conditions on the ability to parallelize pseudospectral multidomain algorithms is investigated. Using the properties of spectral expansions, a novel parallel two domain procedure is generalized to an arbitrary number of domains each of which can be solved on a separate processor. This interface boundary condition considerably simplifies influence matrix techniques.

Gottlieb, David↗

Parallel pseudospectral domain decomposition techniques

The influence of interface boundary conditions on the ability to parallelize pseudospectral multidomain algorithms is investigated. Using the properties of spectral expansions, a novel parallel two domain procedure is generalized to an arbitrary number of domains each of which can be solved on a separate processor. This interface boundary condition considerably simplifies influence matrix techniques.

Gottlieb, David↗

A split band-Cholesky equation solving strategy for finite element analysis of transient field problems

The paper describes the split-Cholesky strategy for banded matrices arising from the large systems of equations in certain fluid mechanics problems. The basic idea is that for a banded matrix the computation can be carried out in pieces, with only a small portion of the matrix residing in core. Mesh considerations are discussed by demonstrating the manner in which the assembly of finite element equations proceeds for linear trial functions on a triangular mesh. The FORTRAN code which implements the out-of-core decomposition strategy for banded symmetric positive definite matrices (mass matrices) of a coupled initial value problem is given.

Cooke, C. H.↗

How to Estimate Attitude from Vector Observations

The most robust estimators minimizing Wahba's loss function are Davenport's q method and the Singular Value Decomposition (SVD) method. The q method is faster than the SVD method with three or more measurements. The other algorithms are less robust since they solve the characteristic polynomial equation to find the maximum eigenvalue of Davenport's K matrix. They are only preferable when speed or processor power is an important consideration. Of these, Fast Optimal Attitude Matrix (FOAM) is the most robust and faster than the q method. Robustness is only an issue for measurements with widely differing accuracies, so the fastest algorithms, Quaternion ESTimator (QUEST), EStimator of the Optimal Quaternion (ESOQ), and ESOQ2, are well suited to star sensor applications.

Markley, F. Landis↗

A finite element method for the thermochemical decomposition of polymeric materials. I - Theory

The governing differential equations are developed to model the thermomechanical behavior of chemically decomposing, polymeric materials. These equations account for thermal and gaseous diffusion through a poroelastic, transversely isotropic solid. The Bubnov-Galerkin finite element method is applied to the governing equations to cast the coupled set into a single matrix equation. A method for solving these equations simultaneously at each time step is discussed.

Sullivan, R. M.↗

Completion of the Viking Labeled Release experiment on Mars

The final Labeled Release (LR) cycle on each Viking lander tested a surface sample that had been stored for several months at approximately 10 C prior to the onset of the active sequence. At each lander site, activity was strongly diminished. This thermal sensitivity of the active agent on the surface of Mars is consistent with a biological explanation of the LR experiment. At the end of one of these cycles, the incubation mixture was heated to 50 C to release any radioactive gas trapped in the sample matrix. The results suggest that more than one carbon substrate is involved in the LR reaction on Mars. The thermal data from the stored samples, coupled with data from previous cycles, have formed the basis for evaluation of the thermal decomposition of the Mars active agent. The slope of the resulting Arrhenius plot has been used to test the fit of other flight data and to calculate the activation energy for thermal decomposition of the Mars agent. The results and their interpretation still leave unresolved the question of whether the Mars LR data were generated by biological or chemical activity.

Levin, G. V.↗

Resistivity of Carbon-Carbon Composites Halved

Carbon-carbon composites have become the material of choice for applications requiring strength and stiffness at very high temperatures (above 2000 C). These composites comprise carbon or graphite fibers embedded in a carbonized or graphitized matrix. In some applications, such as shielding sensitive electronics in very high temperature environments, the performance of these materials would be improved by lowering their electrical resistivity. One method to lower the resistivity of the composites is to lower the resistivity of the graphite fibers, and a proven method to accomplish that is intercalation. Intercalation is the insertion of guest atoms or molecules into a host lattice. In this study the host fibers were highly graphitic pitch-based graphite fibers, or vapor-grown carbon fibers (VGCF), and the intercalate was bromine. Intercalation compounds of graphite are generally thought of as being only metastable, but it has been shown that the residual bromine graphite fiber intercalation compound is remarkably stable, resisting decomposition even at temperatures at least as high as 1000 C. The focus of this work was to fabricate composite preforms, determine whether the fibers they were made from were still intercalated with bromine after processing, and determine the effect on composite resistivity. It was not expected that the resistivity would be lowered as dramatically as with graphite polymer composites because the matrix itself would be much more conductive, but it was hoped that the gains would be substantial enough to warrant its use in high-performance applications. In a collaborative effort supporting a Space Act Agreement between the NASA Glenn Research Center and Applied Sciences, Inc. (Cedarville, OH), laminar preforms were fabricated with pristine and bromine-intercalated pitch-based fibers (P100 and P100-Br) and VGCF (Pyro I and Pyro I-Br). The green preforms were carbonized at 1000 C and then heat treated to 3000 C. To determine whether the fibers in the samples were still intercalated after composite fabrication, they were subjected to X-ray diffraction. The composites containing intercalated graphite fibers showed much higher background scatter than that of pristine fibers, indicating the presence of bromine in the samples. More importantly, faint features indicative of intercalation were visible in the diffraction pattern, showing that the fibers were still intercalated.

Gaier, James R.↗

Structure of correlation tensors in homogeneous anisotropic turbulence

The theory of isotropic tensors, developed by Robertson (1940), Batchelor (1946), and Chandrasekhar (1951), is extended to cover the general case of turbulence with a pseudo-vector-preferred direction, without assuming mirror-reflection invariance. Attention is focused on two-point-correlation functions, and it is shown that the form of the decomposition into proper and pseudo-tensor contributions is restricted by the homogeneity requirement. The somewhat unexpected result that the vector- and pseudo-vector-preferred-direction cases yield different results is presented: A pseudo-vector-preferred direction allows the correlation matrix one more functional degree of freedom than does the proper vector case. An explicit form of the two-point-correlation tensor in the presence of a uniform mean magnetic field is presented which may be appropriate for use in analysis of magnetic fluctuations in plasma containment devices or the interplanetary medium. A procedure for determining the magnetic helicity from experimental data is presented.

Matthaeus, W. H.↗

Optimization of Time-Dependent Particle Tracing Using Tetrahedral Decomposition

An efficient algorithm is presented for computing particle paths, streak lines and time lines in time-dependent flows with moving curvilinear grids. The integration, velocity interpolation and step-size control are all performed in physical space which avoids the need to transform the velocity field into computational space. This leads to higher accuracy because there are no Jacobian matrix approximations or expensive matrix inversions. Integration accuracy is maintained using an adaptive step-size control scheme which is regulated by the path line curvature. The problem of cell-searching, point location and interpolation in physical space is simplified by decomposing hexahedral cells into tetrahedral cells. This enables the point location to be done analytically and substantially faster than with a Newton-Raphson iterative method. Results presented show this algorithm is up to six times faster than particle tracers which operate on hexahedral cells yet produces almost identical particle trajectories.

Kenwright, David↗

Research on output feedback control of systems with ill-conditioned dynamics

The accuracy of a two time scale approach to the output feedback regulator design problem is examined. An approximate quadratic performance index which reflects a two time scale decomposition of the system dynamics is developed. A sequential numerical algorithm is defined which obtains output feedback gains minimizing a broad class of performance indices, including the standard LQ case. A procedure for optimally zeroing selected gain elements in an output feedback gain matrix is developed and demonstrated. A summary of conference and journal publications from this research is also provided.

Calise, A. J.↗

Evaluation of constraint stabilization procedures for multibody dynamical systems

Comparative numerical studies of four constraint treatment techniques for the simulation of general multibody dynamic systems are presented, and results are presented for the example of a classical crank mechanism and for a simplified version of the seven-link manipulator deployment problem. The staggered stabilization technique (Park, 1986) is found to yield improved accuracy and robustness over Baumgarte's (1972) technique, the singular decomposition technique (Walton and Steeves, 1969), and the penalty technique (Lotstedt, 1979). Furthermore, the staggered stabilization technique offers software modularity, and the only data each solution module needs to exchange with the other is a set of vectors plus a common module to generate the gradient matrix of the constraints, B.

Park, K. C.↗

Biobased Polybenzoxazine Derived from Furfurylamine and Piceol with Low Cure Temperature and Advanced Properties for Composite Matrices

A novel benzoxazine made from furfurylamine, paraformaldehyde, and piceol, Bz-FA-HA, is assessed for applications in fiber reinforced (FR) composites. Piceol is a biobased phenolic compound derived from the roots of Norwegian spruce trees and contains a methyl ketone group at the para position. Bz-FA-HA is a liquid at room temperature, has a viscosity of < 1 Pa.s at temperatures above 90 °C, and a Tonset of cure at 132 °C. The carbonyl is found to react with the furan ring, yielding a crosslinking reaction, when cured above 180 °C as indicated by differential scanning calorimetry and thermogravimetric analysis coupled with Fourier transform infrared spectroscopy. Poly(Bz-FA-HA) has a Tg > 350 °C, attributed to the crosslinking reaction. Furthermore, the storage modulus is > 3 GPa, regardless of cure temperature. Poly(Bz-FA-HA) has a char yield at 800 °C of 65.0 % (62.3 % at 1000 °C), and a Tonset of decomposition of 357 °C in nitrogen. The resulting carbon formed during pyrolysis shrinks during the carbonization reaction and scanning electron microscopy imaging shows a cross-section with micro cracks. The high processability, advanced mechanical properties, and exceptional char yield make it a promising candidate as the matrix for FR composites.

Benzoxazine↗

Solving periodic block tridiagonal systems using the Sherman-Morrison-Woodbury formula

Many algorithms for solving the Navier-Stokes equations require the solution of periodic block tridiagonal systems of equations. By applying a splitting to the matrix representing this system of equations, it may first be reduced to a block tridiagonal matrix plus an outer product of two block vectors. The Sherman-Morrison-Woodbury formula is then applied. The algorithm thus reduces a periodic banded system to a non-periodic banded system with additional right-hand sides and is of higher efficiency than standard Thomas algorithm/LU decompositions.

Yarrow, Maurice↗

Inversion of imaging spectrometry data using singular value decomposition

The use of imaging spectrometers, which acquire data that are both spectrally contiguous images and spatially contiguous spectra, for quantitative remote sensing of the earth is addressed. Such data sets cannot be analyzed fully using either existing spectroscopic or image techniques. Singular value decomposition (SVD) is used here for spectral unmixing and determination of the spatial scales of mixing. It is shown that when it is used to invert the mixing endmember library, SVD allows more insight into library characteristics and more control of the inversion process than other commonly used matrix inversion techniques.

Boardman, Joe W.↗