Search NASA⌕ Search

SEARCH · Search NASA

Results for “Eigenvector Continuation”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

On the eigenvalue and eigenvector derivatives of a non-defective matrix

A novel approach is introduced to address the problem of existence of differentiable eigenvectors for a nondefective matrix which may have repeated eigenvalues. The existence of eigenvector derivatives for a unique set of continuous eigenvectors corresponding to a repeated eigenvalue is rigorously established for nondefective and analytic matrices. A numerically implementable method is then developed to compute the differentiable eigenvectors associated with repeated eigenvalues. The solutions of eigenvalue and eigenvector derivatives for repeated eigenvalues are then derived. An example is given to illustrate the validity of formulations developed in this paper.

Juang, Jer-Nan↗
NASA NTRS1993

A Spectral Algorithm for Envelope Reduction of Sparse Matrices

The problem of reordering a sparse symmetric matrix to reduce its envelope size is considered. A new spectral algorithm for computing an envelope-reducing reordering is obtained by associating a Laplacian matrix with the given matrix and then sorting the components of a specified eigenvector of the Laplacian. This Laplacian eigenvector solves a continuous relaxation of a discrete problem related to envelope minimization called the minimum 2-sum problem. The permutation vector computed by the spectral algorithm is a closest permutation vector to the specified Laplacian eigenvector. Numerical results show that the new reordering algorithm usually computes smaller envelope sizes than those obtained from the current standard algorithms such as Gibbs-Poole-Stockmeyer (GPS) or SPARSPAK reverse Cuthill-McKee (RCM), in some cases reducing the envelope by more than a factor of two.

Barnard, Stephen T.↗
NASA NTRS1987

Investigation, development and application of optimal output feedback theory. Vol. 4: Measures of eigenvalue/eigenvector sensitivity to system parameters and unmodeled dynamics

Some measures of eigenvalue and eigenvector sensitivity applicable to both continuous and discrete linear systems are developed and investigated. An infinite series representation is developed for the eigenvalues and eigenvectors of a system. The coefficients of the series are coupled, but can be obtained recursively using a nonlinear coupled vector difference equation. A new sensitivity measure is developed by considering the effects of unmodeled dynamics. It is shown that the sensitivity is high when any unmodeled eigenvalue is near a modeled eigenvalue. Using a simple example where the sensor dynamics have been neglected, it is shown that high feedback gains produce high eigenvalue/eigenvector sensitivity. The smallest singular value of the return difference is shown not to reflect eigenvalue sensitivity since it increases with the feedback gains. Using an upper bound obtained from the infinite series, a procedure to evaluate whether the sensitivity to parameter variations is within given acceptable bounds is developed and demonstrated by an example.

Halyo, Nesim↗
NASA NTRS1992

Evolution of the shock front and turbulence structures in the shock/turbulence interaction

The interaction of a weak shock front with isotropic turbulence has been investigated using Direct Numerical Simulation (DNS). Two problems were considered: the ability of the field equation (the equation for a propagating surface) to model the shock; and a quantitative study of the evolution of turbulence structure using the database generated by Lee et al. Field equation model predictions for front shape have been compared with DNS results; good agreement is found for shock wave interaction with 2D turbulence and for a single steady vorticity wave. In the interaction of 3D isotropic turbulence with a normal shock, strong alignment of vorticity with the intermediate eigenvector of the rate of strain tensor (S(sup *)(sub ij) = S(sub ij) - (1/3)(delta(sub ij))(S(sub kk))) is seen to develop upstream of the shock and to be further amplified on passage through the shock. Vorticity tends to align at 90 deg to the largest eigenvector, but there is no preferred alignment with the smallest eigenvector. Upstream of the shock, the alignments continue to develop even after the velocity derivative skewness saturates. There is a significant tendency, which increases with time throughout the computational domain, for velocity to align with vorticity. The alignment between velocity and vorticity is strongest in eddy regions and weakest in convergence regions.

Kevlahan, N.↗
NASA NTRS1996

Visualization of 3-D tensor fields

