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 307 records · Page 17

Application of Variational Methods to the Thermal Entrance Region of Ducts

A variational method is presented for solving eigenvalue problems which arise in connection with the analysis of convective heat transfer in the thermal entrance region of ducts. Consideration is given, to both situations where the temperature profile depends upon one cross-sectional coordinate (e.g. circular tube) or upon two cross-sectional coordinates (e.g. rectangular duct). The variational method is illustrated and verified by application to laminar heat transfer in a circular tube and a parallel-plate channel, and good agreement with existing numerical solutions is attained. Then, application is made to laminar heat transfer in a square duct as a check, an alternate computation for the square duct is made using a method indicated by Misaps and Pohihausen. The variational method can, in principle, also be applied to problems in turbulent heat transfer.

Sparrow, E. M.↗

On Dynamics of Spinning Structures

This paper provides details of developments pertaining to vibration analysis of gyroscopic systems, that involves a finite element structural discretization followed by the solution of the resulting matrix eigenvalue problem by a progressive, accelerated simultaneous iteration technique. Thus Coriolis, centrifugal and geometrical stiffness matrices are derived for shell and line elements, followed by the eigensolution details as well as solution of representative problems that demonstrates the efficacy of the currently developed numerical procedures and tools.

Rotating bodies↗

Multidisciplinary Design Optimization and Cruise Mach Number Study of Truss-Braced Wing Aircraft

The Subsonic Ultra Green Aircraft Research (SUGAR) Phase III was led by Dr. Rakesh K. Kapania and Dr. Joseph A. Schetz at the Multidisciplinary Analysis and Design Center for Advanced Vehicles, Department of Aerospace and Ocean Engineering, Virginia Tech, Blacksburg VA. The research was performed from December 2014 to December 2015. Three major areas were investigated: Multidisciplinary Design Optimization (MDO) studies of truss braced wing (TBW) and strut braced wing (SBW) vehicles at cruise Mach numbers of 0.7 and 0.8 for a flight mission similar to current market single aisle configurations. The performance and the characteristics of the optimized vehicles were compared to the SUGAR Phase II TBW vehicle. These results were obtained without applying any of the extended transonic aerodynamic and aeroelastic tools that will be discussed later. It was found that the cruise Mach number has a large effect on the best “truss” configuration. At Mach 0.7, an SBW has a better fuel consumption and better take-off gross weight (TOGW). However, at Mach 0.8, the TBW is superior because the jury strut aids in satisfying the flutter constraint; Two-dimensional, steady, transonic aerodynamic analysis of the Boeing Airfoil J (BACJ) airfoil was performed for a range of thickness ratios, Mach numbers and lift coefficients. Reynolds-averaged Navier-Stokes (RANS) equations were solved to obtain the lift-curve slope, wave drag coefficient, the location of the center of pressure and to predict the separation at the trailing edge, which may lead to buffeting. One of the goals was to develop a database of lift-curve slope and the location of center of pressure, which could be used in a transonic aeroelastic analysis. Another goal was to compare the wave drag coefficients to those predicted by Lock’s fourth-power law and also to compare the transonic effects obtained from RANS simulations to those predicted by the Korn equations. A third goal was to develop a buffet boundary, which can be integrated into the MDO framework to prevent the optimized designs from probable buffeting; A state-space transonic aeroelastic analysis tool was developed, which can incorporate the nonlinear transonic effects in the unsteady aerodynamics but is yet computationally cheap when used within the VT MDO framework. The aeroelastic analysis uses Leishman- Beddoes (LB) indicial functions, which generated a state-space representation of the aeroelastic system. The indicial functions allow the incorporation of data for steady lift-curve slope and location of the center of pressure. Thus, the steady transonic effects are included, and the unsteady aerodynamic responses are a linearization about the steady results. The aeroelastic approach discretizes the wing into numerous strips, which results in a large eigenvalue problem as each strip has eight augmented aerodynamic states as per the LB theory. Thus, to reduce the computation expense, a reduced order model (ROM) was developed. The approach was validated using a few examples.

Kapania, Rakesh K.↗

Adjoint-based Sensitivities of Flutter Predictions based on the Linearized Frequency-domain Approach

Flutter is a critical factor in designing and certifying aircraft. The linearized frequency-domain method offers a lower cost alternative to time-marching computational fluid dynamics for high-fidelity flutter analysis. In this work, adjoint-based sensitivities are added to a flutter analysis based on the linearized frequency-domain method to efficiently compute derivatives of flutter cost functions with respect to design variables or uncertain parameters. The derivation of the adjoint equations, which involve complications such as derivatives of a nonlinear generalized eigenvalue problem with complex-valued inputs and derivatives of the linearized Navier-Stokes equations, is provided. The implemented adjoint terms and derivatives are verified before demonstrating the approach for derivatives of flutter dynamic pressure with respect to Mach number for the AGARD 445.6 wing.

