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 631 records · Page 35

Practical Methodology for the Inclusion of Nonlinear Slosh Damping in the Stability Analysis of Liquid-propelled Space Vehicles

One of the challenges of developing flight control systems for liquid-propelled space vehicles is ensuring stability and performance in the presence of parasitic minimally damped slosh dynamics in the liquid propellants. This can be especially difficult when the fundamental frequencies of the slosh motions are in proximity to the frequency used for vehicle control. The challenge is partially alleviated since the energy dissipation and effective damping in the slosh modes increases with amplitude. However, traditional launch vehicle control design methodology is performed with linearized systems using a fixed slosh damping corresponding to a slosh motion amplitude based on heritage values. This papers presents a method for performing the control design and analysis using damping at slosh amplitudes chosen based on the resulting limit cycle amplitude of the vehicle thrust vector system due to a control-slosh interaction under degraded phase and gain margin conditions.

Ottander, John A.↗

Hypersonic Transition Analysis for HIFiRE Experiments

The HIFiRE-1 flight experiment provided a valuable database pertaining to boundary layer transition over a 7-degree half-angle, circular cone model from supersonic to hypersonic Mach numbers, and a range of Reynolds numbers and angles of attack. This paper reports selected findings from the ongoing computational analysis of the measured in-flight transition behavior. Transition during the ascent phase at nearly zero degree angle of attack is dominated by second mode instabilities except in the vicinity of the cone meridian where a roughness element was placed midway along the length of the cone. The growth of first mode instabilities is found to be weak at all trajectory points analyzed from the ascent phase. For times less than approximately 18.5 seconds into the flight, the peak amplification ratio for second mode disturbances is sufficiently small because of the lower Mach numbers at earlier times, so that the transition behavior inferred from the measurements is attributed to an unknown physical mechanism, potentially related to step discontinuities in surface height near the locations of a change in the surface material. Based on the time histories of temperature and/or heat flux at transducer locations within the aft portion of the cone, the onset of transition correlated with a linear N-factor, based on parabolized stability equations, of approximately 13. Due to the large angles of attack during the re-entry phase, crossflow instability may play a significant role in transition. Computations also indicate the presence of pronounced crossflow separation over a significant portion of the trajectory segment that is relevant to transition analysis.

Li, Fei↗

Boundary-layer stability analysis of NLF and LFC experimental data at subsonic and transonic speeds

NASA-Lewis has conducted wind tunnel experiments to ascertain the effectiveness of state-of-the-art in natural laminar flow (NLF) and LFC airfoils for subsonic and transonic speeds, such as the NLF(1)-0414F and the SCLFC(1)-0513F. Attention is given to the effects of Tollmien-Schlichting (TS) and/or crossflow linear mechanisms amplifying small disturbances to generate turbulence. It is found that the incompressible TS transitional n-factors were generally in the 9-12 range, in agreement with earlier correlation studies; the TS instability was the dominant instability mode on a swept-planform LFC airfoil over the entire range of test conditions.

Berry, Scott A.↗

Stability analysis of double-peaked proton distribution functions in the solar wind

Twelve proton velocity distribution functions as measured for different solar wind speeds by the plasma analyzer at the Helios spacecraft are modeled analytically by a superposition of two bi-Maxwellians. The stability of these distributions is investigated by solving the linear dispersion equation for parallel propagation of both left- and right-hand polarized modes. Due to the temperature anisotropy of the main proton component the presence of the resonant left-hand ion cyclotron instability is found to be typical for high-speed streams and is independent of an ion beam. The results indicate that the ratio of the beam density to total particle density and the beam drift velocity are generally close to the threshold values necessary to drive a highly unstable nonresonant right-hand mode. This instability is found for one proton distribution and is a consequence of the beam particles. It is concluded that both the ion cyclotron instability and the ion beam instability act as a regulating mechanism on the shape of solar wind double humped proton distribution functions.

Leubner, Manfred P.↗

Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics

The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations for initial-boundary-value problems of the hyperbolic type. For the numerical experiments, select the one-dimensional inviscid gas-dynamic equations in conservation-law form is selected. A class of implicit schemes based on linear multistep methods for ordinary differential equations is chosen and the use of space or space-time extrapolations as implicit or explicit boundary schemes is emphasized. Some numerical examples with various inflow-outflow conditions highlight the commonly discussed issues: explicit versus implicit boundary schemes, unconditionally stable schemes, and underspecification or overspecification of boundary conditions.

Yee, H. C.↗

Formation Control of the MAXIM L2 Libration Orbit Mission

