Nonlinear rocket motor stability prediction: Limit amplitude, triggering, and mean pressure shift

Gary A. Flandro, Sean R. Fischbach, Joseph Majdalani
2007 Physics of Fluids  
High-amplitude pressure oscillations in solid propellant rocket motor combustion chambers display nonlinear effects including: 1) limit cycle behavior in which the fluctuations may dwell for a considerable period of time near their peak amplitude, 2) elevated mean chamber pressure (DC shift), and 3) a triggering amplitude above which pulsing will cause an apparently stable system to transition to violent oscillations. Along with the obvious undesirable vibrations, these features constitute the
more » ... ost damaging impact of combustion instability on system reliability and structural integrity. The physical mechanisms behind these phenomena and their relationship to motor geometry and physical parameters must, therefore, be fully understood if instability is to be avoided in the design process, or if effective corrective measures must be devised during system development. Predictive algorithms now in use have limited ability to characterize the actual time evolution of the oscillations, and they do not supply the motor designer with information regarding peak amplitudes or the associated critical triggering amplitudes. A pivotal missing element is the ability to predict the mean pressure shift; clearly, the designer requires information regarding the maximum chamber pressure that might be experienced during motor operation. In this paper, a comprehensive nonlinear combustion instability model is described that supplies vital information. The central role played by steep-fronted waves is emphasized. The resulting algorithm provides both detailed physical models of nonlinear instability phenomena and the critically needed predictive capability. In particular, the true origin of the DC shift is revealed. Nomenclature a 0 Mean speed of sound e Oscillatory energy density E Time-averaged oscillatory system energy E m 2 Normalization constant for mode m k m Wave number for axial mode m L Chamber length
doi:10.1063/1.2746042 fatcat:qgv2zlg2yzffbc7y724bebgkdy