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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Buckling of rotating beams

The stability of a beam subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. The differential equations of motion are solved numerically using an integrating matrix method in combination with an eigenanalysis to determine the eigenvalues from which stability is assessed. The results clarify several differences which have been identified in the literature relating to the proper behavior of the critical rotational speed for buckling as the radius of rotation of the clamped end of the beam is reduced.

White, W. F., Jr.↗

NASA NTRS1980

Improvement of analytical dynamic models using modal test data

A method developed to determine maximum changes in analytical mass and stiffness matrices to make them consistent with a set of measured normal modes and natural frequencies is presented. The corrected model will be an improved base for studies of physical changes, boundary condition changes, and for prediction of forced responses. The method features efficient procedures not requiring solutions of the eigenvalue problem, and the ability to have more degrees of freedom than the test data. In addition, modal displacements are obtained for all analytical degrees of freedom, and the frequency dependence of the coordinate transformations is properly treated.

Berman, A.↗

NASA NTRS1986

Numerical simulation of a controlled boundary layer

The problem of interest is the boundary layer over a flat plate. The three standard laminar flow control (LFC) techniques are pressure gradient, suction, and heating. The parameters used to describe the amount of control in the context of the boundary layer equations are introduced. The numerical method required to find the mean flow, the linear eigenvalues of the Orr-Sommerfeld equation, and the full, nonlinear, 3-D solution of the Navier-Stokes equations are outlined. A secondary instability exists for the parallel boundary subject to uniform pressure gradient, suction, or heating. Selective control of the spanwise mode reduces the secondary instability in the parallel boundary layer at low Reynolds number.

Zang, Thomas A.↗

NASA NTRS1992

Calculation of two-dimensional shear-driven cavity flows at high Reynolds numbers

The time-dependent Navier-Stokes equations are numerically integrated for two-dimensional incompressible viscous flow in a shear-driven square cavity. Using a time-splitting method and finite differences on a staggered mesh, the momentum and pressure equations are directly solved by a tensor product method where one finite difference direction is diagonalized by eigenvalue decomposition. The effects of increasing Reynolds number are studied and the developing boundary layer is captured by using a finely clustered mesh. At Re = 30,000 the flow is in a continuously developing unsteady regime. Power spectrum plots indicate that the unsteady flow oscillates with one fundamental frequency and exhibits some characteristics of transition between laminar and turbulent states.

Huser, A.↗

NASA NTRS1990

A study of delamination buckling of laminates

The subject of this paper is the buckling of laminated plates, with a preexisting delamination, subjected to in-plane loading. Each laminate is modelled as an orthotropic Mindlin plate. The analysis is carried out by a combination of the finite element and asymptotic expansion methods. By applying the finite element method, plates with general delamination regions can be studied. The asymptotic expansion method reduces the number of unknown variables of the eigenvalue equation to that of the equation for a single Kirchhoff plate. Numerical results are presented for several examples. The effects of the shape, size, and position of the delamination on the buckling load are studied through these examples.

Mukherjee, Yu-Xie↗

NASA NTRS1991

Identification of linear multivariable systems from a single set of data by identification of observers with assigned real eigenvalues

This paper presents a formulation for identification of linear multivariable systems from a single set of input-output data. The identification method is formulated with the mathematical framework of learning identification, by extension of the repetition domain concept to include shifting time intervals. This contrasts existing learning approaches that require data from multiple experiments. In this method, the system input-output relationship is expressed in terms of an observer, which is made asymptotically stable by an embedded real eigenvalue assignment procedure. Through this relationship, the Markov parameters of the observer are identified. The Markov parameters of the actual system are recovered from those of the observer, and then used to obtain a state space model of the system by standard realization techniques. The basic mathematical formulation is derived, and numerical examples presented to illustrate the proposed method.

Phan, Minh↗

NASA NTRS2008

Generation of Boundary Manikin Anthropometry

The purpose of this study was to develop 3D digital boundary manikins that are representative of the anthropometry of a unique population. These digital manikins can be used by designers to verify and validate that the components of the spacesuit design satisfy the requirements specified in the Human Systems Integration Requirements (HSIR) document. Currently, the HSIR requires the suit to accommodate the 1st percentile American female to the 99th percentile American male. The manikin anthropometry was derived using two methods: Principal Component Analysis (PCA) and Whole Body Posture Based Analysis (WBPBA). PCA is a statistical method for reducing a multidimensional data set by using eigenvectors and eigenvalues. The goal is to create a reduced data set that encapsulates the majority of the variation in the population. WBPBA is a multivariate analytical approach that was developed by the Anthropometry and Biomechanics Facility (ABF) to identify the extremes of the population for a given body posture. WBPBA is a simulation-based method that finds extremes in a population based on anthropometry and posture whereas PCA is based solely on anthropometry. Both methods yield a list of subjects and their anthropometry from the target population; PCA resulted in 20 female and 22 male subjects anthropometry and WBPBA resulted in 7 subjects' anthropometry representing the extreme subjects in the target population. The subjects anthropometry is then used to 'morph' a baseline digital scan of a person with the same body type to create a 3D digital model that can be used as a tool for designers, the details of which will be discussed in subsequent papers.

