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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1992

A numerical method for systems of conservation laws of mixed type admitting hyperbolic flux splitting

The present treatment of elliptic regions via hyperbolic flux-splitting and high order methods proposes a flux splitting in which the corresponding Jacobians have real and positive/negative eigenvalues. While resembling the flux splitting used in hyperbolic systems, the present generalization of such splitting to elliptic regions allows the handling of mixed-type systems in a unified and heuristically stable fashion. The van der Waals fluid-dynamics equation is used. Convergence with good resolution to weak solutions for various Riemann problems are observed.

Shu, Chi-Wang↗

NASA NTRS2000

Low-dimensional Representation of Error Covariance

Ensemble and reduced-rank approaches to prediction and assimilation rely on low-dimensional approximations of the estimation error covariances. Here stability properties of the forecast/analysis cycle for linear, time-independent systems are used to identify factors that cause the steady-state analysis error covariance to admit a low-dimensional representation. A useful measure of forecast/analysis cycle stability is the bound matrix, a function of the dynamics, observation operator and assimilation method. Upper and lower estimates for the steady-state analysis error covariance matrix eigenvalues are derived from the bound matrix. The estimates generalize to time-dependent systems. If much of the steady-state analysis error variance is due to a few dominant modes, the leading eigenvectors of the bound matrix approximate those of the steady-state analysis error covariance matrix. The analytical results are illustrated in two numerical examples where the Kalman filter is carried to steady state. The first example uses the dynamics of a generalized advection equation exhibiting nonmodal transient growth. Failure to observe growing modes leads to increased steady-state analysis error variances. Leading eigenvectors of the steady-state analysis error covariance matrix are well approximated by leading eigenvectors of the bound matrix. The second example uses the dynamics of a damped baroclinic wave model. The leading eigenvectors of a lowest-order approximation of the bound matrix are shown to approximate well the leading eigenvectors of the steady-state analysis error covariance matrix.

Tippett, Michael K.↗

NASA NTRSDate not supplied

Variational Coupled Loads Analysis using the Hybrid Parametric Variation Method

Time-domain coupled loads analysis (CLA)is used to determine the response of a launch vehicle and payload system to transient forces, such as liftoff, engine ignitions and shutdowns, jettison events, and atmospheric flight loads, such as buffet. CLA, using Hurty/Craig-Bampton (HCB)component models, is the accepted method for the establishment of design-level loads for launch systems. However, uncertainty in the component models flows into uncertainty in predicted system results. Uncertainty in the structural responses during launch is a significant concern because small variations in launch vehicle and payload mode shapes and their interactions can result in significant variations in system loads. Uncertainty quantification (UQ)is used to determine statistical bounds on prediction accuracy based on model uncertainty. In this paper uncertainty is treated at the HCB component-model level. In an effort to account for model uncertainties and statistically bound their effect on CLA predictions, this work combines CLA with UQ in a process termed variational coupled loads analysis (VCLA). The modeling of uncertainty using a parametric approach, in which input parameters are represented by random variables, is common, but its major drawback is the resulting uncertainty is limited to the form of the nominal model. Uncertainty in model form is one of the biggest contributors to uncertainty in complex built-up structures. Model-form uncertainty can be represented using a nonparametric approach based on random matrix theory (RMT). In this work, UQ is performed using the hybrid parametric variation (HPV)method, which combines parametric with nonparametric uncertainty at the HCB component model level. The HPV method requires the selection of dispersion values for the HCB fixed-interface (FI)eigenvalues, and the HCB mass and stiffness matrices. The dispersions are based upon component test-analysis modal correlation results. During VCLA, random component models are assembled into an ensemble of random systems using a Monte Carlo (MC)approach. CLA is applied to each of the ensemble members to produce an ensemble of system-level responses for statistical analysis. The proposed methodology is demonstrated through its application to a buffet loads analysis of NASA’s Space Launch System (SLS)during the transonic regime fifty seconds after liftoff. Core stage (CS)section shears and moments are recovered, and statistics are computed.

