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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

The eigenvalues of the pseudospectral Fourier approximation to the operator sin (2x) d/dx

It is shown that the eigenvalues Z sub i of the pseudospectral Fourier approximation to the operator sin(2x) curly d/curly dx satisfy (R sub e) (Z sub i) = + or - 1 or (R sub e)(Z sub I) = 0. Whereas this does not prove stability for the Fourier method, applied to the hyperbolic equation U sub t = sin (2x)(U sub x) - pi x pi; it indicates that the growth in time of the numerical solution is essentially the same as that of the solution to the differential equation.

Tal-Ezer, H.↗
NASA NTRS1987

Calculation of stresses and forces between the skin and stiffener in composite panels

A methodology for computing the stresses at the interface between the skin and stiffener in stiffened composite panels is described. The methodology is based on finite-element analyses and on elasticity solution. The finite-element analyses are standard, while the elasticity solution is based on eigenvalue expansions of the stress function. The eigenvalue expansions are assumed to be valid in the local region where the stiffener flange terminates. The local elasticity solution is coupled to the global finite-element analysis using collocation on the boundary of the local region. Accuracy and convergence of the method are discussed and several examples of its utility are presented.

Hyer, M. W.↗
NASA NTRS1985

Eigenvalue Placement and Stabilization by Constrained Optimization

A pole placement algorithm is proposed which uses constrained nonlinear optimization techniques on a finite dimensional model of a linear n degree of freedom system. Low order feedback control is assumed where r poles may be assigned; r being the rank of the sensor coefficient matrix. It is shown that by combining feedback control theory methods with optimization techniques, one can ensure the stability characteristics of a system, and can alter its transient response.

Decaro, S. M.↗
NASA NTRS1989

Parallel Eigenvalue extraction

A new numerical algorithm for the solution of large-order eigenproblems typically encountered in linear elastic finite element systems is presented. The architecture of parallel processing is utilized in the algorithm to achieve increased speed and efficiency of calculations. The algorithm is based on the frontal technique for the solution of linear simultaneous equations and the modified subspace eigenanalysis method for the solution of the eigenproblem. Assembly, elimination and back-substitution of degrees of freedom are performed concurrently, using a number of fronts. All fronts converge to and diverge from a predefined global front during elimination and back-substitution, respectively. In the meantime, reduction of the stiffness and mass matrices required by the modified subspace method can be completed during the convergence/divergence cycle and an estimate of the required eigenpairs obtained. Successive cycles of convergence and divergence are repeated until the desired accuracy of calculations is achieved. The advantages of this new algorithm in parallel computer architecture are discussed.

Akl, Fred A.↗
NASA NTRS1988

An upwind differencing scheme for the time-accurate incompressible Navier-Stokes equations

The two-dimensional incompressible Navier-Stokes equations are solved in a time-accurate manner in using the method of pseudocompressibility. Using this method, subiterations in pseudotime are required to satisfy the continuity equation at each time step. An upwind differencing scheme based on flux-difference splitting is used to compute the convective terms. The upwind differencing is biased based on the sign of the local eigenvalue of the Jacobian matrix. Third-order or fifth-order spatial accuracy is maintained throughout the interior grid points. The equations are solved using an implicit line-relaxation scheme. This solution scheme is stable and is capable of running at large time steps in pseudotime, leading to fast convergence for each physical time step. A variety of computed results are presented to validate the present scheme. Results for the flow over an oscillating plate are compared with the exact analytic solution, good agreement is seen. Excellent comparison is obtained between the computed solution and the analytical results for inviscid channel flow with an oscillating back pressure. Flow solutions over a circular cylinder with vortex shedding are also presented. Finally, the flow past an airfoil at -90 deg angle-of-attack is also computed.

Rogers, Stuart E.↗
NASA NTRS1988

Attitude determination and parameter estimation using vector observations

Procedures for attitude determination based on Wahba's loss function are generalized to include the estimation of parameters other than the attitude, such as sensor biases. Optimization with respect to the attitude requires either the singular value decomposition of a 3x3 matrix or finding the maximum eigenvalue and corresponding eigenvector of a 4x4 symmetric matrix, but does not require an a priori estimate of the attitude. Optimization with respect to the other parameters employs an iterative approach, which does require an a priori estimate of these parameters. Conventional state estimation methods require a priori estimates of both the parameters and the attitude, while the algorithms presented in this paper always compute the exact optimal attitude for given values of the parameters. The proposed method is shown to give the correct solution of an example problem. An expression for the covariance of the attitude and parameter estimates is derived.

