Search NASA⌕ Search

SEARCH · Search NASA

Results for “eigenvalue problems”

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 451 records · Page 25

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.↗

The inviscid stability of supersonic flow past heated or cooled axisymmetric bodies

The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.

Shaw, Stephen J.↗

Two-Stage Reduction Of Dynamical Models

No longer necessary to solve eigenvalue problems of high order. Component-mode projection-and-assembly model-reduction (COMPARE) method provides approximation of dynamics of vibrations of complicated, multiple flexible bodies by use of mathematical models of reduced order. Incorporates component-mode synthesis (CMS) method and enhanced projection-and-assembly (EP&A) method, described in "Enhanced Method of Reduction of Dynamical Models" (NPO-18402), providing for somewhat simplified two-stage process in which order of applicable mathematical models reduced. Reduced-order models used to design algorithms of control systems to suppress vibrations or otherwise control structure.

Lee, Allan Y.↗

Effect of surface tension on the onset of convection in a double-diffusive layer

The effect of surface tension on the stability of a double-diffusive layer is considered using linear stability analysis. The surface tension is assumed to vary linearly with temperature and solute concentration. The eigenvalue problem is solved by the Galerkin method. Results show that the predicted stability boundary based on Marangoni effects alone is completely altered in the presence of buoyancy effects induced by low gravity levels (about 10 exp -5 g). At reduced gravity levels, salt-finger instability may onset in the overstable mode due to the stabilizing effect of surface tension. Fluid properties in terms of the Prandtl and the Lewis numbers have a profound effect on the stability conditions; opposite stability characteristics are found in salt solutions and in molten metals.

Chen, C. F.↗

On the shear stabilization of capillary break-up of finite liquid bridges

In this paper we consider an isothermal finite liquid bridge under zero-gravity. A sinusoidal interfacial shear stress drives a nonparallel flow in the bridge. The linear stability of this flow to three-dimensional disturbances (which may deform the cylindrical gas-liquid interface) is determined numerically by solving an elliptic eigenvalue problem. Previous results on shear stabilization of capillary break-up of axially unbounded cylindrical interfaces containing a parallel flow are hereby extended to a nonparallel flow. The turning flow regions influence the stability significantly. However, for the particular cases considered, a small area in parameter space remains where the capillary instability is suppressed through interfacial shear. Second, non-axisymmetric oscillatory instabilities are found which originate from an interaction of the interface deformation and the basic flow. These instabilities may be the isothermal limit of the oscillatory instabilities observed in float-zone crystal growth.

Dijkstra, Henk A.↗

Sensitivity analysis of flutter response of a typical section and a wing in transonic flow

A sensitivity analysis of flutter response of a two-degree of freedom airfoil with plunging and pitching degrees of freedom in transonic flow was performed using a state-space representation of the unsteady aerodynamic behavior. The structural equations of motion of the airfoil with bending and torsional degrees of freedom are coupled to the unsteady airloads, and the aeroelastic system so modeled is solved as an eigenvalue problem to determine the stability. The results of eigenanalysis showed good agreement with flutter calculations performed using a time-integration of the aeroelastic equations. The sensitivities of the flutter speed with respect to the mass and stiffness parameters wee computed by both the analytical and finite difference methods, showing excellent agreement.

Kapania, Rakesh K.↗

Effect of curvature on stationary crossflow instability of a three-dimensional boundary layer

An incompressible three-dimensional laminar boundary-layer flow over a swept wing is used as a model to study both the wall-curvature and streamline-curvature effects on the stationary crossflow instability. The basic state is obtained by solving the full Navier-Stokes (N-S) equations numerically. The linear disturbance equations are cast on a fixed, body-intrinsic, curvilinear coordinate system. Those nonparallel terms which contribute mainly to the streamline-curvature effect are retained in the formulation of the disturbance equations and approximated by their local finite difference values. The resulting eigenvalue problem is solved by a Chebyshev collocation method. The present results indicate that the convex wall curvature has a stabilizing effect, whereas the streamline curvature has a destabilizing effect. A validation of these effects with an N-S solution for the linear disturbance flow is provided.

Lin, Ray-Sing↗

A hydrodynamical model for the explosion of a neutron star just below the minimum mass

