Search NASA⌕ Search

SEARCH · Search NASA

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 19 records

A study of the effects of state transition matrix approximations

The effects of using an approximate state transition matrix in orbit estimation are investigated. The approximate state transition matrix results when higher order geopotential terms in the equations of motion are ignored in the formation of the variational equations. Two methods of orbit estimation are considered: the differential correction procedure (DC) and the extended Kalman filter (EKF). The system used for the study is the Research & Development version of the Goddard Trajectory Determination System. The effects of the approximation are analyzed on a number of orbits. These orbits of various inclinations and semimajor axes. Other parameters studied include geopotential models and DC arc length.

May, J. A.↗

Optimal matrix approximants in structural identification

Problems of model correlation and system identification are central in the design, analysis, and control of large space structures. Of the numerous methods that have been proposed, many are based on finding minimal adjustments to a model matrix sufficient to introduce some desirable quality into that matrix. In this work, several of these methods are reviewed, placed in a modern framework, and linked to other previously known ideas in computational linear algebra and optimization. This new framework provides a point of departure for a number of new methods which are introduced here. Significant among these is a method for stiffness matrix adjustment which preserves the sparsity pattern of an original matrix, requires comparatively modest computational resources, and allows robust handling of noisy modal data. Numerical examples are included to illustrate the methods presented herein.

Beattie, C. A.↗

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↗

Low-dimensional Representation of Error Covariance

Ensemble and reduced-rank approaches to prediction and assimilation rely on low-dimensional approximations of the estimation error covariances. Here stability properties of the forecast/analysis cycle for linear, time-independent systems are used to identify factors that cause the steady-state analysis error covariance to admit a low-dimensional representation. A useful measure of forecast/analysis cycle stability is the bound matrix, a function of the dynamics, observation operator and assimilation method. Upper and lower estimates for the steady-state analysis error covariance matrix eigenvalues are derived from the bound matrix. The estimates generalize to time-dependent systems. If much of the steady-state analysis error variance is due to a few dominant modes, the leading eigenvectors of the bound matrix approximate those of the steady-state analysis error covariance matrix. The analytical results are illustrated in two numerical examples where the Kalman filter is carried to steady state. The first example uses the dynamics of a generalized advection equation exhibiting nonmodal transient growth. Failure to observe growing modes leads to increased steady-state analysis error variances. Leading eigenvectors of the steady-state analysis error covariance matrix are well approximated by leading eigenvectors of the bound matrix. The second example uses the dynamics of a damped baroclinic wave model. The leading eigenvectors of a lowest-order approximation of the bound matrix are shown to approximate well the leading eigenvectors of the steady-state analysis error covariance matrix.

Tippett, Michael K.↗

Organic Matter in Extraterrestrial Water-Bearing Salt Crystals

Introduction: Direct samples of early Solar System fluids are present in two thermally-metamorphosed ordinary chondrite regolith breccias (Monahans (1998) [H5] and Zag [H3-6]), which were found to contain brine-bearing halite (NaCl) crystals that have been added to the regolith of an S-type asteroid following asteroidal metamorphism [1, 2]. The brine-bearing halite grains were proposed to be formed on an icy C-type asteroids (possibly Ceres), and transferred to an S-type asteroid via cryovolcanic event(s) [3]. A unique aspect of these halites is that they contain abundant organic rich solid inclusions hosted within the halites alongside the water inclusions. Methods: We analyzed in detail the compositions of the organic solids and the amino acid content of the halite crystals with two-step laser desorption/laser ionization mass spectrometry (L(sup 2) MS), Raman spectroscopy, X-ray absorption near edge structure (XANES), nanoscale secondary ion mass spectrometry (NanoSIMS), and ultra-performance liquid chromatography fluorescence detection and quadrupole time of flight hybrid mass spectrometry (UPLC-FD/QToF-MS). Results and Discussion: The L(sup 2) MS results show signatures of low-mass polyaromatic hydro-carbons (PAHs) indicated by sequences of peaks separated by 14 atomic mass units (amu) due to successive addition of methylene (CH2) groups to the PAH skeletons [4]. Raman spectra of the micron-sized solid inclusions of the halites indicate the presence of abundant and highly variable organic matter that include a mixture of short-chain aliphatic compounds and macromolecular carbon. C-XANES analysis identified C-rich areas with peaks at 285.0 eV (aromatic C=C) and 286.6 eV (vinyl-keto C=O). However, there is no 1s-sigma* exciton peak (291.7 eV) that is indicative of the development of graphene structure [5], which suggests the organics were synthesized cold. Na-noSIMS analyses show C-rich and N-rich areas that exhibit similar isotopic values with that of the IOM in the unweathered CR chondrites and less metamorphosed meteorites [6], and are moderately enriched in N-15 (delta N-15 = 106.1-164.5 per mille). The total amino acid distribution and abundance of the Zag matrix (approximately 1,940 parts per billion [ppb]) is comparable to other ordinary chondrites (60-3,330 ppb) [7, 8]. While the Zag matrix is gamma-ABA and EACA-deficient, the halite is shown to exhibit an opposite trend and is almost depleted in amino acids. The striking difference in the amino acid contents between the halite and matrix indicates their separate synthetic origins. Conclusion: Abundant, primitive, and highly-diverse N-15-rich organic compounds were detected in brine-water bearing halite crystals that were synthesized on a cryovolcanically-active asteroid. Our study suggests that the asteroidal parent body where the halite precipitated, potentially Ceres, is a host to abundance large variety organic precursors. Insoluble organic matter and amino acids can be synthesized from similar organic precusors under hydrous conditions [9].We envision that similar organic synthetic processes could have occurred on Ceres that synthesized organic solids as well as biologically relevant molecules.