Aeroelasticity↗

Aerodynamic Design Optimization for Natural Laminar Flow Airfoils

Natural laminar flow technology is a passive laminar flow control (LFC) strategy that seeks to delay the onset of boundary-layer transition (BLT) through shape optimization to reduce the drag of the aerodynamic vehicle. Adjoint-based design optimization for LFC is proposed in an integrated multidisciplinary framework, which includes the computational fluid dynamics (CFD), geometry and grid deformation, and linear stability analysis (LSA) for transition prediction. In particular, the BLT location is predicted using the dual N-factor method that is based on a linear stability theory (LST) eigenvalue problem. The dual N-factor criterion accounts for the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) boundary-layer instabilities to predict the transition location in three-dimensional boundary-layer flows. The adjoint-based shape optimization procedure is based on an iteratively coupled CFD and LSA methodology to converge the transition location and flow solutions, as well as to calculate the sensitivities of the aerodynamic metrics of interest with respect to the flow and shape design parameters. The RAE 2822 airfoil at 0 and 30 degrees yaw angles, an angle of attack of 0.72 degrees, and subsonic conditions (M∞ = 0.19, Rec = 5.6 × 106 ) are used as baseline configurations for design optimization. The angle of attack and the vertical displacement of free-form-deformation control points are used as design variables to reduce the drag coefficient while reaching a specified lift coefficient. The optimized unswept airfoil designs achieve a 30% drag reduction accompanied by a downstream shift of the transition locations over both suction and pressure sides of the airfoil. The initial design iterations for the swept case also show a favorable trend in the drag reduction with transition delay over both sides.

Transition↗

Eigenvalues and stability problems of rotors

The essential theoretical results of the application of a developed transfer matrix method to the free transverse vibration of a rotor are shown. Gyroscopic and shear effects, rotary inertia, and external and internal damping as well as the influence of sleeve bearings and rotor supports are taken into consideration. The eigenvalues of the motion equations of the rotor are searched by using a modified determinant method.

Walczyk, Z.↗

Arbitrary eigenvalue assignments for linear time-varying multivariable control systems

The problem of eigenvalue assignments for a class of linear time-varying multivariable systems is considered. Using matrix operators and canonical transformations, it is shown that a time-varying system that is 'lexicography-fixedly controllable' can be made via state feedback to be equivalent to a time-invariant system whose eigenvalues are arbitrarily assignable. A simple algorithm for the design of the state feedback is provided.

Nguyen, Charles C.↗

The eigenvalue spectrum of the Orr-Sommerfeld problem

A numerical investigation of the temporal eigenvalue spectrum of the ORR-Sommerfeld equation is presented. Two flow profiles are studied, the plane Poiseuille flow profile and the Blasius boundary layer (parallel): flow profile. In both cases a portion of the complex c-plane bounded by 0 less than or equal to CR sub r 1 and -1 less than or equal to ci sub i 0 is searched and the eigenvalues within it are identified. The spectra for the plane Poiseuille flow at alpha = 1.0 and R = 100, 1000, 6000, and 10000 are determined and compared with existing results where possible. The spectrum for the Blasius boundary layer flow at alpha = 0.308 and R = 998 was found to be infinite and discrete. Other spectra for the Blasius boundary layer at various Reynolds numbers seem to confirm this result. The eigenmodes belonging to these spectra were located and discussed.

Antar, B. N.↗

The stability of numerical methods for second order ordinary differential equations

An important characterization of a numerical method for first order ODE's is the region of absolute stability. If all eigenvalues of the linear problem dy/dt = Ay are inside this region, the numerical method is stable. If the second order system d/dt(dy/dt) = 2Ady/dt - By is solved as a first order system, the same result applies to the eigenvalues of the generalized eigenvalue problem (lambda-squared)I 2(lambda)A + B. No such region exists for general methods for second order equations, but in some cases a region of absolute stability can be defined for methods for the single second order equation d/dt(dy/dt) = 2ady/dt - by. The absence of a region of absolute stability can occur when different members of a system of first order equations are solved by different methods.

Gear, C. W.↗

A study of the response of nonlinear springs