The instability of a neutron star at the minimum mass is investigated using a hydrodynamical model to follow the evolution of the unstable star. A detailed analysis of the linear stability of the equilibrium configurations near the minimum mass is performed by solving the radial eigenvalue problem for the fundamental mode. The minimum mass configuration M(mmc) of 0.196 solar mass is found to be stable to small perturbations. Stability to radial perturbations is lost only at a lower critical mass of 0.16 solar mass, corresponding to about 0.8 M(mmc). It is shown that a simple dynamical model constructed using a 3-polytrope equation of state for hot dense matter reproduces the key dynamical features of the instability in the explosion phase.

Colpi, Monica↗

Variational modal identification of gyroscopic distributed-parameter systems

Variational modal identification of gyroscopic distribution parameter systems is discussed. Topics covered include: equation of motion for gyroscopic distributed-parameter systems; a variational formulation for the eigenvalue problem; a variational formulation for the modal identification of gyroscopic distributed-parameter systems; discretization; and numerical example.

Silverberg, Larry M.↗

On the three-dimensional instability of laminar boundary layers on concave walls

A study is made of the stability of laminar boundary-layer profiles on slightly curved walls relative to small disturbances that result from vortices whose axes are parallel to the principal direction of flow. The result is an eigenvalue problem by which, for a given undisturbed flow at a prescribed wall, the amplification or decay is computed for each Reynolds number and each vortex thickness. For neutral disturbances (zero amplification) a critical Reynolds number is determined for each vortex distribution. The numerical calculation produces amplified disturbances on concave walls only.

FLOW, LAMINAR↗

The effect of initial imperfections on shell stability: An updated review

The development of the Delft Interactive Shell Design Code (DISDECO) is described. The purpose is to make the accumulated theoretical, numerical, and practical knowledge of the last years readily accessible to users interested in the analysis of buckling sensitive structures. With this open ended, hierarchical, interactive computer code, the user can access workstation programs of increasing complexity. As a first update of DISDECO, the computational modules of the initial level were modified so as to handle combined loads consisting of axial compression, internal and external pressure, and clockwise or counter clockwise torque. The applied load is assumed to have a uniform spatial distribution and is divided into a fixed and a variable part. The magnitude of the variable part is allowed to vary in proportion to a load parameter Lambda. This leads to an eigenvalue problem for the critical lambda sub C.

Arbocz, J.↗

Comparison of level-1 and level-2 buckling and postbuckling solutions

The accuracy of the solutions obtained with the Level-1 computational modules of DISDECO (Delft Interactive Shell Design Code) are investigated by comparing them with solutions obtained with ANILISA, a Level-2 computational module. Whereas the Level-1 solutions are based on the use of trigonometric functions to represent the out-of-plane displacement W and involve Galerkin type approximations to reduce the governing partial differential equations to simple algebraic expressions, the Level-2 solutions employ a truncated Fourier expansion in the circumferential direction and solve the resulting ordinary differential equations forming a nonlinear eigenvalue problem by a Stodola like technique. This involves the numerical integration of the set of ordinary differential equations, which in turn makes it possible to satisfy the specified boundary conditions exactly. It is shown that in many applications the Level-1 predictions concerning the critical buckling load of the perfect structure and its imperfection sensitivity are quite accurate and reliable. For those cases where this is not true, it is shown how the pre and postprocessing facilities of DISDECO can be used efficiently to arrive at a reliable prediction. This involves always the use of the Level-2 computational module ANILISA with its rigorous prebuckling solution and its feature which makes it possible to satisfy the specified boundary conditions exactly.

Arbocz, J.↗

Inviscid instability of streamwise corner flow

Linear stability of the incompressible flow along a stream wise corner is studied by solving the two-dimensional eigenvalue problem governed by partial differential equations. It is found that this fully three-dimensional flow is subject to inviscid instability due to the inflectional nature of the stream wise velocity profile. The higher growth rates for the inviscid instability mode, which is symmetric about the corner bisector, as compared to the viscous Tollmien-Schlichting instability operative away from the corner is consistent with the experimental findings that the corner flow transitions to turbulence earlier than the two-dimensional Blasius flow away from the corner.

Balachandar, S.↗

Automatic A-set selection for dynamics analysis

