Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “eigenvalue analysis”

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.

481 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

An investigative study of multispectral data compression for remotely-sensed images using vector quantization and difference-mapped shift-coding

A study is conducted to investigate the effects and advantages of data compression techniques on multispectral imagery data acquired by NASA's airborne scanners at the Stennis Space Center. The first technique used was vector quantization. The vector is defined in the multispectral imagery context as an array of pixels from the same location from each channel. The error obtained in substituting the reconstructed images for the original set is compared for different compression ratios. Also, the eigenvalues of the covariance matrix obtained from the reconstructed data set are compared with the eigenvalues of the original set. The effects of varying the size of the vector codebook on the quality of the compression and on subsequent classification are also presented. The output data from the Vector Quantization algorithm was further compressed by a lossless technique called Difference-mapped Shift-extended Huffman coding. The overall compression for 7 channels of data acquired by the Calibrated Airborne Multispectral Scanner (CAMS), with an RMS error of 15.8 pixels was 195:1 (0.41 bpp) and with an RMS error of 3.6 pixels was 18:1 (.447 bpp). The algorithms were implemented in software and interfaced with the help of dedicated image processing boards to an 80386 PC compatible computer. Modules were developed for the task of image compression and image analysis. Also, supporting software to perform image processing for visual display and interpretation of the compressed/classified images was developed.

Jaggi, S.↗
NASA NTRS1986

Landsat D Thematic Mapper image dimensionality reduction and geometric correction accuracy

To characterize and quantify the performance of the Landsat thematic mapper (TM), techniques for dimensionality reduction by linear transformation have been studied and evaluated and the accuracy of the correction of geometric errors in TM images analyzed. Theoretical evaluations and comparisons for existing methods for the design of linear transformation for dimensionality reduction are presented. These methods include the discrete Karhunen Loeve (KL) expansion, Multiple Discriminant Analysis (MDA), Thematic Mapper (TM)-Tasseled Cap Linear Transformation and Singular Value Decomposition (SVD). A unified approach to these design problems is presented in which each method involves optimizing an objective function with respect to the linear transformation matrix. From these studies, four modified methods are proposed. They are referred to as the Space Variant Linear Transformation, the KL Transform-MDA hybrid method, and the First and Second Version of the Weighted MDA method. The modifications involve the assignment of weights to classes to achieve improvements in the class conditional probability of error for classes with high weights. Experimental evaluations of the existing and proposed methods have been performed using the six reflective bands of the TM data. It is shown that in terms of probability of classification error and the percentage of the cumulative eigenvalues, the six reflective bands of the TM data require only a three dimensional feature space. It is shown experimentally as well that for the proposed methods, the classes with high weights have improvements in class conditional probability of error estimates as expected.

Ford, G. E.↗
NASA NTRS1990

Implicit application of polynomial filters in a k-step Arnoldi method

The Arnoldi process is a well known technique for approximating a few eigenvalues and corresponding eigenvectors of a general square matrix. Numerical difficulties such as loss of orthogonality and assessment of the numerical quality of the approximations as well as a potential for unbounded growth in storage have limited the applicability of the method. These issues are addressed by fixing the number of steps in the Arnoldi process at a prescribed value k and then treating the residual vector as a function of the initial Arnoldi vector. This starting vector is then updated through an iterative scheme that is designed to force convergence of the residual to zero. The iterative scheme is shown to be a truncation of the standard implicitly shifted QR-iteration for dense problems and it avoids the need to explicitly restart the Arnoldi sequence. The main emphasis of this paper is on the derivation and analysis of this scheme. However, there are obvious ways to exploit parallelism through the matrix-vector operations that comprise the majority of the work in the algorithm. Preliminary computational results are given for a few problems on some parallel and vector computers.

Sorensen, D. C.↗
NASA NTRS2003

Algorithmic Enhancements to the VULCAN Navier-Stokes Solver

VULCAN (Viscous Upwind aLgorithm for Complex flow ANalysis) is a cell centered, finite volume code used to solve high speed flows related to hypersonic vehicles. Two algorithms are presented for expanding the range of applications of the current Navier-Stokes solver implemented in VULCAN. The first addition is a highly implicit approach that uses subiterations to enhance block to block connectivity between adjacent subdomains. The addition of this scheme allows more efficient solution of viscous flows on highly-stretched meshes. The second algorithm addresses the shortcomings associated with density-based schemes by the addition of a time-derivative preconditioning strategy. High speed, compressible flows are typically solved with density based schemes, which show a high level of degradation in accuracy and convergence at low Mach numbers (M less than or equal to 0.1). With the addition of preconditioning and associated modifications to the numerical discretization scheme, the eigenvalues will scale with the local velocity, and the above problems will be eliminated. With these additions, VULCAN now has improved convergence behavior for multi-block, highly-stretched meshes and also can solve the Navier-Stokes equations for very low Mach numbers.

