SEARCH · Search NASA

Results for “EIGENVECTOR”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Lifetimes and eigenstates in atmospheric chemistry

The time scales and mode of the atmosphere's response to chemical perturbations are defined by the eigenvalues and eigenvectors of the system. The eigenstates of a simplified one-box CH4-CO-OH system are analyzed. The longest time constant (smallest eigenvalue) always exceeds the lifetime defined by the steady-state loss frequency for CH4, the longest lived gas. Thus, the extent of a CH4 perturbation -- the methane response time -- is always longer than predicted by the steady-state lifetime and is independent of size of the perturbation in the linear limit. This lengthening of the atmospheric recovery time can be diagnosed by how close we are to a chemically unstable troposphere, i.e., how much OH production exceeds that minimum needed to oxidize just the global emissions of CH4, CO, and other hydrocarbons and species.

Prather, Michael J.↗

NASA NTRS1993

Linear static structural and vibration analysis on high-performance computers

Parallel computers offer the oppurtunity to significantly reduce the computation time necessary to analyze large-scale aerospace structures. This paper presents algorithms developed for and implemented on massively-parallel computers hereafter referred to as Scalable High-Performance Computers (SHPC), for the most computationally intensive tasks involved in structural analysis, namely, generation and assembly of system matrices, solution of systems of equations and calculation of the eigenvalues and eigenvectors. Results on SHPC are presented for large-scale structural problems (i.e. models for High-Speed Civil Transport). The goal of this research is to develop a new, efficient technique which extends structural analysis to SHPC and makes large-scale structural analyses tractable.

Baddourah, M. A.↗

NASA NTRS1995

Spectral ordering techniques for incomplete LU preconditoners for CG methods

The effectiveness of an incomplete LU (ILU) factorization as a preconditioner for the conjugate gradient method can be highly dependent on the ordering of the matrix rows during its creation. Detailed justification for two heuristics commonly used in matrix ordering for anisotropic problems is given. The bandwidth reduction and weak connection following heuristics are implemented through an ordering method based on eigenvector computations. This spectral ordering is shown to be a good representation of the heuristics. Analysis and test cases in two and three dimensional diffusion problems demonstrate when ordering is important, and when an ILU decomposition will be ordering insensitive. The applicability of the heuristics is thus evaluated and placed on a more rigorous footing.

Clift, Simon S.↗

NASA NTRS1995

Minimization of the Vibration Energy of Thin-Plate Structures and the Application to the Reduction of Gearbox Vibration

While the vibration analysis of gear systems has been developed, a systematic approach to the reduction of gearbox vibration has been lacking. The technique of reducing vibration by shifting natural frequencies is proposed here for gearboxes and other thin-plate structures using the theories of finite elements, modal analysis, and optimization. A triangular shell element with 18 degrees of freedom is developed for structural and dynamic analysis. To optimize, the overall vibration energy is adopted as the objective function to be minimized at the excitation frequency by varying the design variable (element thickness) under the constraint of overall constant weight. Modal analysis is used to determine the sensitivity of the vibration energy as a function of the eigenvalues and eigenvectors. The optimum design is found by the gradient projection method and a unidimensional search procedure. By applying the computer code to design problems for beams and plates, it was verified that the proposed method is effective in reducing vibration energy. The computer code is also applied to redesign the NASA Lewis gear noise rig test gearbox housing. As one example, only the shape of the top plate is varied, and the vibration energy levels of all the surfaces are reduced, yielding an overall reduction of 1/5 compared to the initial design. As a second example, the shapes of the top and two side plates are varied to yield an overall reduction in vibration energy of 1/30.

Inoue, Katsumi↗

NASA NTRS1994

Structural Health Monitoring of Large Structures

This paper describes a damage detection and health monitoring method that was developed for large space structures using on-orbit modal identification. After evaluating several existing model refinement and model reduction/expansion techniques, a new approach was developed to identify the location and extent of structural damage with a limited number of measurements. A general area of structural damage is first identified and, subsequently, a specific damaged structural component is located. This approach takes advantage of two different model refinement methods (optimal-update and design sensitivity) and two different model size matching methods (model reduction and eigenvector expansion). Performance of the proposed damage detection approach was demonstrated with test data from two different laboratory truss structures. This space technology can also be applied to structural inspection of aircraft, offshore platforms, oil tankers, ridges, and buildings. In addition, its applications to model refinement will improve the design of structural systems such as automobiles and electronic packaging.

Kim, Hyoung M.↗

NASA NTRS1996

Coherent States on Riemann Surfaces as m-Photon States

Coherent states on the m-sheeted sphere (for the SU(2) group) are used to define analytic representations. The corresponding generators create and annihilate clusters of m-photons. Non-linear Hamiltonians that contain these generators are considered and their eigenvectors and eigenvalues are explicitly calculated. The Holstein-Primakoff and Schwinger formalisms in this context are also discussed.

