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 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Nonlinear structural dynamics analysis using a modified modal method

The procedure for predicting the nonlinear dynamic response of structural components subjected to a step loading is presented. The procedure is a modified modal method that involves a change of dependent variables from the unknown nodal degrees of freedom of the finite element model of the structure to a smaller set of generalized coordinates. This change of dependent variables uses a combination of the nonlinear static solution and some selected vibration mode shapes. The vibration mode shapes correspond to the eigenvectors obtained by solving a standard free vibration eigenvalue problem wherein the stiffness matrix is expanded about the nonlinear static solution. A strategy is also presented for determining which and how many vibration mode shapes to include in the transformation. The effect of inaccurate representation of the spatial distribution of the applied load on the nonlinear dynamic response is discussed for two classes of structural behavior. Application of the procedure to structures which exhibit a stiffening behavior and to those with a softening behavior is presented.

Knight, N. F., Jr.↗

NASA NTRS1985

Numerical Techniques for Scattering from Submerged Objects

To represent the final results in terms of matrices, one expands all appropriate physical quantities in terms of partial wave basis states. This includes expansions for the incident and scattered fields and the surface quantities. The method then utilizes the Huygen-Poincare integral representation for both the exterior and interior solutions, leading to the required matrix equations. One thus deals with matrix equations, the complexity of which depends on the nature of the problem. It is shown that in general a transition matrix T can be obtained relating the incident field A with the scattered field f having the form T = PQ(-1), where f = TA. The structure of Q can be quite complicated and can itself be composed of other matrix inversions such as arise from layered objects. Recent improvements in this method appropriate for a variety of physical problems are focused on, and on their implementation. Results are outlined from scattering simulations for very elongated submerged objects and resonance scattering from elastic solids and shells. The final improvement concerns eigenfunction expansions of surface terms, arising from solution of the interior problem, obtained via a preconditioning technique. This effectively reduces the problem to that of obtaining eigenvalues of a Hermitian operator. This formalism is reviewed for scattering from targets that are rigid, sound-soft, acoustic, elastic solids, elastic shells, and elastic layered objects. Two sets of the more interesting results are presented. The first concerns scattering from elongated objects, and the second to thin elastic spheroids.

Werby, M. F.↗

NASA NTRS1976

A FORTRAN program for the analysis of linear continuous and sample-data systems

A FORTRAN digital computer program which performs the general analysis of linearized control systems is described. State variable techniques are used to analyze continuous, discrete, and sampled data systems. Analysis options include the calculation of system eigenvalues, transfer functions, root loci, root contours, frequency responses, power spectra, and transient responses for open- and closed-loop systems. A flexible data input format allows the user to define systems in a variety of representations. Data may be entered by inputing explicit data matrices or matrices constructed in user written subroutines, by specifying transfer function block diagrams, or by using a combination of these methods.

Edwards, J. W.↗

NASA NTRS1978

A simple method for designing structural models with closely spaced modes of vibration

A simple method for designing a mathematical model with closely spaced vibration modes is described. The design process begins with a reference model having specified geometry, continuous inertia and stiffness distributions, and degrees of freedom, all of which remain unchanged. Two natural frequencies of this model are then forced together by means of systematic perturbation of the model's discrete inertia and stiffness parameters. There is only one eigenvalue solution per design cycle, and the gradient vector is calculated directly from the resulting modal quantities. The minimization procedure employed is unconstrained. As applications, a cantilevered plane grid model with five degrees of freedom and a bending-torsion-oscillator with eleven degrees of freedom are treated.

Hallauer, W. L., Jr.↗

NASA NTRS1981

Multigrid solvers on parallel computers

Massively parallel computers, as considered in this investigation, are not yet available. However, a large-scale parallel computer cannot usefully be designed before the hypothetical algorithms which will employ it are studied. Most of the studies of parallel partial differential equations (PDE) solvers are based on solution techniques much slower (on sequential machines) than multigrid methods. Multigrid methods are highly parallelizable. Each of their processes can simultaneously be performed at all grid points. The present investigation is concerned with a preliminary exploration of the potential of multigrid, or, more generally, Multi-Level Adaptive Techniques (MLAT) on computers with many processors. Basic processes are considered, taking into account coarse-grid approximation, relaxation, coarse-grid corrections, full multigrid algorithms, nonlinear problems and eigenvalue problems, fine-to-coarse correction, and chains of problems. Details of parallel multigrid processing are also examined.

Brandt, A.↗

NASA NTRS1989

Sequential design of discrete linear quadratic regulators via optimal root-locus techniques

A sequential method employing classical root-locus techniques has been developed in order to determine the quadratic weighting matrices and discrete linear quadratic regulators of multivariable control systems. At each recursive step, an intermediate unity rank state-weighting matrix that contains some invariant eigenvectors of that open-loop matrix is assigned, and an intermediate characteristic equation of the closed-loop system containing the invariant eigenvalues is created.

