This study presents a comprehensive analysis of the frequency response characteristics in a gas generator cycle liquid rocket engine, employing modular decomposition and linearised frequency-domain modeling to simulate dynamic behaviours under forced oscillations. The engine is dissected into key subsystems, including liquid pipelines, turbopump assembly, valves, flow regulation components, thrust chamber, gas generator and pyrotechnic starter, highlighting features such as centrifugal pump pressurisation, staged combustion and cavitation mitigation via venturis. Three oscillation scenarios are examined: supply system responses to thrust chamber pressure disturbances, combustion component responses to fluid disturbances and combustion component responses to pump speed disturbances. Simulations over 0–2000 Hz reveal acoustic-dominated traits in the thrust chamber with oxidiser pathway dominance, low-frequency emphasis in the gas generator driven by fuel disturbances, and heightened instability risks from pump pulsations. Parametric analyses demonstrate that increased pipeline lengths shift resonant frequencies downward, elevated injector pressure drops enhance stability margins by 1.6% with a 20% pressure drop increase, and chamber structural/gas parameter variations erode system stability. These insights, validated against benchmark models, inform strategies for mitigating combustion instability, optimising design parameters, and improving reliability in high-thrust propulsion applications.

Although the problem of locked-in deep stall is well documented over many years, there currently exists no consistent procedure that can guarantee recovery. Past studies have suggested that it might be possible to rock the aircraft in pitch to destabilise the statically stable deep stall trim point, thereby gaining enough momentum to push the nose down. However, the methods used in these studies are either of preliminary or empirical nature and cannot guarantee recovery. In this paper, we use bifurcation analysis to derive a recovery manoeuvre, specifically by assessing the aircraft’s nonlinear frequency response under an elevator forcing. The ensuing nonlinear Bode plot detects unstable (divergent) solutions near resonance that contribute to a successful deep stall recovery. Moreover, the nonlinear resonant frequency is slightly lower than the result obtained using linear analysis, and time simulation shows that relying on the linear result does not lead to a successful recovery. It is also found that at the high angles of attack associated with deep stall, the frequency separation between the short period and phugoid mode is significantly reduced, leading to only one visible peak in the frequency response. This feature is also reflected in the time-domain step response.