The Micro-Arcsecond X-ray Imaging Mission (MAXIM), a proposed concept for the Structure and Evolution of the Universe (SEU) Black Hole Imager mission, is designed to make a ten million-fold improvement in X-ray image clarity of celestial objects by providing better than 0.1 micro-arcsecond imaging. Currently the mission architecture comprises 25 spacecraft, 24 as optics modules and one as the detector, which will form sparse sub-apertures of a grazing incidence X-ray interferometer covering the 0.3-10 keV bandpass. This formation must allow for long duration continuous science observations and also for reconfiguration that permits re-pointing of the formation. To achieve these mission goals, the formation is required to cooperatively point at desired targets. Once pointed, the individual elements of the MAXIM formation must remain stable, maintaining their relative positions and attitudes below a critical threshold. These pointing and formation stability requirements impact the control and design of the formation. In this paper, we provide analysis of control efforts that are dependent upon the stability and the configuration and dimensions of the MAXIM formation. We emphasize the utilization of natural motions in the Lagrangian regions to minimize the control efforts and we address continuous control via input feedback linearization (IFL). Results provide control cost, configuration options, and capabilities as guidelines for the development of this complex mission.

Folta, David↗

Whirl Motion of a Seal Test Rig with Squeeze-Film Dampers

This paper presents the experimental behavior and dynamic analysis of a high speed test rig with rolling element bearings mounted in squeeze film oil damper bearings. The test rotor is a double overhung configuration with rolling element ball bearings mounted in uncentered squeeze-film oil dampers. The damper design is similar to that employed with various high-speed aircraft HP gas turbines. The dynamic performance of the test rig with the originally installed dampers with an effective damper length of length 0.23-inch was unacceptable. The design speed of 40,000 RPM could not be safely achieved as nonsynchronous whirling at the overhung seal test disk and high amplitude critical speed response at the drive spline section occurred at 32,000 RPM. In addition to the self excited stability and critical speed problems, it was later seen from FFT data analysis, that a region of supersynchronous dead band whirling occurs between 10,000 to 15,000 RPM which can lead to bearing distress and wear. The system was analyzed using both linear and nonlinear techniques. The extended length damper design resulting from the analysis eliminated the rotor subsynchronous whirling, high amplitude critical speed, and the dead band whirling region allowing the system to achieve a speed of 45,000 RPM. However, nonlinear analysis shows that damper lockup could occur with high rotor unbalance at 33,000 RPM, even with the extended squeeze-film dampers. The control of damper lockup will be addressed in a future paper.

Proctor, Margaret P.↗

Artificial neural network implementation of a near-ideal error prediction controller

A theory has been developed at the University of Virginia which explains the effects of including an ideal predictor in the forward loop of a linear error-sampled system. It has been shown that the presence of this ideal predictor tends to stabilize the class of systems considered. A prediction controller is merely a system which anticipates a signal or part of a signal before it actually occurs. It is understood that an exact prediction controller is physically unrealizable. However, in systems where the input tends to be repetitive or limited, (i.e., not random) near ideal prediction is possible. In order for the controller to act as a stability compensator, the predictor must be designed in a way that allows it to learn the expected error response of the system. In this way, an unstable system will become stable by including the predicted error in the system transfer function. Previous and current prediction controller include pattern recognition developments and fast-time simulation which are applicable to the analysis of linear sampled data type systems. The use of pattern recognition techniques, along with a template matching scheme, has been proposed as one realizable type of near-ideal prediction. Since many, if not most, systems are repeatedly subjected to similar inputs, it was proposed that an adaptive mechanism be used to 'learn' the correct predicted error response. Once the system has learned the response of all the expected inputs, it is necessary only to recognize the type of input with a template matching mechanism and then to use the correct predicted error to drive the system. Suggested here is an alternate approach to the realization of a near-ideal error prediction controller, one designed using Neural Networks. Neural Networks are good at recognizing patterns such as system responses, and the back-propagation architecture makes use of a template matching scheme. In using this type of error prediction, it is assumed that the system error responses be known for a particular input and modeled plant. These responses are used in the error prediction controller. An analysis was done on the general dynamic behavior that results from including a digital error predictor in a control loop and these were compared to those including the near-ideal Neural Network error predictor. This analysis was done for a second and third order system.

Mcvey, Eugene S.↗

A split delta-V technique for drift control of geosynchronous spacecraft

A split-delta-V technique is presented for the east-west stationkeeping of geosynchronous spacecraft using the linearized optimal transfer strategies for near-circular coplanar orbits. The maneuvers are carried out by taking advantage of the dynamics of the spacecraft in order to cancel out the effect of the drift accumulated between successive impulsive delta-V's. A stability analysis is then carried out in order to determine the effect of implementing a series of precanned maneuvers at the end of each drift cycle in an autonomous mode by taking into account initial orbit determination uncertainties and execution errors emanating from each delta-V. It is shown that small tolerance deadbands are violated after only a few cycles and the system is unstable.