Shieh, Leang S.↗

NASA NTRS1986

Preconditioning for first-order spectral discretization

Efficient solution of the equations from spectral discretizations is essential if the high-order accuracy of these methods is to be realized. Direct solution of these equations is rarely feasible, thus iterative techniques are required. A preconditioning scheme for first-order Chebyshev collocation operators is proposed herein, in which the central finite difference mesh is finer than the collocation mesh. Details of the proper techniques for transferring information between the meshes are given here, and the scheme is analyzed by examination of the eigenvalue spectra of the preconditioned operators. The effect of artificial viscosity required in the inversion of the finite difference operator is examined. A second preconditioning scheme, involving a high-order upwind finite difference operator of the van Leer type is also analyzed to provide a comparison with the present scheme. Finally, the performance of the present scheme is verified by application to several test problems.

Streett, C. L.↗

NASA NTRS1987

Aeroelastic effects in multirotor vehicles. Part 2: Methods of solution and results illustrating coupled rotor/body aeromechanical stability

This report is a sequel to the earlier report titled, Aeroelastic Effects in Multi-Rotor Vehicles with Application to Hybrid Heavy Lift System, Part 1: Formulation of Equations of Motion (NASA CR-3822). The trim and stability equations are presented for a twin rotor system with a buoyant envelope and an underslung load attached to a flexible supporting structure. These equations are specialized for the case of hovering flight. A stability analysis, for such a vehicle with 31 degrees of freedom, yields a total of 62 eigenvalues. A careful parametric study is performed to identify the various blade and vehicle modes, as well as the coupling between various modes. Finally, it is shown that the coupled rotor/vehicle stability analysis provides information on both the aeroelastic stability as well as complete vehicle dynamic stability. Also presented are the results of an analytical study aimed at predicting the aeromechanical stability of a single rotor helicopter in ground resonance. The theoretical results are found to be in good agreement with the experimental results, thereby validating the analytical model for the dynamics of the coupled rotor/support system.

Venkatesan, C.↗

NASA NTRS1972

Model potential calculations of lithium transitions.

Semi-empirical potentials are constructed that have eigenvalues close in magnitude to the binding energies of the valence electron in lithium. The potentials include the long range polarization force between the electron and the core. The corresponding eigenfunctions are used to calculate dynamic polarizabilities, discrete oscillator strengths, photoionization cross sections and radiative recombination coefficients. A consistent application of the theory imposes a modification on the transition operator, but its effects are small for lithium. The method presented can be regarded as a numerical generalization of the widely used Coulomb approximation.

Caves, T. C.↗

NASA NTRS1992

On development of a finite dynamic element and solution of associated eigenproblem by a block Lanczos procedure

The paper first presents the details of the development of a new six-noded plane triangular finite dynamic element. A block Lanczos algorithm is developed next for the accurate and efficient solution of the quadratic matrix eigenvalue problem associated with the finite dynamic element formulation. The resulting computer program fully exploits matrix sparsity inherent in such a discretization and proves to be most efficient for the extraction of the usually required first few roots and vectors, including repeated ones. Most importantly, the present eigenproblem solution is shown to be comparable to that of the corresponding finite element analysis, thereby rendering the associated dynamic element method rather attractive owing to superior convergence characteristics of such elements, presented herein.

Gupta, K. K.↗

DOE OSTI2025

Stochastic noise can be helpful for variational quantum algorithms

Saddle points constitute a crucial challenge for first-order gradient descent algorithms. In notions of classical machine learning, they are avoided, for example, by means of stochastic gradient descent methods. In this work, we provide evidence that the saddle-points problem can be naturally avoided in variational quantum algorithms by exploiting the presence of stochasticity. We prove convergence guarantees and present practical examples in numerical simulations and on quantum hardware. We argue that the natural stochasticity of variational algorithms can be beneficial for avoiding strict saddle points, i.e., those saddle points with at least one negative Hessian eigenvalue. This insight that some levels of shot noise could help is expected to add a new perspective to notions of near-term variational quantum algorithms. Published by the American Physical Society 2025

Liu, Junyu↗

NASA NTRS1984

Recent advances in methods for numerical solution of O.D.E. initial value problems

In the mathematical modeling of physical systems, it is often necessary to solve an initial value problem (IVP), consisting of a system of ordinary differential equations (ODE). A typical program produces approximate solutions at certain mesh points. Almost all existing codes try to control the local truncation error, while the user is really interested in controlling the true or global error. The present investigation provides a review of recent advances regarding the solution of the IVP, giving particular attention to stiff systems. Stiff phenomena are customarily defined in terms of the eigenvalues of the Jacobian. There are, however, some difficulties connected with this approach. It is pointed out that an estimate of the Lipschitz constant proves to be a very practical way to determine the stiffness of a problem.

