Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

HADY-I, a FORTRAN program for the compressible stability analysis of three-dimensional boundary layers

A computer program HADY-I for calculating the linear incompressible or compressible stability characteristics of the laminar boundary layer on swept and tapered wings is described. The eigenvalue problem and its adjoint arising from the linearized disturbance equations with the appropriate boundary conditions are solved numerically using a combination of Newton-Raphson interative scheme and a variable step size integrator based on the Runge-Kutta-Fehlburh fifth-order formulas. The integrator is used in conjunction with a modified Gram-Schmidt orthonormalization procedure. The computer program HADY-I calculates the growth rates of crossflow or streamwise Tollmien-Schlichting instabilities. It also calculates the group velocities of these disturbances. It is restricted to parallel stability calculations, where the boundary layer (meanflow) is assumed to be parallel. The meanflow solution is an input to the program.

El-Hady, N. M.↗

Iterative optical vector-matrix processors (survey of selected achievable operations)

An iterative optical vector-matrix multiplier with a microprocessor-controlled feedback loop capable of performing a wealth of diverse operations was described. A survey and description of many of its operations demonstrates the versatility and flexibility of this class of optical processor and its use in diverse applications. General operations described include: linear difference and differential equations, linear algebraic equations, matrix equations, matrix inversion, nonlinear matrix equations, deconvolution and eigenvalue and eigenvector computations. Engineering applications being addressed for these different operations and for the IOP are: adaptive phased-array radar, time-dependent system modeling, deconvolution and optimal control.

Casasent, D.↗

Statistical properties of MHD fluctuations associated with high speed streams from HELIOS 2 observations

Helios 2 magnetic data were used to obtain several statistical properties of MHD fluctuations associated with the trailing edge of a given stream served in different solar rotations. Eigenvalues and eigenvectors of the variance matrix, total power and degree of compressibility of the fluctuations were derived and discussed both as a function of distance from the Sun and as a function of the frequency range included in the sample. The results obtained add new information to the picture of MHD turbulence in the solar wind. In particular, a dependence from frequency range of the radial gradients of various statistical quantities is obtained.

Bavassano, B.↗

Application of the Finite Element Method to Rotary Wing Aeroelasticity

A finite element method for the spatial discretization of the dynamic equations of equilibrium governing rotary-wing aeroelastic problems is presented. Formulation of the finite element equations is based on weighted Galerkin residuals. This Galerkin finite element method reduces algebraic manipulative labor significantly, when compared to the application of the global Galerkin method in similar problems. The coupled flap-lag aeroelastic stability boundaries of hingeless helicopter rotor blades in hover are calculated. The linearized dynamic equations are reduced to the standard eigenvalue problem from which the aeroelastic stability boundaries are obtained. The convergence properties of the Galerkin finite element method are studied numerically by refining the discretization process. Results indicate that four or five elements suffice to capture the dynamics of the blade with the same accuracy as the global Galerkin method.

Straub, F. K.↗

Effects of rotation and magnetic field on the onset of convective instability in a liquid layer due to buoyancy and surface tension

Thermocapillary stability characteristics of a horizontal liquid layer heated from below rotating about a vertical axis and subjected to a uniform vertical magnetic field are analyzed under a variety of thermal and electromagnetic boundary conditions. Results based on analytical solutions to the pertinent eigenvalue problems are discussed in the light of earlier work on special cases of the more general problem considered here to show in particular the effects of the heat transfer, nonzero curvature and gravity waves at the two-fluid interface. Although the expected stabilizing action of the Coriolis and Lorentz force fields in this configuration are in evidence the optimal choice of an appropriate range for the relevant parameters is shown to be critically dependent on the interfacial effects mentioned above.

Sarma, G. S. R.↗

Eigenspace techniques for active flutter suppression

Mathematical models to be used in the control system design were developed. A computer program, which takes aerodynamic and structural data for the ARW-2 aircraft and converts these data into state space models suitable for use in modern control synthesis procedures, was developed. Reduced order models of inboard and outboard control surface actuator dynamics and a second order vertical wind gust model were developed. An analysis of the rigid body motion of the ARW-2 was conducted. The deletion of the aerodynamic lag states in the rigid body modes resulted in more accurate values for the eigenvalues associated with the plunge and pitch modes than were obtainable if the lag states were retained.