Uncertainty Quantification↗

NASA NTRS1977

Extraction of the large eddy structure of a turbulent boundary layer

Much attention is currently being focused on turbulence modeling by separation of large-scale and small-scale turbulence. One major difficulty with this method is the problem of separating and defining the large-scale structure. In this paper, a method is presented which does this separation by analyzing time-averaged, two-point double velocity correlation measurements. The basis of the method is a maximization of the inner product of a candidate velocity field with the measured turbulent velocity field. An eigenvalue problem results, whose solution provides a set of orthogonal eigenfunctions associated with specific modes of turbulent motion. The large-scale eddy, containing most of the turbulent energy, is associated with the first (largest) mode. The turbulent velocity field is then represented as a sum of the eigenfunctions. A short description of the equations is given, the methods of solution are discussed, and some results for a two-dimensional fully developed turbulent boundary layer are presented.

Lemmerman, L. A.↗

NASA NTRS1982

Calculation of exact vibration modes for plane grillages by the dynamic stiffness method

A dynamic stiffness method is developed for the calculation of the exact modal parameters for plane grillages which consist of straight and uniform beams with coincident elastic and inertial axes. Elementary bending-torsion beam theory is utilized, and bending translation is restricted to one direction. The exact bending-torsion dynamic stiffness matrix is obtained for a straight and uniform beam element with coincident elastic and inertial axes. The element stiffness matrices are assembled using the standard procedure of the static stiffness method to form the dynamic stiffness matrix of the complete grillage. The exact natural frequencies, mode shapes, and generalized masses of the grillage are then calculated by solving a nonlinear eigenvalue problem based on the dynamic stiffness matrix. The exact modal solutions for an example grillage are calculated and compared with the approximate solutions obtained by using the finite element method.

Hallauer, W. L., Jr.↗

NASA NTRS1984

Vibration and buckling of general periodic lattice structures

A method is presented for vibration and buckling analysis of arbitrary lattice structures having repetitive geometry in any combination of coordinate directions. The approach is based on exact member theory for representing the stiffness of an individual member subject to axial load, and in the case of vibration, undergoing harmonic oscillation. The method is an extension of previous work that was limited to specific geometries. The resulting eigenvalue problem is of the size associated with the repeating element of the structure. A computer program has been developed incorporating the theory and results are given for vibration of rectangular platforms and a large antenna structure having rotational symmetry. Buckling and vibration results for cable-stiffened rings are also given.

Anderson, M. S.↗

NASA NTRS1989

Finite element large-amplitude free and forced vibrations of rectangular thin composite plates

A finite element formulation is presented for determining the large-amplitude free and steady-state forced vibration response of arbitrarily laminated anisotropic composite rectangular thin plates. The nonlinear stiffness and harmonic force matrices of an arbitrarily laminated composite rectangular plate element are developed for nonlinear free and forced vibration analyses. The linearized updated-mode method with nonlinear time function approximation is employed for the solution of the system nonlinear eigenvalue equations. The amplitude-frequency relations for convergence with gridwork refinement, different boundary conditions, aspect ratios, lamination angles and number of plies are presented. The finite element results are compared with available approximate continuum solutions.

Chiang, C. K.↗

NASA NTRS1992

Quantum spatial propagation of squeezed light in a degenerate parametric amplifier

Differential equations which describe the steady state spatial evolution of nonclassical light are established using standard quantum field theoretic techniques. A Schroedinger equation for the state vector of the optical field is derived using the quantum analog of the slowly varying envelope approximation (SVEA). The steady state solutions are those that satisfy the time independent Schroedinger equation. The resulting eigenvalue problem then leads to the spatial propagation equations. For the degenerate parametric amplifier this method shows that the squeezing parameter obey nonlinear differential equations coupled by the amplifier gain and phase mismatch. The solution to these differential equations is equivalent to one obtained from the classical three wave mixing steady state solution to the parametric amplifier with a nondepleted pump.