A method for selecting optimum NASTRAN analysis set degrees of freedom for the dynamic eigenvalue problem is described. Theoretical development of the Guyan reduction procedure on which the method is based is first summarized. The algorithm used to select the analysis set degrees of freedom is then developed. Two example problems are provided to demonstrate the accuracy of the algorithm.

Allen, Tom↗

Separation analysis, a tool for analyzing multigrid algorithms

The separation of vectors by multigrid (MG) algorithms is applied to the study of convergence and to the prediction of the performance of MG algorithms. The separation operator for a two level cycle algorithm is derived. It is used to analyze the efficiency of the cycle when mixing of eigenvectors occurs. In particular cases the separation analysis reduces to Fourier type analysis. The separation operator of a two level cycle for a Schridubger eigenvalue problem, is derived and analyzed in a Fourier basis. Separation analysis gives information on how to choose performance relaxations and inter-level transfers. Separation analysis is a tool for analyzing and designing algorithms, and for optimizing their performance.

Costiner, Sorin↗

Buckling of imperfect, anisotropic, ring-stiffened cylinders under combined loads

The objective of this study is to develop an anlysis to predict buckling loads of ring-stiffened anisotropic cylinders subject to axial compression, torsion, and internal pressure. This structure is modeled as a branched shell. A nonlinear axisymmetric prebuckling equilibrium state is assumed which is amenable to an exact solution within each branch. Axisymmetric geometric imperfections are included. Buckling displacements are represented by a Fourier series in the circumferential coordination and the finite-element method in the axial coordinate. Application of the Trefftz criterion to the second variation of the total potential energy leads to a nonlinear eigenvalue problem for the buckling load and mode. Results are presented for both unstiffened and ring-stiffened cylinders in the form of buckling interaction diagrams. Imperfections can cause an unexpected buckling mode in the ring web which would not occur for the perfect structure, and pressurization diminishes the benefit of adding rings to the unstiffened shell to increase the buckling load. The implementation of the analysis methodology into a structural sizing algorithm is discussed.

Ley, Robert P.↗

Sensitivity Analysis of Wing Aeroelastic Responses

Design for prevention of aeroelastic instability (that is, the critical speeds leading to aeroelastic instability lie outside the operating range) is an integral part of the wing design process. Availability of the sensitivity derivatives of the various critical speeds with respect to shape parameters of the wing could be very useful to a designer in the initial design phase, when several design changes are made and the shape of the final configuration is not yet frozen. These derivatives are also indispensable for a gradient-based optimization with aeroelastic constraints. In this study, flutter characteristic of a typical section in subsonic compressible flow is examined using a state-space unsteady aerodynamic representation. The sensitivity of the flutter speed of the typical section with respect to its mass and stiffness parameters, namely, mass ratio, static unbalance, radius of gyration, bending frequency, and torsional frequency is calculated analytically. A strip theory formulation is newly developed to represent the unsteady aerodynamic forces on a wing. This is coupled with an equivalent plate structural model and solved as an eigenvalue problem to determine the critical speed of the wing. Flutter analysis of the wing is also carried out using a lifting-surface subsonic kernel function aerodynamic theory (FAST) and an equivalent plate structural model. Finite element modeling of the wing is done using NASTRAN so that wing structures made of spars and ribs and top and bottom wing skins could be analyzed. The free vibration modes of the wing obtained from NASTRAN are input into FAST to compute the flutter speed. An equivalent plate model which incorporates first-order shear deformation theory is then examined so it can be used to model thick wings, where shear deformations are important. The sensitivity of natural frequencies to changes in shape parameters is obtained using ADIFOR. A simple optimization effort is made towards obtaining a minimum weight design of the wing, subject to flutter constraints, lift requirement constraints for level flight and side constraints on the planform parameters of the wing using the IMSL subroutine NCONG, which uses successive quadratic programming.

Issac, Jason Cherian↗

Robust Passification via Optimal Sensor Blending and Control Allocation

Robust passification is considered for uncertain linear, time invariant, systems having redundant actuators and sensors. The approach is to obtain optimal sensor blending and control allocation matrices that maximize the region in the parameter space in which the system remains passive. The approach results in a generalized eigenvalue problem consisting of a number of linear matrix inequalities (LMIs). Reduction of the number of LMIs is investigated, and a numerical example is given for demonstrating the approach.

Joshi, S. M.↗