Markley, F. Landis↗
NASA NTRS1996

A Computer Program for the Computation of Running Gear Temperatures Using Green's Function

A new technique has been developed to study two dimensional heat transfer problems in gears. This technique consists of transforming the heat equation into a line integral equation with the use of Green's theorem. The equation is then expressed in terms of eigenfunctions that satisfy the Helmholtz equation, and their corresponding eigenvalues for an arbitrarily shaped region of interest. The eigenfunction are obtalned by solving an intergral equation. Once the eigenfunctions are found, the temperature is expanded in terms of the eigenfunctions with unknown time dependent coefficients that can be solved by using Runge Kutta methods. The time integration is extremely efficient. Therefore, any changes in the time dependent coefficients or source terms in the boundary conditions do not impose a great computational burden on the user. The method is demonstrated by applying it to a sample gear tooth. Temperature histories at representative surface locatons are given.

GREENS FUNCTIONS↗
NASA NTRS1993

Optimal design of imperfect, anisotropic, ring-stiffened cylinders under combined loads

Development of an algorithm to perform the optimal sizing of buckling resistant, imperfect, anisotropic ring-stiffened cylinders subjected to axial compression, torsion, and internal pressure is presented. The enforcement of stability constraints is treated in a way that does not require any eigenvalue analysis. Assumption of criticality of these stability constraints during the optimal sizing of the cylinders produced designs that nevertheless satisfied all of the stress constraints as well as the stability constraints. Case studies performed using a combination of penalty function and feasible direction optimization methods indicate that the presence of the axisymmetric initial imperfection in the cylinder wall can significantly affect the optimal designs. Weight savings associated with the addition of two rings to the unstiffened cylinder and/or the addition of internal pressure is substantial when torsion makes up a significant fraction of the combined load state.

Ley, Robert P.↗
NASA NTRS2013

Mixture-Tuned, Clutter Matched Filter for Remote Detection of Subpixel Spectral Signals

Mapping localized spectral features in large images demands sensitive and robust detection algorithms. Two aspects of large images that can harm matched-filter detection performance are addressed simultaneously. First, multimodal backgrounds may thwart the typical Gaussian model. Second, outlier features can trigger false detections from large projections onto the target vector. Two state-of-the-art approaches are combined that independently address outlier false positives and multimodal backgrounds. The background clustering models multimodal backgrounds, and the mixture tuned matched filter (MT-MF) addresses outliers. Combining the two methods captures significant additional performance benefits. The resulting mixture tuned clutter matched filter (MT-CMF) shows effective performance on simulated and airborne datasets. The classical MNF transform was applied, followed by k-means clustering. Then, each cluster s mean, covariance, and the corresponding eigenvalues were estimated. This yields a cluster-specific matched filter estimate as well as a cluster- specific feasibility score to flag outlier false positives. The technology described is a proof of concept that may be employed in future target detection and mapping applications for remote imaging spectrometers. It is of most direct relevance to JPL proposals for airborne and orbital hyperspectral instruments. Applications include subpixel target detection in hyperspectral scenes for military surveillance. Earth science applications include mineralogical mapping, species discrimination for ecosystem health monitoring, and land use classification.

Thompson, David R.↗
NASA NTRS1987

On the acceleration of charged particles at relativistic shock fronts

The diffusive acceleration of highly relativistic particles at a shock is reconsidered. Using the same physical assumptions as Blandford and Ostriker (1978), but dropping the restriction to nonrelativistic shock velocities, the authors find approximate solutions of the particle kinetic equation by generalizing the diffusion approximation to higher order terms in the anisotropy of the particle distribution. The general solution of the transport equation on either side of the shock is constructed, which involves the solution of an eigenvalue problem. By matching the two solutions at the shock, the spectral index of the resulting power law is found by taking into account a sufficiently large number of eigenfunctions. Low-order truncation corresponds to the standard diffusion approximation and to a somewhat more general method described by Peacock (1981). In addition to the energy spectrum, the method yields the angular distribution of the particles and its spatial dependence.

