SEARCH · Search NASA
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.
Stable limit cycles and triggering limits of first radial mode in nonlinearly unstable liquid propellant rocket motors, using wave equation and Galerkin method
The Cassini spacecraft dynamics-related telemetry during long Main Engine (ME) burns has indicated the presence of stable limit cycles between 0.03-0.04 Hz frequencies. These stable limit cycles cause the spacecraft to possess non-zero oscillating rates for extended periods of time. This indicates that the linear ME guidance control system does not model the complete dynamics of the spacecraft. In this study, we propose that the observed limit cycles in the spacecraft dynamics telemetry appear from a stable interaction between the unmodeled nonlinear elements in the ME guidance control system. Many nonlinearities in the control system emerge from translating the linear engine gimbal actuator (EGA) motion into a spacecraft rotation. One such nonlinearity comes from the gear backlash in the EGA system, which is the focus of this paper. The limit cycle characteristics and behavior can be predicted by modeling this gear backlash nonlinear element via a describing function and studying the interaction of this describing function with the overall dynamics of the spacecraft. The linear ME guidance controller and gear backlash nonlinearity are modeled analytically. The frequency, magnitude, and nature of the limit cycle are obtained from the frequency response of the ME guidance controller and nonlinear element. In addition, the ME guidance controller along with the nonlinearity is simulated. The simulation response contains a limit cycle with similar characterstics as predicted analytically: 0.03-0.04 Hz frequency and stable, sustained oscillations. The analytical and simulated limit cycle responses are compared to the flight telemetry for long burns such as the Saturn Orbit Insertion and Main Engine Orbit Trim Maneuvers. The analytical and simulated limit cycle characteristics compare well with the actual observed limit cycles in the flight telemetry. Both have frequencies between 0.03-0.04 Hz and stable oscillations. This work shows that the stable limit cycles occur due to the interaction between the unmodeled nonlinear elements and linear ME guidance controller.
Attention is given to the physical concepts and mathematical techniques useful in the analysis of the stabilizing effect of friction on aerodynamically unstable rotor stages. Results are presented for three-, four-, and five-bladed disks. In the present multidegree-of-freedom model of an aerodynamically unstable rotor stage, a harmonic steady state solution due to the friction dampers may be either a stability limit, a stable cycle limit, or neither. A criterion is established in the form of an energy function which determines whether the solution is a stability limit. In the event that the initial displacement and velocity exeed those associated with the steady state solution corresponding to a stability limit, the reponse becomes unbounded.
Compressor stall phenomena are analyzed from the point of view of nonlinear control theory, based on bifurcation-catastrophe techniques. This new approach appears promising and offers insight into such well-known compressor instability problems as surge and rotating stall and suggests strategies for recovery. Three interlocking dynamic nonlinear state space models are developed. It is shown that the problem of rotating stall can be viewed as an induced bifurcation of solution of the unstalled model. Hysteresis effects are shown to exist in the stall/recovery process. Surge cycles are observed for some critical parameter values. The oscillatory behavior is seen to be due to development of limit cycles, generated by Hopf bifurcation of solutions. More specifically, it is observed that at certain critical values of parameters, a family of stable limit cycles with growning and then diminishing amplitudes is generated, then giving rise to an unstable family of limit cycles. This unstable family in turn bifurcates into other unstable families. To further illustrate the utility of the methodology, some partial computation of domains is carried out, and parameter sensitivity analysis is performed.
Pulsed oscillations in solid rocket motors are investigated with emphasis on nonlinear combustion response. The study employs a wave equation governing the unsteady motions in a two-phase flow, and a solution technique based on spatial- and time-averaging. A wide class of combustion response functions is studied to second-order in fluctuation amplitude to determine if, when, and how triggered instabilities arise. Conditions for triggering are derived from analysis of limit cycles, and regions of triggering are found in parametric space. Based on the behavior of model dynamical systems, introduction of linear cross-coupling and quadratic self-coupling among the acoustic modes appears to be the manner in which the nonlinear combustion response produces triggering to a stable limit cycle. Regions of initial conditions corresponding to stable pulses were found, suggesting that stability depends on initial phase angle and harmonic content, as well as the composite amplitude, of the pulse.
Non-linear dynamics of multiple infernal Alfvén eigenmodes—a subset of global Alfvén eigenmodes in tokamak plasmas with extended low-shear central core [Marchenko et al., Phys. Plasmas 16, 092502 (2009)]—is studied. The analysis is carried out for a mode triplet with toroidal mode-numbers n=1, 2, 3. It was assumed that the n = 1 mode was linearly unstable due to precession resonance with trapped fast ions, whereas the other modes were linearly damped. The modes were coupled due to a non-linearity in a bounce-averaged drift kinetic equation for the distribution function of fast ions. Nonlinear equations for the mode amplitudes and phases are derived and solved numerically. It is found that the temporal evolution of the amplitudes and the phase (responsible for the frequency chirping) of the modes exhibit Hopf bifurcations to stable limit cycles. This can explain a synchronous cyclic destabilization of multiple modes in Alfvén avalanches (sudden growth of amplitudes of the mode cluster with different n and approximately equal frequency spacing) in NSTX and bursting modes in MAST—events, which resulted in enhanced loss of fast ions.
Most published work relating to dynamically loaded journal bearings are directed to determining the minimum film thickness from the predicted journal trajectories. These do not give any information about the subsynchronous whirl stability of journal bearing systems since they do not consider the equations of motion. It is, however, necessary to know whether the bearing system operation is stable or not under such an operating condition. The stability characteristics of the system are analyzed. A linearized perturbation theory about the equilibrium point can predict the threshold of stability; however it does not indicate postwhirl orbit detail. The linearized method may indicate that a bearing is unstable for a given operating condition whereas the nonlinear analysis may indicate that it forms a stable limit cycle. For this reason, a nonlinear transient analysis of a rigid rotor supported on oil journal bearings under: (1) a unidirectional constant load, (2) a unidirectional periodic load, and (3) variable rotating load are performed. The hydrodynamic forces are calculated after solving the time-dependent Reynolds equation by a finite difference method with a successive overrelaxation scheme. Using these forces, equations of motion are solved by the fourth-order Runge-Kutta method to predict the transient behavior of the rotor. With the aid of a high-speed digital computer and graphics, the journal trajectories are obtained for several different operating conditions.
The general equations describing the nonlinear fluttering oscillations of shallow, curved, heated orthotropic panels have been derived. The formulation takes into account the location of the panel on the surface of a generic hypersonic vehicle, when calculating the aerodynamic loads. It is also shown that third order piston theory produces unsteady aerodynamic loading which is in close agreement with that based upon direct solution of the Euler equations. Results, for simply supported panels, are obtained using Galerkin's method combined with direct numerical integration in time to compute stable limit cycle amplitudes. These results illustrate the sensitivity of the aeroelastic behavior to the unsteady aerodynamic assumptions, temperature, orthotropicity and flow orientation.
Oscillations of supported liquid drops are the subject of wide scientific interest, with applications in areas as diverse as liquid-liquid extraction, synthesis of ceramic powders, growing of pure crystals in low gravity, and measurement of dynamic surface tension. In this research, axisymmetric forced oscillations of arbitrary amplitude of viscous liquid drops of fixed volume which are pendant from or sessile on a rod with a fixed or moving contact line and surrounded by an inviscid ambient gas are induced by moving the rod in the vertical direction sinusiodally in time. In this paper, a preliminary report is made on the computational analysis of the oscillations of supported drops that have 'clean' interfaces and whose contact lines remain fixed throughout their motions. The relative importance of forcing to damping can be increased by either increasing the amplitude of rod motion A or Reynolds number Re. It is shown that as the ratio of forcing to damping rises, for drops starting from an initial rest state a sharp increase in deformation can occur when they are forced to oscillate in the vicinity of their resonance frequencies, indicating the incipience of hysteresis. However, it is also shown that the existence of a second stable limit cycle and the occurrence of hysteresis can be observed if the drop is subjected to a so-called frequency sweep, where the forcing frequency is first increased and then decreased over a suitable range. Because the change in drop deformation response is abrupt in the vicinity of the forcing frequencies where hysteresis occurs, it should be possible to exploit the phenomenon to accurately measure the viscosity and surface tension of the drop liquid.
To facilitate analysis, ALS systems are often assumed to be linear and time invariant, but they usually have important nonlinear and dynamic aspects. Nonlinear dynamic behavior can be caused by time varying inputs, changes in system parameters, nonlinear system functions, closed loop feedback delays, and limits on buffer storage or processing rates. Dynamic models are usually cataloged according to the number of state variables. The simplest dynamic models are linear, using only integration, multiplication, addition, and subtraction of the state variables. A general linear model with only two state variables can produce all the possible dynamic behavior of linear systems with many state variables, including stability, oscillation, or exponential growth and decay. Linear systems can be described using mathematical analysis. Nonlinear dynamics can be fully explored only by computer simulations of models. Unexpected behavior is produced by simple models having only two or three state variables with simple mathematical relations between them. Closed loop feedback delays are a major source of system instability. Exceeding limits on buffer storage or processing rates forces systems to change operating mode. Different equilibrium points may be reached from different initial conditions. Instead of one stable equilibrium point, the system may have several equilibrium points, oscillate at different frequencies, or even behave chaotically, depending on the system inputs and initial conditions. The frequency spectrum of an output oscillation may contain harmonics and the sums and differences of input frequencies, but it may also contain a stable limit cycle oscillation not related to input frequencies. We must investigate the nonlinear dynamic aspects of advanced life support systems to understand and counter undesirable behavior.
The Common Extensive Cryogenic Engine program demonstrated the operation of a deep throttling engine design. The program, spanning five years from August 2005 to July 2010, funded testing through four separate engine demonstration test series. Along with successful completion of multiple objectives, a discrete response of approximately 4000 Hz was discovered and explored throughout the program. The typical low-amplitude acoustic response was evident in the chamber measurement through almost every operating condition; however, at certain off-nominal operating conditions, the response became discrete with higher amplitude. This paper summarizes the data reduction, characterization, and analysis of the 4,000 Hz response for the entire program duration, using the large amount of data collected. Upon first encountering the response, new objectives and instrumentation were incorporated in future test series to specifically collect 4,000 Hz data. The 4,000 Hz response was identified as being related to the first tangential acoustic mode by means of frequency estimation and spatial decomposition. The latter approach showed that the effective node line of the mode was aligned with the manifold propellant inlets with standing waves and quasi-standing waves present at various times. Contour maps that contain instantaneous frequency and amplitude trackings of the response were generated as a significant improvement to historical manual approaches of data reduction presentation. Signal analysis and dynamic data reduction also uncovered several other features of the response including a stable limit cycle, the progressive engagement of subsequent harmonics, the U-shaped time history, an intermittent response near the test-based neutral stability region, other acoustic modes, and indications of modulation with a separate subsynchronous response. Although no engine damage related to the acoustic mode was noted, the peak-to-peak fluctuating pressure amplitude achieved 12.1% of the mean chamber pressure at its highest. The identification of this response in terms of an instability is also discussed.
The nonlinear dynamics of propellant sloshing during orbital ascent are usually neglected in the flight control analysis of large boost vehicles under the assumption that the viscous damping of the fluid is sufficient to suppress nonlinear phenomena and confine the fluid to small, planar free surface displacements. In this case, the sloshing dynamics can be modeled using a spring-mass-damper or linearized pendulum mechanical analog. However, large, smooth-wall tanks without significant internal hardware or ring baffles are still susceptible to nonlinear effects. In particular, rotary sloshing can present a risk to flight control as it involves the formation of a stable limit cycle which can lead to undesirable roll coupling. The underlying phenomena of jump resonance does not manifest in linear models, but can be reproduced using a nonlinear spherical pendulum or the Bauer paraboloid model developed during the Apollo/Saturn program. In this paper, a detailed analysis of the rotary sloshing dynamics of these mechanical analogs is presented, and discussed in the context of flight control stability. High-fidelity simulations of a representative boost vehicle are used to verify the semi-analytical predictions of the nonlinear dynamic response.
Compressor stall phenomena is analyzed from nonlinear control theory viewpoint, based on bifurcation-catastrophe techniques. This new approach appears promising and offers insight into such well known compressor instability problems as surge and rotating stall; furthermore it suggests strategies for recovery from stall. Three interlocking dynamic nonlinear state space models are developed. It is shown that the problem of rotating stall can be viewed as an (induced) bifurcation of solution of the unstalled model. Hysteresis effect is shown to exist in the stall/recovery process. Surge cycles are observed to develop for some critical parameter values. It is shown that the oscillatory behavior is due to development of limit cycles, generated by Hopf bifurcation of solutions. Both stable and unstable limit cycles are observed. To further illustrate the usefulness of the methodology some partial computation of domains of attraction of equilibria is carried out, and parameter sensitivity analysis is performed.
Aircraft that do not possess inherent (aerodynamic) stability must rely on closed-loop control systems for stable operation. Because there are limits on the deflections of an aircraft's control surfaces, the region of stable operation also is bounded. These boundaries are investigated for a lateral-directional example in which vertical fin size is inadequate to provide directional stability and where aileron and rudder deflections are subject to saturation. Fourth-order models are used in this study, with flight control logic based on minimum-control-energy linear-quadratic-regulatory theory. It is found that the stability boundaries can be described by unstable limit cycles surrounding stable equilibrium points. Variations in regions of stability with gain levels and command inputs are illustrated. Current results suggest guidelines for permissible limits on the open-loop instability of an aircraft's lateral-directional modes.
In order to phase lock the flow at the desired shedding cycle, particularly at Phi,best, We designed a feedback compensator. (Even though the open-loop forcing at Wf below Wn can lead to phase-locked limit cycles with a high average lift,) This feedback controller resulted in the phase-locked limit cycles that the open-loop control could not achieve for alpha=30 and 40 Particularly for alpha=40, the feedback was able to stabilize the limit cycle that was not stable with any of the open-loop periodic forcing. This results in stable phase-locked limit cycles for a larger range of forcing frequencies than the open-loop control. Also, it was shown that the feedback achieved the high-lift unsteady flow states that open-loop control could not sustain even after the states have been achieved for a long period of time.
A new analytical technique for the solution of nonlinear longitudinal combustion instability problems in rocket combustors is developed. Using relatively little computation time, this technique is capable of predicting the transient and limit cycle behavior of the combustion instability oscillations as well as the disturbance amplitude required to trigger an instability in a linearly stable motor. The limit cycle waveforms are found to exhibit shock wave characteristics for most unstable engine operating conditions. It is shown that the characteristics of the resulting instability are independent of the nature of the initial disturbance and they depend solely upon the engine operating conditions and the characteristics of the unsteady combustion process.
A number of issues related to the flutter and postflutter of two-dimensional supersonic lifting surfaces are addressed. Among them there are the 1) investigation of the implications of the nonlinear unsteady aerodynamics and structural nonlinearities on the stable/unstable character of the limit cycle and 2) study of the implications of the incorporation of a control capability on both the flutter boundary and the postflutter behavior. To this end, a powerful methodology based on the Lyapunov first quantity is implemented. Such a treatment of the problem enables one to get a better understanding of the various factors involved in the nonlinear aeroelastic problem, including the stable and unstable limit cycle. In addition, it constitutes a first step toward a more general investigation of nonlinear aeroelastic phenomena of three-dimensional lifting surfaces.
Conditions are predicted under which limit cycles may be caused by the friction torque of the control moment gyroscope (CMG) output gimbal used in the Large Space Telescope (LST). The parameter plane analysis technique is described, and three separate models of the CMG output friction nonlinearity are used in three separate analyses. The results suggest that the Dahl model (included in model 3) is the best available friction model for use in dynamic analysis. The analysis and analog simulation indicate the absence of limit cycle behavior due to the CMG output gimbal friction nonlinearity for presently estimated LST numerical parameters and the occurrence of one stable and one unstable limit cycle if substantially larger numerical values are allowed for the CMG gimbal pivot friction parameter and/or the running friction torque.