SEARCH · Engineering Papers

Results for “eigenvalues”

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 1,711 records · Page 95

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Flutter analysis using transversality theory

A new method of calculating flutter boundaries of undamped aeronautical structures is presented. The method is an application of the weak transversality theorem used in catastrophe theory. In the first instance, the flutter problem is cast in matrix form using a frequency domain method, leading to an eigenvalue matrix. The characteristic polynomial resulting from this matrix usually has a smooth dependence on the system's parameters. As these parameters change with operating conditions, certain critical values are reached at which flutter sets in. Our approach is to use the transversality theorem in locating such flutter boundaries using this criterion: at a flutter boundary, the characteristic polynomial does not intersect the axis of the abscissa transversally. Formulas for computing the flutter boundaries and flutter frequencies of structures with two degrees of freedom are presented, and extension to multi-degree of freedom systems is indicated. The formulas have obvious applications in, for instance, problems of panel flutter at supersonic Mach numbers.

Afolabi, D.↗

NASA NTRS1993

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↗

NASA NTRS1993

Application of Cloude's target decomposition theorem to polarimetric imaging radar data

We apply Cloude's decomposition to imaging radar polarimetry. We derive the general expressions for the eigenvalues and eigenvectors for the case of terrain with reflection symmetry, and show in detail how the decomposition results can guide the interpretation of scattering from vegetated areas. For multi-frequency polarimetric radar measurements of a clear-cut area, the decomposition leads us to conclude that the vegetation is probably thin compared to even the C-band radar wavelength of 6 cm. For a forested area, we notice an increased amount of even number of reflection scattering at P-band and L-band, probably the result of penetration through the coniferous canopy resulting in trunk-ground double reflection scattering. The scattering for the forested area is still dominated by scattering from randomly oriented cylinders, however. It is found that these cylinders are thicker than in the case of clear-cut areas, leading us to conclude that scattering from the branches probably dominate in this case.

Vanzyl, Jakob J.↗

NASA NTRS1993

The multigrid preconditioned conjugate gradient method

A multigrid preconditioned conjugate gradient method (MGCG method), which uses the multigrid method as a preconditioner of the PCG method, is proposed. The multigrid method has inherent high parallelism and improves convergence of long wavelength components, which is important in iterative methods. By using this method as a preconditioner of the PCG method, an efficient method with high parallelism and fast convergence is obtained. First, it is considered a necessary condition of the multigrid preconditioner in order to satisfy requirements of a preconditioner of the PCG method. Next numerical experiments show a behavior of the MGCG method and that the MGCG method is superior to both the ICCG method and the multigrid method in point of fast convergence and high parallelism. This fast convergence is understood in terms of the eigenvalue analysis of the preconditioned matrix. From this observation of the multigrid preconditioner, it is realized that the MGCG method converges in very few iterations and the multigrid preconditioner is a desirable preconditioner of the conjugate gradient method.

Tatebe, Osamu↗

NASA NTRS1994

Structural design criteria for filament-wound composite shells

Advanced composite cylinders, manufactured by filament winding, provide a cost effective solution to many present structural applications; however, the compressive performance of filament-wound cylinders is lower than comparable shells fabricated from unidirectional tape. The objective of this study was to determine the cause of this reduction in thin filament-wound cylinders by relating the manufacturing procedures to the quality of the cylinder and to its compressive performance. The experiments on cylinder buckling were complemented by eigenvalue buckling analysis using a detailed geometric model in a finite element analysis. The applicability of classical buckling analyses was also investigated as a design tool.

Hahn, H. T.↗

NASA NTRS1991

Computation of vibration mode elastic-rigid and effective weight coefficients from finite-element computer program output

Post-processing algorithms are given to compute the vibratory elastic-rigid coupling matrices and the modal contributions to the rigid-body mass matrices and to the effective modal inertias and masses. Recomputation of the elastic-rigid coupling matrices for a change in origin is also described. A computational example is included. The algorithms can all be executed by using standard finite-element program eigenvalue analysis output with no changes to existing code or source programs.

Levy, R.↗

NASA NTRS1994

Parabolized stability equations

The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Herbert, Thorwald↗

NASA NTRS1994

Dynamical observer for a flexible beam via finite element approximations

The purpose of this view-graph presentation is a computational investigation of the closed-loop output feedback control of a Euler-Bernoulli beam based on finite element approximation. The observer is part of the classical observer plus state feedback control, but it is finite-dimensional. In the theoretical work on the subject it is assumed (and sometimes proved) that increasing the number of finite elements will improve accuracy of the control. In applications, this may be difficult to achieve because of numerical problems. The main difficulty in computing the observer and simulating its work is the presence of high frequency eigenvalues in the finite-element model and poor numerical conditioning of some of the system matrices (e.g. poor observability properties) when the dimension of the approximating system increases. This work dealt with some of these difficulties.

