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 739 records · Page 41

Applications Study of Aero-Maneuvering Orbit-to-Orbit Shuttle (AMOOS), Executive Summary

Preliminary designs and the supporting analysis for both the Aeromaneuvering Orbit-to-Orbit Shuttle (AMOOS) and the Aeromaneuvering Recovery System (AMRS) are presented. The AMOOS design is shown to yield from twice to almost three times the high energy, round-trip payloads as a purely propulsive vehicle of the same weight. Typically AMOOS can perform a crew rotation mission to equatorial geosynchronous orbit in one space shuttle launch. The flight test program analysis modeled parameters such as maximum dynamic pressure, heating rates, guidance, stability and recovery. Two model flight schedules were developed, one consisting of four flights and the other of two flights. The former is considered a very low risk, high information return program whereas the latter is a minimal cost program consistent with reasonable data returns and chance of success. The AMOOS and AMRS guidance scheme developed using linear regulator theory proved a precise and accurate guidance scheme. Both it and a classical linear system's based scheme were evaluated using 65 simulated trajectories in which the position in the entry corridor and the atmospheric density were varied randomly.

White, J.↗

A genuinely multi-dimensional upwind cell-vertex scheme for the Euler equations

The solution of the two-dimensional Euler equations is based on the two-dimensional linear convection equation and the Euler-equation decomposition developed by Hirsch et al. The scheme is genuinely two-dimensional. At each iteration, the data are locally decomposed into four variables, allowing convection in appropriate directions. This is done via a cell-vertex scheme with a downwind-weighted distribution step. The scheme is conservative, and third-order accurate in space. The derivation and stability analysis of the scheme for the convection equation, and the derivation of the extension to the Euler equations are given. Preconditioning techniques based on local values of the convection speeds are discussed. The scheme for the Euler equations is applied to two channel-flow problems. It is shown to converge rapidly to a solution that agrees well with that of a third-order upwind solver.

Powell, Kenneth G.↗

Stability of a class of interconnected evolution systems

Stability conditions for a class of interconnected systems modeled by linear abstract evolution equations and a memoryless nonlinearity are derived. These conditions are stated in terms of the passivity of each of the subsystems and can be considered as a partial generalization of the hyperstability theorem. A Liapunov function approach is used in the proof without requiring the positive definiteness of the Liapunov function. Application to the robustness analysis of the infinite-dimensional linear quadratic regulator is also discussed.

Wen, John T.↗

Unsteady linearized transonic flow analysis for slender bodies

An unsteady linearized formulation based on Oswatitsch-Keune's parabolic method is developed to analyze transonic flow past oscillating slender bodies. In contrast to the widely used integral transform method, it is shown that all solutions can be derived by a simpler method directly in the physical plane. By various expansion procedures, low-frequency solutions then are derived according to two clearly defined frequency ranges. Adams-Sears' iteration is employed to account for the second-order effects. Stability derivatives are compared with available theories and data. It is found that the derivatives depend more sensitively on thickness than on the reduced frequency. Finally, a critical assessment of the present method is given.

Liu, D. D.↗

Analysis and design of an adaptive multi-loop controlled two winding buck/boost regulator

Small signal low frequency linear average model is derived for a multi-loop controlled two-winding buck/boost converter employing average techniques and the describing function method. The model reveals that a well-designed multi-loop control can provide a second-order zero adaptive to output filter parameter changes due to component tolerances, temperature changes, aging, and the effect of duty cycle modulation. It also can provide stabilization effect by shifting the positive zero to the left-half S-plane. Design guidelines are formulated to optimize regulator-loop dependent characteristics.

Mahmoud, M. F.↗

Variation objective analyses for cyclone studies

The objectives were to: (1) develop an objective analysis technique that will maximize the information content of data available from diverse sources, with particular emphasis on the incorporation of observations from satellites with those from more traditional immersion techniques; and (2) to develop a diagnosis of the state of the synoptic scale atmosphere on a much finer scale over a much broader region than is presently possible to permit studies of the interactions and energy transfers between global, synoptic and regional scale atmospheric processes. The variational objective analysis model consists of the two horizontal momentum equations, the hydrostatic equation, and the integrated continuity equation for a dry hydrostatic atmosphere. Preliminary tests of the model with the SESMAE I data set are underway for 12 GMT 10 April 1979. At this stage of purpose of the analysis is not the diagnosis of atmospheric structures but rather the validation of the model. Model runs for rawinsonde data and with the precision modulus weights set to force most of the adjustment of the wind field to the mass field have produced 90 to 95 percent reductions in the imbalance of the initial data after only 4-cycles through the Euler-Lagrange equations. Sensitivity tests for linear stability of the 11 Euler-Lagrange equations that make up the VASP Model 1 indicate that there will be a lower limit to the scales of motion that can be resolved by this method. Linear stability criteria are violated where there is large horizontal wind shear near the upper tropospheric jet.

