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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

On the stability of the decelerating laminar boundary layer

The stability of a decelerating boundary-layer flow is investigated experimentally and numerically. Experimentally, a flat plate having a Blasius boundary layer is decelerated in an 18 m towing tank. The boundary layer becomes unstable to two-dimensional waves, which break down into three-dimensional patterns, hairpin vortices, and finally turbulent bursts when the vortices lift off the wall. The unsteady boundary-layer equations are solved numerically to generate instantaneous velocity profiles for a range of boundary and initial conditions. A quasi-steady approximation is invoked and the stability of local velocity profiles is determined by solving the Orr-Sommerfeld equation using Chebyshev matrix methods. Comparisons are made between the numerical predictions and the experimentally observed instabilities.

Gad-El-hak, M.↗
NASA NTRS1982

Scout trajectory error propagation computer program

Since 1969, flight experience has been used as the basis for predicting Scout orbital accuracy. The data used for calculating the accuracy consists of errors in the trajectory parameters (altitude, velocity, etc.) at stage burnout as observed on Scout flights. Approximately 50 sets of errors are used in Monte Carlo analysis to generate error statistics in the trajectory parameters. A covariance matrix is formed which may be propagated in time. The mechanization of this process resulted in computer program Scout Trajectory Error Propagation (STEP) and is described herein. Computer program STEP may be used in conjunction with the Statistical Orbital Analysis Routine to generate accuracy in the orbit parameters (apogee, perigee, inclination, etc.) based upon flight experience.

Myler, T. R.↗
NASA NTRS1988

PIA update: Correlation analyses of mass spectra

The PIA instrument aboard the Giotto spacecraft (a time of flight spectrometer) has been presented elsewhere. The mass spectra used in this analysis were decoded and mass numbers assigned according to the presence of carbon and silver, using the global values for these elements in their spectral absence. The results presented here were obtained using a frequency of occurrence based on analysis which correlated how often mass numbers appear in the mass spectra and which mass numbers tend to occur together in the same spectra; no amplitude information is utilized. The data are presented as plots of mass vs coincident mass for different subsets of the PIA data set, with both axes having units of atomic mass. Frequency contours are plotted at approximately five percent contour intervals, relative to the maximum AMU occurrence in that plot. The plots presented are symmetrical about the matrix diagonal, i.e., every mass is coincident with itself in a given spectra.

Mason, L. W.↗
NASA NTRS1990

A review of reduced Navier-Stokes computations for compressible viscous flows

A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.

Rubin, S. G.↗
NASA NTRS1989

Low-order design and high-order simulation of active closed-loop control for aerospace structures under construction

Partially constructed/assembled structures in space are complicated enough but their dynamics will also be operating in closed-loop with feedback controllers. The dynamics of such structures are modeled by large-scale finite element models. The model dimension L is extremely large (approximately 10,000) while the numbers of actuators (M) and sensors (P) are small. The model parameters M(sub m) mass matrix, D(sub o) damping matrix, and K(sub o) stiffness matrix, are all symmetric and sparse (banded). Thus simulation of open-loop structure models of very large dimension can be accomplished by special integration techniques for sparse matrices. The problem of simulation of closed-loop control of such structures is complicated by the addition of controllers. Simulation of closed-loop controlled structures is an essential part of the controller design and evaluation process. Current research in the following areas is presented: high-order simulation of actively controlled aerospace structures; low-order controller design and SCI compensation for unmodeled dynamics; prediction of closed-loop stability using asymptotic eigenvalue series; and flexible robot manipulator control experiment.

Balas, Mark J.↗
NASA NTRS1995

Krylov methods for compressible flows

We investigate the application of Krylov methods to compressible flows, and the effect of implicit boundary conditions on the implicit solution of nonlinear problems. Two defect-correction procedures, namely, approximate factorization (AF) for structured grids and ILU/GMRES for general grids, are considered. Also considered here are Newton-Krylov matrix-free methods that we combined with the use of mixed discretization schemes in the implicitly defined Jacobian and its preconditioner. Numerical experiments that show the performance of our approaches are then presented.

Tidriri, M. D.↗
NASA NTRS2003

The Importance of Chaos and Lenticulae on Europa for the JIMO Mission

