Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear stability analysis”

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 775 records · Page 43

Implementation of a partitioned algorithm for simulation of large CSI problems

The implementation of a partitioned numerical algorithm for determining the dynamic response of coupled structure/controller/estimator finite-dimensional systems is reviewed. The partitioned approach leads to a set of coupled first and second-order linear differential equations which are numerically integrated with extrapolation and implicit step methods. The present software implementation, ACSIS, utilizes parallel processing techniques at various levels to optimize performance on a shared-memory concurrent/vector processing system. A general procedure for the design of controller and filter gains is also implemented, which utilizes the vibration characteristics of the structure to be solved. Also presented are: example problems; a user's guide to the software; the procedures and algorithm scripts; a stability analysis for the algorithm; and the source code for the parallel implementation.

Alvin, Kenneth F.↗

Orbital stability in combined uniform axial and three-dimensional wiggler magnetic fields for free-electron lasers

Zachary Phys. Rev. A 29 (6), 3224 (1984) recently analyzed the instability of relativistic-electron helical trajectories in combined uniform axial and helical wiggler magnetic fields when the radial variation of the wiggler field is taken into account. It is shown here that the type 2 instability comprised of secular terms growing linearly in time, identified by Zachary and earlier by Diament Phys. Rev. A 23 (5), 2537 (1981), is an artifact of simple perturbation theory. A multiple-time-scale perturbation analysis reveals a nonsecular evolution on a slower time scale which accommodates an arbitrary initial perturbation. It is shown that, in the absence of exponential instability, the electron seeks a modified helical orbit more appropriate to its perturbed state and oscillates stably about it. Thus, the perturbed motion is oscillatory but nonsecular, and hence the helical orbits are stable.

Johnston, S.↗

Adaptive control of Space Station with control moment gyros

An adaptive approach to Space Station attitude control is investigated. The main components of the controller are the parameter identification scheme, the control gain calculation, and the control law. The control law is a full-state feedback space station baseline control law. The control gain calculation is based on linear-quadratic regulator theory with eigenvalues placement in a vertical strip. The parameter identification scheme is a recursive extended Kalman filter that estimates the inertias and also provides an estimate of the unmodeled disturbances due to the aerodynamic torques and to the nonlinear effects. An analysis of the inertia estimation problem suggests that it is possible to estimate Space Station inertias accurately during nominal control moment gyro operations. The closed-loop adaptive control law is shown to be capable of stabilizing the Space Station after large inertia changes. Results are presented for the pitch axis.

Bishop, Robert H.↗

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗

Supersonic Coaxial Jets: Noise Predictions and Measurements

The noise from perfectly expanded coaxial jets was measured in an anechoic chamber for different operating conditions with the same total thrust, mass flow, and exit area. The shape of the measured noise spectrum at different angles to the jet axis was found to agree with spectral shapes for single, axisymmetric jets. Based on these spectra, the sound was characterized as being generated by large turbulent structures or fine-scale turbulence. Modeling the large scale structures as instability waves, a stability analysis was conducted for the coaxial jets to identify the growing and decaying instability waves in each shear layer and predict their noise radiation pattern outside the jet. When compared to measured directivity, the analysis identified the region downstream of the outer potential core, where the two shear layers were merging, as the source of the peak radiated noise where instability waves, with their origin in the inner shear layer, reach their maximum amplitude. Numerical computations were also performed using a linearized Euler equation solver. Those results were compared to both the results from the instability wave analysis and to measured data.

Dahl, Milo D.↗

Investigation of temperature and concentration oscillations in the directional solidification of Pb-Sn-Te

Directional solidification of the pseudobinary compound semiconductor material Pb sub 1-x Sn sub x Te by the Bridgman crystal growth process will be studied. Natural convection in the molten sample will be visualized with a novel electrochemical cell technique that employs the solid electrolyte material yttria-stabilized zirconia. Mass transfer by both diffusion and convection will be measured by detecting the motion of oxygen tracer in the liquid. Additional applications for electrochemical cells in semiconductor crystal growth are suggested. Unsteady convection in the melt will also be detected by the appearance of temperature oscillations. The purpose of this study is to experimentally characterize the overstable conditions for a Pb sub 1-x Sn sub x Te melt in the vertical Bridgman crystal growth technique and use a linear analysis to predict the onset of convection for this system.

Anderson, T. J.↗

Guaranteed properties for nonlinear gain scheduled control systems

An analysis is presented of two types of nonlinear gain-scheduled control systems: scheduling on a reference trajectory and scheduling on the plant output. Conditions are given which guarantee stability, robustness, and performance properties of the global gain-scheduled designs. These conditions confirm and formalize popular notions regarding gain-scheduled designs; e.g., that the scheduling variable should vary slowly and capture the plant's nonlinearities. These results extend the work of Shamma and Athans (1987), which addressed the case of linear plants whose dynamics depends on exogenous parameters.

Shamma, Jeff S.↗

Determination of Longitudinal Stability and Control Characteristics from Free-Flight Model Tests with Results at Transonic Speeds for Three Airplane Configurations

