Search NASA⌕ Search

SEARCH · Search NASA

Results for “Eigensolution”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

On computing eigensolution sensitivity data using free vibration solutions

A simplified method of computing eigensolution sensitivity derivatives in structural dynamics is developed. It is shown that if the elements of stiffness and mass matrices associated with a design variable are homogeneous functions of that design variable, then eigenvalue derivatives can be computed from element strain and kinetic energies. Furthermore, if cross-mode energies are known, eigensolution derivatives of modified systems can be computed approximately using assume mode reanalysis formulation. A ten bar truss example is used to illustrate the present formulations.

Wang, B. P.↗
NASA NTRS1990

Efficient eigensolution reanalysis of nonclassically damped structures

Effective methods of approximate eigensolution reanalysis of modified nonclassically damped structures are developed. For structures with passive or active discrete damping devices or with damping treatment, the system becomes nonproportionally damped and the computation of its dynamic responses may require the use of complex modes. For larger systems, the computation of complex modes is very expensive. Thus it is desirable to have approximate reanalysis techniques for the efficient evaluation of the effect of design changes. In recent years, the assumed mode reanalysis method was successfully applied to minimum weight design of undamped structures with natural frequency constraints. The accuracy of the assumed mode reanalysis method can be improved dramatically if the global approximation function includes the normal modes of the original system and their derivatives. This approach was demonstrated to be effective even for a system with shape changes. The approach used by Noor et al. for eigensolution reanalysis of undamped structures is extended to treat a nonclassically damped system.

Wang, Bo Ping↗
NASA NTRS1980

Note on the eigensolution of a homogeneous equation with semi-infinite domain

The 'variation-iteration' method using Green's functions to find the eigenvalues and the corresponding eigenfunctions of a homogeneous Fredholm integral equation is employed for the stability analysis of fluid hydromechanics problems with a semiinfinite (infinite) domain of application. The objective of the study is to develop a suitable numerical approach to the solution of such equations in order to better understand the full set of equations for 'real-world' flow models. The study involves a search for a suitable value of the length of the domain which is a fair finite approximation to infinity, which makes the eigensolution an approximation dependent on the length of the interval chosen. In the examples investigated y = 1 = a seems to be the best approximation of infinity; for y greater than unity this method fails due to the polynomial nature of Green's functions.

Wadia, A. R.↗
NASA NTRS1994

Stiffness-generated rigid-body mode shapes for Lanczos eigensolution with SUPORT DOF by way of a MSC/NASTRAN DMAP alter

When using all MSC/NASTRAN eigensolution methods except Lanczos, the analyst can replace the coupled system rigid-body modes calculated within DMAP module READ with mass orthogonalized and normalized rigid-body modes generated from the system stiffness. This option is invoked by defining MSC/NASTRAN r-set degrees of freedom via the SUPORT bulk data card. The newly calculated modes are required if the rigid-body modes calculated by the eigensolver are not 'clean' due to numerical roundoffs in the solution. When performing transient structural dynamic load analysis, the numerical roundoffs can result in inaccurate rigid-body accelerations which affect steady-state responses. Unfortunately, when using the Lanczos method and defining r-set degrees of freedom, the rigid-body modes calculated within DMAP module REIGL are retained. To overcome this limitation and to allow MSC/NASTRAN to handle SUPORT degrees of freedom identically for all eigensolvers, a DMAP Alter has been written which replaces Lanczos-calculated rigid-body modes with stiffness-generated rigid-body modes. The newly generated rigid-body modes are normalized with respect to the system mass and orthogonalized using the Gram-Schmidt technique. This algorithm has been implemented as an enhancement to an existing coupled loads methodology.

Abdallah, Ayman A.↗
NASA NTRS1981

On eigensolutions for discontinuous liners in a duct containing uniform mean flow

Sound attenuation in a rectangular acoustically lined duct containing uniform mean flow is analytically investigated using the generalized Wiener-Hopf technique. Uniqueness of the solution is enforced for lined sections of the finite axial extent by imposing edge conditions at the liner interface. Possible edge conditions are considered, including the Kutta condition, and the causal solution corresponding to edge conditions is considered the best choice. Solution methods such as the mode matching and singularity methods imply differing edge conditions, and results show that power attenuation is insensitive to the imposed edge conditions, although significant differences are observed for the reflection coefficient. The amplitude of the exponentially increasing instability mode in the lined section must be set to zero as a first approximation to the nonlinear situation, and results indicate that measurements of the reflection factor can be used to make a more definite decision about physically appropriate edge conditions.

Koch, W.↗
NASA NTRS1989

Eigensolution of finite element problems in a completely connected parallel architecture