Kechichian, F. A.↗

Multivariable control of the Space Shuttle remote manipulator system using H2 and H(infinity) optimization

Three linear controllers are desiged to regulate the end effector of the Space Shuttle Remote Manipulator System (SRMS) operating in Position Hold Mode. In this mode of operation, jet firings of the Orbiter can be treated as disturbances while the controller tries to keep the end effector stationary in an orbiter-fixed reference frame. The three design techniques used include: the Linear Quadratic Regulator (LQR), H2 optimization, and H-infinity optimization. The nonlinear SRMS is linearized by modelling the effects of the significant nonlinearities as uncertain parameters. Each regulator design is evaluated for robust stability in light of the parametric uncertanties using both the small gain theorem with an H-infinity norm and the less conservative micro-analysis test. All three regulator designs offer significant improvement over the current system on the nominal plant. Unfortunately, even after dropping performance requirements and designing exclusively for robust stability, robust stability cannot be achieved. The SRMS suffers from lightly damped poles with real parametric uncertainties. Such a system renders the micro-analysis test, which allows for complex peturbations, too conservative.

Prakash, OM, II↗

Coronal closed structures. IV - Hydrodynamical stability and response to heating perturbations

The response of magnetically confined atmospheres to perturbations in the temperature and density distribution, and the local heating rate by means of a one-dimensional time-dependent hydrodynamical code, which incorporates the full energy, momentum and mass conservation equations is studied. These studies extend the linear instability analysis of Habbal and Rosner (1979) into the finite-amplitude regime, and generalize the confined atmosphere models of Serio et al., to the time-dependent domain. The results show that closed coronal atmospheres are stable against finite-amplitude perturbations if the chromospheric response is taken into account; and observed correlated increases in coronal density and temperature can only be achieved under quiescent conditions by increasing the heat deposition rate relatively more in the chromosphere than in the corona.

Peres, G.↗

Conic Sector Analysis of Hybrid Control Systems

A hybrid control system contains an analog plant and a hybrid (or sampled-data) compensator. In this thesis a new conic sector is determined which is constructive and can be used to: (1) determine closed loop stability, (2) analyze robustness with respect to modelling uncertainties, (3) analyze steady state response to commands, and (4) select the sample rate. The use of conic sectors allows the designer to treat hybrid control systems as though they were analog control systems. The center of the conic sector can be used as a rigorous linear time invariant approximation of the hybrid control system, and the radius places a bound on the errors of this approximation. The hybrid feedback system can be multivariable, and the sampler is assumed to be synchronous. Algorithms to compute the conic sector are presented. Several examples demonstrate how the conic sector analysis techniques are applied. Extensions to single loop multirate hybrid feedback systems are presented. Further extensions are proposed for multiloop multirate hybrid feedback system and for single rate systems with asynchronous sampling.

Thompson, P. M.↗

Venus gravity - A high-resolution map

The Doppler data from the radio tracking of the Pioneer Venus Orbiter (PVO) have been used in a two-stage analysis to develop a high-resolution map of the gravitational potential of Venus, represented by a central mass and a surface mass density. The two-stage procedure invokes a Kalman filter-smoother to determine the orbit of the spacecraft, and a stabilized linear inverter to estimate the surface mass density. The resultant gravity map is highly correlated with the topographic map derived from the PVO radar altimeter data. However, the magnitudes of the gravity variations are smaller than would be expected if the topography were uncompensated, indicating that at least partial compensation has taken place.

Reasenberg, R. D.↗

A note on the effects of viscosity on the stability of a trailing-line vortex

The linear stability of the Batchelor (1964) vortex is examined with emphasis on new viscous modes recently found numerically by Khorrami (1991). Unlike the previously reported inviscid modes of instability, these modes are destabilized by viscosity and exhibit small growth rates at large Reynolds numbers. The analysis presented here uses a combination of asymptotic and numerical techniques. The results confirm the existence of the additional modes of instability due to viscosity.

Duck, Peter W.↗

Multidisciplinary aeroelastic analysis of a generic hypersonic vehicle