Litton, D. K.↗
NASA NTRS1987

A general method for dynamic analysis of structures overview

The presented research deals with the development of a dynamic analysis method for structural systems. The modeling approach is essentially a finite element method in the sense that the structure is divided into n elements. An element is defined as any structural unit whose degree of freedom (dofs) can be categorized as either interface or non-interface dofs. An element could be a fundamental unit such as a rod, a beam, a plate etc., or it could be an entire structural component. Furthermore, the parameters for the element could be distributed or lumped. The choice of elements is totally arbitrary and is a matter of user convenience. In particular, issues of accuracy and convergence do not enter on the level of example that bookkeeping is reduced to a minimum. Each element is modeled using a set of interface constraint modes (ICM) combined with a set of interface restrained normal models (IRNM). The next step is the solution of the system eigenvalue problem. The procedure calls for the sequential solution of a number of small eigenvalue problems based on a truncation principle for IRNM. In addition, the form of these eigenvalue problems is very simple such that an escalator type of eigenvalue problem solver can be used which is extremely cost-effective and fast.

Engels, Remi C.↗
NASA NTRS1976

On a numerical solution of the supersonic panel flutter eigenproblem

An automated digital computer procedure is presented in this paper which enables efficient solution of the eigenvalue problem associated with the supersonic panel flutter phenomena. The step-by-step incremental solution procedure is based on an inverse iteration technique which effectively utilizes solution results from the previous step in determining such results during the current solution step. Also, the computations are limited to the determination of a few specific roots only, which are expected to contain the flutter mode, and this is achieved at each step without having to compute any other root. The structural discretization achieved by the finite-element method yields highly banded stiffness, mass, and aerodynamic matrices; the aerodynamic matrix evaluated by the linearized piston theory is real but unsymmetric in nature. The solution algorithm presented in this paper fully exploits the banded form of the associated matrices, and the resulting computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most efficient and economical when compared to existing procedures of such analysis. Numerical results are presented for a two-dimensional panel flutter problem.

Gupta, K. K.↗
NASA NTRSDate not supplied

High-Fidelity Flight Dynamic Analysis of Transonic Truss-Braced Wing

This paper presents a high-fidelity flight dynamic analysis of the Mach 0.8 Transonic Truss-Braced Wing (TTBW). Unsteady RANS CFD simulations of the Mach 0.8 TTBW in pitch, plunge, roll, and yaw oscillations are conducted in FUN3D. The time-domain data are transformed into the frequency-domain data by Fourier series. Transfer functions of the dynamic stability derivatives are then estimated by a frequency-domain re- gression. The dynamic stability derivatives with respect to the angle of attack are determined by the regression of the unsteady aerodynamic coefficients for the plunge motion. The dynamic stability derivatives with re- spect to the pitch rate are determined by the regression of the differential unsteady aerodynamic coefficients for the pitch motion upon the removal of the angle of attack contribution by the plunge motion. Similarly, dynamic stability derivatives with respect to the angle of sideslip, roll rate, and yaw rates are determined from the frequency domain regression. The longitudinal and lateral-direction flight dynamic models of the Mach 0.8 TTBW are constructed from these dynamic stability derivatives. The eigenvalues of the aircraft modes are analyzed to determine the vehicle stability.

Aircraft Stability and Control↗
NASA NTRSDate not supplied

Stability Analysis of the Flow over a Swept Forward-Facing Step using PIV Base Flowsin a Non-Orthogonal Coordinate System

Understanding the flow of a crossflow-vortex-dominated boundary layer over a forward-facing step excrescence is necessary in order to mitigate an early transition scenario on the wings and tail of commercial aircraft. By performing BiGlobal stability analysis on the flowfield measured with high-resolution, stereographic Particle Image Velocimetry (PIV), previous work has shown that a family of unstable disturbances exists in the direct downstream vicinity of a supercritical, forward-facing step. The goal of the present paper is to improve the previously used stability approach on two fronts: 1) the stability problem is formulated in a non-orthogonal coordinate system and 2) the contribution of the out-of-plane base-flow derivatives is accounted for locally. The resulting eigen solutions feature significantly stronger growth rates and eigenfunctions that are localized above regions of reverse flow. The changes in the problem formulation furthermore establish a significant improvement in the comparison of the stability solutions with the results from the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content. A large destabilizing effect by the out-of-plane base-flow derivatives is determined to be justified, despite the assumption that these derivatives are small, because the effect can be reconstructed upon using an eigenvalue-correction formula that assumes the responsible terms are infinitesimally small. The unstable perturbation mechanisms are demonstrated to have a convective nature by assessing the relation between the frequency and the out-of-plane wavenumber, which indicates that their group speed does not approach zero. The last 2 facts, that 1) the perturbations are convective and 2) that the out-of-plane base-flow-derivative terms are small, remove any qualitative suspicions that the perturbation problem is not conducive to parabolization. This opens the path to analyzing this problem with a plane-marching approach.

