Research Article  Open Access
Weixuan Li, Dan Zhao, Xiong Chen, Liang Zhu, Siliang Ni, "Numerical Investigation of Inlet Thermodynamic Conditions on Solid Fuel Ramjet Performances", International Journal of Aerospace Engineering, vol. 2021, Article ID 8868288, 20 pages, 2021. https://doi.org/10.1155/2021/8868288
Numerical Investigation of Inlet Thermodynamic Conditions on Solid Fuel Ramjet Performances
Abstract
In this work, 2D numerical RANS (Reynolds Average NavierStokes) simulations were carried out to investigate the thermodynamic performance of a solid fuel ramjet (SFRJ) with different inlet conditions. This is achieved by using an inhouse FORTRAN code to simulate a 2D turbulent, reacting, unsteady flow in the ramjet engine. The inlet conditions are characterized by three key parameters: (1) swirl number (), (2) mass flow rate (), and (3) inlet temperature (). With the code numerically validated by benchmarking with a number of computed cases, it is applied to perform systematic studies on the turbulent flow recirculation, combustion, and heat transfer characteristics. It is found that increasing , , or can dramatically enhance the combustion heat release rate, regression rate, and combustor average temperature. Furthermore, the analysis on the chemical reaction intermediate (CO) reveals that the chemical reaction is more sufficient with increased , but . In addition, a secondary vortex is generated at the corner of the backward facing step in the presence of a swirl flow resulting from the instability of the shear layer. Finally, the nonlinear correlations between the heat transfer, combustion characteristics, and flow field characteristics and the corresponding inlet thermodynamic parameters are identified.
1. Introduction
The ramjet is one of the most popular airbreathing propulsion systems due to its high flight speed and high reliability. The flight of a ramjet is achieved by ram air compression [1]. The compressed incoming air flows through the combustor and is mixed with fuel, and then it is ignited by the hightemperature gas.
Ramjets could be divided into three main types in terms of the type of fuel, namely, the gaseous fuel ramjet (GFRJ), the liquid fuel ramjet (LFRJ), and the solid fuel ramjet (SFRJ). As for the ramjet with gaseous and liquid fuels, a fuel supply system is required for the application of the ramjet, which will lead to an increase in cost and complexity of the propulsion system. Thus, comparing with GFRJ and LFRJ, the ramjet with solid fuel has the simplest structure and the highest reliability, because it does not require a fuel supply system and contains no moving parts as well.
The performance of SFRJ could be mainly determined by the regression rate of the solid fuel. However, due to the lack of a fuel supply system, the solid fuel production rate (regression rate) could not be controlled actively like gaseous and liquid fuels. A previous study [2] suggested that the regression rate was sensitive to the geometries [3–5], flight conditions [6], and the fuel types [7–9]. The thermodynamic performance of SFRJ, scramjet, and hybrid rockets were extensively investigated and reported in the literature [10–20]. A systematic investigation on a smallscale solid fuel ramjet and the combustion characteristic of boronfueled SFRJ were conducted in the Israel Institute of Technology, via both experimental and numerical approaches [21–24]. Research of geometry effects on the performance of SFRJ was investigated in Ref. [4]. It is reported that the local regression rate is closely related to the heat flux. And as described in Ref. [25], the lobe geometry also plays an important role in the heat release and characteristics of a scramjet combustor. Liou et al. [26] carried out a series of numerical studies on the combustion behaviors in SFRJ. In their studies, the correlation between the engine geometry, mass flow rate, and sustained combustion is obtained. Gany [27, 28] conducted a theoretical examination of ideal ramjet performance. The results showed that the specific thrust is increased with an increase of fuel/air ratio up to the stoichiometric ratio. Li et al. [29] preformed an investigation on the lowtemperature ignition characteristics of a supersonic combustor with and without struts, and the results showed that flame stability can be sensitive to the addition of the struts. Li et al. [30] performed an investigation on the combustion characteristics of SFRJ. It was found that variation in engine geometry could lead to a different combustion behavior.
The foregoing studies confirmed that the sensitivity between the SFRJ performance and the geometry, fuel type, and flight conditions is owed mainly to the variation of the thermodynamic parameters. A previous study [31] reported the performance of a solid fuel ramjet with different types of fuel. The results indicated that the thermodynamic constants of the solid fuel could be sensitive to the fuel type. Li et al. [32] conducted a research on fuel reactivitycontrolled selfstarting and propulsion performance of a scramjet, and the results showed that the lowest flight Mach number for selfstarting of a hydrogenfueled scramjet is reduced from 6.2 to 5.1 when the activation energy is decreased by 50%. Li et al. [33] conducted a research of the effects of pressure and oxygen concentration on the combustion characteristics of Al/Mg fuelrich propellants, and the result indicated that the heat conduction and heat feedback to the burning surface could be enhanced by higher environmental pressure, and the combustion performance could be sensitive to the ambient oxygen concentration. Hashim [34] investigated the effects of Ti and Mg particles on the combustion characteristics of boron–HTPBbased solid fuels in a hybrid rocket, and according to their paper, the fuel regression rate is found to be closely related to the thermal conductivity. Reference [35] conducted an investigation on thermodynamic analysis of the specific thrust of a hydrocarbonfueled scramjet, and the results showed that the specific thrust initially can be enhanced with the fuel equivalence ratio, then it can reach a maximum, and finally it can reduce rapidly for a given flight Mach number. Tian et al. [36] numerically investigated the combustion performance of a hybrid rocket motor with segmented grain, and they found that the heat transfer could be sensitive to the flow conditions of the engine. The regression rate was shown to be closely related to the heat transfer of the engine internal flow field. However, the literature of thermal analysis is mainly focused on hybrid rockets and fuels with metal particles. The detailed thermal analysis of a solid fuel ramjet with HDPE is quite limited. This partially motivated the present research.
In this work, a detailed analysis of the effect of the inlet condition on heat transfer and combustion characteristics of SFRJ is investigated. The governing equations and the numerical model are described in Section 2. A CFD code is developed to simulate the twodimensional, turbulent, reacting, unsteady flow in SFRJ. The effects of the swirl number, mass flow rate, and inlet temperature on SFRJ’s thermodynamic performance are discussed in Section 3. Key findings are summarized in Section 4.
2. Numerical Method
2.1. Governing Equations of Fluid Domain
Equation (1) describes the governing equations of the fluid domain [37]: where , , , , , , , and are expressed as follows:
In equation (3), the shear stress and heat fluxes ( and ) in and directions can be calculated by in which is calculated by
The laminar and turbulent viscosities ( and ) of the mixture are calculated using the following [38]:
The equation shown below is applied to calculate the thermal conductivity .
The pressure is given by equation (7), and the Newton iteration method (equation (8)) is used to calculate the temperature () of the flow field: in which that is used in equation (8) is defined by equation (9), and the parameter () is obtained from the chemical kinetics package [39]:
2.2. Governing Equation of Solid Domain
The governing equation of the solid domain can be written as follows [40]:
where denotes the solid phase energy change caused by the chemical reaction [41].
2.3. Numerical Solution Method
To simulate the 2dimensional, turbulent, reacting flow in SFRJ, the RANS equations of the fluid domain are discretized by the finitevolume method. The thirdorder monotone upstreamcentered scheme for conservation laws (MUSCL) and advection upstream splitting method by pressurebased eight function (AUSMPW+) scheme [42] are utilized to compute convective flux quantities. To improve the accuracy and computation efficiency of the unsteady cases, the lower upper symmetric GaussSeidel (LU_SGS) [43] implicit method is used for temporal discretization. Meanwhile, the SST (shearstress transport) turbulence model proposed by Menter [44] is used to simulate the turbulent flow in this research.
At the first 0.5 s of the ignition process, the ignition gas of 0.3 kg/s of mass flow rate, at 2500 K of total temperature, flows through the engine together with incoming air and ignites the fuel. After ignition, sustained combustion is achieved, and the simulation will terminate when a stable result is obtained (at the global time of 1.35 s).
2.4. Chemical Reaction Model
In our previous research, two chemical reaction models (eddydissipation and finite rate models) were compared in the combustion process of swirl intake SFRJ [45]. The result shows little difference between the two models. Thus, in order to save computational resources, the finite rate model is used in this work. However, the chemical reaction of highdensity polyethylene (HDPE) in SFRJ is very complicated, and it is very difficult to describe accurately. Thus, to simplify the combustion, according to Ref. [46], the pyrolysis products of polyethylene can be approximated as C_{2}H_{4} only, and the main products of the combustion are O_{2}, CO_{2}, H_{2}O, CO, and N_{2}. Therefore, a simplified kinetic of chemical reaction [46, 47] shown in Table 1 is applied to describe the chemical reaction of HDPE. In this research, the chemical reaction process is described by the finite rate model, and the Arrhenius formula is utilized to determine reaction rate constants.