Vourdas, A.↗

NASA NTRS1996

Predictions of Control Inputs, Periodic Responses and Damping Levels of an Isolated Experimental Rotor in Trimmed Flight

Since the early 1990s the Aeroflightdynamics Directorate at the Ames Research Center has been conducting tests on isolated hingeless rotors in hover and forward flight. The primary objective is to generate a database on aeroelastic stability in trimmed flight for torsionally soft rotors at realistic tip speeds. The rotor test model has four soft inplane blades of NACA 0012 airfoil section with low torsional stiffness. The collective pitch and shaft tilt are set prior to each test run, and then the rotor is trimmed in the following sense: the longitudinal and lateral cyclic pitch controls are adjusted through a swashplate to minimize the 1/rev flapping moment at the 12 percent radial station. In hover, the database comprises lag regressive-mode damping with pitch variations. In forward flight the database comprises cyclic pitch controls, root flap moment and lag regressive-mode damping with advance ratio, shaft angle and pitch variations. This report presents the predictions and their correlation with the database. A modal analysis is used, in which nonrotating modes in flap bending, lag bending and torsion are computed from the measured blade mass and stiffness distributions. The airfoil aerodynamics is represented by the ONERA dynamic stall models of lift, drag and pitching moment, and the wake dynamics is represented by a state-space wake model. The trim analysis of finding, the cyclic controls and the corresponding, periodic responses is based on periodic shooting with damped Newton iteration; the Floquet transition matrix (FTM) comes out as a byproduct. The stabillty analysis of finding the frequencies and damping levels is based on the eigenvalue-eigenvector analysis of the FTM. All the structural and aerodynamic states are included from modeling to trim analysis. A major finding is that dynamic wake dramatically improves the correlation for the lateral cyclic pitch control. Overall, the correlation is fairly good.

Gaonkar, G. H.↗

NASA NTRS1990

A Shifted Block Lanczos Algorithm 1: The Block Recurrence

In this paper we describe a block Lanczos algorithm that is used as the key building block of a software package for the extraction of eigenvalues and eigenvectors of large sparse symmetric generalized eigenproblems. The software package comprises: a version of the block Lanczos algorithm specialized for spectrally transformed eigenproblems; an adaptive strategy for choosing shifts, and efficient codes for factoring large sparse symmetric indefinite matrices. This paper describes the algorithmic details of our block Lanczos recurrence. This uses a novel combination of block generalizations of several features that have only been investigated independently in the past. In particular new forms of partial reorthogonalization, selective reorthogonalization and local reorthogonalization are used, as is a new algorithm for obtaining the M-orthogonal factorization of a matrix. The heuristic shifting strategy, the integration with sparse linear equation solvers and numerical experience with the code are described in a companion paper.

Grimes, Roger G.↗

NASA NTRS1997

HARP: A Dynamic Inertial Spectral Partitioner

Partitioning unstructured graphs is central to the parallel solution of computational science and engineering problems. Spectral partitioners, such recursive spectral bisection (RSB), have proven effecfive in generating high-quality partitions of realistically-sized meshes. The major problem which hindered their wide-spread use was their long execution times. This paper presents a new inertial spectral partitioner, called HARP. The main objective of the proposed approach is to quickly partition the meshes at runtime in a manner that works efficiently for real applications in the context of distributed-memory machines. The underlying principle of HARP is to find the eigenvectors of the unpartitioned vertices and then project them onto the eigerivectors of the original mesh. Results for various meshes ranging in size from 1000 to 100,000 vertices indicate that HARP can indeed partition meshes rapidly at runtime. Experimental results show that our largest mesh can be partitioned sequentially in only a few seconds on an SP2 which is several times faster than other spectral partitioners while maintaining the solution quality of the proven RSB method. A parallel WI version of HARP has also been implemented on IBM SP2 and Cray T3E. Parallel HARP, running on 64 processors SP2 and T3E, can partition a mesh containing more than 100,000 vertices into 64 subgrids in about half a second. These results indicate that graph partitioning can now be truly embedded in dynamically-changing real-world applications.

Simon, Horst D.↗

NASA NTRS1997

Improvements in Block-Krylov Ritz Vectors and the Boundary Flexibility Method of Component Synthesis

