SEARCH · Engineering Papers

Results for “Eigenvalues and eigenvectors”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

NASA NTRS1998

Model Refinement Using Eigensystem Assignment

This paper presented a 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, which 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 self-conjugate pair of closed-loop eigenvalues at each step using symmetric output feedback gain matrices, and the closed-loop 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 in order to guarantee the stability of the refined model. A numerical example, involving finite element model refinement for a structural testbed at NASA Langley (CSI Evolutionary Model) is presented to demonstrate the feasibility of the proposed approach.

Maghami, Pieman G.↗

NASA NTRS1974

Analytical and experimental investigation of a 1/8-scale dynamic model of the shuttle orbiter. Volume 1: Summary report

A 1/8-scale structural dynamics model of the space shuttle orbiter was analyzed using the NASA Structural Analysis System (NASTRAN). Comparison of the calculated eigenvalues with preliminary test data for the unrestrained condition indicate that the analytical model was consistently stiffer, being about 20% higher in the first mode. The eigenvectors show reasonably good agreement with test data. A series of analytical and experimental investigations undertaken to resolve the discrepancy are described. Modifications in the NASTRAN model based upon these investigations resulted in close agreement for both eigenvalues and eigenvectors.

Mason, P. W.↗

DOE OSTI2023

High-Impedance Non-Linear Fault Detection via Eigenvalue Analysis with low PMU Sampling Rates

This work presents a hybrid data-driven and physics-based framework for high-impedance fault detection in power systems. An innovative method based on eigenvalue analysis is expanded and validated. Phasor Measurement Unit data is used to estimate eigenvalues corresponding to the powerlines being monitored. The projection and drift of these eigenvalues is then tracked and evaluated. Faults are detected as they drive eigenvalues outside of their normal zones. Eigenvectors are leveraged to support and validate the decisions made by the main algorithm. This technique holds several advantages over contemporary techniques in that it utilizes technology that is already deployed in the field, it offers a significant degree of generality, and so far it has displayed a very high-level of sensitivity without sacrificing accuracy. Validation takes place in the form of simulations in the IEEE 13 Node System considering a popular high-impedance non-linear fault model. Test results are encouraging indicating potential for real-life applications.

Paramo, Gian↗

DOE OSTI2022

Eigenmode Analysis of Pulsed Neutron Transport Simulations

We discuss the time dependent behavior of some simple pulsed neutron simulations of subcritical problems in slab geometry. Our intent is to investigate the eigenvalue structure of the discretized neutron transport equation and to show under some reasonable assumptions that a dominant time eigenvalue exists that has a nonnegative eigenvector that determines the long time dependent behavior of the solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1984

Design of an optimal output feedback control system with modal insensitivity

This paper deals with the design of an output feedback controller which results in selected modal insensitivity, and at the same time optimizes a quadratic performance index representative of desired system performance for nominal plant parameter values. The approach taken here is to characterize the class of attainable eigenvectors for a given set of eigenvalues (distinct or non-distinct) which lie in a subspace called the 'Modal Insensitivity Subspace'. A constraint is established on the feedback matrix which results in modal insensitivity. Necessary conditions for optimality subject to the constraint on the feedback matrix are given. This forms the basis for a numerical algorithm to compute the optimal feedback gain which analyzed for convergence. To illustrate the procedure, a design is carried out using the lateral dynamics of an L-1011 aircraft.

Raman, K. V.↗

NASA NTRS1986

Eigenfunction methods in magnetospheric radial-diffusion theory

Complete sets of orthonormal basis functions constructed according to a generalization of the quantum-mechanical WKB approximation can be used to generate a nearly-diagonal matrix representation of the radial-transport operator for ring-current ions in the presence of radial diffusion and charge exchange. The resulting eigenfunctions (constructed by weighting the basis functions in proportion to the respective components of the eigenvectors of the matrix representation) and eigenvalues provide a spatial and temporal description of the evolving phase-space density during and following a magnetospheric disturbance (e.g., a magnetic storm). A linear superposition of the basis functions can also be used to eliminate any discrepancy between the steady-state solution of the transport equation and the appropriate WKB approximation of this steady-state solution.

Schulz, Michael↗

NASA NTRS1987

Modal decoupling insensitivity

This note deals with the design of an output feedback control system such that selected closed-loop eigenvalues are insensitive, while the corresponding eigenvectors preserve a prescribed decoupling under parameter variation. This results in restricting specified elements of an eigenvector to be insensitive, rather than the entire eigenvector. This requirement is henceforth referred to as modal decoupling insensitivity. Sufficient conditions for modal decoupling insensitivity are derived. These conditions offer the advantage of increased flexibility in eigenvector assignment, but the constraint on the feedback matrix is more restrictive. A numerical example is presented to illustrate the design procedure.

Raman, K. V.↗

NASA NTRSDate not supplied

