SEARCH · Engineering Papers
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.
Amplification matrices, numerical kernels, stable, and exponentially stable numerical solutions are examined. The various techniques involved in these concepts are applied to certain systems that have Jordan forms, which are nondiagonal, with particular interest in the case of imaginary or zero eigenvalues.
Analysis and techniques are developed for the multivariable root locus and the multivariable optimal root locus. The generalized eigenvalue problem is used to compute angles and sensitivities for both types of loci, and an algorithm is presented that determines the asymptotic properties of the optimal root locus.
Analytical and empirical studies of a finite difference method for the solution of the transonic flow about harmonically oscillating wings and airfoils are presented. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady equations for small disturbances. The steady velocity potential is obtained first from the well-known nonlinear equation for steady transonic flow. The unsteady velocity potential is then obtained from a linear differential equation in complex form with spatially varying coefficients. Since sinusoidal motion is assumed, the unsteady equation is independent of time. An out-of-core direct solution procedure was developed and applied to two-dimensional sections. Results are presented for a section of vanishing thickness in subsonic flow and an NACA 64A006 airfoil in supersonic flow. Good correlation is obtained in the first case at values of Mach number and reduced frequency of direct interest in flutter analyses. Reasonable results are obtained in the second case. Comparisons of two-dimensional finite difference solutions with exact analytic solutions indicate that the accuracy of the difference solution is dependent on the boundary conditions used on the outer boundaries. Homogeneous boundary conditions on the mesh edges that yield complex eigenvalues give the most accurate finite difference solutions. The plane outgoing wave boundary conditions meet these requirements.
The mathematical theory for decoupling mth-order matrix differential equations is presented. It is shown that the decoupling precedure can be developed from the algebraic theory of matrix polynomials. The role of eigenprojectors and latent projectors in the decoupling process is discussed and the mathematical relationships between eigenvalues, eigenvectors, latent roots, and latent vectors are developed. It is shown that the eigenvectors of the companion form of a matrix contains the latent vectors as a subset. The spectral decomposition of a matrix and the application to differential equations is given.
Special modeling techniques were developed to simulate an accurate mathematical model of the transporter/canister/payload system during ground transport of the Inertial Upper Stage/Tracking and Data Relay Satellite (IUS/TDRS) payload. The three finite element models - the transporter, the canister, and the IUS/TDRS payload - were merged into one model and used along with the NASTRAN normal mode analysis. Deficiencies were found in the NASTRAN program that make a total analysis using modal transient response impractical. It was also discovered that inaccuracies may exist for NASTRAN rigid body modes on large models when Given's method for eigenvalue extraction is employed. The deficiencies as well as recommendations for improving the NASTRAN program are discussed.
The neutral g-modes of a degenerate fluid at zero temperature are analyzed. The g-modes of a degenerate fluid at finite but small temperatures are then expanded in terms of those of the zero temperature fluid. For nonrelativistic degenerate fluids it is found that (1) the g-eigenvalues are proportional to T mu(6)sub e mu(-1)sub i, where T is the internal temperature of the fluid, mu sub e and mu sub i are the mean molecular weights of electrons and ions, respectively; (2) the ion pressure is solely responsible for driving the g-modes. For white dwarfs of about a solar mass, the periods of the g-oscillations are in the range of a few hundredths of seconds.
An analytical technique for the prediction of fan blade flutter was evaluated by utilizing first stage fan flutter data from tests on an advanced high performance engine. The formulation includes both aerodynamic and mechanical coupling among all the blades of the assembly. Mistuning is accounted for in the analysis so that individual blade inertias, frequencies, or damping can be considered. Airfoil stability was predicted by calculating a flutter determinant, the eigenvalues of which indicate the extent of susceptibility to flutter. When blade to blade differences in frequencies are considered, a stable system is predicted for the test points examined. For a tuned system, it was found that torsional flutter can be predicted at a limited number of interblade phase angles. Examination of these phase angles indicated that they were "close" to the condition of acoustic resonance. For the range of Mach numbers and reduced frequencies considered, the so called subcritical flutter cannot be predicted. The essential influence of mechanical coupling among the blades is to change the frequencies of the system with little or no change in damping; however, aerodynamic coupling together with mechanical coupling could change not only frequencies, but also damping in the system, with a trend toward instability.
Noise generated by air flow through a coaxial obstruction in a long, straight pipe was investigated with concentration on the modal characteristics of the noise field inside the pipe and downstream of the restriction. Two measurement techniques were developed for separation of the noise into the acoustic duct modes. The instantaneous mode separation technique uses four microphones, equally spaced in the circumferential direction, at the same axial location. The time-averaged mode separation technique uses three microphones mounted at the same axial location. A matrix operation on time-averaged data produces the modal pressure levels. This technique requires the restrictive assumption that the acoustic modes are uncorrelated with each other. The measured modal pressure spectra were converted to modal power spectra and integrated over the frequency range 200-6000 Hz. The acoustic efficiency levels (acoustic power normalized by jet kinetic energy flow), when plotted vs. jet Mach number, showed a strong dependence on the ratio of restriction diameter to pipe diameter. The acoustic energy flow analyses based on the thermodynamic energy equation and on the results of Mohring both resulted in orthogonality properties for the eigenfunctions of the radial mode shape equation. These orthogonality relationships involve the eigenvalues and derivatives of the radial mode shape functions.
A mixed formulation for calculating static equilibrium and stability eigenvalues of nonuniform rotor blades in hover is presented. The static equilibrium equations are nonlinear and are solved by an accurate and efficient collocation method. The linearized perturbation equations are solved by a one step, second order integration scheme. The numerical results correlate very well with published results from a nearly identical stability analysis based on a displacement formulation. Slight differences in the results are traced to terms in the equations that relate moments to derivatives of rotations. With the present ordering scheme, in which terms of the order of squares of rotations are neglected with respect to unity, it is not possible to achieve completely equivalent models based on mixed and displacement formulations. The one step methods reveal that a second order Taylor expansion is necessary to achieve good convergence for nonuniform rotating blades. Numerical results for a hypothetical nonuniform blade, including the nonlinear static equilibrium solution, were obtained with no more effort or computer time than that required for a uniform blade.
A method for selecting a multivariable state feedback controller that simultaneously achieves an a priori specification on closed loop eigenvalues and good mode mixing is presented. The problem is solved by projecting a desired modal matrix onto a constraint set containing the null space of the closed-loop state matrix, while assuring that the projection is in the null space. The feedback matrix follows immediately in the formulation. An example involving a helicopter hover controller is presented.
The paper considers principal component analysis of solar flares in the soft X-ray flux, a technique for extracting the salient features from a mass of data. The method applies particularly to the analysis of nonstationary ensembles, and its computations require the evaluation of eigenvalues of matrices. The Eispack matrix eigen system routines were used to analyze full-disk proportional-counter data from the X-ray event analyzer which was part of the Skylab experiment. Empirical orthogonal functions were derived for events in the soft X-ray spectrum between 2.5 and 20 A during different time periods, indicating that about 90% of the cumulative power of each analyzed flare is contained in the largest eigenvector. The first two largest eigenvectors are sufficient for an empirical curve fit through the raw data and a characterization of solar flares in the soft X-ray flux, and power spectra of two largest eigenvectors reveal a reported periodicity of about 5 min.
The stability of a beam subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. The differential equations of motion are solved numerically using an integrating matrix method in combination with an eigenanalysis to determine the eigenvalues from which stability is assessed. The results clarify several differences which have been identified in the literature relating to the proper behavior of the critical rotational speed for buckling as the radius of rotation of the clamped end of the beam is reduced.
This brief paper summarizes the approach the authors will take in designing a feedback controller for the F-100 turbofan engine. The technique to be utilized simultaneously realizes dominant closed-loop eigenvalues, approximates specified modal behavior, and achieves low eigensystem sensitivity with respect to certain plant parameter variations.
An analysis is presented of the parallel neutral stability of three-dimensional incompressible, isothermal boundary-layer flows. A Taylor-series expansion of the dispersion relation is used to derive the general eigenvalues. These equations are functions of the complex group velocity. These relations are verified by numerical results obtained for two- and three-dimensional disturbances in two- and three-dimensional flows.
This paper is an expanded version of a talk given at the 1979 T.I.C.O.M. conference. It is a self-contained introduction, for applied mathematicians and numerical analysts, to quantum mechanics and quantum field theory. It also contains a brief description of the authors' numerical approach to the problems of quantum field theory, which may best be summarized by the question; Can we compute the eigenvalues and eigenfunctions of Schrodinger operators in infinitely many variables.
Equations are developed for the buckling of a general lattice structure that has repetitive geometry. Equilibrium at a typical node is expressed using finite element techniques, and the only assumption is that the response is periodic. By basing the stiffness matrix on the exact solution of the beam column equation, accurate results are obtained for complex buckling behavior that would require a very large system of equations using conventional techniques. The present method requires the eigenvalues of only a 6x6 determinant. The results are used to study the buckling of isogrid cylinders, three-element truss columns and polygonal rings. Details of the analysis including expressions for all terms in the governing stability determinant are given.
A method developed to determine maximum changes in analytical mass and stiffness matrices to make them consistent with a set of measured normal modes and natural frequencies is presented. The corrected model will be an improved base for studies of physical changes, boundary condition changes, and for prediction of forced responses. The method features efficient procedures not requiring solutions of the eigenvalue problem, and the ability to have more degrees of freedom than the test data. In addition, modal displacements are obtained for all analytical degrees of freedom, and the frequency dependence of the coordinate transformations is properly treated.
This paper presents a new method of order reduction based on recent work of a similar nature applicable to system dynamics and control. This method provides a practical computational procedure for producing a condensed model which exactly preserves the slowest n1 modes (where n1 is the number of 'slow' eigenvalues) of the total n modes of the original model for almost any set of n1 degrees of freedom retained in the condensed model. The method can also be used to compute the eigensolutions corresponding to the n1 slowest modes of the original structural dynamics problem. The accuracy of the condensed model and the speed/accuracy performance of the eigensolver are compared with standard methods for a 90 DOF cantilevered plate.