where , , and could be obtained from Table 1. The value of for reactions 1, 2, and 3 are 0, 2, and 1, respectively.
Moreover, to calculate the pyrolysis rate of the solid fuel more accurately, the zeroth order () Arrhenius equation [48] is employed:
in which the value of and are 8750 m/s and 130 kJ/mol, respectively [49].
2.5. Physical Model and Boundary Conditions
Figure 1(a) describes the physical model of the SFRJ used in this investigation. The boundary conditions can be summarized as follows: (1) mass flow inlet condition at the inlet, (2) pressure outlet condition at the engine outlet, (3) axisymmetric boundary condition at the axis, and (4) application of noslip adiabatic wall condition at the walls.
(a)
(b)
(c)
The heat transfer between the solid and fluid domains is considered, and the coupled wall temperature is calculated by equation (13). Moreover, for the cases with 0.6 swirl number, the radial profiles of the tangential and axial components of the inlet momentum are defined by the experimental data of Dellenback et al. [50] (see Figures 1(b) and 1(c)).
2.6. Case Description
A series of simulations are conducted with an inhouse code written in FORTRAN, and the investigations concentrate on the different inlet conditions of SFRJ as shown in Table 2. The simulations are performed with different swirl intensities, which are defined by equation (14) (see group A), various inlet mass flow rates (group B), and different inlet temperatures (group C). The model of SFRJ has 40 mm of inlet diameter, 70 mm of port diameter, 500 mm of the fuel length, and 28.5 mm of throat diameter. where is the inlet diameter, and and represent the axial and tangential velocities, respectively.

2.7. Mesh Independency and Model Validation
To save computing resources and guarantee the accuracy of the simulation, the meshindependent studies are conducted at first. The simulations with three different mesh numbers (161580, 113566 and 78569) of case A2 are conducted, and the mesh is clustered near the wall boundary and the solid fuel surface, as shown in Figure 1. The of the cases with different mesh numbers is shown in Figure 2. It could be seen from Figure 2 that thevalues for the fine mesh and the medium mesh are controlled under the value of 5, while the of the coarse mesh is higher than 5. The result of the local regression rate is shown in Figure 3, and the regression rate for the fine mesh and the medium mesh shows little difference with each other but shows significant difference with the case of the coarse mesh. The deviations between the calculated and experimental average regression rates are shown in Table 3; as shown in Table 3, the deviation between the coarse mesh and the experiment is much higher than that of the other two cases. Thus, considering saving the computational resources and improving the accuracy of the simulation, the case with the medium mesh number (113566) is used in this research.

