1. Introduction
Liquid rocket engines (LRE) have been widely used in spacecraft launch and recovery due to their following advantages: high specific impulse, high thrust, the ability to start repeatedly, variable operating time, adjustable thrust, and repeated usage [
1]. Combustion instability is widely found in rocket engines, aero engines, large air heaters and so on, which will lead to uncontrollable, huge pressure pulsation in the combustors. This pulsation can cause backfire, flameout, and overheating of the combustor wall, which may lead to a series of problems, such as shortening the engine life and even destroying the engine [
2]. Predicting and suppressing the combustion instability of the combustor in LRE is of great significance to the development of aerospace industry.
In order to reduce the occurrence of combustion instability in experiments, theoretical and numerical simulation methods are generally used to predict the combustion instability. Prediction methods of combustion instability are divided into two categories [
3].
The first kind is directly solving the coupled combustion instability system by computational fluid mechanics and compressible solver, which can obtain the frequency of combustion instability. Leonardi et al. [
4] implemented a specific module based on the double time lag model and the coupling of the combustion chamber, and feed line oscillations were investigated by using a complete set of nonlinear equations. This method generally uses large-eddy simulation (LES) to obtain the information of the coupled combustion instability system. Chen et al. [
5] studied self-exciting combustion instability in twin scroll GTMC (gas turbine model combustors) and its interaction with air-fuel mixing by the LES technique. The detached eddy simulation (DES) method was developed in recent years and has been used to the prediction of combustion instability because its advantages in the computational efficiency of Reynolds-averaged Navier–Stokes (RANS) and high accuracy of LES. Yuan et al. [
6] captured self-exciting high-frequency combustion instability in a rocket engine combustion chamber by the two-dimensional DES technique and verified the existence of self-exciting high-frequency combustion instability.
The second kind decouples the combustion system and acoustic system and calculates the flame response and acoustic response with the incoming flow rate disturbance, respectively. Generally, the acoustic response is calculated by a low-order acoustic model, and the flame response is represented as the transfer relationship between the combustion heat release rate and incoming flow rate disturbance, such as the flame transfer function. Urbano et al. [
7] simulated a full-scale combustor by the LES technique and studied the combustion instability of BKD with the acoustic modal identification method. The study showed that the instability of first-order tangential and first-order longitudinal plays an important role in the combustion instability of the chamber. Li et al. [
8] used the LES technique to simulate the combustor and substituted the transfer function into OSCILOS, which is a combustion instability prediction and simulation program, and finally obtained accurate prediction values of thermoacoustic oscillation frequency and amplitude.
Influence factors of combustion instability were studied by experimental and numerical simulation methods to suppress the occurrence of combustion instability. There are many factors that affect combustion instability, such as Reynolds number (Re), combustion chamber geometry [
9], mixing ratio, type of fuel, degree of premix and turbulence intensity [
10]. By affecting the heat release rate and its coupling with acoustic waves, these factors cause different combustion instability phenomena in combustors.
Experimental research studies on combustion instability can be traced back to the 1950s. A lot of research has been carried out on the combustion instability of LRE. The F-1 rocket engine used in the Apollo moon landing program has undergone abundant experiments so as to solve the combustion instability problem. By costing billions of dollars and testing thousands of combustion chamber geometries, the program accumulated plenty of experience in combustion instability research [
11]. Bazarov et al. [
12] found the injector structure influence on combustion instability. The dynamic characteristics of injectors were analyzed theoretically and experimentally. At present, experimental methods are mainly based on two types of advanced optical testing methods: the particle image velocimetry (PIV) and planar laser-induced Fluorescence (PLIF), which can measure the process data and then research on combustion instability. Soller et al. in Germany [
13,
14] studied the oxygen/kerosene coaxial swirl combustor. They observed the longitudinal high-frequency combustion instability phenomenon during the experiment and concluded that the injector structure is the key factor on combustion dynamics. Wang et al. [
15] conducted experimental studies on a single-injector engine. The length of injector recess influenced combustion instability, but the influence of the combustion chamber length on the longitudinal high-frequency combustion instability was more obvious. Xue et al. [
16] carried out experiments on a single-injector rectangular combustion chamber, which showed that there is a relatively optimal recess ratio to make the combustion more stable under supercritical conditions. Bai et al. [
17] conducted an experimental study on a combustor with liquid-centered coaxial swirling injectors, and their research showed that when the self-exciting oscillation frequency of the injectors coupled with the natural acoustic frequency of first order longitudinal mode (1L) of the combustor, the combustor pressure would oscillate at the same frequency. This phenomenon suggested that self-exciting oscillation may be a key factor in the combustion instability of LRE. Armbruster et al. [
18] conducted experiments on combustion instability on the BKD combustion chamber. The research showed that the self-exciting oscillation of injectors makes the combustion chamber transform from the first-order tangential mode, with a larger amplitude, to the first-order longitudinal mode, which is lower frequency and couples with the self-exciting oscillation frequency. Stefan et al. [
19] conducted an in-depth study on the coupling of self-exciting oscillation and combustion instability, which showed that the flow oscillation of injectors has a greater impact on the heat release rate change than sound pressure. Laera et al. [
20] proposed an experimental device that can detect combustion stability and combined the self-exciting oscillation experiment with numerical simulation to study the effect of a single-injector combustion chamber length on the longitudinal combustion instability [
21,
22]. Hardi et al. [
23] took the BKH combustion chamber as an example to study the interaction between the sound field and combustion field, experimentally. The study showed that the oscillation of the sound velocity seriously affects jet breakup when it gets close to the jet velocity.
In numerical simulation, Fureby [
24] simulated combustion instability by the LES technique, and the simulation results showed that vortex shedding is an important factor causing unsteady heat release. Nie [
25], Cheng [
26] and Feng [
27] conducted a lot of simulation research on the injector structure parameters and physical parameters of hypergolic propellant engines, hydrogen oxygen engines and hydrocarbon engines, and concluded that structures such as baffles can suppress high-frequency combustion instability. Qin et al. [
28] found that, compared with the short and thick combustion chamber, the elongated chamber is beneficial to suppress tangential mode combustion instability which is the most harmful factor to the engine. Fang et al. [
29] conducted a series of studies on combustion chamber with liquid oxygen/methane injectors and studied the influence of injector structure parameters on the atomization angle and combustion performance experimentally and numerically. Kraus et al. [
30,
31] predicted combustion instability considering wall heat transfer by the LES technique. The results showed that the heat transfer of wall has a certain impact on combustion field simulation. Wu et al. [
32] simulated the engine with UDMH/N
2O
4 and compared with the experimental results of the hydrogen oxygen engine. The results showed that the two engines have similar variation trends of temperature and Mach number, but the tail flame core temperature of UDMH/N
2O
4 engine is relatively low.
In this paper, the decoupling thermoacoustic system method is adopted. The flame transfer function is obtained by computational fluid dynamics, while the acoustic mode is solved by Helmholtz equation, and the combustion instability of a combustion chamber is predicted. The influence of the flow rate amplitude and flow rate distribution on combustion instability is explored by studying different conditions.
5. Conclusions
In this paper, the combustion instability of 10 different schemes are predicted, and the following conclusions are drawn:
The increase in the disturbance amplitude will increase the absolute value of the modal growth rate. The stable combustion will be more stable when increasing the flow rate amplitude, while the unstable combustion will be more unstable when increasing the flow rate amplitude.
The circumferential unevenness of the flow rate distribution is unfavorable to the stability of the combustion, but the influence of the circumferential unevenness of the flow rate distribution on stability is not monotonic.
For this combustion chamber structure, within a certain range, the smaller the oxidant radial flow rate gradient is, the more favorable it is to the stability of the combustion.