A test technique and data analysis method has been developed for determining the longitudinal aerodynamic characteristics from free-flight tests of rocket-propelled models. The technique makes use of accelerometers and an angle-of-attack indicator to permit instantaneous measurements of lift, drag, and pitching moments. The data, obtained during transient oscillations resulting from control-surface disturbances, are analyzed by essentially nonlinear direct methods (such as cross plots of the variation of lift coefficient with angle of attack) and by linear indirect methods by using the equations of motion for a transient oscillation. The analysis procedure has been set forth in some detail and the feasibility of the method has been demonstrated by data measured through the transonic speed range on several airplane configurations. It was shown that the flight conditions and dynamic similitude factors for the tests described were reasonably close to typical full-scale airplane conditions.

Gillis, Clarence L↗

Aeroelastic stability and response of horizontal axis wind turbine blades

The coupled flap-lag-torsion equations of motion of an isolated horizontal axis wind turbine blade are formulated. Quasi-steady blade-element strip theory was applied to derive the aerodynamic operator which includes boundary layer type gradient winds. The final equations which have periodic coefficients were solved in order to obtain the aeroelastic response and stability of large horizontal axis wind turbine blade. A new method of generating an appropriate time-dependent equilibrium position (required for the stability analysis) has been implemented. Representative steady-state responses and stability boundaries, applicable mainly to an existing blade design (NASA/-ERDA MOD-0), are presented. The results indicate that the MOD-0 configuration is a basically stable design and that blade stability is not sensitive to offsets between blade elastic axis and aerodynamic center. Blade stability appears to be sensitive to precone. The tower shadow (or wake) has a considerable effect on the flap response but leaves blade stability unchanged. Finally, it was found that non linear terms in the equations of motion can significantly affect the linearized stability boundaries, however, these terms have a negligible effect on blade response at operating conditions.

Kottapalli, S. B. R.↗

A modal aeroelastic analysis scheme for turbomachinery blading

An aeroelastic analysis is developed which has general application to all types of axial-flow turbomachinery blades. The approach is based on linear modal analysis, where the blade's dynamic response is represented as a linear combination of contributions from each of its in-vacuum free vibrational modes. A compressible linearized unsteady potential theory is used to model the flow over the oscillating blades. The two-dimensional unsteady flow is evaluated along several stacked axisymmetric strips along the span of the airfoil. The unsteady pressures at the blade surface are integrated to result in the generalized force acting on the blade due to simple harmonic motions. The unsteady aerodynamic forces are coupled to the blade normal modes in the frequency domain using modal analysis. An iterative eigenvalue problem is solved to determine the stability of the blade when the unsteady aerodynamic forces are included in the analysis. The approach is demonstrated by applying it to a high-energy subsonic turbine blade from a rocket engine turbopump power turbine. The results indicate that this turbine could undergo flutter in an edgewise mode of vibration.

Smith, Todd E.↗

Computation of Steady and Unsteady Laminar Flames: Theory

In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.

Hagstrom, Thomas↗

Quantifying Fractional Ground Cover on the Climate Sensitive High Plains Using AVIRIS and Landsat TM Data

The High Plains is an economically important and climatologically sensitive region of the United States and Canada. The High Plains contain 100,000 sq km of Holocene sand dunes and sand sheets that are currently stabilized by natural vegetation. Droughts and the larger threat of global warming are climate phenomena that could cause depletion of natural vegetation and make this region susceptible to sand dune reactivation. This thesis is part of a larger study that is assessing the effect of climate variability on the natural vegetation that covers the High Plains using Landsat 5 and Landsat 7 data. The question this thesis addresses is how can fractional vegetation cover be mapped with the Landsat instruments using linear spectral mixture analysis and to what accuracy. The method discussed in this thesis made use of a high spatial and spectral resolution sensor called AVIRIS (Airborne Visible and Infrared Imaging Spectrometer) and field measurements to test vegetation mapping in three Landsat 7 sub-scenes. Near-simultaneous AVIRIS images near Ft. Morgan, Colorado and near Logan, New Mexico were acquired on July 10, 1999 and September 30, 1999, respectively. The AVIRIS flights preceded Landsat 7 overpasses by approximately one hour. These data provided the opportunity to test spectral mixture algorithms with AVIRIS and to use these data to constrain the multispectral mixed pixels of Landsat 7. The comparisons of mixture analysis between the two instruments showed that AVIRIS endmembers can be used to unmix Landsat 7 data with good estimates of soil cover, and reasonable estimates of non-photosynthetic vegetation and green vegetation. Landsat 7 derived image endmembers correlate with AVIRIS fractions, but the error is relatively large and does not give a precise estimate of cover.

Warner, Amanda Susan↗

Supersonic Traveling Crossflow Wave Characteristics in Ground and Flight Tests