Young, Karen S.↗

NASA NTRS1993

The multigrid preconditioned conjugate gradient method

A multigrid preconditioned conjugate gradient method (MGCG method), which uses the multigrid method as a preconditioner of the PCG method, is proposed. The multigrid method has inherent high parallelism and improves convergence of long wavelength components, which is important in iterative methods. By using this method as a preconditioner of the PCG method, an efficient method with high parallelism and fast convergence is obtained. First, it is considered a necessary condition of the multigrid preconditioner in order to satisfy requirements of a preconditioner of the PCG method. Next numerical experiments show a behavior of the MGCG method and that the MGCG method is superior to both the ICCG method and the multigrid method in point of fast convergence and high parallelism. This fast convergence is understood in terms of the eigenvalue analysis of the preconditioned matrix. From this observation of the multigrid preconditioner, it is realized that the MGCG method converges in very few iterations and the multigrid preconditioner is a desirable preconditioner of the conjugate gradient method.

Tatebe, Osamu↗

NASA NTRS2002

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one- dimensional convection-diffusion-reaction (CDR) equations. Accuracy, stability, conservation, and dense-output are first considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, (N = 2), additive Runge-Kutta (ARK(sub 2)) methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms of the partitioned method are of equal order to those of the elemental methods. Derived ARK(sub 2) methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z['] yields -infinity, and retain high stability efficiency in the absence of stiffness, z['] yield 0. Extrapolation-type stage- value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK(sub 2) error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK(sub 2) methods.

Kennedy, Christopher A.↗

NASA NTRS1988

Vibration and Control of Flexible Rotor Supported by Magnetic Bearings

Active vibration control of flexible rotors supported by magnetic bearings is discussed. Using a finite-element method for a mathematical model of the flexible rotor, the eigenvalue problem is formulated taking into account the interaction between a mechanical system of the flexible rotor and an electrical system of the magnetic bearings and the controller. However, for the sake of simplicity, gyroscopic effects are disregarded. It is possible to adapt this formulation to a general flexible rotor-magnetic bearing system. Controllability with and without collocation sensors and actuators located at the same distance along the rotor axis is discussed for the higher order flexible modes of the test rig. In conclusion, it is proposed that it is necessary to add new active control loops for the higher flexible modes even in the case of collocation. Then it is possible to stabilize for the case of uncollocation by means of this method.

Nonami, Kenzou↗

NASA NTRS1981

Acoustic transmission in non-uniform ducts with mean flow. I - The method of weighted residuals. II - The finite element method

The problem of acoustic transmission through nonuniform ducts containing a high-speed subsonic flow is studied by means of the method of weighted residuals in the form of a modified Galerkin method and a Galerkin formulation of the finite element method. The method of weighted residuals is shown to employ the basis functions generated from eigenvalue calculations for the case of no flow, and is verified by comparison with exact eigenvalue calculations in the uniform duct case and numerical solutions of the one-dimensional form of the equations in the nonuniform duct case. The finite element scheme based on both the Galerkin method and the residual least squares method and employing eight-noded isoparametric elements is presented and used to investigate multimodal propagation by the coupling of the solution in the duct nonuniform section to modal expansions in uniform sections. Comparison of the results of the two methods reveals them to be in substantial agreement, and predicts the importance of multimodal interactions at high Mach numbers.

Eversman, W.↗

NASA NTRS1993

On the cross-stream spectral method for the Orr-Sommerfeld equation

Cross-stream models are defined as solutions to the Orr-Sommerfeld equation which are propagating normal to the flow direction. These models are utilized as a basis for a Hilbert space to approximate the spectrum of the Orr-Sommerfeld equation with plane Poiseuille flow. The cross-stream basis leads to a standard eigenvalue problem for the frequencies of Poiseuille flow instability waves. The coefficient matrix in the eigenvalue problem is shown to be the sum of a real matrix and a negative-imaginary diagonal matrix which represents the frequencies of the cross-stream modes. The real coefficient matrix is shown to approach a Toeplitz matrix when the row and column indices are large. The Toeplitz matrix is diagonally dominant, and the diagonal elements vary inversely in magnitude with diagonal position. The Poiseuille flow eigenvalues are shown to lie within Gersgorin disks with radii bounded by the product of the average flow speed and the axial wavenumber. It is shown that the eigenvalues approach the Gersgorin disk centers when the mode index is large, so that the method may be used to compute spectra with an essentially unlimited number of elements. When the mode index is large, the real part of the eigenvalue is the product of the axial wavenumber and the average flow speed, and the imaginary part of the eigen value is identical to the corresponding cross-stream mode frequency. The cross-stream method is numerically well-conditioned in comparison to Chebyshev based methods, providing equivalent accuracy for small mode indices and superior accuracy for large indices.

