Search NASA⌕ Search

SEARCH · Search NASA

Results for “eigenvalue method”

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

Rotordynamic Analysis of the SSME Turbopumps Using Reduced Models

Alternative methods for the rotor-dynamic and sensitivity analysis of large rotor systems are examined. The methods are assessed for their ability to utilize accurate models of reduced size along with effective procedures for describing the dynamic behavior of the systems. Frequency response-based techniques are developed for determining the steady state response to imbalance of the space shuttle main engine (SSME) turbopumps and the related eigenvalue problem. The rotor and housing are represented by reduced receptances associated with their coupling points. The housing may be described by all of its normal modes within a frequency range of interest. The effects of truncated higher and lower modes are accounted for in an approximate manner. A procedure is described for determining the sensitivity of the coupling forces to changes in the coupling elements and rotor speed of the turbopump systems. In addition, an eigenvalue sensitivity analysis technique is adopted for application to the systems. Computer programs were developed for the numerical implementation of the impedance and eigenvalue sensitivity formulated in this study.

Noah, S. T.↗
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 NTRS1988

Multiple boundary condition test (MBCT) - Identification with mode shapes

The multiple boundary condition test (MBCT) approach is a ground test method to test a class of large flexible structures which cannot be ground tested by state-of-the-art test methods due to the adverse terrestrial environment. The ultimate objective of a ground test is considered to be the validation and update of a mathematical model of the structure. The research to date has indicated the MBCT does work on numerical simulations and on experimental laboratory hardware. To date only the eigenvalue has been used in the model correlation/update by inclusion of the information in the nonlinear terms resulting from the difference between the analytical and measured eigenvectors. This paper presents the results of utilizing additional information, namely the difference in the analytical and the test eigenvectors, in the validation and update of the mathematical model.

Kuo, C. P.↗
NASA NTRS1990

Efficient numerical simulation of electron states in quantum wires

A new algorithm is presented for the numerical simulation of electrons in a quantum wire as described by a two-dimensional eigenvalue problem for Schroedinger's equation coupled with Poisson's equation. Initially, the algorithm employs an underrelaxed fixed point iteration to generate an approximation which is reasonably close to the solution. Subsequently, this approximate solution is employed as an initial guess for a Jacobian-free implementation of an approximate Newton method. In this manner the nonlinearity in the model is dealt with effectively. The effectiveness of this approach is demonstrated in a set of numerical experiments which study the electron states on the cross section of a quantum wire structure based on III-V semiconductors at 4.2 and 77 K.

Kerkhoven, Thomas↗
NASA NTRS1991

A time-accurate implicit method for chemically reacting flows at all Mach numbers

The objective of this work is to develop a unified solution algorithm capable of treating time-accurate chemically reacting flows at all Mach numbers, ranging from molecular diffusion velocities to supersonic speeds. A rescaled pressure term is used in the momentum equation to circumvent the singular behavior of pressure at low Mach numbers. A dual time-stepping integration procedure is established. The system eigenvalues become well behaved and have the same order of magnitude, even in the very low Mach number regime. The computational efficiency for moderate and high speed flow is competitive with the conventional density-based scheme. The capabilities of the algorithm are demonstrated by applying it to selected model problems including nozzle flows and flame dynamics.

Withington, J. P.↗
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 NTRS1971

Numerical solution of stiff systems of ordinary differential equations with applications to electronic circuits

Systems of ordinary differential equations in which the magnitudes of the eigenvalues (or time constants) vary greatly are commonly called stiff. Such systems of equations arise in nuclear reactor kinetics, the flow of chemically reacting gas, dynamics, control theory, circuit analysis and other fields. The research reported develops an A-stable numerical integration technique for solving stiff systems of ordinary differential equations. The method, which is called the generalized trapezoidal rule, is a modification of the trapezoidal rule. However, the method is computationally more efficient than the trapezoidal rule when the solution of the almost-discontinuous segments is being calculated.

Rosenbaum, J. S.↗
NASA NTRS1987

A new method to real-normalize measured complex modes