This paper presents details of a flutter and stability analysis of aerospace structures such as hypersonic vehicles. Both structural and aerodynamic domains are discretized by the common finite element technique. A vibration analysis is first performed by the STARS code employing a block Lanczos solution scheme. This is followed by the generation of a linear aerodynamic grid for subsequent linear flutter analysis within subsonic and supersonic regimes of the flight envelope; the doublet lattice and constant pressure techniques are employed to generate the unsteady aerodynamic forces. Flutter analysis is then performed for several representative flight points. The nonlinear flutter solution is effected by first implementing a CFD solution of the entire vehicle. Thus, a 3-D unstructured grid for the entire flow domain is generated by a moving front technique. A finite element Euler solution is then implemented employing a quasi-implicit as well as an explicit solution scheme. A novel multidisciplinary analysis is next effected that employs modal and aerodynamic data to yield aerodynamic damping characteristics. Such analyses are performed for a number of flight points to yield a large set of pertinent data that define flight flutter characteristics of the vehicle. This paper outlines the finite-element-based integrated analysis procedures in detail, which is followed by the results of numerical analyses of flight flutter simulation.

Gupta, K. K.↗

A linear analysis of the transition curve for the baroclinic annulus

A primitive equation linear model is shown to be capable of quantitatively predicting the neutral curve for a baroclinic annulus. Flows unstable to small nonaxisymmetric perturbations are assumed to lie on the nonaxisymmetric side of a transition curve, while the stable flows are on the axisymmetric side. The stability in the presence of perturbations is tested with linearized hydrostatic primitive equations. Attention is given to the fastest growing eigenmodes, demonstrating that the linear and nonlinear wave have similar structures at a particular point in the nonaxisymmetric regime. All the waves are similar to the Eddy wave except at the extreme to the transition curve, where the waves maintain little structure with height. The eddy kinetic energy is baroclinic with a potential energy suppresed by the large static stability of the basic state. Eddy dissipation is noted to be significant near all boundaries.

Miller, T. L.↗

Non-Linear Analysis of the NASA Super Pressure Balloons: Whole Flight Simulations

Tensys have a long-established background in the shape generation and load analysis of architectural stressed membrane structures. Founded upon their inTENS finite element analysis suite, these activities have broadened to encompass 'lighter than air' structures such as aerostats, airships, hybrid air vehicles and stratospheric balloons. Tensys have acted as consultants to the NASA Super Pressure Balloon (SPB) Program since 2004. Previous papers have focussed upon the application of inTENS to the overall structural and stability analysis of pumpkin type balloons as used in the SPB Program. Particular emphasis has been placed upon the ability to study both stress and stability at all stages of a flight. As the program has developed, increasing modelling fidelity has been introduced as design refinement has moved emphasis from the overall shape to the performance of individual details. Examples include the introduction of a contact capability into inTENS to better represent the separate cap layer of film. Localised investigations have considered the consequence of local geometric anomalies introduced during fabrication and the effects of debonding of the PBO tendon within its sleeve. Analysis to date has used a material model for the polyethylene shell film developed by Dr Rand of Winzen Engineering, supported by a program of fine resolution material tests at the Balloon Research and Development Laboratory (BRDL) at the GSFC Wallops Island Facility. Based upon work first reported by Schapery, this model has been incorporated into so-called 'snapshot' analyses within inTENS. For a given elapsed time point and temperature distribution, non-linear material properties are iteratively updated for the current stress state in each individual element until changes in those material properties are insignificant. This process has to be iterative as the film stresses and material properties are interdependent. Attempts to derive a time-stepping incremental viscoelastic capability using the Schapery Rand model encountered problems when dealing with the varying temperatures associated with balloon deployment, pressurisation and diurnal behaviour. An alternative approach has been developed by Pellegrino et al at Caltech, again with support from the NASA Balloon Program Office. This is a large strain non-linear viscoelastic model that includes film out-of-plane mechanical and thermal effects. Working from a unified base, this model can be utilised for time-stepping analyses in both a modulus or compliance mode, or a combination of both. This new capability is intended to enhance the current NASA balloon design process. This paper presents the implementation of the Caltech model into a specialist finite element suite, inTENS and its application to whole flight simulations.

Wakefield, David↗

Aeroelastic stability of complete rotors with application to a teetering rotor in forward flight

The derivation of a set of nonlinear coupled flap-lag-torsion equations of motion for moderately large deflections of an elastic, two bladed teetering helicopter rotor in forward flight is concisely outlined. The following degrees of freedom are included in the mathematical model: rigid body flapping, rigid body lead lag, elastic bending in flap and lead-lag blade root torsion and shaft torsion. Quasi-steady aerodynamic loads are considered and the effects of reversed flow are included. The aeroelastic stability of the complete rotor is investigated using a linearized system of equations of motion. The equilibrium position about which the equations are linearized is obtained by considering the trim state of the helicopter, in true or simulated forward flight conditions. The sensitivity of the aeroelastic stability boundaries to interblade structural and mechanical coupling is illustrated by comparing the complete rotor stability boundaries with those obtained from a single blade analysis for a number of hover and forward flight cases.

Shamie, J.↗