Zorumski, William E.↗

NASA NTRS1975

Relaxation methods in fluid mechanics

The present work considers the iterative solution of a coupled set of difference equations and examines methods that carry successive approximates to a state that is invariant with further iteration and independent of the initial guess. Methods are studied with regard to their efficiency and economy of computer resources. The basic principles of classical relaxation are set forth, with attention confined to linear elliptic equations. This discussion involves the evaluation of the spectral radius that is the magnitude of the eigenvalue with largest modulus. The subject of relaxation is then related to the study of ordinary differential equations and hyperbolic partial differential equations. Problems that occur when linearly dependent eigenvectors appear in the relaxation matrix are discussed, leading to multiply connected eigenvalues in the Jordan canonical form. Finally, a brief survey of relaxation methods used in aerodynamics is given.

Lomax, H.↗

NASA NTRS1988

Investigating a hybrid perturbation-Galerkin technique using computer algebra

A two-step hybrid perturbation-Galerkin method is presented for the solution of a variety of differential equations type problems which involve a scalar parameter. The resulting (approximate) solution has the form of a sum where each term consists of the product of two functions. The first function is a function of the independent field variable(s) x, and the second is a function of the parameter lambda. In step one the functions of x are determined by forming a perturbation expansion in lambda. In step two the functions of lambda are determined through the use of the classical Bubnov-Gelerkin method. The resulting hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Bubnov-Galerkin methods applied separately, while combining some of the good features of each. In particular, the results can be useful well beyond the radius of convergence associated with the perturbation expansion. The hybrid method is applied with the aid of computer algebra to a simple two-point boundary value problem where the radius of convergence is finite and to a quantum eigenvalue problem where the radius of convergence is zero. For both problems the hybrid method apparently converges for an infinite range of the parameter lambda. The results obtained from the hybrid method are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Andersen, Carl M.↗

NASA NTRS1993

On eigensystem assignment with dissipativity constraints

Stability of the closed-loop system is of vital importance when controller synthesis is being performed using a reduced order model of the open-loop plant. For passive systems such as flexible space structures, closed-loop stability is often ensured by low-authority dissipative controllers using collocated and compatible sensors and actuators. This paper considers the problem of eigensystem assignment for flexible structures using constant-gain output feedback, under the dissipativity constraints requiring that the gain matrices be nonnegative definite. A novel method is developed in which a sequential procedure is used to assign one complex-conjugate pair of closed-loop eigenvalues in each step. At each step, an iterative technique is employed for generating a nonnegative definite gain matrix that assigns a new pair while keeping the previously assigned pairs unchanged. This method can assign up to m closed-loop system eigenvalues, where m is the number of control inputs. A numerical example of damping enhancement for a large flexible structure is presented to demonstrate the approach.

Maghami, Peiman G.↗

NASA NTRS1991

Role of artificial viscosity in Euler and Navier-Stokes solvers

A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small time-dependent perturbation about a steady inviscid flow over an NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 deg. The method has been applied to inviscid flows here, but as discussed is also applicable to viscous flows. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundaries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time-marching scheme are also presented. The eigenvalue procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of artificial viscosity.

Mahajan, Aparajit J.↗

NASA NTRS1992

Optimum design of a gearbox for low vibration

A computer program was developed for designing a low vibration gearbox. The code is based on a finite element shell analysis, a modal analysis, and a structural optimization method. In the finite element analysis, a triangular shell element with 18 degrees-of-freedom is used. In the optimization method, the overall vibration energy of the gearbox is used as the objective function and is minimized at the exciting frequency by varying the finite element thickness. Modal analysis is used to derive the sensitivity of the vibration energy with respect to the design variable. The sensitivity is representative of both eigenvalues and eigenvectors. The optimum value is computed by the gradient projection method and a unidimensional search procedure under the constraint condition of constant weight. The computer code is applied to a design problem derived from an experimental gearbox in use at the NASA Lewis Research Center. The top plate and two side plates of the gearbox are redesigned and the contribution of each surface to the total vibration is determined. Results show that optimization of the top plate alone is effective in reducing total gearbox vibration.

Inoue, Katsumi↗

NASA NTRS2001

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one-dimensional convection-diffusion-reaction (CDR) equations. First, accuracy, stability, conservation, and dense output are considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, N = 2, additive Runge-Kutta ARK2 methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms are of equal order to those of the elemental methods. Derived ARK2 methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z(exp [I]) goes to infinity, and retain high stability efficiency in the absence of stiffness, z(exp [I]) goes to zero. Extrapolation-type stage-value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK2 error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK2 methods while offering error control and stage-value predictors.

Kennedy, Christopher A.↗