Second-order tensor fields have applications in many different areas of physics, such as general relativity and fluid mechanics. The wealth of multivariate information in tensor fields makes them more complex and abstract than scalar and vector fields. Visualization is a good technique for scientists to gain new insights from them. Visualizing a 3-D continuous tensor field is equivalent to simultaneously visualizing its three eigenvector fields. In the past, research has been conducted in the area of two-dimensional tensor fields. It was shown that degenerate points, defined as points where eigenvalues are equal to each other, are the basic singularities underlying the topology of tensor fields. Moreover, it was shown that eigenvectors never cross each other except at degenerate points. Since we live in a three-dimensional world, it is important for us to understand the underlying physics of this world. In this report, we describe a new method for locating degenerate points along with the conditions for classifying them in three-dimensional space. Finally, we discuss some topological features of three-dimensional tensor fields, and interpret topological patterns in terms of physical properties.

Hesselink, L.↗
NASA NTRS1978

On the append and continue features in NASTRAN

The APPEND feature is described which is applicable in the case of real eigenvalue analysis. This feature permitted the addition of new eigenvalues and eigenvectors to those already computed in a previously checkpointed run without re-executing the entire problem. The next feature was the CONTINUE feature which was applicable in the case of transient analysis of coupled equations. This feature enabled the integration of coupled equations to be continued beyond the last output time for which the solution was obtained in a previously checkpointed run (without re-executing the entire problem). These two features are illustrated by suitable examples.

Pamidi, P. R.↗
NASA NTRS1998

Using Solar Radiation Pressure to Control L2 Orbits

The main perturbations at the Sun-Earth Lagrange points L1 and L2 are from solar radiation pressure (SRP), the Moon and the planets. Traditional approaches to trajectory design for Lagrange-point orbits use maneuvers every few months to correct for these perturbations. The gravitational effects of the Moon and the planets are small and periodic. However, they cannot be neglected because small perturbations in the direction of the unstable eigenvector are enough to cause exponential growth within a few months. The main effect of a constant SRP is to shift the center of the orbit by a small distance. For spacecraft with large sun-shields like the Microwave Anisotropy Probe (MAP) and the Next Generation Space Telescope (NGST), the SRP effect is larger than all other perturbations and depends mostly on spacecraft attitude. Small variations in the spacecraft attitude are large enough to excite or control the exponential eigenvector. A closed-loop linear controller based on the SRP variations would eliminate one of the largest errors to the orbit and provide a continuous acceleration for use in controlling other disturbances. It is possible to design reference trajectories that account for the periodic lunar and planetary perturbations and still satisfy mission requirements. When such trajectories are used the acceleration required to control the unstable eigenvector is well within the capabilities of a continuous linear controller. Initial estimates show that by using attitude control it should be possible to minimize and even eliminate thruster maneuvers for station keeping.

Tene, Noam↗
NASA NTRS1994

Multigrid techniques for nonlinear eigenvalue probems: Solutions of a nonlinear Schroedinger eigenvalue problem in 2D and 3D

This paper presents multigrid (MG) techniques for nonlinear eigenvalue problems (EP) and emphasizes an MG algorithm for a nonlinear Schrodinger EP. The algorithm overcomes the mentioned difficulties combining the following techniques: an MG projection coupled with backrotations for separation of solutions and treatment of difficulties related to clusters of close and equal eigenvalues; MG subspace continuation techniques for treatment of the nonlinearity; an MG simultaneous treatment of the eigenvectors at the same time with the nonlinearity and with the global constraints. The simultaneous MG techniques reduce the large number of self consistent iterations to only a few or one MG simultaneous iteration and keep the solutions in a right neighborhood where the algorithm converges fast.

Costiner, Sorin↗
NASA NTRS2012

Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

Orr, Jeb S.↗
NASA NTRS1992

Efficient use of high performance computers for integrated controls and structures design

Simply transporting design codes from sequential-scalar computers to parallel-vector computers does not fully utilize the computational benefits offered by high performance computers. By performing integrated controls and structures design on an experimental truss platform with both sequential-scalar and parallel-vector design codes, conclusive results are presented to substantiate this claim. The efficiency of a Cholesky factorization scheme in conjunction with a variable-band row data structure is presented. In addition, the Lanczos eigensolution algorithm has been incorporated in the design code for both parallel and vector computations. Comparisons of computational efficiency between the initial design code and the parallel-vector design code are presented. It is shown that the Lanczos algorithm with the Cholesky factorization scheme is far superior to the sub-space iteration method of eigensolution when substantial numbers of eigenvectors are required for control design and/or performance optimization. Integrated design results show the need for continued efficiency studies in the area of element computations and matrix assembly.