A time domain subspace iteration technique is presented to compute a set of normal modes from the measured complex modes. By using the proposed method, a large number of physical coordinates are reduced to a smaller number of model or principal coordinates. Subspace free decay time responses are computed using properly scaled complex modal vectors. Companion matrix for the general case of nonproportional damping is then derived in the selected vector subspace. Subspace normal modes are obtained through eigenvalue solution of the (M sub N) sup -1 (K sub N) matrix and transformed back to the physical coordinates to get a set of normal modes. A numerical example is presented to demonstrate the outlined theory.

Wei, Max L.↗
NASA NTRS1988

Dynamic analysis of the joint dominated beam

A method for determining the vibrational modes of the joint dominated beam is presented with reference to space erectable structures. Mathematical models with and without joint imperfection are developed. Joint imperfection is represented by a simple flexibility, k, in the axial direction of the chord members only. An interactive computer program is used to obtain eigenvalue characteristics and mode shapes of the system. Results are presented on frequencies of the given truss beam with and without joint imperfections for three-panel and 10-panel models. Also presented are results of parametric studies of the beam truss models of three and 10 panels with k = 0.05, 0.10, 0.50, and 1.0 x 10 to the -7th.

Abu-Saba, Elias G.↗
NASA NTRS2007

An Investigation of Synchrony in Transport Networks

The cumulative degree distributions of transport networks, such as air transportation networks and respiratory neuronal networks, follow power laws. The significance of power laws with respect to other network performance measures, such as throughput and synchronization, remains an open question. Evolving methods for the analysis and design of air transportation networks must address network performance in the face of increasing demands and the need to contain and control local network disturbances, such as congestion. Toward this end, we investigate functional relationships that govern the performance of transport networks; for example, the links between the first nontrivial eigenvalue of a network's Laplacian matrix - a quantitative measure of network synchronizability - and other global network parameters. In particular, among networks with a fixed degree distribution and fixed network assortativity (a measure of a network's preference to attach nodes based on a similarity or difference), those with the small eigenvalue are shown to be poor synchronizers, to have much longer shortest paths and to have greater clustering in comparison to those with large. A simulation of a respiratory network adds data to our investigation. This study is a beginning step in developing metrics and design variables for the analysis and active design of air transport networks.

Kincaid, Rex K.↗
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 NTRS1967

Matrix Analysis of Longitudinal and Torsional Vibrations in Nonuniform Multibranch Beams

Since longitudinal modes and frequencies provide basic data for dynamic analyses of arbitrary beam-like structures and since closed-form solutions for the modes are generally not feasible to obtain, an approximate method is developed for computing the natural frequencies and the corresponding mode shapes for a variable-section, unconstrained multibranch beam. A lumped mass analogy employing influence coefficients is used to represent the beam. The simultaneous equations of motion for the lumped mass system are derived in matrix form and algebraically manipulated to yield a classical eigenvalue equation solvable by standard procedures. The orthogonality relationship of the natural modes is derived and used to form the basis of an orthogonal sweeping process for determination of modes above the fundamental. Numerical examples including an application to a solid-fuel launch system are presented. Also, a detailed discussion is devoted to the theoretical verifications of the approximate modes and frequencies.

Wingate, Robert T.↗
NASA NTRS1998

Design of Constant Gain Dissipative Controllers for Eigensystem Assignment in Passive Systems

Partial eigensystem assignment with output feedback can lead to an unstable closed-loop system. However, output feedback with passive linear time-invariant systems, such as flexible space structures, is guaranteed to be stable if the controller is dissipative. This paper presents a novel approach for synthesis of dissipative output feedback gain matrices for assigning a selected number of closed-loop poles. Dissipativity of a gain matrix is known to be equivalent to positive semidefiniteness of the symmetric part of the matrix. A sequential procedure is presented to assign one self-conjugate pair of closed-loop eigenvalues at each step using dissipative output feedback gain matrices, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The problem of assigning one closed-loop pair is reduced to a constrained solution of a system of quadratic equations, and necessary and sufficient conditions for the existence of a solution are presented. A minimax approach is presented for determining parameters which satisfy these conditions. This method can assign as many closed-loop system poles as the number of control inputs. A numerical example of damping enhancement for a flexible structure is presented to demonstrate the approach.