Deutsch, Ivan H.↗

NASA NTRS1993

Nonlinear vibrations of thin arbitrarily laminated composite plates subjected to harmonic excitations using DKT elements

A finite element formulation is presented for determining the large-amplitude free and steady-state forced vibration response of arbitrarily laminated anisotropic composite thin plates using the Discrete Kirchhoff Theory (DKT) triangular elements. The nonlinear stiffness and harmonic force matrices of an arbitrarily laminated composite triangular plate element are developed for nonlinear free and forced vibration analyses. The linearized updated-mode method with nonlinear time function approximation is employed for the solution of the system nonlinear eigenvalue equations. The amplitude-frequency relations for convergence with gridwork refinement, triangular plates, different boundary conditions, lamination angles, number of plies, and uniform versus concentrated loads are presented.

Chiang, C. K.↗

DOE OSTI2026

Laplace Transform–Based Quantum Eigenvalue Transformation via Linear Combination of Hamiltonian Simulation

Eigenvalue transformations, which include solving time-dependent differential equations as a special case, have a wide range of applications in scientific and engineering computation. While quantum algorithms for singular value transformations are well studied, eigenvalue transformations are distinct, especially for nonnormal matrices. Here, we propose an efficient quantum algorithm for performing a class of eigenvalue transformations that can be expressed as a certain type of matrix Laplace transformation. This allows us to significantly extend the recently developed linear combination of Hamiltonian simulation method [D. An, J.-P. Liu, and L. Lin, Phys. Rev. Lett., 131 (2023), 150603; D. An, A. M. Childs, and L. Lin, Commun. Math. Phys. 407, 19 (2026)] to represent a wider class of eigenvalue transformations, such as powers of the matrix inverse, 𝐴 −𝑘 , and the exponential of the matrix inverse, 𝑒 −𝐴 −1 . The latter can be interpreted as the solution of a mass-matrix differential equation of the form form 𝐴⁢𝑢′⁡⁡(𝑡) =−𝑢⁡(𝑡). We demonstrate that our eigenvalue transformation approach can solve this problem without explicitly inverting 𝐴, thereby reducing the computational complexity.

Laplace transform↗

NASA NTRS1988

Linear quadratic regulators with eigenvalue placement in a specified region

A linear optimal quadratic regulator is developed for optimally placing the closed-loop poles of multivariable continuous-time systems within the common region of an open sector, bounded by lines inclined at + or - pi/2k (k = 2 or 3) from the negative real axis with a sector angle of pi/2 or less, and the left-hand side of a line parallel to the imaginary axis in the complex s-plane. The design method is mainly based on the solution of a linear matrix Liapunov equation, and the resultant closed-loop system with its eigenvalues in the desired region is optimal with respect to a quadratic performance index.

Shieh, Leang S.↗

NASA NTRS1982

Integrating matrix formulations for vibrations of rotating beams including the effects of concentrated masses

By expressing partial differential equations of motion in matrix notation, utilizing the integrating matrix as a spatial operator, and applying the boundary conditions, the resulting ordinary differential equations can be cast into standard eigenvalue form upon assumption of the usual time dependence. As originally developed, the technique was limited to beams having continuous mass and stiffness properties along their lengths. Integrating matrix methods are extended to treat the differential equations governing the flap, lag, or axial vibrations of rotating beams having concentrated masses. Inclusion of concentrated masses is shown to lead to the same kind of standard eigenvalue problem as before, but with slightly modified matrices.

Lakin, W. D.↗

NASA NTRS1979

The application of finite element techniques to acoustic transmission in lined ducts with flow