Belvin, W. K.↗
NASA NTRS1999

Aerodynamic and Aeroelastic Insights using Eigenanalysis

This paper presents novel analytical results for eigenvalues and eigenvectors produced using discrete time aerodynamic and aeroelastic models. An unsteady, incompressible vortex lattice aerodynamic model is formulated in discrete time; the importance of several modeling parameters is examined. A detailed study is made of the behavior of the aerodynamic eigenvalues both in discrete and continuous time. The aerodynamic model is then incorporated into aeroelastic equations of motion. Eigenanalyses of the coupled equations produce stability results and modal characteristics which are valid for critical and non-critical velocities. Insight into the modeling and physics associated with aeroelastic system behavior is gained by examining both the eigenvalues and the eigenvectors. Potential pitfalls in discrete time model construction and analysis are examined.

Heeg, Jennifer↗
NASA NTRS2004

Aerodynamic and Aeroelastic Insights using Eigenanalysis

This paper presents novel analytical results for eigenvalues and eigenvectors produced using discrete time aerodynamic and aeroelastic models. An unsteady, incompressible vortex lattice aerodynamic model is formulated in discrete time; the importance of several modeling parameters is examined. A detailed study is made of the behavior of the aerodynamic eigenvalues both in discrete and continuous time. The aerodynamic model is then incorporated into aeroelastic equations of motion. Eigenanalyses of the coupled equations produce stability results and modal characteristics which are valid for critical and non-critical velocities. Insight into the modeling and physics associated with aeroelastic system behavior is gained by examining both the eigenvalues and the eigenvectors. Potential pitfalls in discrete time model construction and analysis are examined.

Heeg, Jennifer↗
NASA NTRS1972

Spectral factorization in periodically time-varying systems and application to navigation problems.

Spectral factorization has been used previously to derive the steady-state solution of Kalman filtering equations without iteration for constant coefficient systems. The present work extends the spectral factorization algorithm to time-varying systems having periodic coefficient matrices for cases of both discrete and continuous systems. Time-consuming, expensive iterations of sequential covariance equations are not required to reach the final solution since this is an algebraic algorithm employing existing eigenvalue, eigenvector subroutines. The computer program incorporating the algorithm is suitable for sensitivity studies in formulating navigation and guidance strategies of low-thrust interplanetary missions. The determination of an optimum tracking pattern from an earth station is examined as an example.

Nishimura, T.↗
NASA NTRS1980

The finite element duct eigenvalue problem - An improved formulation with Hermitian elements and no-flow condensation

Hermitian elements are used in a finite element solution for the eigenvalue problem in lined ducts with flow. These elements give significantly greater accuracy for reduced dimensionality when compared with Lagrangian elements. Spurious mode generation associated with the Lagrangian formulation is eliminated. A dramatic improvement in the ratio of the number of reliable eigenvalues to the total number of computed eigenvalues is effected by the use of a condensation scheme based on the no-flow eigenvectors. Results are presented for two dimensional and axisymmetric ducts. In the axisymmetric case good resolution is obtained even for high order, high frequency modes by the use of continuously graded meshes.

Astley, R. J.↗
NASA NTRS1995

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗
NASA NTRS1979

A finite element formulation of the eigenvalue problem in lined ducts with flow

A finite element method is used to formulate the eigenvalue problem for a lined duct with flow. Either two dimensional or axially symmetric ducts with sheared flows can be studied, although the examples cited are two dimensional. The primitive variables of velocity and pressure are used with quadratic shape functions in each element. Results indicate that a useful level of accuracy can be achieved with a modest number of elements. Spurious eigenvalues, easily identified by obviously inconsistent eigenvectors, occur in certain instances. If the number of elements is not too small, these spurious modes are among the higher order eigenvalues of otherwise questionable accuracy. The possibility of using higher order elements which have slope continuity is proposed for future investigations to eliminate the spurious results.

Astley, R. J.↗