The Galileo Solid State Imaging (SSI) experiment provided high-resolution images of Europa's surface allowing identification of surface features barely distinguishable at Voyager's resolution. SSI revealed the visible pitting on Europa's surface to be due to large disrupted features, chaos, and smaller sub-circular patches, lenticulae. Chaos features contain a hummocky matrix material and commonly contain dislocated blocks of ridged plains. Lenticulae are morphologically interrelated and can be divided into three classes: domes, spots, and micro-chaos. Domes are broad, upwarped features that generally do not disrupt the texture of the ridged plains. Spots are areas of low albedo that are generally smooth in texture compared to other units. Micro-chaos are disrupted features with a hummocky matrix material, resembling that observed within chaos regions. Chaos and lenticulae are ubiquitous in the SSI regional map observations, which average approximately 200 meters per pixel (m/pxl) in resolution, and appear in several of the ultra-high resolution, i.e., better than 50 m/pxl, images of Europa as well. SSI also provided a number of multi-spectral observations of chaos and lenticulae. Using this dataset we have undertaken a thorough study of the morphology, size, spacing, stratigraphy, and color of chaos and lenticulae to determine their properties and evaluate models of their formation. Geological mapping indicates that chaos and micro-chaos have a similar internal morphology of in-situ degradation suggesting that a similar process was operating during their formation. The size distribution denotes a dominant size of 4-8 km in diameter for features containing hummocky material (i.e., chaos and micro-chaos). Results indicate a dominant spacing of 15 - 36 km apart. Chaos and lenticulae are generally among the youngest features stratigraphically observed on the surface, suggesting a recent change in resurfacing style. Also, the reddish non-icy materials on Europa's surface have high concentrations in many chaos and lenticulae features. Nonetheless, a complete global map of the distribution of chaos and lenticulae is not possible with the SSI dataset. Only <20% of the surface has been imaged at 200 m/pxl or better resolution, mostly of the near-equatorial regions. Color and ultra-high-res images have much less surface coverage. Thus we suggest that full global imaging of Europa at 200 m/pxl or better resolution, preferably in multi-spectral wavelengths, should be a high priority for the JIMO mission.

Spaun, Nicole A.↗
NASA NTRS2014

Direct Demonstration of the Concept of Unrestricted Effective-Medium Approximation

The modified unrestricted effective-medium refractive index is defined as one that yields accurate values of a representative set of far-field scattering characteristics (including the scattering matrix) for an object made of randomly heterogeneous materials. We validate the concept of the modified unrestricted effective-medium refractive index by comparing numerically exact superposition T-matrix results for a spherical host randomly filled with a large number of identical small inclusions and Lorenz-Mie results for a homogeneous spherical counterpart. A remarkable quantitative agreement between the superposition T-matrix and Lorenz-Mie scattering matrices over the entire range of scattering angles demonstrates unequivocally that the modified unrestricted effective-medium refractive index is a sound (albeit still phenomenological) concept provided that the size parameter of the inclusions is sufficiently small and their number is sufficiently large. Furthermore, it appears that in cases when the concept of the modified unrestricted effective-medium refractive index works, its actual value is close to that predicted by the Maxwell-Garnett mixing rule.

Scattering↗
NASA NTRS1987

Experimental evaluation of longitudinal splitting in unidirectional composites

An experimental study was conducted to determine the fracture behavior of centernotched, unidirectional graphite/epoxy laminates when subjected to tensile loading. The actual behavior is compared to the behavior predicted by a mathematical model based on classical shear-lag assumptions. Experimental results show the fracture behavior to consist of two modes of longitudinal matrix splitting. After split initiation, a region of slow, stable split growth occurs. This is followed by a distinct transition into a rapid split growth region. However, the shear-lag model predicts a single region of rapid longitudinal split growth behavior with unbounded splitting occurring after approximately a ten percent increase in the applied stress over the split initiation stress. It is therefore concluded that the classical shear-lag model does not accurately predict the complete longitudinal split growth behavior in a brittle matrix composite. This shortcoming in the model and some results obtained from related elasticity solutions are discussed and it is suggested that the slow, stable split growth region is likely dominated by transverse matrix normal stresses which are not accounted for in the classical shear-lag model.

Wolla, Jeffrey M.↗
NASA NTRS2013

Quantifying Effects of Voids in Woven Ceramic Matrix Composites

Randomness in woven ceramic matrix composite architecture has been found to cause large variability in stiffness and strength. The inherent voids are an aspect of the architecture that may cause a significant portion of the variability. A study is undertaken to investigate the effects of many voids of random sizes and distributions. Response surface approximations were formulated based on void parameters such as area and length fractions to provide an estimate of the effective stiffness. Obtaining quantitative relationships between the properties of the voids and their effects on stiffness of ceramic matrix composites are of ultimate interest, but the exploratory study presented here starts by first modeling the effects of voids on an isotropic material. Several cases with varying void parameters were modeled which resulted in a large amount of variability of the transverse stiffness and out-of-plane shear stiffness. An investigation into a physical explanation for the stiffness degradation led to the observation that the voids need to be treated as an entity that reduces load bearing capabilities in a space larger than what the void directly occupies through a corrected length fraction or area fraction. This provides explanation as to why void volume fraction is not the only important factor to consider when computing loss of stiffness.