Garrard, W. L.↗

Application of the reduced basis method to nonlinear transient thermal analysis

An effort to apply the reduced basis method to nonlinear transient thermal analysis is described. The method combines the classical Rayleigh-Ritz and modal superposition techniques with contemporary finite element methods to retain modeling versatility as the degrees of freedom in a problem are reduced. The essence of the method is to use a few thermal modes from eigenvalue analyses as basis vectors to represent the temperature response for a given thermal problem similar to the use of vibration modes to represent displacements in a dynamic response problem. Approximate temperature distributions were obtained using the reduced basis method for a small section of the Shuttle Orbiter lower wing undergoing reentry heating. Good agreement was obtained between the reduced basis method solutions and full system solutions with reductions in the degrees of freedom of up to a factor of four. The good agreement indicates the reduced basis method has the potential for significant reduction in computing effort for thermal analysis; however, considerable work remains to determine techniques for selecting the type and number of basis vectors needed for approximate solutions to more complex transient thermal problems.

Shore, C. P.↗

Constrained multibody system dynamics: An automated approach

The governing equations for constrained multibody systems are formulated in a manner suitable for their automated, numerical development and solution. The closed loop problem of multibody chain systems is addressed. The governing equations are developed by modifying dynamical equations obtained from Lagrange's form of d'Alembert's principle. The modifications is based upon a solution of the constraint equations obtained through a zero eigenvalues theorem, is a contraction of the dynamical equations. For a system with n-generalized coordinates and m-constraint equations, the coefficients in the constraint equations may be viewed as constraint vectors in n-dimensional space. In this setting the system itself is free to move in the n-m directions which are orthogonal to the constraint vectors.

Kamman, J. W.↗

Modeling of human operator dynamics in simple manual control utilizing time series analysis

Time series analysis is applied to model human operator dynamics in pursuit and compensatory tracking modes. The normalized residual criterion is used as a one-step analytical tool to encompass the processes of identification, estimation, and diagnostic checking. A parameter constraining technique is introduced to develop more reliable models of human operator dynamics. The human operator is adequately modeled by a second order dynamic system both in pursuit and compensatory tracking modes. In comparing the data sampling rates, 100 msec between samples is adequate and is shown to provide better results than 200 msec sampling. The residual power spectrum and eigenvalue analysis show that the human operator is not a generator of periodic characteristics.

Agarwal, G. C.↗

A quadratic weight selection algorithm

A new numerical algorithm is presented which determines a positive semi-definite state weighting matrix in the linear-quadratic optimal control design problem. The algorithm chooses the weighting matrix by placing closed-loop eigenvalues and eigenvectors near desired locations using optimal feedback gains. A simplified flight control design example is used to illustrate the algorithms capabilities.

Broussard, J. R.↗

Eikonal approximation for charge transfer from a multielectron atom to fast projectiles

The eikonal approach developed previously for calculating electron-capture cross sections for bare projectiles colliding with hydrogenic targets is extended here to allow for multielectron targets. Both the impact and wave pictures are employed and their equivalence is discussed. As a first approximation, each atomic orbital is specified by the three hydrogenic quantum numbers, an effective nuclear charge Z sub t, and an energy eigenvalue in the impact picture, or ionization potential in the wave picture. The Z sub t prime appearing in the eikonal phase factor is left undetermined because of incomplete information on the many-body target. However, analytic expressions are derived for the theoretical cross sections, and numerical values are calculated for simple choices of Z sub t prime. Those results are compared with existing experimental data for C, Ne, Ar, N2, O2, and He targets.

Ho, T. S.↗

An integration scheme to determine the dynamic response of a launch vehicle with several payloads

The coupled equations of motion of a system composed of a launch vehicle and multiple payloads are derived. In the process, it is shown how superfluous interface degrees of freedom on the booster side can be accommodated in the formulation. The discrete system equations are directly integrated, avoiding the solution of an expensive system eigenvalue problem. A modified Newmark-Chan-Beta numerical integration scheme is used to obtain the response. The unique form of the equations of motion allow for a quick and accurate solution both for the response and the internal loads. Although completely general, the presented technique is most effective when the number of interface degrees of freedom is relatively small compared to the overall number of degrees of freedom of the system. The technique is applied to the case of the S.T.S.-S.T.-OMS Kit system. Finally, an evaluation of the presented method is included.