Before the inhouse code can be used for this investigation, the numerical models should be validated against the benchmark cases. An experiment performed by Nejad [51] is utilized to validate the swirl flow model, and a dump combustor model with the swirl number of 0.5 is used in the validation. The result is described in Figure 4, and the agreement between the predicted results of the normalized axial (/) and tangential (/) velocities and experimental results are shown to be excellent.
(a)
(b)
The heat diffusion model of the solid domain is validated by using a semiinfinite plate case. In this validation, the semiinfinite plate is represented by a plate with the size of , density of 7840 kg/m^{3}, specific heat of 465 J/(kg·K), and 49.8 W/(m·K) of thermal conductivity. In addition, the simulated heat transfer is calculated by equation (8), and two different boundary conditions are included at the heated wall surface: (a) heat flux and (b) temperature . Equations (15) and (16) are utilized for the analytical solution of the two cases.
The error function is represented by , and the thermal diffusivity of the material is denoted by . The results are displayed in Figure 5, and the excellent agreements are achieved.
(a)
(b)
A shockinduced combustion benchmark case performed by Lehr [52] is used for the validation of the finite rate chemical reaction model. This validation is performed with a 15 mm diameter spherical projectile moving through a stoichiometric mixture (H_{2} and air) under the pressure of 42662 Pa, temperature of 250 K, and speed of 1685 m/s. The chemical reaction mechanism involving seven species and eightelement reactions is employed. The results are shown in Figure 6, and the experimental data are well reproduced by the results obtained via the inhouse code.
(a)
(b)
A firing experiment on case A2 is performed with a connected pipe facility to validate the accuracy of the inhouse code. In this experiment, the engine (fueled by HDPE) with 70 mm of port diameter, 40 mm of inlet, and 500 mm of fuel length is employed. The comparison between the regression rate of the experiment (measured by a 3D scanner) and that of the simulation is shown in Figure 7. It could be seen from Figure 7 that the trends of the experiment and the simulation are the same, but the location of the reattachment point (the location with the maximum regression rate) has a slight difference in axial direction. The reason behind this is the complicated flow conditions and higher velocity decay near the backward facing step in the experiment. In general, a good agreement is obtained for the predicted results and the experimental result, which indicates that the SFRJ combustion process can be well predicted by the inhouse code.
3. Results and Discussion
3.1. Flow Field Characteristics
Figure 8 illustrates the temperature contours of different cases. It can be seen from Figure 8 that the flow field is divided into two zones by the flame surface (the regions with high temperature over 2200 K), where the core flow area with low temperature is an airrich zone and the regions with high temperature can be regarded as a fuelrich zone. The phenomenon mentioned above is also observed experimentally in Ref. [8]. In the case with a swirl number of 0.6, compared to the nonswirl flow case, the area of the airrich zone is dramatically decreased. However, when the swirl number is increased to 0.74, the airrich zone will be further reduced, but with an insignificant decrease. In addition, under the effect of a swirling inlet flow, the flame surface is pushed forward to the solid domain where the temperatures around the reattachment point and the shear layer are significantly promoted.
(a)
(b)
(c)
For the purpose of better describing the temperature change of the flame surface in the engine, the temperature profile of the flame surface is obtained (see Figure 9). It can be seen from Figure 9 that the introduction of swirling flow into the engine can increase the temperature of the shear layer and the reattachment point. This corresponds with the flame surface temperature distribution shown in the contour; however, the temperature is not enhanced downstream of the flow field. The reason behind this may be due to the attenuation of the swirl intensity at the end of the engine. The temperature profiles of cases B1B3 indicate that with the increase of the inlet mass flow rate, the temperature is also enhanced but with little difference. As for the cases with different (cases C1C3), the temperatures in the shear layer are almost the same. However, in the reattachment point and the recirculation zone, the temperatures are proportional to the inlet temperature. It is interesting that a temperature fluctuation exists in the shear layer, which can indicate that the instability behavior appears in the shear layer.
Figure 10 shows the streamlines of the recirculation zone. The length of the recirculation zone is decreased with the increase of (see Figure 10(a)). The reason behind this is mainly due to the attenuated axial velocity and the enhancement of the tangential velocity caused by the increase of swirl intensity. However, the effect of and on the area of the recirculation zone is not significant.
(a)
(b)
(c)
3.2. Heat Transfer Characteristics
Figure 11 describes the regression rates of the cases with various working conditions. As shown in Figure 11, the local regression rates are all climbing rapidly in the recirculation zone and reach the maximum value at the reattachment point; then, they decrease gradually at the turbulent redevelopment zone. It could be seen from Figure 11 that the regression rate can be enhanced with the increase of , , and . However, the degree of the enhancement on the regression rate with the different working conditions is not the same.
(a)
(b)
(c)
The heat flux of each case is shown in Figure 12, corresponding to the profile of the regression rate; the heat flux has the same trend, and both increase rapidly in the recirculation zone and reach their maximum values at the reattachment point. Moreover, the heat flux can be dramatically improved with a higher swirl number and a higher inlet mass flow rate. Combining with the analysis of the regression rate, it is found that the regression can significantly depend on the surface heat flux.
(a)
(b)
(c)
The time trace of the heat transfer coefficients is shown in Figure 13. In Figure 13, / represents the normalized axial distance, which is calculated by equations (17) and (18). As described in equations (17) and (18), a negative value means that the fuel is in an endothermic state; on the contrary, positive heat transfer coefficients mean that the fuel is in an exothermic state. As shown in Figure 13, at the period of 00.5 s, ignition occurs, and for the cases with a swirling inlet flow, the absolute value of is decreased with the increase of /. When the ignition gas is quenched at 0.5 s, the heat transfer coefficient changes rapidly to positive and reaches its a peak. At the period of 0.50.6 s, is decreased with the increase of /. After 0.7 s, the heat transfer coefficient atis decreased gradually until it is lower than the value at . However, for the case with a nonswirl flow, during the ignition process, the absolute value of the heat transfer coefficient reaches the maximum value at the position of , and the same trend is observed when the selfsustained combustion is achieved.
(a)
(b)
(c)
Based on the streamlines shown in Figure 10, it is known that the position of is within the recirculation zone, and the location of is near the reattachment point. As described in Ref. [53], the heat transfer coefficient can be promoted with the effect of vortices. Thus, the aforementioned phenomenon could be explained as follows. For the swirl flow cases, the tangential momentum is generated with the effect of the swirl flow, which leads to the instability of the shear layer near the backward facing step. At the time of 0.5 s, the quench of the ignition gas enhances the instability of the shear layer. The unstable shear layer then causes the separated boundary layer to roll up into the secondary vortex [54]. As shown in Figure 14(b), at the time of 0.5 s, a secondary vortex appears at the position near for the cases with a swirl flow, and the heat transfer coefficient may be enhanced by the vortices. After the time of 0.6 s, as the combustion process goes by, the secondary vortices vanish. However, as described in Figure 14(a), for the case with a nonswirl flow, the secondary vortex does not exist from the start. Thus, it is believed that the existence of the secondary vortex may lead to the heat transfer coefficient behavior shown in Figure 13.
(a) Case A1
(b) Case A2
Moreover, it could be concluded from Figure 13 that the overall heat transfer coefficient is increased with the increase of , , and .
3.3. Combustion Chrematistics
Figure 15 introduces the time trace of the mass fraction of C_{2}H_{4} (C_{2}H_{4%}) near the fuel surface. It could be seen from Figure 15(a), for the case with a nonswirl flow during the ignition process (00.5 s), that the mass fraction of C_{2}H_{4} increases rapidly over the ignition time and reaches the value of 33% at the location of (recirculation zone) and 0.14 (reattachment point); however, at the redevelopment zone, the concentration of C_{2}H_{4} climbs gradually up to the value of 35%. As the ignition stops, C_{2}H_{4%} drops suddenly and increases gradually until the combustion reaches a stable state, and the maximum value occurs at the reattachment point. For the case with a 0.6 swirl number, after the ignition, the mass fraction of C_{2}H_{4} grows up with the combustion time; then, C_{2}H_{4%} at the recirculation zone starts to be higher than that at the reattachment point. When the swirl number further increases to 0.74, after the ignition, C_{2}H_{4%} at the recirculation zone maintains the maximum value (followed by the C_{2}H_{4} concentration at the reattachment point) among all of the positions until the end of the combustion. The same phenomenon could be observed in Figures 15(b) and 15(c).
(a)
(b)
(c)
Figure 16 reflects the contour of the mass fraction of CO_{2}, O_{2}, and CO, and it could be seen from Figure 16 that CO_{2} is mainly concentrated in the fuelrich zone and O_{2} is concentrated in the airrich zone, and these two zones are divided by the flame surface (see Figure 8). As analyzed in Ref. [3], the combustion in the recirculation zone is taken as a chemical process controlled by the longer residence time of the fuel, and as for the reaction in the shear layer and redevelopment zone, the chemical reaction is subject to a diffusion process. As shown in Figure 16(c), the higher mass fraction of the intermediate of the chemical reaction (CO) appears in the recirculation zone, and it is believed that this is because the chemical reaction in the recirculation zone is incomplete. And this situation occurs in other cases as well.
(a)
(b)
(c)
In order to further analyze the distribution characteristics of the aforementioned components in the combustor, the radial profiles of the mass fraction of CO_{2} (CO_{2%}), O_{2} (O_{2%}), and CO (CO_{%}) are shown in Figures 17–19. As shown in Figure 17, the maximum value of CO_{2%} occurs around the flame surface and fuelrich zone (see Figure 16), and the flame surface is moving towards the fuel surface with the increase of the swirling number, mainly owing to the centrifugal force generated by the swirl flow; however, at the outlet of the combustor, no obvious difference is detected for the position of the flame surface between the cases with and without swirl flow; the reason behind this is the attenuation of the swirl intensity at the end of the combustor. And in the core flow area (airrich zone), for the case of a nonswirl flow, the mass fraction of CO_{2} is 0, and as shown in Figure 16(c), no CO is detected in the core flow area, which indicates that the combustion mainly takes place at the fuelrich zone near the fuel surface. As for the distribution of O_{2} concentration, in the core flow area, the mass fraction of O_{2} is increased with a decreasing swirl number, and it is believed that introducing the swirl flow can increase the mixing degree of fuel and incoming air. According to the chemical mechanism mentioned in Table 1, the concentration of the intermediate could be used to determine whether the chemical reaction is sufficient. As shown in Figure 16(c), the CO_{%} has the highest value in the recirculation zone, and gradually reduces downstream of the flow field. Moreover, the mass fraction of CO is promoted with the enhancement of the swirl number at the flame surface and fuelrich zone, which indicates that the combustion is more sufficient for the case with a nonswirl flow, and the combustion controlled by the chemical process in the recirculation zone is incomplete.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
Moreover, as shown in Figure 18, the distribution of CO_{2%} and O_{2%} has an insignificant difference among the three cases. However, the mass fraction of CO in the fuelrich zone can be promoted with a lower inlet mass flow rate. As described in Figure 19, the combustion affected by the inlet temperature is not obvious.
The heat release rate () is shown in Figure 20, and it could be seen that the flame surface is the place where the chemical reaction mainly occurs. This corresponds to the hightemperature zone (shown in Figure 8). However, a gap is detected in the region of heat release around the mixing plate between the combustor and the aftburning chamber, and it is thought that the compression structure in the flow field of SFRJ may inhibit the chemical reactions. In addition, the maximum value of the heat release rate occurs in the shear layer (see Figure 21); thus, it is considered that the combustion may be enhanced with the expansion structure in SFRJ.
(a)
(b)
(c)
As shown in Figure 21, the heat release rate in the shear layer is increased with a higher swirl number, and in the reattachment point and redevelopment zone, no significant difference is detected. Moreover, as illustrated in Figure 20(b), with a higher inlet mass flow rate, the value of is promoted. It could be seen from Figure 21, for the cases with various inlet temperatures, that the heat release rate is enhanced in the shear layer.
The parameters of ramjet performance and the combustion characteristics of the test cases are shown in Table 4, in which , , , and are calculated by the equations below. It could be seen from Table 4 that the thrust, the characteristic velocity, and the average regression rate can all be enhanced with the increase of the inlet swirl number, mass flow rate, and temperature. However, the specific impulse could be reduced with a higher swirl number and temperature. The reason behind this may be the enhanced regression rate with a relatively lower thrust.