A parallel algorithm is presented for the solution of the generalized eigenproblem in linear elastic finite element analysis. The algorithm is based on a completely connected parallel architecture in which each processor is allowed to communicate with all other processors. The algorithm is successfully implemented on a tightly coupled MIMD parallel processor. A finite element model is divided into m domains each of which is assumed to process n elements. Each domain is then assigned to a processor or to a logical processor (task) if the number of domains exceeds the number of physical processors. The effect of the number of domains, the number of degrees-of-freedom located along the global fronts, and the dimension of the subspace on the performance of the algorithm is investigated. For a 64-element rectangular plate, speed-ups of 1.86, 3.13, 3.18, and 3.61 are achieved on two, four, six, and eight processors, respectively.

Akl, F.↗
NASA NTRS1989

Eigensolution of finite element problems in a completely connected parallel architecture

A parallel algorithm for the solution of the generalized eigenproblem in linear elastic finite element analysis, (K)(phi)=(M)(phi)(omega), where (K) and (M) are of order N, and (omega) is of order q is presented. The parallel algorithm is based on a completely connected parallel architecture in which each processor is allowed to communicate with all other processors. The algorithm has been successfully implemented on a tightly coupled multiple-instruction-multiple-data (MIMD) parallel processing computer, Cray X-MP. A finite element model is divided into m domains each of which is assumed to process n elements. Each domain is then assigned to a processor, or to a logical processor (task) if the number of domains exceeds the number of physical processors. The macro-tasking library routines are used in mapping each domain to a user task. Computational speed-up and efficiency are used to determine the effectiveness of the algorithm. The effect of the number of domains, the number of degrees-of-freedom located along the global fronts and the dimension of the subspace on the performance of the algorithm are investigated. For a 64-element rectangular plate, speed-ups of 1.86, 3.13, 3.18 and 3.61 are achieved on two, four, six and eight processors, respectively.

Akl, Fred A.↗
NASA NTRS1991

Lanczos eigensolution method for high-performance computers

The theory, computational analysis, and applications are presented of a Lanczos algorithm on high performance computers. The computationally intensive steps of the algorithm are identified as: the matrix factorization, the forward/backward equation solution, and the matrix vector multiples. These computational steps are optimized to exploit the vector and parallel capabilities of high performance computers. The savings in computational time from applying optimization techniques such as: variable band and sparse data storage and access, loop unrolling, use of local memory, and compiler directives are presented. Two large scale structural analysis applications are described: the buckling of a composite blade stiffened panel with a cutout, and the vibration analysis of a high speed civil transport. The sequential computational time for the panel problem executed on a CONVEX computer of 181.6 seconds was decreased to 14.1 seconds with the optimized vector algorithm. The best computational time of 23 seconds for the transport problem with 17,000 degs of freedom was on the the Cray-YMP using an average of 3.63 processors.

Bostic, Susan W.↗
NASA NTRS1991

The unsteady laminar boundary layer on an axisymmetric body subject to small-amplitude fluctuations in the free-stream velocity

The effect of small-amplitude, time-periodic, free-stream disturbances on an otherwise steady axisymmetric boundary layer on a circular cylinder is considered. Numerical solutions to the problem are presented, and an asymptotic solution to the flow, valid far downstream along the axis of the cylinder is detailed. Particular emphasis is placed on the unsteady eigensolutions that occur far downstream, which turn out to be very different from the analogous planar eigensolutions. These axisymmetric eigensolutions are computed numerically and also are described by asymptotic analyses valid for low and high frequencies of oscillation.

Duck, Peter W.↗
NASA NTRS1988

Modal parameter identification in space structures

The objective of this project is to identify modal properties such as the eigenvalues and eigenfunctions of structures. The formal means for accomplishing this task, Structural Identification, is viewed as a two step procedure: (1) identify the eigensolution; and (2) using the identified eigensolution, identify the mass and stiffness. The eigensolution is identified as a correction on a postulated model based on erroneous parameters.

Baruh, H.↗
NASA NTRS1973

Fast modal extraction in NASTRAN via the FEER computer program

A new eigensolution routine, FEER (Fast Eigensolution Extraction Routine), used in conjunction with NASTRAN at Israel Aircraft Industries is described. The FEER program is based on an automatic matrix reduction scheme whereby the lower modes of structures with many degrees of freedom can be accurately extracted from a tridiagonal eigenvalue problem whose size is of the same order of magnitude as the number of required modes. The process is effected without arbitrary lumping of masses at selected node points or selection of nodes to be retained in the analysis set. The results of computational efficiency studies are presented, showing major arithmetic operation counts and actual computer run times of FEER as compared to other methods of eigenvalue extraction, including those available in the NASTRAN READ module. It is concluded that the tridiagonal reduction method used in FEER would serve as a valuable addition to NASTRAN for highly increased efficiency in obtaining structural vibration modes.

Newman, M. B.↗
NASA NTRS1977

On the theory of solitary Rossby waves