The various phases to developing a methodology for studying the response of a spring-reinforced arch subjected to a point load are discussed. The arch is simply supported at its ends with both the spring and the point load assumed to be at midspan. The spring is present to off-set the typical snap through behavior normally associated with arches, and to provide a structure that responds with constant resistance over a finite displacement. The various phases discussed consist of the following: (1) development of the closed-form solution for the shallow arch case; (2) development of a finite difference analysis to study (shallow) arches; and (3) development of a finite element analysis for studying more general shallow and nonshallow arches. The two numerical analyses rely on a continuation scheme to move the solution past limit points, and to move onto bifurcated paths, both characteristics being common to the arch problem. An eigenvalue method is used for a continuation scheme. The finite difference analysis is based on a mixed formulation (force and displacement variables) of the governing equations. The governing equations for the mixed formulation are in first order form, making the finite difference implementation convenient. However, the mixed formulation is not well-suited for the eigenvalue continuation scheme. This provided the motivation for the displacement based finite element analysis. Both the finite difference and the finite element analyses are compared with the closed form shallow arch solution. Agreement is excellent, except for the potential problems with the finite difference analysis and the continuation scheme. Agreement between the finite element analysis and another investigator's numerical analysis for deep arches is also good.

Hyer, M. W.↗

The application of the Routh approximation method to turbofan engine models

The Routh approximation technique is applied in the frequency domain to a 16th order state variable turbofan engine model. The results obtained motivate the extension of the frequency domain formulation of the Routh method to the time domain to handle the state variable formulation directly. The time domain formulation is derived, and a characterization, which specifies all possible Routh similarity transformations, is given. The characterization is computed by the solution of two eigenvalue eigenvector problems. The application of the time domain Routh technique to the state variable engine model is described, and some results are given.

Merrill, W. C.↗

Eigenvalue/eigenvector assignment using output feedback

The problem of eigenvalue assignment in a linear time-invariant system using output feedback is considered. New sufficient conditions are derived to assign an almost arbitrary set of minimum (n,m + r - 1) distinct eigenvalues where n, m, and r are the number of states, inputs, and outputs, respectively. These conditions precisely identify the class of systems where such an assignment is impossible. The synthesis technique also highlights the freedom in selection of closed-loop eigenvectors under output feedback. The utility of eigenvalue/eigenvector assignment in transient response shaping is illustrated by the design of a controller for the lateral dynamics of an aircraft.

Srinathkumar, S.↗

Applied Routh approximation

The Routh approximation technique for reducing the complexity of system models was applied in the frequency domain to a 16th order, state variable model of the F100 engine and to a 43d order, transfer function model of a launch vehicle boost pump pressure regulator. The results motivate extending the frequency domain formulation of the Routh method to the time domain in order to handle the state variable formulation directly. The time domain formulation was derived and a characterization that specifies all possible Routh similarity transformations was given. The characterization was computed by solving two eigenvalue-eigenvector problems. The application of the time domain Routh technique to the state variable engine model is described, and some results are given. Additional computational problems are discussed, including an optimization procedure that can improve the approximation accuracy by taking advantage of the transformation characterization.

Merrill, W. C.↗

Stability of constant gain systems with vector feedback

The state space, the controllability, and the observability concepts are discussed in connection with the proposed stability analysis which permits drastic dimensional reductions for a vector feedback problem. Any constant gain system's stability can thus be analyzed in the frequency domain with a single Nyquist plot. The analysis considers the total system with all loops closed, a disturbance vector as input, and the feedback vector as output. All constant gain systems are shown to be decomposable into stable subsystems where the degree of the decomposition determines the dimensions. The maximum decomposition results in the state-space approach which is the limit case. The method is demonstrated with the stability analysis of the pogo phenomenon, an oscillatory interaction between the propulsion and the structure of a space vehicle. This problem, with eigenvalues over a hundred, was drastically but rigorously reduced to a stability analysis of a 4x4 matrix.

Vonpragenau, G. L.↗

The application of the Routh approximation method to turbofan engine models

The Routh approximation technique is applied in the frequency domain to a 16th order state variable turbofan engine model. The results obtained motivate the extension of the frequency domain formulation of the Routh method to the time domain to handle the state variable formulation directly. The time domain formulation is derived and a new characterization, which specifies all possible Routh similarity transformations, is given. This characterization is computed by the solution of two eigenvalue-eigenvector problems. The application of the time domain Routh technique to the state variable engine model is described and some results are given. Additional computational problems are discussed including an optimization procedure which can improve the approximation accuracy by taking advantage of the transformation characterization.

Merrill, W. C.↗