Dynamic Mode Decomposition of Unsteady Pressure-Sensitive Paint Measurements for the NASA Unitary Plan Wind Tunnel Tests

This paper describes the Dynamic Mode Decomposition (DMD) of the pressures on the scale model of the Space Launch System (SLS) Block 1 cargo vehicle with the Unsteady Pressure-Sensitive Paint (uPSP) measurements, which were collected in the Ascent Transient Aerodynamics Tests with the Unitary Plan Wind Tunnel 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. The conventional DMD algorithm is based on the Singular Value Decomposition (SVD) of the data matrix. For the matrix of the uPSP measurements of the SLS ATAT, the number of rows is equal to the number of nodes in the grid of the scale model, and the number of columns is equal to the number of frames in the videos taken with 4 Phantom high-speed cameras. In this paper, it is verified that, for the time series with zero mean value, the DMD is equivalent to the decomposition with the Discrete Fourier Transform (DFT). Considering the uPSP is mainly used in the assessment of the unsteady, aerodynamic phenomena, the DMD of the uPSP measurements can be implemented in two steps: (1) subtract the mean value from the uPSP measurement on each of the grid nodes; (2) apply the Fast Fourier Transform (FFT) on the resulting zero-mean time series. The DMD of the uPSP measurements with FFT has two advantages: (1) the computational complexity of FFT is O(N*logN), where N is the length of the time series; (2) compared to the SVD-based DMD algorithm, the DMD with FFT can be easily implemented in parallel processing. A sample matrix of uPSP measurements, at the size of 341 grid nodes and 128 frames, is generated. Figures 1 and 2 show the eigenvalues and the ratios of the eigenvectors, respectively, of the sample matrix, without and with the mean value removed on each of the grid nodes, computed with the SVD-based DMD and the FFT. The figures demonstrate the equivalence of the SVD-based DMD and the decomposition with DFT/FFT for the time series with zero mean value. The results of DMD of the uPSP measurements of the SLS ATAT in September 2019 are presented in the paper. The DMD modes at different frequencies are shown, the aerodynamic phenomena (e.g. shockwave and vortex shedding) are demonstrated and the correlation of the DMD modes with the test configuration parameter (e.g., the Mach Number) is discussed. Figure 3 shows a software tool to visualize the DMD modes. The code to implement the algorithm described in this paper was written in C, with libraries of FFTW for FFT and MPI/OpenMP for parallel processing, and executed on the NASA Pleiades supercomputer. Funding for this research was provided by the NASA Aerosciences Evaluation and Test Capabilities Project.

Pressure-Sensitive Paint↗

NASA NTRS1981

Eigenvalues of the Rayleigh-Benard and Marangoni problems

The eigenvalues of the linear Benard-Marangoni stability problem are discussed. Pearson and Nield boundary conditions, which correspond to a rigid, isothermal lower boundary and a stress-free conducting upper boundary are considered. It is shown that although a critical value of the Marangoni number can be determined, the number is not, strictly speaking, an eigenvalue and cannot be used as an eigenvalue parameter for the determination of an eigenvector set.

Rosenblat, S.↗

NASA NTRS1984

Eigenspace design techniques for active flutter suppression

The application of eigenspace design techniques to an active flutter suppression system for the DAST ARW-2 research drone is examined. Eigenspace design techniques allow the control system designer to determine feedback gains which place controllable eigenvalues in specified configurations and which shape eigenvectors to achieve desired dynamic response. Eigenspace techniques were applied to the control of lateral and longitudinal dynamic response of aircraft. However, little was published on the application of eigenspace techniques to aeroelastic control problems. This discussion will focus primarily on methodology for design of full-state and limited-state (output) feedback controllers. Most of the states in aeroelastic control problems are not directly measurable, and some type of dynamic compensator is necessary to convert sensor outputs to control inputs. Compensator design are accomplished by use of a Kalman filter modified if necessary by the Doyle-Stein procedure for full-state loop transfer function recovery, by some other type of observer, or by transfer function matching.

Garrard, W. L.↗

NASA NTRS1983

An eigenvalue/eigenvector assignment algorithm using output feedback

An eigenvalue/eigenvector assignment algorithm using stationary output feedback is presented. The algorithm permits assignment of min (n, m + r - 1) eigenvalues and max (m-1, r-1) eigenvectors, where n, m, r refer to the system state, input and output dimensions, respectively. An example is given to illustrate the design procedures.

Mielke, R. R.↗

NASA NTRS1986

Eigensystem synthesis for active flutter suppression on an oblique-wing aircraft

The application of the eigensystem synthesis technique to place the closed-loop eigenvalues and shape the closed-loop eigenvectors has not been practical for active flutter suppression, primarily because of the availability of only one control surface (aileron) for flutter suppression. The oblique-wing aircraft, because of its configuration, provides two independent surfaces (left and right ailerons), making the application of eigensystem synthesis practical. This paper presents the application of eigensystem synthesis using output feedback for the design of an active flutter suppression system for an oblique-wing aircraft. The results obtained are compared with those obtained by linear quadratic Gaussian techniques.