Boundary- Layer Transition↗
NASA NTRS1995

The Physics of Turbulence in the Boundary Layer

The geometry of the velocity field in a numerically simulated incompressible turbulent boundary layer over a flat plate at Re theta=670 has been studied using the invariants of the velocity gradient tensor. These invariants are computed at every grid point in the flow and used to form the discriminant. Of primary interest are those regions in the flow where the discriminant is positive; regions where, according to the characteristic equation, the eigenvalues of the velocity gradient tensor are complex. An observer moving with a frame of reference which is attached to a fluid particle lying within such a region would see a local flow pattern of the type stable-focus-stretching or unstable-focus-compressing. When the flow is visualized this way, continuous, connected, large-scale structures are revealed that extend from the point just below the buffer layer out to the beginning of the wake region. These structures are aligned with the mean shear close to the wall and arch in the cross-stream direction away from the wall. In some cases the structures observed are very similar to to the hairpin eddy vision of boundary layer structure proposed by Theodorsen. That the structure of the flow is revealed more effectively by the discriminant rather than by the vorticity is important and adds support to recent observations of the discriminant in a channel flow simulation. Of particular importance is the fact that the procedure does not require the use of an arbitrary threshold in the discriminant. Further analysis using computer flow visualization shows a high degree of spatial correlation between regions of positive discriminant, extreme negative pressure fluctuations and large instantaneous values of Reynolds shear stress.

Kline, Stephen↗
NASA NTRS1995

Temporal variation of meandering intensity and domain-wide lateral oscillations of the Gulf Stream

The path of the Gulf Stream exhibits two modes of variability: wavelike spatial meanders associated with instability processes and large-sale lateral shifts of the path presumably due to atmospheric forcing. The objectives of this study are to examine the temporal variation of the intensity of spatial meandering in the stream, to characterize large-scale lateral oscillations in the stream's path, and to study the correlation betwen these two dynamically distinct modes of variability. The data used for this analysis are path displacemets ofthe Gulf Stream between 75 deg and 60 deg W obtained from AVHRR-derived (Advanced Very High Resolution Radiometer) infrared images for the period April 1982 through December 1989. Meandering intensity, measured by the spatial root-mean-sqaure displacement of the stream path, displays a 9-month dominant periodicity which is persistent through the study period. The 9-month fluctuation in meandering intensity may be related to the interaction of Rosseby waves with the stream. Interannual variation of meandering intensity is also found to be significant, with meandering being mich more intense during 1985 than it was in 1987. Annual variation, however,is weak and not well-defined.The spatially averaged position of the stream, which reflects nonmeandering large-scale lateral oscillations of the stream path, is dominated by an annual cycle. On average, the mean position is farthest north in November and farthest south in April. The first empirical orthogonal function mode of the space-time path displacements represents lateral oscillatins that are in-phase over the space-time domain. Interannual oscillations are also observed and are found to be weaker than the annual oscillation. The eigenvalue of the first mode indicates that about 21.5% of the total space-time variability of the stream path can be attibuted to domain-wide lateral oscillation. The correlation between meandering intensity and domain-wide lateral oscillations is very weak.

Lee, Tong↗
NASA NTRS2001

Cross-Axis Proportional Gains Used to Control Gyroscopic Effects in a Magnetic- Bearing-Supported Flywheel