Goldsmith, Marlana B.↗
NASA NTRS1986

Preconditioning matrices for Chebyshev derivative operators

The problem of preconditioning the matrices arising from pseudo-spectral Chebyshev approximations of first order operators is considered in both one and two dimensions. In one dimension a preconditioner represented by a full matrix which leads to preconditioned eigenvalues that are real, positive, and lie between 1 and pi/2, is already available. Since there are cases in which it is not computationally convenient to work with such a preconditioner, a large number of preconditioners were studied which were more sparse (in particular three and four diagonal matrices). The eigenvalues of such preconditioned matrices are compared. The results were applied to the problem of finding the steady state solution to an equation of the type u sub t = u sub x + f, where the Chebyshev collocation is used for the spatial variable and time discretization is performed by the Richardson method. In two dimensions different preconditioners are proposed for the matrix which arises from the pseudo-spectral discretization of the steady state problem. Results are given for the CPU time and the number of iterations using a Richardson iteration method for the unpreconditioned and preconditioned cases.

Rothman, Ernest E.↗
NASA NTRS1993

Control design variable linking for optimization of structural/control systems

A method is presented to integrate the design space of structural/control system optimization problems in the case of linear state feedback control. Conventional structural sizing variables and elements of the feedback gain matrix are both treated as strictly independent design variables in optimization by extending design variable linking concepts to the control gains. Several approximation concepts including new control design variable linking schemes are used to formulate the integrated structural/control optimization problem as a sequence of explicit nonlinear mathematical programming problems. Examples which involve a variety of behavior constraints, including constraints on dynamic stability, damped frequencies, control effort, peak transient displacement, acceleration, and control force limits, are effectively solved by using the method presented.

Jin, Ik Min↗
NASA NTRS2021

Thermal Inspection of a Composite Fuselage Section Using theMethod of Proper Generalized Decomposition

Proper Generalized Decomposition (PGD) is a reduced order modeling technique for the simulation of physical systems whose governing equations depend on boundary conditions, initial conditions, material properties, and geometric parameters. It uses separated representations of system covariates combined with an iterative approximation method known as successive enrichment in order to compute an accurate parameter-dependent approximation to the full governing equations. PGD can also be used as an alternative to the Singular Value Decomposition (SVD) of a matrix and therefore as an alternative to PCA thermography. In this paper PGD was used to analyze data derived from the inspection of a composite fuselage forward section using flash thermography, and the results were compared against the standard PCA approach.

Nondestructive Evaluation↗
NASA NTRS1991

CRRES microelectronic test chip

The JPL CRRES chip was designed and fabricated in 1985 and included in the CRRES MEP. MOSFET Matrix results show the effect of shielding on radiation-induced MOSFET threshold voltage shifts and channel mobility degradation. Shielded (middle board) MOSFETs have a threshold-voltage damage factor that is approximately three orders of magnitude smaller than would be estimated from Co-60 ground tests. Unshielded (outer board) MOSFETs have a threshold-voltage damage factor that would be estimated from Co-60 ground tests. Temperature swings as large as 23 C with a 22.5 orbit periodicity affected the MOSFET data and were removed from the data in order to reveal the radiation effects. This experiment demonstrated the feasibility of characterizing MOSFETs in a matrix, thus reducing the complexity and mass of the experiment.

Lin, Y.-S.↗
NASA NTRS1978

On least squares approximations to indefinite problems of the mixed type

A least squares method is presented for computing approximate solutions of indefinite partial differential equations of the mixed type such as those that arise in connection with transonic flutter analysis. The method retains the advantages of finite difference schemes namely simplicity and sparsity of the resulting matrix system. However, it offers some great advantages over finite difference schemes. First, the method is insensitive to the value of the forcing frequency, i.e., the resulting matrix system is always symmetric and positive definite. As a result, iterative methods may be successfully employed to solve the matrix system, thus taking full advantage of the sparsity. Furthermore, the method is insensitive to the type of the partial differential equation, i.e., the computational algorithm is the same in elliptic and hyperbolic regions. In this work the method is formulated and numerical results for model problems are presented. Some theoretical aspects of least squares approximations are also discussed.

Fix, G. J.↗
NASA NTRS1993