Maghami, Peiman G.↗
NASA NTRS1980

Numerical methods for studying anharmonic oscillator approximations to the phi super 4 sub 2 quantum field theory

This paper is an expanded version of a talk given at the 1979 T.I.C.O.M. conference. It is a self-contained introduction, for applied mathematicians and numerical analysts, to quantum mechanics and quantum field theory. It also contains a brief description of the authors' numerical approach to the problems of quantum field theory, which may best be summarized by the question; Can we compute the eigenvalues and eigenfunctions of Schrodinger operators in infinitely many variables.

Isaacson, D.↗
NASA NTRS1974

On the flame spreading over a polymer surface

The characteristics of flame propagation along the interface between a solid fuel and a gaseous oxidizer are considered. It is found that difficulties in a study of this problem are related to the eigenvalue character of the relative velocity of the flame with respect to the oxidizer flow. In the proposed model of flame spreading the difficulties are overcome by retaining the ellipticity of the governing equations and by taking the surface pyrolysis into account. The solution procedure makes use of a sink-source approach and the Fourier-transform method. Flame propagation velocity data computed on the basis of the described approach are compared with experimental results.

Ohki, Y.↗
NASA NTRS1991

Aeroelastic modal characteristics of mistuned blade assemblies: Mode localization and loss of eigenstructure

An investigation of the effects of small mistuning on the aeroelastic modes of bladed disk assemblies with aerodynamic coupling between blades is presented. The cornerstone of the approach is the use and development of perturbation methods that exhibit the crucial role of the interblade coupling and yield general findings regarding mistuning effects. It is shown that blade assemblies with weak aerodynamic interblade coupling are highly sensitive to small blade mistuning, and that their dynamics is quantitatively altered in the following ways: the regular pattern that characterizes the root locus of the tuned aeroelastic eigenvalues in the complex plane is totally lost; the aeroelastic mode shapes becomes severely localized to only a few blades of the assembly and lose their constant interblade phase angle feature; and curve veering phenomena take place when the eigenvalues are plotted versus a mistuning parameter.

Pierre, Christophe↗
NASA NTRS1991

Aeroelastic modal characteristics of mistuned blade assemblies - Mode localization and loss of eigenstructure

An investigation of the effects of small mistuning on the aeroelastic modes of bladed-disk assemblies with aerodynamic coupling between blades is presented. The cornerstone of the approach is the use and development of perturbation methods that exhibit the crucial role of the interblade coupling and yield general findings regarding mistuning effects. It is shown that blade assemblies with weak aerodynamic interblade coupling are highly sensitive to small blade mistuning, and that their dynamics is qualitatively altered in the following ways: the regular pattern that characterizes the root locus of the tuned aeroelastic eigenvalues in the complex plane is totally lost; the aeroelastic mode shapes become severely localized to only a few blades of the assembly and lose their constant interblade phase angle feature; curve veering phenomena take place when the eigenvalues are plotted versus a mistuning parameter.

Pierre, Christophe↗
NASA NTRS1983

Flux vector splitting and approximate Newton methods

In the present investigation, the basic approach is employed to view an iterative scheme as Newton's method or as a modified Newton's method. Attention is given to various modified Newton methods which can arise from differencing schemes for the Euler equations. Flux vector splitting is considered as the basic spatial differencing technique. This technique is based on the partition of a flux vector into groups which have certain properties. The Euler equations fluxes can be split into two groups, the first group having a flux Jacobian with all positive eigenvalues, and the second group having a flux Jacobian with all negative eigenvalues. Flux vector splitting based on a velocity-sound speed split is considered along with the use of numerical techniques to analyze nonlinear systems, and the steady Euler equations for quasi-one-dimensional flow in a nozzle. Results are given for steady flows with shocks.

Jespersen, D. C.↗