Alag, G. S.↗

NASA NTRS1986

Eigensystem synthesis for active flutter suppression on an oblique-wing aircraft

Alag, G. S.↗

NASA NTRS1989

Flutter suppression using eigenspace freedoms to meet requirements

A constrained optimization methodology has been developed which allows specific use of eigensystem freedoms to meet design requirements. A subset of the available eigenvector freedoms was employed. The eigenvector freedoms associated with a particular closed-loop eigenvalue are coefficients of basis vectors which span the subspace in which that closed-loop vector must lie. Design requirements are included as a vector of inequality constraints. The procedure was successfully applied to develop an unscheduled controller which stabilizes symmetric flutter of an aeroelastic vehicle to a dynamic pressure 44 percent above the open-loop flutter point. The design process proceeded from full-state feedback to the inclusion of a full-order observer to the selection of an eighth-order controller which preserved the full-state sensitivity characteristics. Only a subset of the design freedoms was utilized (i.e., assuming full-state feedback only four out of 26 eigenvectors were used, and no variations were made in the closed-loop eigenvalues). Utilization of additional eigensystem freedoms could further improve the controller.

Adams, William M., Jr.↗

NASA NTRS1980

Preliminary demonstration of a robust controller design method

Alternative computational procedures for obtaining a feedback control law which yields a control signal based on measurable quantitites are evaluated. The three methods evaluated are: (1) the standard linear quadratic regulator design model; (2) minimization of the norm of the feedback matrix, k via nonlinear programming subject to the constraint that the closed loop eigenvalues be in a specified domain in the complex plane; and (3) maximize the angles between the closed loop eigenvectors in combination with minimizing the norm of K also via the constrained nonlinear programming. The third or robust design method was chosen to yield a closed loop system whose eigenvalues are insensitive to small changes in the A and B matrices. The relationship between orthogonality of closed loop eigenvectors and the sensitivity of closed loop eigenvalues is described. Computer programs are described.

Anderson, L. R.↗

NASA NTRS1991

Monte Carlo simulation of errors in the anisotropy of magnetic susceptibility - A second-rank symmetric tensor

Monte Carlo perturbations of synthetic tensors to evaluate the Hext/Jelinek elliptical confidence regions for anisotropy of magnetic susceptibility (AMS) eigenvectors are used. When the perturbations are 33 percent of the minimum anisotropy, both the shapes and probability densities of the resulting eigenvector distributions agree with the elliptical distributions predicted by the Hext/Jelinek equations. When the perturbation size is increased to 100 percent of the minimum eigenvalue difference, the major axis of the 95 percent confidence ellipse underestimates the observed eigenvector dispersion by about 10 deg. The observed distributions of the principal susceptibilities (eigenvalues) are close to being normal, with standard errors that agree well with the calculated Hext/Jelinek errors. The Hext/Jelinek ellipses are also able to describe the AMS dispersions due to instrumental noise and provide reasonable limits for the AMS dispersions observed in two Hawaiian basaltic dikes. It is concluded that the Hext/Jelinek method provides a satisfactory description of the errors in AMS data and should be a standard part of any AMS data analysis.

Lienert, Barry R.↗

NASA NTRS1994

Gain weighted eigenspace assignment

This report presents the development of the gain weighted eigenspace assignment methodology. This provides a designer with a systematic methodology for trading off eigenvector placement versus gain magnitudes, while still maintaining desired closed-loop eigenvalue locations. This is accomplished by forming a cost function composed of a scalar measure of error between desired and achievable eigenvectors and a scalar measure of gain magnitude, determining analytical expressions for the gradients, and solving for the optimal solution by numerical iteration. For this development the scalar measure of gain magnitude is chosen to be a weighted sum of the squares of all the individual elements of the feedback gain matrix. An example is presented to demonstrate the method. In this example, solutions yielding achievable eigenvectors close to the desired eigenvectors are obtained with significant reductions in gain magnitude compared to a solution obtained using a previously developed eigenspace (eigenstructure) assignment method.

Davidson, John B.↗

NASA NTRS1974

Results on the two population feature selection problem using probability of correct classification as a criterion

Variational equations are presented for maximizing the probability of correct classification as a function of a 1xn feature selection matrix B for the two-population problem. For the special case of equal covariance matrices the optimal B is unique up to scalar multiples and rank one sufficient. For equal population means, the best 1xn B is an eigenvector corresponding either to the largest or smallest eigenvalue of sigma sub 2 to the minus 1 power sigma sub 1 where sigma sub 1 and sigma sub 2 are the nxn covariance matrices of the two populations. The transformed probability of correct classification depends only on the eigenvalue. Finally, a procedure is proposed for constructing an optimal or nearly optimal kxn matrix of rank k without solving the k-dimensional variational equation.

Peters, B. C.↗