Bui, T. D.↗

NASA NTRS1984

A comparison of measured spacecraft modal damping values

The prediction and the measurement of modal damping values of spacecrafts for use in the determination of design loads were summarized. Realistic transient forcing functions were used. The solution of complex eigenvalue response solutions resulting from modal synthesis of subsystems in which each subsystem is assigned modal damping values is discussed. Damping data on both the Voyager and the Galileo spacecrafts using up to eight different types of techniques to analyze modal data are presented. It is concluded that it is difficult to estimate the true damping of a structure and the measured damping values which are dependent on the analysis method of the data.

Wada, B. K.↗

NASA NTRS1991

Flexible body dynamic stability for high performance aircraft

Dynamic equations which include the effects of unsteady aerodynamic forces and a flexible body structure were developed for a free flying high performance fighter aircraft. The linear and angular deformations are assumed to be small in the body reference frame, allowing the equations to be linearized in the deformation variables. Equations for total body dynamics and flexible body dynamics are formulated using the hybrid coordinate method and integrated in a state space format. A detailed finite element model of a generic high performance fighter aircraft is used to generate the mass and stiffness matrices. Unsteady aerodynamics are represented by a rational function approximation of the doublet lattice matrices. The equations simplify for the case of constant angular rate of the body reference frame, allowing the effect of roll rate to be studied by computing the eigenvalues of the system. It is found that the rigid body modes of the aircraft are greatly affected by introducing a constant roll rate, while the effect on the flexible modes is minimal for this configuration.

Goforth, E. A.↗

NASA NTRS1990

Effects of streamwise curvature variations on Goertler vortices

The influence of wall curvature oscillations on Goertler vortices is studied using an asymptotic method. A large number of oscillations leads to a stabilization of the Goertler vortices. In addition, basic velocity profile of a flow over a wall with a large number of oscillations is similar to that over a wall without oscillation. The curvature stays positive while oscillating and the wall is concave. On a wall with curvature decreasing to zero, the Goertler vortices move outside of the boundary layer; no discrete eigenvalue of Goertler instability is found on a wall with zero curvature. In this case the discrete eigenvalues merge into a continuous spectrum.

Jallade, Sophie↗

NASA NTRS2008

Toward Optimal Transport Networks

Strictly evolutionary approaches to improving the air transport system a highly complex network of interacting systems no longer suffice in the face of demand that is projected to double or triple in the near future. Thus evolutionary approaches should be augmented with active design methods. The ability to actively design, optimize and control a system presupposes the existence of predictive modeling and reasonably well-defined functional dependences between the controllable variables of the system and objective and constraint functions for optimization. Following recent advances in the studies of the effects of network topology structure on dynamics, we investigate the performance of dynamic processes on transport networks as a function of the first nontrivial eigenvalue of the network's Laplacian, which, in turn, is a function of the network s connectivity and modularity. The last two characteristics can be controlled and tuned via optimization. We consider design optimization problem formulations. We have developed a flexible simulation of network topology coupled with flows on the network for use as a platform for computational experiments.

Alexandrov, Natalia↗

NASA NTRS2001

The Effect of Stochastic Perturbation of Fuel Distribution on the Criticality of a One Speed Reactor and the Development of Multi-Material Multinomial Line Statistics

The effect of random fuel redistribution on the eigenvalue of a one-speed reactor is investigated. An ensemble of such reactors that are identical to a homogeneous reference critical reactor except for the fissile isotope density distribution is constructed such that it meets a set of well-posed redistribution requirements. The average eigenvalue, , is evaluated when the total fissile loading per ensemble element, or realization, is conserved. The perturbation is proven to increase the reactor criticality on average when it is uniformly distributed. The various causes of the change in reactivity, and their relative effects are identified and ranked. From this, a path towards identifying the causes. and relative effects of reactivity fluctuations for the energy dependent problem is pointed to. The perturbation method of using multinomial distributions for representing the perturbed reactor is developed. This method has some advantages that can be of use in other stochastic problems. Finally, some of the features of this perturbation problem are related to other techniques that have been used for addressing similar problems.

Jahshan, S. N.↗

NASA NTRS1972

A fast numerical method for determining the optimum SNR of an array subject to a Q factor constraint.

The problem of finding the optimum current excitations for a point-detector array which extracts a signal from background noise using the criterion of maximizing the SNR subject to a constraint on the supergain ratio (or Q factor) has recently been solved. The formulation reduced the problem to one of solving for the roots of a polynomial whose coefficients are, in general, complex. We demonstrate that by making a certain transformation, it becomes possible to employ a state variable technique which changes the problem into one of finding the eigenvalues of a real matrix which is considerably faster to do.

Winkler, L. P.↗