In the design of solid-propellant rocket motors, the ability to understand and predict the expected behaviour of a given motor under unsteady conditions is important. Research towards predicting, quantifying, and ultimately suppressing undesirable strong transient axial combustion instability symptoms necessitates a comprehensive numerical model for internal ballistic simulation under dynamic flow and combustion conditions. An updated numerical model incorporating recent developments in predicting negative and positive erosive burning, and transient, frequency-dependent combustion response, in conjunction with pressuredependent and acceleration-dependent burning, is applied to the investigation of instabilityrelated behaviour in a small cylindrical-grain motor. Pertinent key factors, like the initial pressure disturbance magnitude and the propellant's net surface heat release, are evaluated with respect to their influence on the production of instability symptoms. Two traditional suppression techniques, axial transitions in grain geometry and inert particle loading, are in turn evaluated with respect to suppressing these axial instability symptoms.
NOMENCLATURE
INTRODUCTIONOver the last fifty years or so, there has been a number of research efforts directed towards understanding the physical mechanisms, or at least the surrounding factors, behind the appearance of symptoms associated with nonlinear axial combustion instability [1-2] in solid rocket motors (SRMs). Traditionally, these symptoms are a sustained axial pressure wave presence in the combustion chamber, sometimes accompanied by a substantial rise in the base chamber pressure [dc shift]. The motivation for these past and present studies was and is of course to bring this better understanding to bear in more precisely suppressing, if not eliminating, these symptoms. Studies of nonlinear axial combustion instability have ranged from numerous experimental test firing series on the one hand [3], and linear/nonlinear acoustic theory modelling on the other [2,[4][5][6]. Largely, the acoustic analysis produces frequency-based standing wave solutions for a given chamber geometry, but typically without some useful quantitative information. On occasion, researchers have employed a numerical modelling approach, to work towards a more comprehensive quantitative understanding of the physics involved [7][8]. The numerical model would typically produce a travelling wave solution to a limit pressure wave amplitude and a corresponding small or larger dc shift in chamber pressure, a time-based result evolving from an initial pulse disturbance introduced into the chamber flow. Available computational power and associated result turnaround times commonly forced some simplifications in the given numerical model. A comprehensive numerical model for simulation of nonlinear dynamic flow and combustion conditions is ultimately essential in the quest for the ability to predict and quantify axial combustion instability symptoms in SRMs. An effective model combines the ef...