The effects of long-duration space exposure on the mechanical properties of some carbon-reinforced resin matrix composites

Long Duration Exposure Facility (LDEF) Experiment A0175 involved the non-instrumented exposure of seven carbon-fiber reinforced resin-matrix advanced composite panels contained in two trays - A7 and A1. These two trays were located, respectively, on the leading and trailing faces of LDEF, obliquely oriented to the RAM (Row 9) and WAKE (Row 3) directions. The identity and location of the seven panels, which consisted of six flat laminates of the following material systems are shown: carbon/epoxy (T300/934), carbon/bismaleimide (T300/F178), and carbon/polyimide (C6000/LARC-160 and C6000/PMR-15), plus one bonded honeycomb sandwich panel (T300/934 face sheets and Nomex core) patterned after the Space Shuttle payload bay door construction. These material systems were selected to represent a range of then-available matrix resins which, by virtue of their differing polymer chemistry, could conceivably exhibit differing susceptibility to the low-earth orbit (LEO) environment. The principal exposure conditions of the LDEF environment at these tray locations are shown. Noteworthy to some of the observations discussed is the four-orders-of magnitude difference in the atomic oxygen (AO) fluence, which made a shallow incidence angle (approximately 22 deg) to Tray A7, while Tray A1 on the trailing face was essentially shielded from AO exposure. This evaluation focused on determining the individual and relative suitability of a variety of resin-matrix composite systems for long-term space structural applications. This was accomplished primarily by measuring and comparing a range of engineering mechanical properties on over 300 test coupons sectioned from the flight panels and from identical control panels, and tested at ambient and elevated temperatures. This testing was supported by limited physical characterization, involving visual examination of flight panel surface features, measurements of weight loss and warpage, and examination for changes in internal integrity (micro cracking, delamination) by ultrasonic c-scan and polished cross-sections.

Vyhnal, Richard F.↗
NASA NTRS1992

Rational approximations from power series of vector-valued meromorphic functions

Let F(z) be a vector-valued function, F: C yields C(sup N), which is analytic at z = 0 and meromorphic in a neighborhood of z = 0, and let its Maclaurin series be given. In this work we developed vector-valued rational approximation procedures for F(z) by applying vector extrapolation methods to the sequence of partial sums of its Maclaurin series. We analyzed some of the algebraic and analytic properties of the rational approximations thus obtained, and showed that they were akin to Pade approximations. In particular, we proved a Koenig type theorem concerning their poles and a de Montessus type theorem concerning their uniform convergence. We showed how optical approximations to multiple poles and to Laurent expansions about these poles can be constructed. Extensions of the procedures above and the accompanying theoretical results to functions defined in arbitrary linear spaces was also considered. One of the most interesting and immediate applications of the results of this work is to the matrix eigenvalue problem. In a forthcoming paper we exploited the developments of the present work to devise bona fide generalizations of the classical power method that are especially suitable for very large and sparse matrices. These generalizations can be used to approximate simultaneously several of the largest distinct eigenvalues and corresponding eigenvectors and invariant subspaces of arbitrary matrices which may or may not be diagonalizable, and are very closely related with known Krylov subspace methods.

Sidi, Avram↗
NASA NTRS2017

Relaxation Matrix for Symmetric Tops with Inversion Symmetry: Line Coupling and Line Mixing Effects on NH3 Lines in the V4 Band

Line shape parameters including the half-widths and the off-diagonal elements of the relaxation matrix have been calculated for self-broadened NH3 lines in the perpendicular v4 band. As in the pure rotational and the parallel v1 bands, the small inversion splitting in this band causes a complete failure of the isolated line approximation. As a result, one has to use formalisms not relying on this approximation. However, due to differences between parallel and perpendicular bands of NH3, the applicability of the formalism used in our previous studies of the v1 band and other parallel bands must be carefully verified. We have found that, as long as potential models only contain components with K1 equals K2 equals 0, whose matrix elements require the selection rule delta k equals 0, the formalism is applicable for the v4 band with some minor adjustments. Based on both theoretical considerations and results from numerical calculations, the non-diagonality of the relaxation matrices in all the PP, RP, PQ, RQ, PR, and RR branches is discussed. Theoretically calculated self-broadened half-widths are compared with measurements and the values listed in HITRAN 2012. With respect to line coupling effects, we have compared our calculated intra-doublet off-diagonal elements of the relaxation matrix with reliable measurements carried out in the PP branch where the spectral environment is favorable. The agreement is rather good since our results do well reproduce the observed k and j dependences of these elements, thus validating our formalism.

Absorption coefficient↗