This paper continues analysis of data acquired in a series of supersonic crossflow transition experiments on a 65° swept- wing model tested in both wind tunnel and flight. Flow visualization data are summarized to provide context for the swept-wing instability mechanisms of interest. Unsteady pressure measurements obtained near the wing model leading edge (LE) are studied using 2nd- and 3rd-order statistical methods to extract traveling crossflow wave characteristics from these data. Comparisons between tunnel and flight data with boundary-layer (BL) Linear Stability Theory (LST) predictions provide insights into BL transition phenomena similarities and differences observed in tests using the same swept-wing model configurations. Comparisons of unsteady pressure data results for wing LE configurations with and without distributed-roughness-element (DRE) patterns installed, suggest traveling crossflow disturbance growth is reduced when this type of flow control is applied. A brief stud y of different background surface roughness levels on the wing LE suggests the potential for another flow control approach to promote laminar flow on swept wings using streamwise-biased surface finishes.

Lewis R Owens↗

Analysis of Aircraft, Radiosonde and Radar Observations in Cirrus Clouds Observed During FIRE II: The Interactions Between Environmental Structure, Turbulence and Cloud Microphysical Properties

Ways to determine the turbulence intensity and the horizontal variability in cirrus clouds have been investigated using FIRE-II aircraft, radiosonde and radar data. Higher turbulence intensities were found within some, but not all, of the neutrally stratified layers. It was also demonstrated that the stability of cirrus layers with high extinction values decrease in time, possibly as a result of radiative destabilization. However, these features could not be directly related to each other in any simple manner. A simple linear relationship was observed between the amount of horizontal variability in the ice water content and its average value. This was also true for the extinction and ice crystal number concentrations. A relationship was also suggested between the variability in cloud depth and the environmental stability across the depth of the cloud layer, which requires further investigation.

Smith, Samantha A.↗

Model Validation for Propulsion - On the TFNS and LES Subgrid Models for a Bluff Body Stabilized Flame

With advances in computational power and availability of distributed computers, the use of even the most complex of turbulent chemical interaction models in combustors and coupled analysis of combustors and turbines is now possible and more and more affordable for realistic geometries. Recent more stringent emission standards have enticed the development of more fuel-efficient and low-emission combustion system for aircraft gas turbine applications. It is known that the NOx emissions tend to increase dramatically with increasing flame temperature. It is well known that the major difficulty, when modeling the turbulence-chemistry interaction, lies in the high non-linearity of the reaction rate expressed in terms of the temperature and species mass fractions. The transport filtered density function (FDF) model and the linear eddy model (LEM), which both use local instantaneous values of the temperature and mass fractions, have been shown to often provide more accurate results of turbulent combustion. In the present, the time-filtered Navier-Stokes (TFNS) approach capable of capturing unsteady flow structures important for turbulent mixing in the combustion chamber and two different subgrid models, LEM-like and EUPDF-like, capable of emulating the major processes occurring in the turbulence-chemistry interaction will be used to perform reacting flow simulations of a selected test case. The selected test case from the Volvo Validation Rig was documented by Sjunnesson.

bluff body↗

Steady vibrations of wing of circular plan form

This paper treats the problem of determining the lift, moment, and induced drag of a thin wing of circular plan form in uniform incompressible flow on the basis of linearized theory. As contrasted to a similar paper by Kinner, in which the acceleration potential method was used, the present paper utilizes the concept of the velocity potential. Calculations of the lift and moment are presented for several deformed shapes. It is shown that considerable deviations exist between the strip theory analysis and the more exact theory. The lift, moment, and induced drag are also determined for a harmonically oscillatory circular plan form wing. As contrasted to a similar paper by Schade, in which the acceleration potential method was used, the present paper utilizes the concept of the velocity potential. Expressions for lift, moment, and induced drag are given and finally specialized to the case of a slowly oscillating circular wing.

FLUTTER - WINGS↗

Handbook of structural stability part III : buckling of curved plates and shells

Available theories and test data on buckling of curved plates and shells are reviewed. For torsion and external-pressure loadings, the test data are correlated in terms of linear buckling theories for both the elastic and inelastic ranges. The cases which exhibit a marked disagreement between linear theory and test data include those of curved plates and cylinders under axial compression, cylinders under bending, and spherical plates under external pressure. These cases have been analyzed by a unified semiempirical approach for both the elastic and inelastic ranges which is satisfactory for analysis and design purposes. The effects of internal pressure on buckling of elements under uniaxial loads are discussed and data on various combined loadings are presented in interaction form. (author)

Gerard, George↗

Attitude control of a spinning flexible spacecraft

The dynamics of rotational motion of a spinning orbiting spacecraft consisting of two rigid bodies connected by a flexible joint and arbitrary number of flexible appendages (two of which are flexible massless booms having masses on their tips) is analyzed. Active attitude control is provided by momentum exchange devices (e.g. control moment gyroscopes) or a mass expulsion system. The linearized equations of motion describing the vehicle are presented, and a large scale digital simulation that has been developed at the Marshall Space Flight Center is presented. A simplified model of the geometrically complex vehicle is selected to make it analytically tractable. The simplified model consists of a single rigid core body with two attached flexible massless booms having tip masses. The states of the vehicle are defined as small perturbations about its steady-state spin. An analysis is performed to determine the domain of stability.

Seltzer, S. M.↗