Manitius, Andre↗

NASA NTRS1991

Alleviation of whirl-flutter on a joined-wing tilt-rotor aircraft configuration using active controls

The feasibility of using active controls to delay the onset of whirl-flutter on a joined-wing tilt rotor aircraft was investigated. The CAMRAD/JA code was used to obtain a set of linear differential equations which describe the motion of the joined-wing tilt-rotor aircraft. The hub motions due to wing/body motion is a standard input to CAMRAD/JA and were obtained from a structural dynamics model of a representative joined-wing tilt-rotor aircraft. The CAMRAD/JA output, consisting of the open-loop system matrices, and the airframe free vibration motion were input to a separate program which performed the closed-loop, active control calculations. An eigenvalue analysis was performed to determine the flutter stability of both open- and closed-loop systems. Sensor models, based upon the feedback of pure state variables and based upon hub-mounted sensors, providing physically measurable accelerations, were evaluated. It was shown that the onset of tilt-rotor whirl-flutter could be delayed from 240 to above 270 knots by feeding back vertical and span-wise accelerations, measured at the rotor hub, to the longitudinal cyclic pitch. Time response calculations at a 270-knot cruise condition showed an active cyclic pitch control level of 0.009 deg, which equates to a very acceptable 9 pound active-control force applied at the rotor hub.

Vanaken, Johannes M.↗

NASA NTRS1994

Eigensystem realization algorithm user's guide forVAX/VMS computers: Version 931216

The eigensystem realization algorithm (ERA) is a multiple-input, multiple-output, time domain technique for structural modal identification and minimum-order system realization. Modal identification is the process of calculating structural eigenvalues and eigenvectors (natural vibration frequencies, damping, mode shapes, and modal masses) from experimental data. System realization is the process of constructing state-space dynamic models for modern control design. This user's guide documents VAX/VMS-based FORTRAN software developed by the author since 1984 in conjunction with many applications. It consists of a main ERA program and 66 pre- and post-processors. The software provides complete modal identification capabilities and most system realization capabilities.

Pappa, Richard S.↗

NASA NTRS1994

Gain weighted eigenspace assignment

This report presents the development of the gain weighted eigenspace assignment methodology. This provides a designer with a systematic methodology for trading off eigenvector placement versus gain magnitudes, while still maintaining desired closed-loop eigenvalue locations. This is accomplished by forming a cost function composed of a scalar measure of error between desired and achievable eigenvectors and a scalar measure of gain magnitude, determining analytical expressions for the gradients, and solving for the optimal solution by numerical iteration. For this development the scalar measure of gain magnitude is chosen to be a weighted sum of the squares of all the individual elements of the feedback gain matrix. An example is presented to demonstrate the method. In this example, solutions yielding achievable eigenvectors close to the desired eigenvectors are obtained with significant reductions in gain magnitude compared to a solution obtained using a previously developed eigenspace (eigenstructure) assignment method.

Davidson, John B.↗

NASA NTRS1994

Prediction of Film Cooling on Gas Turbine Airfoils

A three-dimensional Navier-Stokes analysis tool has been developed in order to study the effect of film cooling on the flow and heat transfer characteristics of actual turbine airfoils. An existing code (Arnone et al., 1991) has been modified for the purpose. The code is an explicit, multigrid, cell-centered, finite volume code with an algebraic turbulence model. Eigenvalue scaled artificial dissipation and variable-coefficient implicit residual smoothing are used with a full-multigrid technique. Moreover, Mayle's transition criterion (Mayle, 1991) is used. The effects of film cooling have been incorporated into the code in the form of appropriate boundary conditions at the hole locations on the airfoil surface. Each hole exit is represented by several control volumes, thus providing an ability to study the effect of hole shape on the film-cooling characteristics. Comparison is fair with near mid-span experimental data for four and nine rows of cooling holes, five on the shower head, and two rows each on the pressure and suction surfaces. The computations, however, show a strong spanwise variation of the heat transfer coefficient on the airfoil surface, specially with shower-head cooling.

Garg, Vijay K.↗

NASA NTRS1994

Coupled Marangoni-Benard/Rayleigh-Benard Instability with Temperature Dependent Viscosity