The air equivalence ratio is also calculated based on equation (22); the combustion is in a highly oxidized state with the value of exceeding one, and the fuelrich status on the contrary, if the value of equals one, which means the combustion occurs in stoichiometry. As Table 4 describes, the combustion occurs in a fuelrich state with a higher swirl number and inlet temperature. In addition, based on the evaluation of all the cases with different inlet conditions, case C3, with a 0.5 kg/s inlet mass flow rate, has a relatively higher thrust, specific impulse, and characteristic velocity, and the air equivalence ratio is approximately equal to 1, which indicate that the engine has optimal performance and the combustion is close to the stoichiometric ratio.
4. Conclusions
In this paper, 2D numerical studies were conducted by using an inhouse FORTRAN code to simulate the axisymmetric, unsteady, turbulent, reacting, swirling flow in SFRJ (solid fuel ramjet). The effects of inlet thermodynamic conditions on SFRJ performance have been investigated by varying , , and . The main conclusions are summarized as follows: (i)The turbulent flow field is separated into two zones by the flame surface: (1) the fuelrich zone near the fuel surface where the chemical reaction mainly occurs and (2) the airrich zone in the core flow area(ii)A secondary vortex at the corner of the backward facing step is generated with the introduction of the swirl flow owing to the instability of the shear layer(iii)Two chemical combustion mechanisms are detected in SFRJ: the chemical process in the recirculation zone and the diffusion process in the redevelopment zone and shear layer. According to the analysis of the CO distribution, the reaction controlled by the diffusion process is more sufficient(iv)Increasing , , and is shown to increase the combustor average temperature and the regression rate. In addition, at the reattachment point, the local heat transfer is enhanced, and so the overall heat transfer can be promoted with the increase in the parameters mentioned above(v)The chemical process of SFRJ is analyzed in detail. The pyrolysis behavior of solid fuel depends strongly on the inlet conditions. Based on the analysis of the distribution of the chemical reaction intermediate (CO), the reaction is more sufficient with and higher (vi)The local heat release rate achieves the maximum value in the shear layer, which corresponds to the hightemperature zone. Moreover, the overall heat release can be enhanced with a higher , , and (vii)Based on the elevation of the engine performance, it could be concluded that the best overall performance is achieved when kg/s (i.e., case B3)
Nomenclature
:  Chemical reaction preexponential factor (cm^{3}·mol^{1}·s^{1}) 
:  Area of outlet (m^{2}) 
:  Pyrolytic reaction preexponential factor (m/s) 
:  Area of throat (m^{2}) 
:  Airfuel ratio 
:  Constant () 
:  Characteristic velocity (m/s) 
:  Mass fraction (component ) 
:  Specific heat at constant pressure (J/(kg·k)) 
:  Specific heat at constant pressure of component (J/(kg·k)) 
:  Heat capacity of unit mass of solid fuel (J/(kg·k)) 
:  Diffusion coefficient for component (m^{2}/s) 
:  Convective flux vectors 
:  Energy (J) 
:  Chemical reaction activation energy (J/mol) 
:  Pyrolytic reaction activation energy (kJ/mol) 
:  Vectors of viscous flux 
:  Vectors of convective flux 
:  Thrust of SFRJ (N) 
:  Vectors of viscous flux 
:  Vectors of axisymmetric source terms (convective flux) 
:  Heat transfer coefficient 
:  Enthalpy of unit mass for component (J/kg) 
:  Energy of unit mass solid fuel (J/kg) 
:  Formation heat (component ) (J/kg) 
:  Vectors of axisymmetric source terms (viscous flux) 
:  Specific impulse of SFRJ (N·s/kg) 
:  Turbulent kinetic energy 
:  Length of the solid fuel (mm) 
:  Inlet Mach number 
:  Molecular weight (kg/mol) 
:  Mass flow rate (kg/s) 
:  Temperature exponent 
:  Pressure (Pa) 
:  Pressure of afterburning chamber (Mpa) 
:  Laminar Prandtl numbers 
:  Conservative vectors 
:  Heat release rate (w/m^{2}) 
:  Heat fluxes (w/m^{2}) 
:  Universal gas constant (J·mol^{1}·k^{1}) 
:  Regression rate (mm/s) 
:  Average regression rate (mm/s) 
:  Pyrolysis rate (m/s) 
:  Chemical reaction source term 
:  Swirl number 
:  Laminar Schmidt number 
:  Turbulent Schmidt number 
:  Temperature (K) 
:  Initial temperature (K) 
:  Inlet temperature (K) 
:  Recovery temperature (K) 
:  Axial velocity (m/s) 
:  Radial velocity (m/s) 
:  Tangential velocity (m/s) 
:  Mass source for component 
:  Mole fraction of component 
:  Mole fraction of component 
:  Distance to the wall surface (m) 
:  Specific heat ratio 
:  Air equivalence ratio 
:  Thermal conductivity (w/(m·k)) 
:  Solid fuel thermal conductivity (w/(m·k)) 
:  Viscosity (Pa·s) 
:  Total density (kg/m^{3}) 
:  Density of solid fuel (kg/m^{3}) 
:  Density for component (kg/m^{3}) 
:  Density of solid fuel (kg/m^{3}) 
:  Shear stress (Pa) 
:  Turbulent dissipation rate 
:  Absolute value of vorticity (s^{1}). 
:  Average 
:  Fuel 
:  Fluid domain 
:  Component 
:  Engine inlet 
:  Component 
:  Laminar flow 
:  Outlet surface 
:  Stoichiometry 
:  Turbulent flow 
:  Wall 
:  X direction 
:  Y direction 
:  direction 
%:  Mass fraction of the components. 
HDPE:  Highdensity polyethylene 
SFRJ:  Solid fuel ramjet. 
Data Availability
Some or all 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 first author would like to acknowledge the financial support from the China Scholarship Council which enables him to visit the University of Canterbury and conduct research in Dan Zhao’s group in New Zealand. The financial support from the National Natural Science Foundation of China (11661141020) is also gratefully acknowledged.
References
 R. S. Fry, “A century of ramjet propulsion technology evolution,” Journal of Propulsion and Power, vol. 20, no. 1, pp. 27–58, 2004. View at: Publisher Site  Google Scholar
 O. Musa, C. Xiong, and Z. Changsheng, “Combustion characteristics and turbulence modeling of swirling reacting flow in solid fuel ramjet,” Acta Astronautica, vol. 139, pp. 1–17, 2017. View at: Publisher Site  Google Scholar
 L. Gong, X. Chen, O. Musa, H. Yang, and C. Zhou, “Numerical and experimental investigation of the effect of geometry on combustion characteristics of solidfuel ramjet,” Acta Astronautica, vol. 141, pp. 110–122, 2017. View at: Publisher Site  Google Scholar
 R. Zvuloni, A. Gany, and Y. Levy, “Geometric effects on the combustion in solid fuel ramjets,” Journal of Propulsion and Power, vol. 5, no. 1, pp. 32–37, 1989. View at: Publisher Site  Google Scholar
 L. Yang, G. Yonggang, S. Lei, C. Zexin, and Y. Xiaojing, “Preliminary experimental study on solid rocket fuel gas scramjet,” Acta Astronautica, vol. 153, pp. 146–153, 2018. View at: Publisher Site  Google Scholar
 D. PelosiPinhas and A. Gany, “Bypassregulated solid fuel ramjet combustor in variable flight conditions,” Journal of Propulsion and Power, vol. 19, no. 1, pp. 73–80, 2003. View at: Publisher Site  Google Scholar
 G. Schulte, “Fuel regression and flame stabilization studies of solidfuel ramjets,” Journal of Propulsion and Power, vol. 2, no. 4, pp. 301–304, 1986. View at: Publisher Site  Google Scholar
 G. Schulte, R. Pein, and A. Hogl, “Temperature and concentration measurements in a solid fuel ramjet combustion chamber,” Journal of Propulsion and Power, vol. 3, no. 2, pp. 114–120, 1987. View at: Publisher Site  Google Scholar
 M. Zhu, X. Chen, C.s. Zhou, and J.s. Xu, “Experimental and numerical study of a fullsize directconnect dualinlet DRE with a fuelrich metalized solid propellant,” International Journal of Heat and Mass Transfer, vol. 140, pp. 965–977, 2019. View at: Publisher Site  Google Scholar
 S. Krishnan and Philmon George, “Solid fuel ramjet combustor design,” Progress in aerospace sciences, vol. 34, no. 34, pp. 219–256, 1998. View at: Google Scholar
 Z.b. Du, W. Huang, and L. Yan, “Parametric study on mixing augmentation mechanism induced by air injection in a shockinduced combustion ramjet engine,” Energy, vol. 186, article 115895, 2019. View at: Publisher Site  Google Scholar
 J. Li, Z. Liao, G. Jiao, and W. Song, “The mode transition characteristics in a dualmode combustor at different total temperatures,” Energy, vol. 188, article 116017, 2019. View at: Publisher Site  Google Scholar
 L. Yan, L. Liao, Y.s. Meng, S.b. Li, and W. Huang, “Investigation on the mode transition of a typical threedimensional scramjet combustor equipped with a strut,” Energy, vol. 208, article 118419, 2020. View at: Publisher Site  Google Scholar
 P. J. M. Elands, P. A. O. G. Korting, T. Wijchers, and F. Dijkstra, “Comparison of combustion experiments and theory in polyethylene solid fuel ramjets,” Journal of Propulsion and Power, vol. 6, no. 6, pp. 732–739, 1990. View at: Publisher Site  Google Scholar
 J. GobboFerreira, M. G. Silva, and J. A. Carvalho Jr., “Performance of an experimental polyethylene solid fuel ramjet,” Acta Astronautica, vol. 45, no. 1, pp. 11–18, 1999. View at: Publisher Site  Google Scholar
 M. Ali and A. K. M. S. Islam, “Study on main flow and fuel injector configurations for scramjet applications,” International Journal of Heat and Mass Transfer, vol. 49, no. 1920, pp. 3634–3644, 2006. View at: Publisher Site  Google Scholar
 H. Zhang, N. Wang, Z. Wu, W. Han, and R. Du, “A new model of regression rate for solid fuel scramjet,” International Journal of Heat and Mass Transfer, vol. 144, article 118645, 2019. View at: Publisher Site  Google Scholar
 A. S. Vishnu, G. P. Aravind, M. Deepu, and R. Sadanandan, “Effect of heat transfer on an angled cavity placed in supersonic flow,” International Journal of Heat and Mass Transfer, vol. 141, pp. 1140–1151, 2019. View at: Publisher Site  Google Scholar
 W. J. You, H. J. Moon, S. P. Jang, and J. K. Kim, “Thermal characteristics of an N_{2}O catalytic igniter with metal foam for hybrid rocket motors,” International Journal of Heat and Mass Transfer, vol. 66, pp. 101–110, 2013. View at: Publisher Site  Google Scholar
 Y. Liu, Y. Gao, Z. Chai, Z. Dong, C. Hu, and X. Yu, “Mixing and heat release characteristics in the combustor of solidfuel rocket scramjet based on DES,” Aerospace Science and Technology, vol. 94, article 105391, 2019. View at: Publisher Site  Google Scholar
 Y. Levy, A. Gany, and R. Zvuloni, “Investigation of a small solid fuel ramjet combustor,” Journal of Propulsion and Power, vol. 5, no. 3, pp. 269–275, 1989. View at: Publisher Site  Google Scholar
 R. BenArosh and A. Gany, “Similarity and scale effects in solidfuel ramjet combustors,” Journal of Propulsion and Power, vol. 8, no. 3, pp. 615–623, 1992. View at: Publisher Site  Google Scholar
 B. Natan and A. Gany, “Combustion characteristics of a boronfueled solid fuel ramjet with aftburner,” Journal of Propulsion and Power, vol. 9, no. 5, pp. 694–701, 1993. View at: Publisher Site  Google Scholar
 B. Natan and A. Gany, “Effects of bypass air on boron combustion in solid fuel ramjets,” Journal of Propulsion and Power, vol. 9, no. 1, pp. 155–157, 1993. View at: Publisher Site  Google Scholar
 G. Yonggang, L. Yang, C. Zexin, L. Xiaocong, H. Chunbo, and Y. Xiaojing, “Influence of lobe geometry on mixing and heat release characteristics of solid fuel rocket scramjet combustor,” Acta Astronautica, vol. 164, pp. 212–229, 2019. View at: Publisher Site  Google Scholar
 T.M. Liou, W.Y. Lien, and P.W. Hwang, “Flammability limits and probability density functions in simulated solidfuel ramjet combustors,” Journal of Propulsion and Power, vol. 13, no. 5, pp. 643–650, 1997. View at: Publisher Site  Google Scholar
 A. Gany, “Effect of fuel properties on the specific thrust of a ramjet engine,” Defence Science Journal, vol. 56, no. 3, pp. 321–328, 2006. View at: Publisher Site  Google Scholar
 A. Gany, “Parametric Analysis of the Ideal Ramjet Performance,” 16th International Symposium on Airbreathing engines, vol. 31, 2003. View at: Google Scholar
 J. Li, J. Li, K. Wang, G. Jiao, and Z. Liao, “Study of lowtemperature ignition characteristics in a supersonic combustor,” Energy, vol. 195, article 117060, 2020. View at: Publisher Site  Google Scholar
 W. Li, X. Chen, O. Musa, L. Gong, and L. Zhu, “Investigation of the effect of geometry of combustor on combustion characteristics of solidfuel ramjet with swirl flow,” Applied Thermal Engineering, vol. 145, pp. 229–244, 2018. View at: Publisher Site  Google Scholar
 B. McDonald and J. Rice, “Solid fuel ramjet fuel optimization for maximum impulsedensity with respect to air to fuel ratio and relative fuel regression rates derived from thermogravimetric analysis,” Aerospace Science and Technology, vol. 86, pp. 478–486, 2019. View at: Publisher Site  Google Scholar
 X. Li, X. Huang, H. Liu, and J. Du, “Fuel reactivity controlled selfstarting and propulsion performance of a scramjet: a model investigation,” Energy, vol. 195, article 116920, 2020. View at: Publisher Site  Google Scholar
 L.b. Li, X. Chen, C.s. Zhou, M. Zhu, and O. Musa, “Experimental and numerical investigations of the effect of pressure and oxygen concentration on combustion characteristics of Al/Mg fuelrich propellants,” Applied Thermal Engineering, vol. 167, article 114695, 2020. View at: Publisher Site  Google Scholar
 S. A. Hashim, S. Karmakar, and A. Roy, “Effects of Ti and Mg particles on combustion characteristics of boron–HTPBbased solid fuels for hybrid gas generator in ducted rocket applications,” Acta Astronautica, vol. 160, pp. 125–137, 2019. View at: Publisher Site  Google Scholar
 Q. Yang, J. Chang, and W. Bao, “Thermodynamic analysis on specific thrust of the hydrocarbon fueled scramjet,” Energy, vol. 76, pp. 552–558, 2014. View at: Publisher Site  Google Scholar
 H. Tian, Y. Duan, and H. Zhu, “Threedimensional numerical analysis on combustion performance and flow of hybrid rocket motor with multisegmented grain,” Chinese Journal of Aeronautics, vol. 33, no. 4, pp. 1181–1191, 2020. View at: Publisher Site  Google Scholar
 K. K. Kuo, Principles of Combustion. No. TJ254. 5 K85 2005, John Wiley & Sons, Inc., Hoboken, New Jersey, 2005.
 J. Chen, Combustion Fundamental of Solid Rocket, Nanjing University of Science and Technology Press, Nanjing, 2011.
 J. Kee, F. M. Rupley, E. Meeks, and J. A. Miller, CHEMKINIII: A FORTRAN Chemical Kinetics Package for the Analysis of GasPhase Chemical and Plasma Kinetics, Technical Report, Sandia National Labs, 1996.
 S. I. Stoliarov and R. N. Walters, “Determination of the heats of gasification of polymers using differential scanning calorimetry,” Polymer Degradation & Stability, vol. 93, no. 2, pp. 422–427, 2008. View at: Publisher Site  Google Scholar
 D. Bianchi, F. Nasuti, and M. Onofri, “Radius of curvature effects on throat thermochemical erosion in solid rocket motors,” in 50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, pp. 320–330, Nashville, Tennessee, 2014. View at: Publisher Site  Google Scholar
 K. H. Kim, K. Chongam, and O. H. Rho, “Methods for the accurate computations of hypersonic flows,” Journal of Computational Physics, vol. 174, no. 1, pp. 38–80, 2001. View at: Publisher Site  Google Scholar
 L. P. Zhang and Z. J. Wang, “A block LUSGS implicit dual timestepping algorithm for hybrid dynamic meshes,” Computers and Fluids, vol. 33, no. 7, pp. 891–916, 2004. View at: Publisher Site  Google Scholar
 F. R. Menter, “Twoequation eddyviscosity turbulence models for engineering applications,” AIAA Journal, vol. 32, no. 8, pp. 1598–1605, 1994. View at: Publisher Site  Google Scholar
 O. Musa, C. Xiong, Z. Changsheng, and Z. Min, “Combustion modeling of unsteady reacting swirling flow in solid fuel ramjet,” in 2017 International Conference on Mechanical, System and Control Engineering (ICMSC), St. Petersburg, Russia, March 2017. View at: Publisher Site  Google Scholar
 M. Mawid and B. Sekar, “Kinetic modeling of ethylene oxidation in high speed reacting flows,” in 33rd Joint Propulsion Conference and Exhibit, p. 3269, Seattle, WA, USA, 1997. View at: Publisher Site  Google Scholar
 R. Baurle, T. Mathur, M. Gruber, and K. Jackson, “A numerical and experimental investigation of a scramjet combustor for hypersonic missile applications,” in 34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, vol. 3121, Cleveland, OH, USA, 1998. View at: Publisher Site  Google Scholar
 M. J. Chiaverini, G. C. Harting, Y. C. Lu et al., “Pyrolysis behavior of hybridrocket solid fuels under rapid heating conditions,” Journal of Propulsion and Power, vol. 15, no. 6, pp. 888–895, 1999. View at: Publisher Site  Google Scholar
 L. Gong, X. Chen, C. S. Zhou, Y. K. Li, and M. Zhu, “Numerical investigation on effect of solid fuel ramjet geometry on solid fuel regression rate,” Acta Armamentarii, vol. 37, pp. 798–807, 2016. View at: Google Scholar
 P. Dellenback, D. Metzger, and G. Neitzel, “Measurements in turbulent swirling flow through an abrupt axisymmetric expansion,” AIAA Journal, vol. 26, no. 6, pp. 669–681, 1988. View at: Publisher Site  Google Scholar
 A. S. Nejad, S. P. Vanka, S. C. Favaloro, M. Samimy, and C. Langenfeld, “Application of laser velocimetry for characterization of confined swirling flow,” Journal of Engineering for Gas Turbines and Power, vol. 111, pp. 36–45, 1989. View at: Publisher Site  Google Scholar
 H. F. Lehr, “Experiments on shockinduced combustion,” Astronautica Acta, vol. 17, pp. 589–597, 1972. View at: Google Scholar
 W. Li, X. Chen, W. Cai, and O. Musa, “Numerical investigation of the effect of sudden expansion ratio of solid fuel ramjet combustor with swirling turbulent reacting flow,” Energies, vol. 12, no. 9, p. 1784, 2019. View at: Publisher Site  Google Scholar
 T. Wei and C. R. Smith, “Secondary vortices in the wake of circular cylinders,” Journal of Fluid Mechanics, vol. 169, pp. 513–533, 1986. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2021 Weixuan Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.