Achtemeier, G. L.↗

The Elasto-Plastic Stability of Plates

This article explains results developed from the following research: 'The Stability of Plates and Shells beyond the Elastic Limit.' A significant improvement is found in the derivation of the relations between the stress factors and the strains resulting from the instability of plates and shells. In a strict analysis, the problem reduces to the solution of two simultaneous nonlinear partial differential equations of the fourth order in the deflection and stress function, and in the approximate analysis to a single linear equation of the Bryan type. Solutions are given for the special cases of a rectangular plate buckling into a cylindrical form, and of an arbitrarily shaped plate under uniform compression. These solutions indicate that the accuracy obtained by the approximate method is satisfactory.

Ilyushin, A. A.↗

Stability Metrics for Simulation and Flight-Software Assessment and Monitoring of Adaptive Control Assist Compensators

Due to a need for improved reliability and performance in aerospace systems, there is increased interest in the use of adaptive control or other nonlinear, time-varying control designs in aerospace vehicles. While such techniques are built on Lyapunov stability theory, they lack an accompanying set of metrics for the assessment of stability margins such as the classical gain and phase margins used in linear time-invariant systems. Such metrics must both be physically meaningful and permit the user to draw conclusions in a straightforward fashion. We present in this paper a roadmap to the development of metrics appropriate to nonlinear, time-varying systems. We also present two case studies in which frozen-time gain and phase margins incorrectly predict stability or instability. We then present a multi-resolution analysis approach that permits on-line real-time stability assessment of nonlinear systems.

Hodel, A. S.↗

Utilization of the Deep Space Atomic Clock for Europa Gravitational Tide Recovery

Estimation of Europa's gravitational tide can provide strong evidence of the existence of a subsurface liquid ocean. Due to limited close approach tracking data, a Europa flyby mission suffers strong coupling between the gravity solution quality and tracking data quantity and quality. This work explores utilizing Low Gain Antennas with the Deep Space Atomic Clock (DSAC) to provide abundant high accuracy uplink-only radiometric tracking data. DSAC's performance, expected to exhibit an Allan Deviation of less than 3e-15 at one day, provides long-term stability and accuracy on par with the Deep Space Network ground clocks, enabling one-way radiometric tracking data with accuracy equivalent to that of its two-way counterpart. The feasibility of uplink-only Doppler tracking via the coupling of LGAs and DSAC and the expected Doppler data quality are presented. Violations of the Kalman filter's linearization assumptions when state perturbations are included in the flyby analysis results in poor determination of the Europa gravitational tide parameters. B-plane targeting constraints are statistically determined, and a solution to the linearization issues via pre-flyby approach orbit determination is proposed and demonstrated.

Seubert, Jill↗

Aeroelastic stability and response of horizontal axis wind turbine blades

Coupled flap-lag-torsion equations of motion of an isolated horizontal axis wind turbine (HAWT) blade have been formulated. The analysis neglects blade-tower coupling. The final nonlinear equations have periodic coefficients. A new and convenient method of generating an appropriate time-dependent equilibrium position, required for the stability analysis, has been implemented and found to be computationally efficient. Steady-state response and stability boundaries for an existing (typical) HAWT blade are presented. Such stability boundaries have never been published in the literature. The results show that the isolated blade under study is basically stable. The tower shadow (wake) has a considerable effect on the out-of-plane response but leaves blade stability unchanged. Nonlinear terms can significantly affect linearized stability boundaries; however, they have a negligible effect on response, thus implying that a time-dependent equilibrium position (or steady-state response), based completely on the linear system, is appropriate for the type of HAWT blades under study.

Kottapalli, S. B. R.↗

Static and dynamic characteristics of parallel-grooved seals