A method of dynamic substructuring is presented which utilizes a set of static Ritz vectors as a replacement for normal eigenvectors in component mode synthesis. This set of Ritz vectors is generated in a recurrence relationship, proposed by Wilson, which has the form of a block-Krylov subspace. The initial seed to the recurrence algorithm is based upon the boundary flexibility vectors of the component. Improvements have been made in the formulation of the initial seed to the Krylov sequence, through the use of block-filtering. A method to shift the Krylov sequence to create Ritz vectors that will represent the dynamic behavior of the component at target frequencies, the target frequency being determined by the applied forcing functions, has been developed. A method to terminate the Krylov sequence has also been developed. Various orthonormalization schemes have been developed and evaluated, including the Cholesky/QR method. Several auxiliary theorems and proofs which illustrate issues in component mode synthesis and loss of orthogonality in the Krylov sequence have also been presented. The resulting methodology is applicable to both fixed and free- interface boundary components, and results in a general component model appropriate for any type of dynamic analysis. The accuracy is found to be comparable to that of component synthesis based upon normal modes, using fewer generalized coordinates. In addition, the block-Krylov recurrence algorithm is a series of static solutions and so requires significantly less computation than solving the normal eigenspace problem. The requirement for less vectors to form the component, coupled with the lower computational expense of calculating these Ritz vectors, combine to create a method more efficient than traditional component mode synthesis.

Carney, Kelly Scott↗

NASA NTRS1981

Free-Vibration Analysis of Rotating Beams by a Variable-Order Finite-Element Method

The free vibration of rotating beams is analyzed by means of a finite-element method of variable order. This method entails displacement functions that are a complete power series of a variable number of terms. The terms are arranged so that the generalized coordinates are composed of displacements and slopes at the element extremities and, additionally, displacements at certain points within the element. The displacement is assumed to be analytic within an element and thus can be approximated to any degree of accuracy desired by a complete power series. Numerical results are presented for uniform beams with zero and nonzero hub radii, tapered beams, and a nonuniform beam with discontinuities. Since the present method reduces to a conventional beam finite-element method for a cubic displacement function, the results are compared and found to be superior to the conventional results in terms of accuracy for a given number of degrees of freedom. Indeed, essentially exact eigenvalues and eigenvectors are obtained with this technique, which is far more rapidly convergent than other approaches in the literature.

Hodges, Dewey H.↗

NASA NTRS1998

Development of a Probabilistic Dynamic Synthesis Method for the Analysis of Nondeterministic Structures