The evolution of long, finite amplitude Rossby waves in a horizontally sheared zonal current is studied. The wave evolution is described by the Korteweg-de Vries equation or the modified Korteweg-de Vries equation depending on the atmospheric stratification. In either case, the cross-stream modal structure of these waves is given by the long-wave limit of the neutral eigensolutions of the barotropic stability equation. Both non-singular and singular eigensolutions are considered and the appropriate analysis is developed to yield a uniformly valid description of the motion in the critical-layer region where the wave speed matches the flow velocity. The analysis demonstrates that coherent, propagating, eddy structures can exist in stable shear flows and that these eddies have peculiar interaction properties quite distinct from the traditional views of turbulent motion.

Redekopp, L. G.↗
NASA NTRS1981

A study of the viscous dissipation and surface loading on a vibrating surface

The energy dissipated by viscosity at the edge of a vibrating flat plate is calculated and compared to the radiated acoustic energy. A correction to the Kirchhoff integral estimate of the noise is derived. For Helmholtz number of order unity and smaller the dissipation can be comparable to or greater than the acoustic energy. A viscous compressible theory of the load distribution on a vibrating two dimensional body is developed. First it is shown that load calculations based on potential theory and the Newmann uniqueness condition (continuity of potential or pressure on the surface) are not in agreement with experiment or the more correct viscous theory. For a flat plate airfoil the eigensolution of potential theory is indeterminant while viscous theory yields a unique solution that has square root singularities at the edges. It is also shown that for compact surfaces the far field acoustics depend only on the magnitude of the eigensolutions of potential theoy and so will be uniquely determined by the viscous theory. It is suggested that the general viscous theory of vibrating surfaces with cross sectional geometry will lead to results in agreement with expermentally measured load distributions.

Yates, J. E.↗
NASA NTRS1983

The evolution of Tollmien-Schlichting waves near a leading edge. II - Numerical determination of amplitudes

In the first part of this investigation, Goldstein (1983) has shown that the amplitude of the spatially growing Tollmien-Schlichting wave generated by a time-harmonic free-stream disturbance is related to the coefficient multiplying the lowest-order asymptotic eigensolution of the unsteady boundary-layer equation. In the present study, a numerical solution of the unsteady boundary-layer equation is used to relate the amplitude of the asymptotic eigensolution, and consequently of the Tollmien-Schlichting wave, to that of the imposed free-stream disturbance for the special case of a uniformly pulsating stream. It is pointed out that the ideas of this study can be extended to other, more complex bodies and free-stream oscillations.

Goldstein, M. E.↗
NASA NTRS1985

Numerical solution of the vertical structure equation in the normal mode method

In the present model of multilayered stability stratification, aimed at obtaining the analytic eigensolutions of the vertical structure equation, each layer is characterized by its own static stability value. By requiring continuity of pressure and vertical velocity across each interface level, and by imposing suitable upper and lower boundary conditions, matching eigensolutions are obtained in terms of the Bessel functions. Attention is given to an explicit example of a double-layered stratified atmosphere which demonstrates the mathematical manipulations involved; the resultant vertical structure functions are used to check the accuracy of the numerical solutions by the finite difference and finite element methods.

Sasaki, Y. K.↗
NASA NTRS1985

Parameter identification in distributed systems

This paper describes a method for the identification of the parameters entering into the equations of motion of distributed systems. Because the motion of distributed systems is described in terms of partial differential equations, these parameters are in general continuous functions of the spatial variables. For vibrating systems, these parameters ordinarily represent the mass, stiffness and damping distributions. In this paper, these distributions are expanded in terms of finite series of known functions of the spatial variables multiplied by undetermined coefficients. It is assumed that the nature of the equations of motion is known and that a limited number of eigenvalues and eigenfunctions is identified in advance. Use is then made of the least squares method, in conjunction with the eigenfunctions' orthogonality, to compute the undetermined coefficients, thus identifying the system distributed parameters. A method for the identification of the eigensolution is also presented. The procedure for the identification of the eigensolution and of the system parameters is demonstrated via a numerical example.

Baruh, H.↗
NASA NTRS1987

Generation of Tollmien-Schlichting waves on interactive marginally separated flows

This paper is concerned with the interaction of very-long-wavelength free-stream disturbances with the small but abrupt changes in the mean flow that occur near the minimum-skin-friction point in an interactive marginally separated boundary layer. The source frequency is chosen so that the eigensolutions with that frequency have an 'interactive' structure in the region of marginal separation. The eigensolution wavelength scale must then differ from the lengthscale of the marginal separation, and a composite expansion technique has to be used to obtain the solution. The initial instability wave amplitude turns out to be exponentially small, but eventually dominates the original disturbance owing to its exponential growth. It then begins to decay but ultimately turns into a standard spatially growing Tollmien-Schlichting wave much further downstream.

Goldstein, M. E.↗