Engels, R. C.↗

Accurate numerical solution of compressible, linear stability equations

The present investigation is concerned with a fourth order accurate finite difference method and its application to the study of the temporal and spatial stability of the three-dimensional compressible boundary layer flow on a swept wing. This method belongs to the class of compact two-point difference schemes discussed by White (1974) and Keller (1974). The method was apparently first used for solving the two-dimensional boundary layer equations. Attention is given to the governing equations, the solution technique, and the search for eigenvalues. A general purpose subroutine is employed for solving a block tridiagonal system of equations. The computer time can be reduced significantly by exploiting the special structure of two matrices.

Malik, M. R.↗

On the modelling and simulation of the dynamics and control of large flexible orbiting systems

This paper attempts to review the steps involved in the development of mathematical models that can be used to simulate the in-orbit dynamics of large flexible systems. The use of graph theoretic techniques can often be used to reduce the computational effort involved for calculating the eigenvalues of large ordered systems. Computer generated graphical techniques may provide additional insight into the understanding of elastic modal shape functions of complex systems. Finally the numerical techniques commonly used to develop shape and attitude control laws will be briefly reviewed.

Bainum, P. M.↗

Symmetric baroclinic instability of a Hadley cell

A symmetric baroclinic instability is examined in terms of a Boussinesq fluid contained between two horizontal plates to determine the effects of the Ekman and thermal layers. Governing equations are written for a rotating reference frame, taking into account the Rossby, Ekman, and Prandtl numbers. Equations are defined for the perturbation functions, treated as an eigenvalue problem, and a numerical integration of the full eighth order differential system is performed by a shooting technique. An instability is found to occur in the Hadley cell containing both Ekman and thermal boundary layers when the Richardson number is close to unity. If the Prandtl number is fixed the critical Richardson number decreases with an increasing Ekman number until the Ekman number reaches a certain value, at which time the fluid is stable.

Antar, B. N.↗

Design and flight testing of a digital optimal control general aviation autopilot

This paper presents the designs of Proportional-Integral-Filter (PIF) autopilots for a General Aviation (NAVION) aircraft. The PIF autopilots use modern control theory to determine heading select and altitude select and hold autopilot modes. The PIF control law uses typical General Aviation sensors for state feedback; command error integration for command tracking; digital complimentary filtering and analog prefiltering for sensor noise suppression; a control filter for computation delay accommodation; and the incremental form to eliminate trim values in implementation. Theoretical developments for the control law are described which combine the sampled-data regulator with command generator tracking for use as a digital flight control system. The digital PIF autopilots are evaluated using closed-loop eigenvalues and simulations. Successful flight test results for the PIF autopilots are presented for different turbulence conditions and quadratic weights.

Broussard, J. R.↗

Sampling errors in the estimation of empirical orthogonal functions

Empirical Orthogonal Functions (EOF's), eigenvectors of the spatial cross-covariance matrix of a meteorological field, are reviewed with special attention given to the necessary weighting factors for gridded data and the sampling errors incurred when too small a sample is available. The geographical shape of an EOF shows large intersample variability when its associated eigenvalue is 'close' to a neighboring one. A rule of thumb indicating when an EOF is likely to be subject to large sampling fluctuations is presented. An explicit example, based on the statistics of the 500 mb geopotential height field, displays large intersample variability in the EOF's for sample sizes of a few hundred independent realizations, a size seldom exceeded by meteorological data sets.

North, G. R.↗

Stability of mean flows over an infinite flat plate

The present investigation is mainly concerned with a stability analysis for the linearized Navier-Stokes equations for parallel and nonparallel mean flows over an infinite flat plate. The system of equations for parallel flows is presented. The system is viewed as a generalized Orr-Sommerfeld equation. Attention is given to an explicit criterion characterizing the case when the stability of all physically reasonable solutions is determined by the eigenvalues. The proof given in the investigation is applicable to both the generalized Orr-Somerfeld equations and the modified equations for nonparallel flow. The fact that the criterion is independent of the completeness or incompleteness of eigenfunctions is contrary to some expectations.

Miklavcic, M.↗