The finite element method (FEM) is used to analyze the propagation of sound in two-dimensional nonuniform ducts carrying a compressible subsonic mean flow. Galerkin and residual least squares (RLS) methods with natural and forced boundary conditions are considered. The accuracy of FEM results for the eigenvalue and transmission problems is assessed by comparison with alternative numerical schemes for nonuniform ducts. The results presented and those from associated investigations indicate that modal coupling is a significant feature of the acoustic field, especially at high Mach numbers. A multimodal model therefore appears to be essential if any reliable conclusions are to be drawn in the context of turbofan inlet regions. Improvements to the eigenvalue scheme following the implementation of higher-order Hermitian elements indicate a similar modification for the transmission problem.

Astley, R. J.↗

NASA NTRS1979

Principal component analysis of solar flares in the soft X-ray flux

The paper considers principal component analysis of solar flares in the soft X-ray flux, a technique for extracting the salient features from a mass of data. The method applies particularly to the analysis of nonstationary ensembles, and its computations require the evaluation of eigenvalues of matrices. The Eispack matrix eigen system routines were used to analyze full-disk proportional-counter data from the X-ray event analyzer which was part of the Skylab experiment. Empirical orthogonal functions were derived for events in the soft X-ray spectrum between 2.5 and 20 A during different time periods, indicating that about 90% of the cumulative power of each analyzed flare is contained in the largest eigenvector. The first two largest eigenvectors are sufficient for an empirical curve fit through the raw data and a characterization of solar flares in the soft X-ray flux, and power spectra of two largest eigenvectors reveal a reported periodicity of about 5 min.

Teuber, D. L.↗

NASA NTRS1989

On the role of artificial viscosity in 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 the airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small perturbation about a steady inviscid flow over a NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 degrees. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time marching scheme are also presented. This procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of terms for artificial viscosity.

Mahajan, Aparajit J.↗

NASA NTRS1991

Cascade flutter analysis with transient response aerodynamics

Two methods for calculating linear frequency domain aerodynamic coefficients from a time marching Full Potential cascade solver are developed and verified. In the first method, the Influence Coefficient, solutions to elemental problems are superposed to obtain the solutions for a cascade in which all blades are vibrating with a constant interblade phase angle. The elemental problem consists of a single blade in the cascade oscillating while the other blades remain stationary. In the second method, the Pulse Response, the response to the transient motion of a blade is used to calculate influence coefficients. This is done by calculating the Fourier Transforms of the blade motion and the response. Both methods are validated by comparison with the Harmonic Oscillation method and give accurate results. The aerodynamic coefficients obtained from these methods are used for frequency domain flutter calculations involving a typical section blade structural model. An eigenvalue problem is solved for each interblade phase angle mode and the eigenvalues are used to determine aeroelastic stability. Flutter calculations are performed for two examples over a range of subsonic Mach numbers.

Bakhle, Milind A.↗

NASA NTRS1991

Cascade flutter analysis with transient response aerodynamics

Two methods for calculating linear frequency domain aerodynamic coefficients from a time-marching Full-Potential cascade solver are developed and verified. In the first method, the Influence Coefficient method, solutions to elemental problems are superposed to obtain the solutions for a cascade in which all blades are vibrating with a constant interblade phase angle. The elemental problem consists of a single blade in the cascade oscillating while the other blades remain stationary. In the second method, the Pulse Response method, the response to the transient motion of a blade is used to calculate influence coefficients. This is done by calculating the Fourier transforms of the blade motion and the response. Both methods are validated by comparison with the Harmonic Oscillation method and give accurate results. The aerodynamic coefficients obtained from these methods are used for frequency domain flutter calculations involving a typical section blade structural model. An eigenvalue problem is solved for each interblade phase angle mode and the eigenvalues are used to determine aeroelastic stability. Flutter calculations are performed for two examples over a range of subsonic Mach numbers using both flat plates and actual airfoils.

Bakhle, Milind A.↗

NASA NTRS1991

Cascade flutter analysis with transient response aerodynamics

Bakhle, M. A.↗