Presented is an analytical method to determine static and dynamic characteristics of annular parallel-grooved seals. The governing equations were derived by using the turbulent lubrication theory based on the law of fluid friction. Linear zero- and first-order perturbation equations of the governing equations were developed, and these equations were analytically investigated to obtain the reaction force of the seals. An analysis is presented that calculates the leakage flow rate, the torque loss, and the rotordynamic coefficients for parallel-grooved seals. To demonstrate this analysis, we show the effect of changing number of stages, land and groove width, and inlet swirl on stability of the boiler feed water pump seals. Generally, as the number of stages increased or the grooves became wider, the leakage flow rate and rotor-dynamic coefficients decreased and the torque loss increased.

Iwatsubo, Takuzo↗

Aeroelasticity and structural optimization of composite helicopter rotor blades with swept tips

This report describes the development of an aeroelastic analysis capability for composite helicopter rotor blades with straight and swept tips, and its application to the simulation of helicopter vibration reduction through structural optimization. A new aeroelastic model is developed in this study which is suitable for composite rotor blades with swept tips in hover and in forward flight. The hingeless blade is modeled by beam type finite elements. A single finite element is used to model the swept tip. Arbitrary cross-sectional shape, generally anisotropic material behavior, transverse shears and out-of-plane warping are included in the blade model. The nonlinear equations of motion, derived using Hamilton's principle, are based on a moderate deflection theory. Composite blade cross-sectbnal properties are calculated by a separate linear, two-dimensional cross section analysis. The aerodynamic loads are obtained from quasi-steady, incompressible aerodynamics, based on an implicit formulation. The trim and steady state blade aeroelastic response are solved in a fully coupled manner. In forward flight, where the blade equations of motion are periodic, the coupled trim-aeroelastic response solution is obtained from the harmonic balance method. Subsequently, the periodic system is linearized about the steady state response, and its stability is determined from Floquet theory.

Yuan, K. A.↗

Anisotropic Alfven-ballooning modes in Earth's magnetosphere

We have carried out a theoretical analysis of the stability and parallel structure of coupled shear Alfven and slow magnetosonic waves in Earth's inner magnetopause (i.e., at equatorial distances between about five and ten Earth radii) including effects of finite anisotropic Grad-Shafranov equation yields an approximate self-consistent magnetohydrodynamic (MHD) equilibrium. This MHD equilibrium is used in the numerical solution of a set of eigenmode equations which describe the field line eigenfrequency, linear stability, and parallel eigenmode structure. We call these modes anisotropic Alfven-ballooning modes. The main results are: (1) The field line eigenfrequency can be significantly lowered by finite pressure effects. (2) The parallel mode structure of the transverse wave components is fairly insensitive to changes in the plasma pressure, but the compressional magnetic component can become highly peaked near the magnetic equator as a result of increased pressure, especially when P(sub perpendicular to) is greater than P(sub parallel) (here P(sub perpendicular to) and P(sub parallel) are the perpendicular and parallel plasma pressure). (3) For the isotropic (P(sub parallel) = P(sub perpendicular to) = P) case ballooning instability can occur when the ratio of the plasma presure to the magnetic pressure, beta = P/(B squared/8 pi), exceeds a critical value beta(sup B)(sub 0) is approximately equal to 3.5 at the equator. (4) Compared to the isotropic case the critical beta value is lowered by anisotropy, either due to decreased field line bending stabilization when P(sub parallel) is greater than P(sub perpendicular to) or due to increased ballooning-mirror destabilization when P(sub perpendicular to) is greater than P(sub parallel). (5) We use a beta-delta stability diagram to display the regions of instability with respect to the equatorial values of the parameters bar beta and delta, where bar beta = (1/3)(beta(sub parallel) + 2 beta(sub perpendicular to)) is an average beta value and delta = 1 - P(sub parallel)/P(sub perpendicular to) is a measure of the plasma anisotropy. The diagram is divided into regions corresponding to the firehose, mirror and ballooning instabilities. It appears that observed values of the plasma pressure are below the critical value for the isotropic ballooning instability but it may be possible to approach a ballooning-mirror instability when P(sub perpendicular to)/P(sub parallel) is greater than or approximately 2.

Chan, Anthony A.↗

Nonlinear stability of solar type III radio bursts. II - Application to observations near 1 AU