Kirk, J. G.↗
NASA NTRS1985

The exact eigenfunctions and eigenvalues of a two-dimensional rigid rotor obtained using Gaussian wave packet dynamics

Exact eigenfunctions for a two-dimensional rigid rotor are obtained using Gaussian wave packet dynamics. The wave functions are obtained by propagating, without approximation, an infinite set of Gaussian wave packets that collectively have the correct periodicity, being coherent states appropriate to this rotational problem. This result leads to a numerical method for the semiclassical calculation of rovibrational, molecular eigenstates. Also, a simple, almost classical, approximation to full wave packet dynamics is shown to give exact results: this leads to an a posteriori justification of the De Leon-Heller spectral quantization method.

Reimers, J. R.↗
NASA NTRS2014

Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes a hybrid approach for determining physical displacements by augmenting the original quadratically constrained least squares (LSQI) algorithm with Direct Shape Mapping (DSM) and modifying the energy constraints. The approach presented is applicable to simulation of the elastic behavior of launch vehicles and other structures that utilize discrete LTI finite element model (FEM) derived mode sets (eigenvalues and eigenvectors) that are propagated throughout time. The time invariant nature of the elastic data presents a problem of how to properly transition elastic states from the prior to the new model while preserving motion across the transition and ensuring there is no truncation or excitation of the system. A previous approach utilizes a LSQI algorithm with an energy constraint to effect smooth transitions between eigenvector sets with no requirement that the models be of similar dimension or have any correlation. This approach assumes energy is conserved across the transition, which results in significant non-physical transients due to changing quasi-steady state energy between mode sets, a phenomenon seen when utilizing a truncated mode set. The computational burden of simulating a full mode set is significant so a subset of modes is often selected to reduce run time. As a result of this truncation, energy between mode sets may not be constant and solutions across transitions could produce non-physical transients. In an effort to abate these transients an improved methodology was developed based on the aforementioned approach, but this new approach can handle significant changes in energy across mode set transitions. It is proposed that physical velocities due to elastic behavior be solved for using the LSQI algorithm, but solve for displacements using a two-step process that independently addresses the quasi-steady-state and non-steady-state contributions to the elastic displacement. For structures subject to large external forces, such as thrust or atmospheric drag, it is imperative to capture these forces when solving for elastic displacement. To simplify the mathematical formulation, assumptions are made regarding mass matrix normalization, constant external forcing, and constant viscous damping. These simplifications allow for direct solutions to the quasi-steady-state displacements through a process titled Direct Shape Mapping. DSM solves for the displacements using the eigenvalues of the elastic modes and the external forcing and returns a set of elastic displacements dictated by the eigenvectors of the post-transition mode set. For the non-steady-state contributions to displacement we formulate a LSQI problem that is constrained by energy of the non-steady state terms. The contributions from the quasi-steady-state and non-steady state solutions are then combined to obtain the physical displacements associated with the new set of eigenvectors. Results for the LSQI-DSM approach show significant reduction/complete removal of transients across mode set transitions while maintaining elastic motion from the prior state. For time propagation applications employing discrete elastic models that need to be transitioned in time and where running with full a full mode set is not feasible, the method developed offers a practical solution to simulating vehicle elasticity.

Hannan, Mike R.↗
NASA NTRS1973

Cluster Seeking Techniques in Pattern Recognition

A cluster seeking technique is defined as a method of dividing data into subsets, called clusters. These clusters contain data points that are similar to each other and different from the elements of other clusters. Various cluster seeking techniques were broken down into seven categories: (1) probabilistic, (2) signal detection, (3) clustering, (4) clumping, (5) eigenvalue, (6) minimal mode seeking, and (7) miscellaneous. Each category is described and one or more algorithms of that type are presented.

Barr, B. J.↗
NASA NTRS1986

A class of fast staircase algorithms for generalized state-space systems