For magnetic-bearing-supported high-speed machines with significant gyroscopic effects, it is necessary to stabilize both forward and backward tilt whirling modes. Instability or the low damping of these modes can prevent the attainment of desired shaft speeds. Previous work elsewhere showed that cross-axis derivative gain in the magnetic bearing control law can improve the stability of the forward whirl mode, but it is commonly recognized that derivative gains amplify high-frequency noise and increase the required control effort. At the NASA Glenn Research Center, it has been shown previously that a simple cross-axis proportional gain can add stability (without adding noise) to either forward whirl or backward whirl, depending on the sign of the gain, but that such a gain destabilizes the other mode. It has been predicted by Glenn analysis that both modes can be stabilized by cross-axis proportional gains by utilizing the large-frequency separation of the two modes at speeds where the gyroscopic effects are significant. We use a modal controller that decouples the tilt and center-of-mass-translation modes. Only the tilt modes exhibit speed-dependent gyroscopic effects. The key to controlling them by the present method is to stabilize the backward whirl tilt mode with the appropriate sign of cross-axis proportional gain in the control law, but to include a low-pass filter on that gain term to restrict its effect only to the low-frequency backward-whirl mode. A second cross-axis term with the opposite sign and a high-pass filter stabilizes the forward whirl, which can have a frequency one or two orders of magnitude higher than the backward whirl, permitting very independent action of the two terms. Because the physical gyroscopic torques are proportional to the spinning speed of the shaft, it is convenient to gain-schedule the cross-axis control terms by making them proportional to shaft speed. This has the added benefit of avoiding a somewhat awkward zero-speed splitting of the tilt-mode eigenvalues.

Brown, Gerald V.↗
NASA NTRS2010

Fast Solution in Sparse LDA for Binary Classification

An algorithm that performs sparse linear discriminant analysis (Sparse-LDA) finds near-optimal solutions in far less time than the prior art when specialized to binary classification (of 2 classes). Sparse-LDA is a type of feature- or variable- selection problem with numerous applications in statistics, machine learning, computer vision, computational finance, operations research, and bio-informatics. Because of its combinatorial nature, feature- or variable-selection problems are NP-hard or computationally intractable in cases involving more than 30 variables or features. Therefore, one typically seeks approximate solutions by means of greedy search algorithms. The prior Sparse-LDA algorithm was a greedy algorithm that considered the best variable or feature to add/ delete to/ from its subsets in order to maximally discriminate between multiple classes of data. The present algorithm is designed for the special but prevalent case of 2-class or binary classification (e.g. 1 vs. 0, functioning vs. malfunctioning, or change versus no change). The present algorithm provides near-optimal solutions on large real-world datasets having hundreds or even thousands of variables or features (e.g. selecting the fewest wavelength bands in a hyperspectral sensor to do terrain classification) and does so in typical computation times of minutes as compared to days or weeks as taken by the prior art. Sparse LDA requires solving generalized eigenvalue problems for a large number of variable subsets (represented by the submatrices of the input within-class and between-class covariance matrices). In the general (fullrank) case, the amount of computation scales at least cubically with the number of variables and thus the size of the problems that can be solved is limited accordingly. However, in binary classification, the principal eigenvalues can be found using a special analytic formula, without resorting to costly iterative techniques. The present algorithm exploits this analytic form along with the inherent sequential nature of greedy search itself. Together this enables the use of highly-efficient partitioned-matrix-inverse techniques that result in large speedups of computation in both the forward-selection and backward-elimination stages of greedy algorithms in general.

Moghaddam, Baback↗
NASA NTRSDate not supplied

Dynamic Aeroelastic Flight Dynamic Modeling of Mach 0.745 Transonic Truss-Braced Wing

This paper presents a modeling approach for dynamic aeroelastic flight dynamic analysis of the Mach0.745 Transonic Truss-Braced Wing. The modeling approach is based on a transonic correction method to correct the Theodorsen’s theory for transonic flow. CFD unsteady Reynolds-averaged Navier-Stokes equations (RANS) simulations are conducted using FUN3D solver for a series of wing sections from the Mach 0.745Transonic Truss-Braced Wing in pitch and plunge oscillations. Unsteady lift and pitching moment coefficients are obtained and used to develop the correction terms in the Theodorsen’s theory to account for transonic aerodynamics. The unsteady lift and pitching moment derivatives with respect to the unsteady angle of attack are obtained as functions of the reduced frequency. These derivatives are used to compute the unsteady lift and pitching moment contributions by the angle of attack, pitch rate, and roll rate. They are then approximated using a frequency-domain regression to obtain the dynamic stability derivatives for the Mach 0.745 Transonic Truss-Braced Wing. The structural dynamic mode shapes of the Mach 0.745 TTBW are extracted from a NAS-TRAN finite-element model. These mode shapes are used to compute the generalized unsteady aerodynamic forces. The aerodynamic mass, damping, and stiffness and the aerodynamic lag states are constructed to couple the dynamic aeroelastic contribution to the flight dynamic model of the Mach 0.745 Transonic Truss-Braced Wing. The coupled dynamic aeroelastic flight dynamic equations of motion are formulated. The eigenvalues of the coupled system are computed. All the flight dynamic modes and structural dynamic modes are stable at Match 0.745. The effect of transonic aerodynamics generally causes all the dynamic modes to have lower damping values

Flight Dynamics↗