Abstract
A hybrid powder-solid ramjet (HPSR) that combined the advantages of a solid rocket ramjet (SRJ) and a powder ramjet engine was investigated in this research. To improve the combustion efficiency and optimize the inner wall thermal protection of the afterburner, the effects of the air inlet angle on the combustion and wall ablation environment were studied. The standard - model, the eddy-dissipation model (EDM), and the boron particle ignition and combustion model were adopted to simulate the two-phase flow in the afterburners with different air inlet angles (45°, 60°, 75°, and 90°). The results showed that the global flow field and the distribution of the vortexes in the afterburner that had a significant influence on the ablation environment of the inner wall and the combustion efficiency were determined by the impact effect and the squeezing effect of the ram air on the primary fuel gas, which was affected by the air inlet angle. As the air inlet angle increased, the total combustion efficiency of the four cases first increased and then decreased, reaching 80.38%, 81.64%, 84.34%, and 83.26% for angles of 45°, 60°, 75°, and 90°, respectively. At the same time, the inner wall ablation became more severe because both the erosion effect of the condensed phase particles and the gas-flow scouring effect were enhanced, and a large temperature gradient was generated on the inner wall. The study results can provide a reference for designing the air inlet angle of an HPSR.
1. Introduction
Ducted rockets with important military significance have been developed for several decades [1–4]. To improve the comprehensive performance of ducted rockets, researchers usually apply a high-energy propellant with metallic particles or optimize the structure of a ramjet engine [5–9].
High-energy density fuel formulations containing polymer with embedded metal particles, including aluminum, magnesium, and boron, have been adopted in ducted rockets for decades [10–14]. Boron has a higher volume and calorific value than other metals, and its combustion products are in the gas phase in a high temperature environment, which will not cause thrust loss, so boron has been considered the most attractive fuel for ducted rockets. The combustion of boron particles in an afterburner is a well-explored subject [15–19]. Odawara et al. [15] researched the ignition delay time of boron particles in a combustion chamber and analyzed the ignition and combustion characteristics of gas generators containing boron. Natan and Gany [16] investigated the burning behavior of a single boron particle and found that there was a particle size and injection velocity range conducive to maintaining the combustion. Gany [17] studied the thermodynamic conditions associated with highly boron-loaded ramjet combustors and predicted the conditions for blocking reactions within a certain range of room pressures and temperatures. Liu et al. [18] developed a new experimental system to study the combustion of boron particles in a ducted rocket, and they provided some suggestions to optimize the structure of an afterburner to obtain a higher completion rate of boron combustion. Hashim et al. [19] used a low-cost screening instrument called an opposed flow burner (OFB) to investigate the effects of the size of boron particles on combustion by testing pure hydroxyl-terminated polybutadiene (HTPB) and HTPB loaded with boron nanoparticles at various concentrations (5%–40% by wt.). In short, previous scholars [15–19] have performed comprehensive studies of boron particles and given detailed descriptions of the ignition and combustion process of boron particles in ramjets. Therefore, in the current work, HTPB propellant containing boron was adopted.
However, it is difficult to adjust the content of boron particles and the thrust of a solid rocket ramjet. Therefore, a new ramjet, called a hybrid powder-solid ramjet (HPSR), was investigated. The HPSR combined the advantages of a powder ramjet and the solid ramjet to adjust the powder flow rate and achieve the self-sustaining burn of metal particles easily [20]. In addition, many studies have focused on optimizing the structure of a ramjet to improve its comprehensive performance [20–26]. Zvuloni et al. [21] investigated the effects of the geometry and size on the combustion in a solid fuel ramjet (SFRJ) through experiments and found that the mean average regression rate decreased with the increasing gas port diameter. Vigot et al. [22] studied a ducted rocket with a double gas inlet with a theoretical analysis, demonstrating that the double gas inlet divided the inner flow field into two parts with a significant difference in the temperature and oxygen concentrations. Pein and Vinnemeier [23] carried out an experiment to study an SFRJ with a swirler and showed that the swirl flow strongly increased the combustion efficiency, until a maximum was reached at about 20%. Furthermore, Liu et al. [24] conducted a study on a turbocharged solid ramjet with an afterburner that combined the head region inlet and the lateral inlet and found that the radial and tangential momentum exchange of the air-fuel streams was enhanced by driving the compressed air to deflect inward. Musa et al. [25] performed an experimental investigation on an SFRJ with a new design swirler. It was revealed that increasing the swirl number and the combustor length increased the thrust and the characteristic velocity but decreased the air-fuel ratio and the specific impulse, while the increase of the port diameter was diametrically opposite. Xu et al. [20] investigated the effect of the tangential swirl air inlet angle in an HPSR and found that the ignition distance was decreased, and the combustion efficiency was increased with the increase of the angle. Li et al. [26] investigated an SFRJ that included a primary gas generator with a swirler and found that increasing the swirl number could enhance the combustion heat release rate, regression rate, and combustor average temperature.
In summary, although many studies on the structural factors of an afterburner have been carried out, few studies have focused on the effect of the air inlet angle [20–26]. Additionally, most studies have focused on combustion efficiency, while few have focused on the inner wall ablation of the combustion chamber. However, when an aircraft with a ramjet is performing a high-speed maneuverer, it has to withstand high overloads and a two-phase flow carrying a high concentration of condensed phase particles forms in the combustion chamber, worsening the working environment of the thermal insulation layer. In severe cases, this leads to the failure of the thermal insulation protection and even causes the ramjet to burn through and explode [27, 28]. Therefore, to ensure the reliable thermal protection of a ramjet, it is important to investigate the inner wall ablation environment of the afterburner.
In this work, an HPSR afterburner with symmetrical air inlets was established, and the effects of various air inlet angles (45°, 60°, 75°, and 90°) on the combustion efficiency and inner wall ablation environment were analyzed. The primary focus was on the following aspects: (i) the effects of the air inlet angle on the inner flow field of the afterburners, (ii) the effects of the air inlet angle on the combustion efficiency, and (iii) the effects of the air inlet angle on the inner wall ablation environment. Due to limited research on the effects of various air inlet angles on the combustion efficiency and the inner wall ablation environment, this paper can provide a reference to improve the comprehensive performance of an HPSR from the perspective of designing the air inlet angle.
2. Engine Model
Figure 1 shows the structure of the afterburner of the HPSR. The total length of the afterburner was 900 mm. The head length was 250 mm. The main body radius was 88 mm. The primary gas generator nozzle radius was 20 mm. The radius of the nozzle throat was 65 mm. The angle of the air inlet was denoted using . To study the comprehensive effect of the air inlet angle on the afterburner, four cases were selected as case 1, case 2, case 3, and case 4, with inlet angles of 45°, 60°, 75°, and 90°, respectively.
To simplify the model, the following assumptions were made for the internal flow field of the afterburner. (1)The gas jet in the afterburner had quasisteady flow, regardless of the heat exchange between the gas and the wall of the afterburner(2)The gas in the afterburner was an ideal gas that obeyed the ideal gas state equation
3. Numerical Setup
3.1. Two-Phase Flow Governing Equation
According to the above assumptions, the Reynolds-Averaged Navier-Stokes equation was adopted as the gas-phase governing equation. The general forms of the continuum equation, the momentum equation, the energy equation, and the composition equation could be written as where is time, is the direction, is the direction, and is the direction. is the conservative vector, is the vector of convective flux, is the vector of viscous flux, and is the vector of source terms.
The particle phase governing equation employed the Lagrangian particle trajectory model [29, 30]. The force balance equation could be written as:
The track equation of the particle in the direction could be written as: where is the instantaneous velocity of particles, is the component of gravitational acceleration in the direction, and is the density of the particles. and are the velocity and the density of the fluid, respectively. is the external force of the particles. is the drag force per unit mass.
To simulate the three-dimensional multiphase turbulent and reacting flow in the afterburner of the HPSR, the finite-volume method was adopted to discretize the RANS equations with the software program FLUENT. The theoretical formulation was solved using a numerical scheme based on the semi-implicit method for pressure linked equations (SIMPLE). The second order upwind scheme was utilized to compute the convective flux quantities. The standard - model with stable convergence, appropriate accuracy, and a small amount of calculation was used to address the turbulence in the afterburner [31, 32].
3.2. Combustion Model
3.2.1. Gas-Phase Combustion Model
The gas-phase combustion model employed the single-step reaction eddy-dissipation model (EDM) [33, 34], in which it was assumed that the combustion in the afterburner was a simple chemical reaction system with only two main reactions occurring, as follows:
3.2.2. Boron Particle Ignition and Combustion Model
The typical ignition and combustion models of boron particles include King’s model [35–37] and Williams’s model [38]. The former has a more definite physical meaning and is easy to use to establish a mathematical model, which was chosen as the boron particle ignition and combustion model in this work. The mathematical model could be written as the following formula: where is the diameter of the boron particle, is the temperature of the boron particle, is the energy change during ignition, is the oxide layer thickness of the boron particle, is the evaporation rate of the boron oxide layer, is the rate at which the boron oxide layer reacts with water vapor, and are the consumption rate and the diameter of boron in the oxide layer, respectively, and are the molar mass and the density of B2O3, and and are the specific heat capacity at a constant pressure of B and B2O3, respectively.
The burning rate of the boron particles can be listed as follows: where is the partial pressure of oxygen in the reaction environment.
3.3. Boundary Conditions
The simulated flight altitude is 20,000 meters above sea level, and the Mach number is 3.
Table 1 shows the boundary conditions in the study, which can be summarized as follows: (1) the mass flow inlet condition at the inlet and the entrance of the air inlet, (2) the pressure outlet at the outlet of the afterburner, and (3) the no-slip adiabatic wall condition at the inner wall.
3.4. Characterization of Combustion Efficiency
The component combustion completion rate, with clear physical meaning and simple calculation, was employed to calculate the combustion efficiency. The combustion efficiency of the gas phase component in any cross-section could be expressed by the following formula: where is the integral section. , , , and are the combustion efficiency, density, viscosity coefficient, and mass fraction of the gas phase component in the section, respectively. is the total mass flow rate of component .
The total combustion efficiency of the gas phase can be written as follows: where QCO is the combustion heat of CO. is the combustion heat of H2. λCO is the mass percentage of CO, and is the mass percentage of H2 in the gas phase components.
The combustion efficiency of the boron particles in any cross section can be written as follows: where is the mass flow rate of boron particles remaining on section A. is the total mass flow rate of the boron particles.
The total combustion efficiency in any cross section can be written as follows: where is the mass percentage of the gas phase components and is the combustion heat of boron.
3.5. Grid Independency
A multiblock grid approach was used to generate the computational domain [39]. The computational domain and grid were generated using the commercial software ICEM. The main part of the afterburner used a double O-grid. The grid was clustered near the wall boundary, the primary combustion inlet, and the powder inlet (see Figure 2). To save computing resources and guarantee the accuracy of the calculation, grid independence validation was conducted. As shown in Figure 3, the simulations with different grid numbers for case 2 were conducted. When the grid number was greater than 1.2 million, the grid number had little effect on the film cooling effectiveness of the wall, so the case with a grid number of 1.2 million (1236080) was used in this work. The height of the first row of cells was set at a distance of 0.1 mm for the inner walls. The of case 2 with a grid number of 1.2 million is shown in Figure 4. It can be seen from the figure that the values for the grid were controlled to be under the value of 100, which met the requirements of the standard - model.
3.6. Model Validation
Several experiments were carried out to validate the numerical models with the conditions of case 2 on a direct-connect test bed [40]. The experiment results, including the pressure on line 1 and the temperature at point 1 (as shown in Figure 5), were collected and compared with the simulation results. The result of the pressure is illustrated in Figure 6. The agreement between the simulation results and the experimental result from A-1 was found to be excellent. Table 2 shows the comparison of the temperature at point 1 between the four experimental results and the simulations. From the table, it can be seen that the absolute value of the temperature difference between the experimental results and the simulation result was controlled at about 100 K, and the deviation between the average temperature of the four experimental results and the simulation temperature was 5.36%. In summary, a good agreement was obtained for the simulation results and the experimental results, indicating that the simulation results in this research were accurate, and the calculation model was reliable.
4. Results and Discussion
4.1. Temperature Distribution
4.1.1. Global Temperature Distribution
To explore the effects of different air inlet angles on the combustion, Figure 5 shows the temperature distributions in each cross section of the afterburners. It can be seen from the figure that the global temperature of the internal flow field was axially symmetric. At the head of the afterburners, the temperature was high and even. From case 1 to case 4, the high temperature regions and the low temperature regions gradually shrank, and between the temperature regions, the transition tended to be smooth. This was likely because as the air inlet angle increased, the radial velocity of the ram air was improved, and the impact effect of the ram air on the primary gas was enhanced, causing the ram air at a low temperature and the primary gas at a high temperature to fully mix, so the temperature distribution in the afterburners was more even.
Figure 7 shows the oxygen concentration on the two marking lines on the inner wall of the afterburners. Line 1 passed through the high temperature region, and line 2 mainly passed through the low temperature region, as shown in Figure 5. It can be seen from the figure that the oxygen concentration on line 2 was three orders of magnitude greater than that on line 1, indicating that oxygen was mainly concentrated in the low temperature regions. From case 1 to case 4, the oxygen concentration on line 2 gradually decreased from 19.1% in case 1 to 16.8% in case 4 behind the air inlet exits (), while the oxygen concentration on line 1 increased. Therefore, it could be seen that the oxygen concentration decreased in the low temperature regions and increased in the high temperature regions, showing that the distribution of oxygen in the afterburners tended to be even as the air inlet angle increased.
(a) Oxygen concentration on line 1
(b) Oxygen concentration on line 2
The streamlines of the primary gas and the air are shown in Figure 8 to explore the reasons for the above temperature distribution and the oxygen distribution. It is shown in the figure that behind the air inlet exits, the primary gas was concentrated near the axis of the afterburners, while the air was mainly concentrated near the wall. From case 1 to case 4, the air streamlines and the primary gas streamlines were more interleaved, and the primary gas streamlines tended to diverge. All of these factors indicated that as the air inlet angle increased, the impact effect of the ram air on the primary gas was enhanced, causing the primary gas and the ram air to be fully mixed, so the distribution of the primary gas and air in the afterburner was more even, which was why the temperature and the oxygen in the afterburners showed the above distribution.
4.1.2. Local Temperature Distribution near the Outlet of the Air Inlets
Figure 5 shows that the local temperature distribution near the air inlet exits was also significantly affected by different air inlet angles (see the temperature distribution in the cross section of ). To further explore this phenomenon, Figure 9 shows the temperature distribution on the axial symmetrical sections of the afterburners, from which it can be seen that as the air inlet angle increased, larger local high temperature areas appeared near the air inlet exits, and the maximum temperatures were improved. For example, the maximum temperatures reached 1400 K, 1700 K, 2700 K, and 2900 K from case 1 to case 4.
For the sake of the analysis of the above local temperature distribution, Figure 10 shows the velocity vector on the front sides of exits. As the air inlet angle increased, the number of local vortices increased from one to two, and their range also extended. The section behind the air inlet exits was selected. Its velocity vector is shown in Figure 11. From case 1 to case 4, it can be seen that the local single vortex extended to double vortices with larger ranges. This was because as the air inlet angle increased, the radial velocity of the ram air was accelerated, enhancing the impact between the ram air and the combustion gas in regions near the air inlet exits, so the local vortexes become more intense, and the mixing combustion was strengthened. Therefore, it could be concluded that as the air inlet angle increases, a larger range and a greater intensity of the local vortices were generated on both sides of the air inlet exits. Therefore, the mixed combustion of primary gas and oxygen at those locations was enhanced, resulting in the above local high temperature area.
4.2. Combustion Efficiency
To explore the effect of different air inlet angles on the combustion efficiency of the afterburners, the gas-phase combustion efficiency, the boron particle combustion efficiency, and the total combustion efficiency in each cross section were calculated in this work. The changes in the gas-phase combustion efficiency in each cross section for different cases can be seen in Figure 12. The gas-phase combustion efficiency in the nozzle outlet section () could reach more than 99% for the four cases, indicating that for this kind of afterburner with symmetrical air inlets, the combustion of the primary gas was completed, and different air inlet angles had little effect on the combustion efficiency of the gas-phase components.
Figure 13 shows the changes in the combustion efficiency of the boron particles in each cross section in different cases. The boron particle combustion efficiency order in the four cases from high to low was , where , , , and . The total combustion efficiency order of each cross section in the four cases can be seen in Figure 14, which from high to low was , where , , , and . In comparison, the changing trend of the total combustion efficiency was similar to that of the boron particle combustion efficiency. Thus, it could be inferred that as the air inlet angle changed, the total combustion efficiency was mainly determined by the boron particle combustion efficiency, both of which first increased (from 45° to 75°) and then decreased (from 75° to 90°).
A further study of the boron particle combustion efficiency was carried out. This combustion efficiency is usually determined by factors such as temperature, oxygen concentration, and boron particle residence time in the afterburner [41, 42].
Figure 15 shows the trajectories of the boron particles in the afterburner, which were closely related to the residence time of the boron particles. It can be seen from the figure that the boron particles were mainly concentrated on the head regions of the afterburners and the axial planes of the symmetry of afterburners. In combination with the analysis of the oxygen distribution described in Section 4.1.1, the oxygen concentration in these regions was low where the boron particles were dense. Therefore, it could be considered that the change in the boron particle combustion efficiency was not mainly caused by the change in the oxygen concentration. From case 1 to case 4, the trajectories of the boron particles at the head of the afterburners and near the air inlet exits tended to be dense, and the average residence time was also longer. Combined with the analysis described in Section 4.1.2, it could be inferred that as the air inlet angle increased, the vortexes in the regions of the afterburner heads and near the air inlet exits were enhanced, making more boron particles gather in these regions and move with the vortexes, thus prolonging the residence time of the boron particles in the afterburners.
In addition to the boron particle residence time, the effect of the environment temperature on the boron particle combustion efficiency was investigated in this work, as shown in Figure 16. It can be seen from the figure that from case 1 to case 4, the high temperature region greater than 2600 K in the afterburner first expanded (from 45° to 60°) and then shrank (from 75° to 90°). This was because as the air inlet angle increased, the radial velocity of the ram air increased, which strengthened the mixing between the ram air and the primary gas, enhancing the combustion and rising global temperature in the afterburners. However, when the air inlet increased to a critical value (in this work, the increase was from 60° to 75°), the ram air with a cold temperature squeezed the high temperature combustion gas in the afterburners, and this effect was stronger than the mixing effect caused by ram air, leading to the concentration of the high temperature regions in a narrow range near Line 1 behind the air inlet exits, causing the reduction of the high temperature regions. Combined with the change tendency of the boron particle combustion efficiency in four cases, we inferred that the boron particle combustion efficiency was increased because of the larger high temperature region. When the high temperature region shrank, it led to the reduction of the combustion efficiency.
In summary, increasing the air inlet angle had little effect on the gas-phase combustion efficiency, but it could lead to the boron particle combustion efficiency first increasing and then decreasing. Consequently, the residence time of the boron particles increased, and the high temperature region expanded first as the air inlet angle increased, both of which promoted the combustion efficiency of the boron particles. However, when the angle reached 75°, although the boron particle residence time continued to be prolonged, the high temperature region shrank, which became the dominant factor and reduced the boron particle combustion efficiency. Since the total combustion efficiency was mainly determined by the boron particle combustion efficiency, the total combustion efficiency also had a similar tendency to first increase and then decrease. When the air inlet angle was 75°, the total combustion efficiency was the highest, reaching 84.34%.
4.3. Inner Wall Flow Field
The insulation layer of an afterburner is generally made of composite materials such as silicon rubber, undergoing intense thermal decomposition at a high temperature and serious oxidation, which is closely related to the environment temperature, so it is important to analyze the temperature distribution of the afterburner inner wall. In addition, in a real work environment, the insulation layer is also resistant to severe particle erosion and airflow scouring, which causes the main damage to the afterburner insulation layer [43]. Therefore, in this work, the ablation conditions of the afterburner insulation layer were studied by analyzing the temperature distribution, oxygen concentration, airflow scouring, and erosion of the condensed phase particles of the inner wall [44].
According to Figure 5, the whole afterburner inner wall could be divided into a high temperature wall and a low temperature wall, and line 1 passed through the whole high temperature inner wall, while line 2 went through the main part of the low temperature inner wall. The characteristics of the inner wall were analyzed by studying the parameter distribution on the two making lines. Figure 17 shows the temperature change on the two marking lines. It can be seen from the figure that the temperature on line 1 for the four cases varied from 2300 K to 2900 K. The value range change was smaller than that on line 2, and the temperature difference at the same position of the coordinate for different cases was small. Combined with Figure 7, the oxygen concentration on line 1 had only a slight change, from 0.005% in case 1 to 0.020% in case 4. This indicated that the change of the air inlet angle had little influence on the high-temperature inner wall. However, for the temperature distribution on line 2, there was a maximum temperature after , and it gradually rose from case 1 to case 4, at 1447 K in case 1, 1734 K in case 2, 2763 K in case 3, and 2851 K in case 4. Combined with the analysis described in Section 4.1.2, it could be inferred that as the air inlet angle increased, the radial velocity of the ram air increased, leading to a more severe impact between the ram air and the combustion gas, generating a larger range and more intense vortexes behind the air inlet exits (from to ), so the mixing combustion there was enhanced and the temperature increased. Therefore, as the air inlet angle increased, some influences affected the low temperature inner wall because the temperature of the local inner wall behind the air inlet exits rose, gradually causing a great temperature gradient. Due to the different thermal strain rates of the thermal protection materials, a concentration of the thermal stress was caused, likely leading to crack formation in the insulation layer, which was not conducive to thermal protection.
(a) Temperature distribution on line 1
(b) Temperature distribution on line 2
To explore the airflow scouring effect on the afterburner inner wall, the velocity distributions on two marking lines were calculated, as seen in Figure 18. It can be seen from the figure that behind the air inlet exits (from to ), the airflow velocities on line 1 and line 2 all increased from case 1 to case 4, and the maximum velocities on line 1 and line 2 for case 4 reached 380 m/s and 220 m/s, respectively. The squeezing effect of the ram air on the combustion gas was enhanced because of the increase of the air inlet angle, which let the combustion gas concentrate in the region near the inner wall so that the velocity on the two lines increased. Some researchers believe that for the high temperature jet in a rocket engine combustion chamber, the shearing effect of the high temperature jet on the thermal insulation material surface is the main reason why the material is ablated and damaged, and the effect is severe as the velocity of jet increases [45, 46]. Therefore, it could be determined that as the air inlet angle increased, the airflow velocity on the afterburner inner wall behind the air inlet exits increased, which intensified the airflow scouring on the inner wall.
(a) Velocity distribution on line 1
(b) Velocity distribution on line 2
The density distribution on the two lines was also calculated to study the erosion effect of the condensed phase particles on the inner wall in this work (see Figure 19). It can be seen from the figure that behind the air inlet exits (from to ), the density on line 1 increased from case 1 to case 4, with increases of 0.365 kg/m3 in case 1, 0.358 kg/m3 in case 2, 0.374 kg/m3 in case 3, and 0.381 kg/m3 in case 4, indicating that a larger number of condensed phase particles were collected on the high temperature wall surface by the flow. This was because the stronger squeezing effect of the ram air on the primary gas with boron particles made the particles concentrate in the regions near the high temperature wall. On line 2, although the density on line 2 was about an order of magnitude larger than that on line 1, as seen in Figures 15 and 16, the temperature on the low temperature wall mostly did not exceed 1200 K, and the velocity was low, approximately 100 m/s. Therefore, the condensed phase particles had little erosion effect on the inner wall [47]. In conclusion, with the increase of the air inlet angle, it could be concluded that the erosion effect of the condensed phase particles on the inner wall increased on the high temperature wall, but there was little change on the low temperature wall.
(a) Density distribution on line 1
(b) Density distribution on line 2
In summary, as the air inlet angle increased, the overall jet environment of the inner wall tended to be harsh. For the high temperature wall, the erosion effect of the condensed particles on the inner wall behind the air inlet exits increased; for the low temperature wall, the large temperature gradient was generated on the inner wall near air inlet exits, which meant that the thermal insulation layer tended to crack. The whole inner wall sustained a notable airflow scouring effect because of the accelerated combustion gas flow.
5. Conclusions
In this study, a hybrid powder-solid ramjet (HPSR) with symmetrical air inlets was investigated. The two-phase flow and the combustion in the afterburners were numerically simulated based on King’s boron particle ignition and combustion model, a standard - model, and an eddy-dissipation model (EDM). The effects of different air inlet angles on the combustion efficiency and the inner wall ablation were evaluated. According to the results, as the air inlet angle increased, the following conclusions were drawn: (i)The impact of the ram air on the primary fuel gas was enhanced, and the mixing between the primary gas and the ram air was more thorough, which resulted in the distributions of the temperature and the oxygen tending to be more even(ii)Vortices with higher intensity and a larger range were gradually generated near the air inlet exits, which enhanced the mixing combustion of the primary gas and oxygen there, resulting in local high temperature regions(iii)The boron particle combustion efficiency and the total combustion efficiency first increased and then decreased as the air inlet angle changed. There was an optimal air inlet angle to maximize the total combustion efficiency from 60° to 90° in this work. The total combustion efficiencies of the afterburner were 80.38%, 81.64%, 84.34%, and 83.26% for air inlet angles of 45°, 60°, 75°, and 90°, respectively(iv)The global jet environment in the afterburners was more severe. For the high temperature wall, the erosion of condensed phase particles was enhanced, intensifying the comprehensive damage to the wall; for the low temperature wall, a large temperature gradient was generated near the air inlet exits, in that the thermal insulation layer tended to be cracked there. The gas-flow scouring effect on the inner wall was enhanced in the whole afterburner inner wall(v)When designing the air inlet angle, the total combustion efficiency and the inner wall thermal protection had to be considered comprehensively. It was better to set the air inlet angle between 60° and 75°, which could improve the total combustion efficiency while reducing the difficulty of the thermal protection of the inner wall
Nomenclature
| : | direction |
| : | direction |
| : | direction |
| : | Time (s) |
| %: | Mass fraction of the components |
| : | Air inlet angle |
| : | Conservative vectors |
| : | Vectors of convective flux |
| : | Vectors of viscous flux |
| : | Vectors of source terms |
| : | Velocity of particles (m/s) |
| : | Density of particle (kg/m3) |
| : | Velocity of fluid (m/s) |
| : | Density of fluid (kg/m3) |
| : | Particle diameter (m) |
| : | Oxide layer thickness of particle (m) |
| : | Particle temperature (K) |
| : | Heat fluxes (w/m2) |
| : | Reaction rate (kg/s) |
| : | Molar mass |
| : | Specific heat (J/(kg·K)) |
| : | Partial pressure of component (Pa) |
| : | Combustion efficiency of component |
| : | Density of component (kg/m3) |
| : | Mass fraction of component |
| : | Viscosity of component (Pa·s) |
| : | Burning rate of particle (kg/s) |
| : | Total mass flow rate of component (kg/s) |
| : | The combustion efficiency of gas phase |
| : | Combustion heat of component (J/mol) |
| : | Mass fraction of component in the gas phase. |
Data Availability
Some or all of the data, models, or code generated or used during the study are available from the corresponding author by request.
Conflicts of Interest
The authors declare that they have no conflict of interest.
Acknowledgments
The financial support from the National Natural Science Foundation of China (52003296), the Nature Science Basic Research Program of Shaanxi Province (2020JQ-488), and the Key Laboratory Foundation of Shaanxi Province (2019SZS-09) is gratefully acknowledged. We thank LetPub (http://www.letpub.com) for its linguistic assistance during the preparation of this manuscript.