Several methods are presented for defining generalized state space models (GSSM), with emphasis on 'fast' techniques for transforming GSSMs to condensed state space models (SSM). The fast forms are configured to yield invariant transfer functions. Details of the decomposition process are summarized in terms of the determination of eigenvalues which are separated within staircase matrices. Applications of the techniques are illustrated through implementation of a unimodular transformation which does not affect the finite eigenvalues, by demonstrating deadbeat control of a GSSM, which the definition of a reduced observer of a GSSM, and by embedding a polynomial matrix into a unimodular matrix.

Beelen, T.↗
NASA NTRS2008

Modeling State-Space Aeroelastic Systems Using a Simple Matrix Polynomial Approach for the Unsteady Aerodynamics

A simple matrix polynomial approach is introduced for approximating unsteady aerodynamics in the s-plane and ultimately, after combining matrix polynomial coefficients with matrices defining the structure, a matrix polynomial of the flutter equations of motion (EOM) is formed. A technique of recasting the matrix-polynomial form of the flutter EOM into a first order form is also presented that can be used to determine the eigenvalues near the origin and everywhere on the complex plane. An aeroservoelastic (ASE) EOM have been generalized to include the gust terms on the right-hand side. The reasons for developing the new matrix polynomial approach are also presented, which are the following: first, the "workhorse" methods such as the NASTRAN flutter analysis lack the capability to consistently find roots near the origin, along the real axis or accurately find roots farther away from the imaginary axis of the complex plane; and, second, the existing s-plane methods, such as the Roger s s-plane approximation method as implemented in ISAC, do not always give suitable fits of some tabular data of the unsteady aerodynamics. A method available in MATLAB is introduced that will accurately fit generalized aerodynamic force (GAF) coefficients in a tabular data form into the coefficients of a matrix polynomial form. The root-locus results from the NASTRAN pknl flutter analysis, the ISAC-Roger's s-plane method and the present matrix polynomial method are presented and compared for accuracy and for the number and locations of roots.

Pototzky, Anthony S.↗
NASA NTRS1972

Sound propagation in duct shear layers.

Consideration of the propagation of sound in a two-dimensional inviscid shear layer for given initial sound pressure profiles. The eigenvalue problem resulting from an assumption of separable solutions in a form first obtained by Pridmore-Brown (1958) is solved numerically. It is shown that although the resulting eigenfunctions cannot be proven to be orthogonal or complete, they can be combined by a least total error squared method to give a good representation of the initial pressure profile. The acoustic pressure in the duct is then easily calculated. The results verify all the predictions made in an earlier perturbation calculation. Moreover, they show that the refraction effect gets saturated at high subsonic Mach numbers and at high frequencies. The proposed technique may be used, if necessary, with impedance boundary conditions at the duct walls.

Shankar, P. N.↗
NASA NTRS1973

Studies in tilt rotor VTOL aircraft aeroelasticity, volume 2

Two methods for natural mode vibration analysis are discussed. The first consists of a direct approach based on a finite element representation of the complete structure as an entity. The mass and stiffness matrices for the complete structure are assembled by properly combining the mass and stiffness matrices of the individual elements into which the structure has been divided. The second approach is that of component mode synthesis. This method is based on the concept of synthesizing the natural modes of the complete structure from modes of conveniently difined substructures, or components, into which the structure has been partitioned. In this way the expedient of reducing the system degrees of freedom, and thus the size of the eigenvalue problem, can be introduced by partial modal synthesis.

Kvaternik, R. G.↗
NASA NTRS1972

A critical discussion of theories of flame spread across solid and liquid fuels.

Mathematical descriptions of flames spreading over liquid and solid fuels are obtained, using basic assumptions derived from observations or physical reasoning. A review of existing theories shows that they are incomplete in that they either treat an uncoupled problem of the condensed phase where the spreading rate and heat flux at the surface are given, or they merely determine the spreading rate in terms of a new vaguely defined eigenvalue. An important difference between the liquid and solid cases, due to convection, is pointed out, and it is shown that solid-fuel flame-spread theories which claim to apply to the liquid case, in reality do not apply to it. A mathematical formulation and a method of solution are presented for the phenomenon of flame spread over solid fuels with forward heat conduction in both the solid and the gas. The method uses an energy integral over the field to determine the spreading rate in terms of the basic properties of the fuel and air.

Sirignano, W. A.↗