Chan, Q. H. S.↗

Studies of mineralization in tissue culture: optimal conditions for cartilage calcification

The optimal conditions for obtaining a calcified cartilage matrix approximating that which exists in situ were established in a differentiating chick limb bud mesenchymal cell culture system. Using cells from stage 21-24 embryos in a micro-mass culture, at an optimal density of 0.5 million cells/20 microliters spot, the deposition of small crystals of hydroxyapatite on a collagenous matrix and matrix vesicles was detected by day 21 using X-ray diffraction, FT-IR microscopy, and electron microscopy. Optimal media, containing 1.1 mM Ca, 4 mM P, 25 micrograms/ml vitamin C, 0.3 mg/ml glutamine, no Hepes buffer, and 10% fetal bovine serum, produced matrix resembling the calcifying cartilage matrix of fetal chick long bones. Interestingly, higher concentrations of fetal bovine serum had an inhibitory effect on calcification. The cartilage phenotype was confirmed based on the cellular expression of cartilage collagen and proteoglycan mRNAs, the presence of type II and type X collagen, and cartilage type proteoglycan at the light microscopic level, and the presence of chondrocytes and matrix vesicles at the EM level. The system is proposed as a model for evaluating the events in cell mediated cartilage calcification.

Non-NASA Center↗

An approximate loop transfer recovery method for designing fixed-order compensators

A method for designing fixed-order dynamic compensators for multivariable time invariant linear systems is presented which is based on the minimization of a linear quadratic performance index. The present formulation is performed in an output feedback setting which uses an observer cononical form to represent the compensator dynamics. Techniques for penalizing the plant and compensator states and for selecting the distribution on initial conditions such that the loop transfer matrix approximates that of a full-state feedback design have been developed. The effectiveness of the method is demonstrated using the examples of the pointing of a flexible structure and a helicopter flight control problem.

Calise, Anthony J.↗

The recovery of polarized apparent temperature distributions of flat scenes from antenna temperature measurements

The reception theory for quasi-monochromatic partially polarized emissions is employed to describe the antenna temperature when observing a flat scene from within or above an absorbing/emitting atmosphere. When the polarization aspects of the antenna and scene are carefully observed, the resulting integral expression for the antenna temperature differs significantly from that reported for polarization invariant temperature distributions. It is observed that the antenna couples to both emitted surface polarizations, and the phase character of the antenna introduces a term involving the product of the polarized and cross-polarized antenna patterns. A matrix approximation of the integral expression is suggested and developed as a means for inverting polarized antenna temperature observations for the underlying partially polarized apparent temperature distributions.

Claassen, J. P.↗

Robust fixed order dynamic compensation for large space structure control

A simple formulation for designing fixed order dynamic compensators which are robust to both uncertainty at the plant input and structured uncertainty in the plant dynamics is presented. The emphasis is on designing low order compensators for systems of high order. The formulation is done in an output feedback setting which exploits an observer canonical form to represent the compensator dynamics. The formulation also precludes the use of direct feedback of the plant output. The main contribution lies in defining a method for penalizing the states of the plant and of the compensator, and for choosing the distribution on initial conditions so that the loop transfer matrix approximates that of a full state design. To improve robustness to parameter uncertainty, the formulation avoids the introduction of sensitivity states, which has led to complex formulations in earlier studies where only structured uncertainty has been considered.

Calise, Anthony J.↗

Application of vector-valued rational approximations to the matrix eigenvalue problem and connections with Krylov subspace methods

Let F(z) be a vectored-valued function F: C approaches C sup N, which is analytic at z=0 and meromorphic in a neighborhood of z=0, and let its Maclaurin series be given. We use vector-valued rational approximation procedures for F(z) that are based on its Maclaurin series in conjunction with power iterations to develop bona fide generalizations of the power method for an arbitrary N X N matrix that may be diagonalizable or not. These generalizations can be used to obtain simultaneously several of the largest distinct eigenvalues and the corresponding invariant subspaces, and present a detailed convergence theory for them. In addition, it is shown that the generalized power methods of this work are equivalent to some Krylov subspace methods, among them the methods of Arnoldi and Lanczos. Thus, the theory provides a set of completely new results and constructions for these Krylov subspace methods. This theory suggests at the same time a new mode of usage for these Krylov subspace methods that were observed to possess computational advantages over their common mode of usage.

Sidi, Avram↗