The onset of convection induced by coupled surface tension gradient and buoyancy forces is investigated with temperature dependent viscosity. Both surface tension and viscosity are assumed to vary linearly with temperature. The limiting case, Ma=0, is the buoyancy driven convection problem typically referred to as the Rayleigh-Benard stability problem. The other limiting case, Ra=0, is the surface tension (gradient) driven flow problem referred to as the Marangoni-Benard problem. The equations and boundary conditions obtained from the linear analysis are solved numerically as a generalized eigenvalue problem. Neutral stability curves for different viscosity slopes have been generated for the Marangoni-Benard and Rayleigh-Benard problems. It is shown that the curves can be collapsed to a single curve by appropriately scaling the results for each of the limiting cases. The critical Marangoni number is determined as a function of the slope of the viscosity temperature variation, epsilon, for different values of the Rayleigh number. When the viscosity decreases linearly with temperature, the coupled buoyancy-surface tension problem, including the limiting cases of Ra=0 and Ma=0, is found to be more stable than the constant viscosity case.

Skarda, J. Raymond Lee↗

NASA NTRS1993

Stability of the Halley cometosheath with resistivity and plasma motion

The MHD stability of the cometary inner sheath determined by the balance between the inward Lorentz body force and the outward ion-neutral drag force is investigated by numerically solving the wave equations which include resistivity, plasma motion and plasma pressure with the help of two-point boundary value method. The eigenvalues and the eigenfunctions are obtained numerically by treating the cometary inner sheath as a layer of finite thickness, bounded by the contact surface, that is, the diamagnetic cavity boundary. To gain insight into the problem, certain limiting cases of the wave equations are also discussed. The diamagnetic cavity boundary and the adjacent layer of about 100-km thickness of Comet Halley is found to be unstable. The effects of finite plasma pressure, dissociative recombination, mass loading due to photoionization, resistivity, and plasma motion are found to be stabilizing but are unable to quench the instability completely. Motion of the Halley ionopause has been confirmed by observations.

Srivastava, Krishna M.↗

NASA NTRS1992

On a finite-state inflow application to flap-lag-torsion damping in hover

An aerodynamic model with a coupled set of generalized dynamic wake equations and hybrid equations of motion for an elastic blade are applied here to a two-blade untwisted stiff in-plane hingeless small-scale model rotor with torsionally soft blades. Blade root offset, precone, blade droop, pitch control stiffness, and blade pitch angle are included in the model rotor. Numerical results show that 3D tip relief effects within the nonuniform steady-state inflow are significant to predict steady-state aerodynamic loads and blade deflections. Eigenvalue results confirm the importance of unsteady 3D aerodynamics in predicting lead-lag damping and frequency. Eigenvector analysis correlations reinforced qualitative and quantitative shortcomings associated with quasi-steady 2D aerodynamic theory for aeroelastic applications in hover.

De Andrade, Donizeti↗

NASA NTRS1994

Quantum statistics of Raman scattering model with Stokes mode generation

The model describing three coupled quantum oscillators with decay of Rayleigh mode into the Stokes and vibration (phonon) modes is examined. Due to the Manley-Rowe relations the problem of exact eigenvalues and eigenstates is reduced to the calculation of new orthogonal polynomials defined both by the difference and differential equations. The quantum statistical properties are examined in the case when initially: the Stokes mode is in the vacuum state; the Rayleigh mode is in the number state; and the vibration mode is in the number of or squeezed states. The collapses and revivals are obtained for different initial conditions as well as the change in time the sub-Poisson distribution by the super-Poisson distribution and vice versa.

Tanatar, Bilal↗

NASA NTRS1994

On Bi-Grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations

A procedure is presented for utilizing a bi-grid stability analysis as a practical tool for predicting multigrid performance in a range of numerical methods for solving Euler and Navier-Stokes equations. Model problems based on the convection, diffusion and Burger's equation are used to illustrate the superiority of the bi-grid analysis as a predictive tool for multigrid performance in comparison to the smoothing factor derived from conventional von Neumann analysis. For the Euler equations, bi-grid analysis is presented for three upwind difference based factorizations, namely Spatial, Eigenvalue and Combination splits, and two central difference based factorizations, namely LU and ADI methods. In the former, both the Steger-Warming and van Leer flux-vector splitting methods are considered. For the Navier-Stokes equations, only the Beam-Warming (ADI) central difference scheme is considered. In each case, estimates of multigrid convergence rates from the bi-grid analysis are compared to smoothing factors obtained from single-grid stability analysis. Effects of grid aspect ratio and flow skewness are examined. Both predictions are compared with practical multigrid convergence rates for 2-D Euler and Navier-Stokes solutions based on the Beam-Warming central scheme.

Ibraheem, S. O.↗

NASA NTRS1995

Covariant deformed oscillator algebras

The general form and associativity conditions of deformed oscillator algebras are reviewed. It is shown how the latter can be fulfilled in terms of a solution of the Yang-Baxter equation when this solution has three distinct eigenvalues and satisfies a Birman-Wenzl-Murakami condition. As an example, an SU(sub q)(n) x SU(sub q)(m)-covariant q-bosonic algebra is discussed in some detail.

Quesne, Christiane↗