Accounting for the statistical geometric and material variability of structures in analysis has been a topic of considerable research for the last 30 years. The determination of quantifiable measures of statistical probability of a desired response variable, such as natural frequency, maximum displacement, or stress, to replace experience-based "safety factors" has been a primary goal of these studies. There are, however, several problems associated with their satisfactory application to realistic structures, such as bladed disks in turbomachinery. These include the accurate definition of the input random variables (rv's), the large size of the finite element models frequently used to simulate these structures, which makes even a single deterministic analysis expensive, and accurate generation of the cumulative distribution function (CDF) necessary to obtain the probability of the desired response variables. The research presented here applies a methodology called probabilistic dynamic synthesis (PDS) to solve these problems. The PDS method uses dynamic characteristics of substructures measured from modal test as the input rv's, rather than "primitive" rv's such as material or geometric uncertainties. These dynamic characteristics, which are the free-free eigenvalues, eigenvectors, and residual flexibility (RF), are readily measured and for many substructures, a reasonable sample set of these measurements can be obtained. The statistics for these rv's accurately account for the entire random character of the substructure. Using the RF method of component mode synthesis, these dynamic characteristics are used to generate reduced-size sample models of the substructures, which are then coupled to form system models. These sample models are used to obtain the CDF of the response variable by either applying Monte Carlo simulation or by generating data points for use in the response surface reliability method, which can perform the probabilistic analysis with an order of magnitude less computational effort. Both free- and forced-response analyses have been performed, and the results indicate that, while there is considerable room for improvement, the method produces usable and more representative solutions for the design of realistic structures with a substantial savings in computer time.

Brown, A. M.↗

NASA NTRS1998

A Control Law Design Method Facilitating Control Power, Robustness, Agility, and Flying Qualities Tradeoffs: CRAFT

A multi-input, multi-output control law design methodology, named "CRAFT", is presented. CRAFT stands for the design objectives addressed, namely, Control power, Robustness, Agility, and Flying Qualities Tradeoffs. The methodology makes use of control law design metrics from each of the four design objective areas. It combines eigenspace assignment, which allows for direct specification of eigenvalues and eigenvectors, with a graphical approach for representing the metrics that captures numerous design goals in one composite illustration. Sensitivity of the metrics to eigenspace choice is clearly displayed, enabling the designer to assess the cost of design tradeoffs. This approach enhances the designer's ability to make informed design tradeoffs and to reach effective final designs. An example of the CRAFT methodology applied to an advanced experimental fighter and discussion of associated design issues are provided.

Murphy, Patrick C.↗

NASA NTRS1999

Technical Report Series on Global Modeling and Data Assimilation: Filtering Techniques on a Stretched Grid General Circulation Model - Volume 16

This report documents the techniques used to filter quantities on a stretched grid general circulation model. Standard high-latitude filtering techniques (e.g., using an FFT (Fast Fourier Transformations) to decompose and filter unstable harmonics at selected latitudes) applied on a stretched grid are shown to produce significant distortions of the prognostic state when used to control instabilities near the pole. A new filtering technique is developed which accurately accounts for the non-uniform grid by computing the eigenvectors and eigenfrequencies associated with the stretching. A filter function, constructed to selectively damp those modes whose associated eigenfrequencies exceed some critical value, is used to construct a set of grid-spaced weights which are shown to effectively filter without distortion. Both offline and GCM (General Circulation Model) experiments are shown using the new filtering technique. Finally, a brief examination is also made on the impact of applying the Shapiro filter on the stretched grid.

Takacs, Lawrence L.↗

NASA NTRS1994

The Topology of Three-Dimensional Symmetric Tensor Fields

We study the topology of 3-D symmetric tensor fields. The goal is to represent their complex structure by a simple set of carefully chosen points and lines analogous to vector field topology. The basic constituents of tensor topology are the degenerate points, or points where eigenvalues are equal to each other. First, we introduce a new method for locating 3-D degenerate points. We then extract the topological skeletons of the eigenvector fields and use them for a compact, comprehensive description of the tensor field. Finally, we demonstrate the use of tensor field topology for the interpretation of the two-force Boussinesq problem.

Lavin, Yingmei↗

NASA NTRS2000

Intrinsic Coupled Ocean-Atmosphere Modes of the Asian Summer Monsoon: A Re-assessment of Monsoon-ENSO Relationships

Using global rainfall and sea surface temperature (SST) data for the past two decades (1979-1998), we have investigated the intrinsic modes of Asian summer monsoon (ASM) and ENSO co-variability. Three recurring ASM rainfall-SST coupled modes were identified. The first is a basin scale mode that features SST and rainfall variability over the entire tropics (including the ASM region), identifiable with those occurring during El Nino or La Nina. This mode is further characterized by a pronounced biennial variation in ASM rainfall and SST associated with fluctuations of the anomalous Walker circulation that occur during El Nino/La Nina transitions. The second mode comprises mixed regional and basin-scale rainfall and SST signals, with pronounced intraseasonal and interannual variabilities. This mode features a SST pattern associated with a developing La Nina, with a pronounced low level anticyclone in the subtropics of the western Pacific off the coast of East Asia. The third mode depicts an east-west rainfall and SST dipole across the southern equatorial Indian Ocean, most likely stemming from coupled ocean-atmosphere processes within the ASM region. This mode also possesses a decadal time scale and a linear trend, which are not associated with El Nino/La Nina variability. Possible causes of year-to-year rainfall variability over the ASM and sub-regions have been evaluated from a reconstruction of the observed rainfall from singular eigenvectors of the coupled modes. It is found that while basin-scale SST can account for portions of ASM rainfall variability during ENSO events (up to 60% in 1998), regional processes can accounts up to 20-25% of the rainfall variability in typical non-ENSO years. Stronger monsoon-ENSO relationship tends to occur in the boreal summer immediately preceding a pronounced La Nina, i.e., 1998, 1988 and 1983. Based on these results, we discuss the possible impacts of the ASM on ENSO variability via the west Pacific anticyclone and articulate a hypothesis that anomalous wind forcings derived from the anticyclone may be instrumental in inducing a strong biennial modulation to natural ENSO cycles.

Lau, K.-M.↗

NASA NTRS2000

Model Refinement Using Eigensystem Assignment

IA novel approach for the refinement of finite-element-based analytical models of flexible structures is presented. The proposed approach models the possible refinements in the mass, damping, and stiffness matrices of the finite element model in the form of a constant gain feedback with acceleration, velocity, and displacement measurements, respectively. Once the free elements of the structural matrices have been defined, the problem of model refinement reduces to obtaining position, velocity, and acceleration gain matrices with appropriate sparsity that reassign a desired subset of the eigenvalues of the model, along with partial mode shapes, from their baseline values to those obtained from system identification test data. A sequential procedure is used to assign one conjugate pair of eigenvalues at each step using symmetric output feedback gain matrices, and the eigenvectors are partially assigned, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The procedure can also impose that gain matrices be dissipative to guarantee the stability of the refined model. A numerical example, involving finite element model refinement for a structural testbed at NASA Langley Research Center (Controls-Structures-Interaction Evolutionary model) is presented to demonstrate the feasibility of the proposed approach.

Maghami, Peiman G.↗

NASA NTRS2001

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