Basis Function Approximation of Transonic Aerodynamic Influence Coefficient Matrix

A technique for approximating the modal aerodynamic influence coefficients [AIC] matrices by using basis functions has been developed and validated. An application of the resulting approximated modal AIC matrix for a flutter analysis in transonic speed regime has been demonstrated. This methodology can be applied to the unsteady subsonic, transonic and supersonic aerodynamics. The method requires the unsteady aerodynamics in frequency-domain. The flutter solution can be found by the classic methods, such as rational function approximation, k, p-k, p, root-locus et cetera. The unsteady aeroelastic analysis for design optimization using unsteady transonic aerodynamic approximation is being demonstrated using the ZAERO(TradeMark) flutter solver (ZONA Technology Incorporated, Scottsdale, Arizona). The technique presented has been shown to offer consistent flutter speed prediction on an aerostructures test wing [ATW] 2 configuration with negligible loss in precision in transonic speed regime. These results may have practical significance in the analysis of aircraft aeroelastic calculation and could lead to a more efficient design optimization cycle

Li, Wesley Waisang↗

State-space formulations for flutter analysis

Two state-space methods for true-damping flutter analysis are presented here. The aerodynamic matrix is assumed to be function of the complex reduced frequency p (off-the-imaginary axis aerodynamics) as the ones provided for instance by the computer code SOUSSA (Steady, Oscillatory and Unsteady, Subsonic and Supersonic Aerodynamics). This matrix is approximated by a third order polynomial or by a Pade's matrix approximant. Then the equations are cast into a suitable state space format including the Laplace transform of the space vector, a frequency-independent matrix, and the flight speed. The solution is then obtained by standard root-locus methods. Results are in good agreement with existing ones in a suitable range of the flutter reduced frequency. A discussion on spurious roots is also included.

Weiss, S. J.↗

Basis Function Approximation of Transonic Aerodynamic Influence Coefficient Matrix

A technique for approximating the modal aerodynamic influence coefficients matrices by using basis functions has been developed and validated. An application of the resulting approximated modal aerodynamic influence coefficients matrix for a flutter analysis in transonic speed regime has been demonstrated. This methodology can be applied to the unsteady subsonic, transonic, and supersonic aerodynamics. The method requires the unsteady aerodynamics in frequency-domain. The flutter solution can be found by the classic methods, such as rational function approximation, k, p-k, p, root-locus et cetera. The unsteady aeroelastic analysis for design optimization using unsteady transonic aerodynamic approximation is being demonstrated using the ZAERO flutter solver (ZONA Technology Incorporated, Scottsdale, Arizona). The technique presented has been shown to offer consistent flutter speed prediction on an aerostructures test wing 2 configuration with negligible loss in precision in transonic speed regime. These results may have practical significance in the analysis of aircraft aeroelastic calculation and could lead to a more efficient design optimization cycle.

Li, Wesley W.↗

Approximate Cartesian state transition matrix

Based on a Taylor series expansion, an easily-computed approximation to the Cartesian state transition matrix is presented for a general velocity-independent force field. Suitable for the short time intervals encountered in onboard navigation applications of the extended Kalman filter, it provides approximately five decimal digits of accuracy for earth orbiting spacecraft with update intervals of one minute, and better accuracy for shorter intervals.

Markley, F. L.↗

Model reduction of nonsquare linear MIMO systems using multipoint matrix continued-fraction expansions

This paper deals with the multipoint Cauer matrix continued-fraction expansion (MCFE) for model reduction of linear multi-input multi-output (MIMO) systems with various numbers of inputs and outputs. A salient feature of the proposed MCFE approach to model reduction of MIMO systems with square transfer matrices is its equivalence to the matrix Pade approximation approach. The Cauer second form of the ordinary MCFE for a square transfer function matrix is generalized in this paper to a multipoint and nonsquare-matrix version. An interesting connection of the multipoint Cauer MCFE method to the multipoint matrix Pade approximation method is established. Also, algorithms for obtaining the reduced-degree matrix-fraction descriptions and reduced-dimensional state-space models from a transfer function matrix via the multipoint Cauer MCFE algorithm are presented. Practical advantages of using the multipoint Cauer MCFE are discussed and a numerical example is provided to illustrate the algorithms.

Guo, Tong-Yi↗

An analytical technique for approximating unsteady aerodynamics in the time domain

An analytical technique is presented for approximating unsteady aerodynamic forces in the time domain. The order of elements of a matrix Pade approximation was postulated, and the resulting polynomial coefficients were determined through a combination of least squares estimates for the numerator coefficients and a constrained gradient search for the denominator coefficients which insures stable approximating functions. The number of differential equations required to represent the aerodynamic forces to a given accuracy tends to be smaller than that employed in certain existing techniques where the denominator coefficients are chosen a priori. Results are shown for an aeroelastic, cantilevered, semispan wing which indicate a good fit to the aerodynamic forces for oscillatory motion can be achieved with a matrix Pade approximation having fourth order numerator and second order denominator polynomials.

Dunn, H. J.↗