A set of rate equations including strong turbulence effects and anomalous resistivity are solved using parameters which model several solar type III bursts. Analysis of these bursts has led to quantitative comparisons between several of the observed phenomena and the theory. Through use of an analytic model for the time evolution of the energetic electron exciter, it is found that the exciter distributions observed at 1 AU are unstable to the excitation of the linear bump-in-tail instability, amplifying Langmuir waves above the threshold for the oscillating two-stream instability (OTSI). The OTSI and the attendant anomalous resistivity produce a rapid spectral transfer of Langmuir waves to short wavelengths, out of resonance with the electron exciter. In addition, the various parameters needed to model the bursts are extrapolated inside 1 AU with similar results. Finally, reabsorption of the Langmuir waves by the beam is shown to be unimportant in all cases, even at 0.1 AU.

Goldstein, M. L.↗

Decomposition-aggregation stability analysis

This report presents the development and description of the decomposition aggregation approach to stability investigations of high dimension mathematical models of dynamic systems. The high dimension vector differential equation describing a large dynamic system is decomposed into a number of lower dimension vector differential equations which represent interconnected subsystems. Then a method is described by which the stability properties of each subsystem are aggregated into a single vector Liapunov function, representing the aggregate system model, consisting of subsystem Liapunov functions as components. A linear vector differential inequality is then formed in terms of the vector Liapunov function. The matrix of the model, which reflects the stability properties of the subsystems and the nature of their interconnections, is analyzed to conclude over-all system stability characteristics. The technique is applied in detail to investigate the stability characteristics of a dynamic model of a hypothetical spinning Skylab.

Siljak, D. D.↗

Three dimensional unsteady aerodynamics and aeroelastic response of advanced turboprops

A method for the prediction of steady and unsteady aerodynamic loads and aeroelastic response of advanced turboprops is presented. The aerodynamic analysis uses three dimensional unsteady linearized compressible flow theory to compute the blade pressure distribution. The aeroelastic analysis is based on a normal mode representation of the structure. The method is applicable to both conventional and advanced turbo-prop configurations, provided that blade stall and transonic shock waves are not important factors. Aerodynamic results are presented which validate the model in various limits by comparisons to alternative theories and experimental data. Finally, results of a stability analysis of an advanced turboprop are given, with comparisons to measurements made at NASA Lewis Research Center.

Williams, M. H.↗

Design of a wing shape for study of hypersonic crossflow transition in flight

Computational fluid dynamics methods were used in the design of a wing shape for study of hypersonic crossflow transition in flight. The flight experiment is to be performed on the delta wing of the first stage of a Pegasus launch vehicle as a piggy-back experiment to support boundary-layer stability code development and validation. The design goal is to obtain crossflow-induced transition at 20-40 percent of the chord for a flight Mach number of approximately six. The present paper describes the design and analysis process utilized to obtain desired glove shape. A variety of schemes were used in the design, ranging from simple empirical crossflow correlations to three-dimensional Navier-Stokes codes in conjunction with linear stability/N-factor computations. The sensitivity to various parameters, such as trajectory variations, allowable wing thickness, leading-edge radius and surface temperature, is also discussed.

Godil, A.↗

TPSAS-NF1676L-10863-DND

One important aspect of successful climate monitoring is satellite calibration. There is at least a 25-year of operational geostationary (GEO) record available for climate studies, if properly calibrated. The Global Space-based Inter-Calibration System (GSICS) project is another effort to provide the climate community consistent calibration across multiple platforms. Also, the Clouds and the Earth's Radiant Energy System (CERES) project, provides consistent GEO retrieved cloud properties and derived broadband fluxes across multiple platforms. This year GSICS is focusing on geostationary visible channel calibration. The GSICS strategy is to use multiple approaches, such as cross-calibrating the GEO radiances against MODIS as a reference, deep convective clouds as bright stable targets, deserts, sun-glint, stars, moon and other model based calibration methods. Each of these methods is independent of each other and can assess the stability of the GEO visible instrument with differing uncertainties. This study focuses on the inter-calibration of the GEO radiances against MODIS as a reference by using collocated, and coincident ray-matched radiances. The advantage of the ray-matching technique is that it provides a full dynamic range of radiances to check the linearity of the GEO visible instrument. This presentation will focus on the technique, some error analysis, comparisons with other methods, and an approach to inter-calibrate historical GEOs before MODIS.

D R Doelling↗