Research Article

Journal of Propulsion and Energy. 30 November 2025. 17-26
https://doi.org/10.6108/JPNE.2025.5.2.017

ABSTRACT


MAIN

  • 1. Introduction

  • 2. Numerical Method

  •   2.1 Pintle Injector and Combustor Design

  •   2.2 Numerical Modeling

  •   2.3 Grid Independence Test

  •   2.4 Validation

  • 3. Results and Discussion

  •   3.1 Effects of Momentum Ratio

  •   3.2 Effects of Blockage Factor

  •   3.3 Effects of Mixture Ratio

  • 4. Conclusion

1. Introduction

Following the successful landing of the Apollo lunar module on the moon surface in 1969 [1], pintle injectors are receiving renewed attention after their applications in SpaceX’s Falcon 9 and the Chinese lunar lander Chang’E [2,3]. A pintle injector consists of a radial orifice on the side of the pintle and an annular axial orifice around the pintle, which inject oxidizer and fuel, respectively [4]. This unique structure is known to provide controllable propellant injection area, simplicity, and inherent combustion stability. Pintle injectors, originally developed to improve the mixing performance of hypergolic liquid-liquid propellants, are currently also applied to a variety of propellant combinations for gas-liquid injection, such as LOx/H2[5], LOx/CH4[6], and C2H6O/H2O2[7,8].

Most of the 20th century pintle injector research was conducted by TRW [9], and it was reported that the highest combustion efficiency was achieved at an interaction momentum ratio (IMR), which is expressed by Fig. 1 and Eq. (1), around 1 [10,11]. These findings have been repeatedly cited by previous studies as a solid guideline for pintle injector design [12,13] and experimentally reaffirmed by the liquid-liquid pintle injector hot firing tests of Sakaki et al. [14,15].

However, for gas-liquid pintle injectors, some recent works have reported that the peak combustion efficiency was observed at momentum ratios far from 1. Kim et al. experimentally investigated the combustion characteristics of a LOx/GCH4 pintle injector and reported that the combustion efficiency reached its maximum at a momentum ratio of around 0.1 [16]. Vasques et al. reported that the maximum combustion efficiency was achieved at a momentum ratio between 0.1 and 0.2, through hot firing tests of a LOx/GCH4 pintle injector [17]. These results raise questions about TRW’s previous findings and necessitate further exploration into combustion efficiency across a broader range of momentum ratios in gas-liquid pintle injectors, along with other influencing factors.

https://cdn.apub.kr/journalsite/sites/jpne/2025-005-02/N0640050202/images/jpne_2025_52_017_F1.jpg
Fig. 1.

Schematics for interaction momentum ratio, IMR.

Therefore, in this study, numerical analysis was conducted over a wide range of momentum ratio to investigate how the combustion efficiency of a gas-liquid pintle injector varies with the momentum ratio. Furthermore, the effects of blockage factor (BF) and mixture ratio (O/F) on the combustion efficiency were also discussed.

(1)
IMR=ρrVr2WLρaVa2(GW+2GLC), where C=LDαVaBF,D=πDpN-W,BF=NWπDp,α200m/s

2. Numerical Method

2.1 Pintle Injector and Combustor Design

Fig. 2 shows the schematics of a multi-slit type gas-liquid pintle injector used in this study [18]. The injector was designed for a 400 N-class LOx/CH4 rocket engine with a combustion pressure of 10 bar, and it injects LOx through the radial orifices and GCH4 through the axial orifice. Table 1 lists the specifications of the pintle injector and combustor.

Table 1.

Specifications of the pintle injector and combustor.

Parameter (Unit) Value
Pintle diameter, Dp(mm) 11.0
Skip length, Ls (mm) 11.0
Sleeve thickness (mm) 1.0
Annular axial orifice gap, G (mm) 0.45
Chamber diameter, Dc (mm) 46.0
Nozzle throat diameter, Dt (mm) 15.8
Chamber length, L (mm) 136.3
Characteristic length, L* (mm) 1,000.0

https://cdn.apub.kr/journalsite/sites/jpne/2025-005-02/N0640050202/images/jpne_2025_52_017_F2.jpg
Fig. 2.

Schematics of the pintle injector.

To investigate the combustion efficiency across momentum ratios above and below unity under a fixed propellant mass flow rate, three types of pintle tips shown in Fig. 3 were used [19]. Additionally, five mixture ratios ranging from 2.72 to 4.08 were used. As a result, numerical analysis was performed for 15 cases. Table 2 lists the specifications of the pintle tips and calculated IMRs depending on pintle tip types and mixture ratios.

https://cdn.apub.kr/journalsite/sites/jpne/2025-005-02/N0640050202/images/jpne_2025_52_017_F3.jpg
Fig. 3.

Schematic of the pintle tips.

Table 2.

Specifications of the pintle tips and IMRs for each case.

Parameter (Unit) Type #1 Type #2 Type #3
Slit height, H (mm) 1.10 1.10 1.10
Slit width, W (mm) 0.50 0.50 0.50
Slit end radius, R (mm) 0.25 0.25 0.25
Number of slits, N 12 18 24
Blockage factor, BF 0.174 0.260 0.347
IMRO/F = 2.72 2.058 0.744 0.319
O/F = 3.06 2.569 0.916 0.386
O/F = 3.40 3.129 1.100 0.458
O/F = 3.74 3.736 1.297 0.533
O/F = 4.08 4.388 1.504 0.611

2.2 Numerical Modeling

The three-dimensional steady-state RANS numerical simulations were performed using ANSYS Fluent 2024 R2 [20]. The numerical modeling of the pintle injector combustion is as follows. The computational domain was set to a region corresponding to 1/3 of the combustor, as shown in Fig. 4. It consisted of the inlets of the oxidizer and fuel, the outlet corresponding to the nozzle, the inner walls of the combustor, and the two symmetry planes. For the turbulence model, the k-ε standard model was used. To represent fast chemical reactions in the combustor with both accuracy and simplicity, the eddy dissipation concept model (EDC) was used for the species model [21]. The DRM-22 mechanism was employed to express the combustion between oxygen and methane [22].

LOx at 90 K was injected as particles at each inlet corresponding to the radial flow orifices, while GCH4 at 298.15 K was injected as a continuum. Additionally, the secondary breakup of the injected LOx particles was represented by the WAVE model, with the recommended breakup coefficients (B0 = 0.61, B1 = 1.73) [23,24].

https://cdn.apub.kr/journalsite/sites/jpne/2025-005-02/N0640050202/images/jpne_2025_52_017_F4.jpg
Fig. 4.

Schematics of the computational domain.

2.3 Grid Independence Test

Four levels of grids were used to ensure grid independence of the numerical model: coarse, moderate, fine, and finest grids. The grid independence test was conducted for pintle tip Type #1 at a mixture ratio of 3.40. The characteristic velocity (c*), calculated using Eq. (2), and the volume-averaged chamber temperature were compared for each grid. As a result, the moderate grid exhibited sufficiently converged values, suggesting that the minimum cell size of this grid level is acceptable for all numerical simulation cases. Table 3 shows the conditions and results of the test.

(2)
c*=PcAtm˙
Table 3.

Conditions and results of the grid independence test.

Grid type Minimum cell size
(mm)
No. of cells Characteristic velocity (m/s) Volume-averaged chamber temperature (K)
Coarse 0.75 295,611 1307.0 1861.4
Moderate 0.60 542,419 1376.5 1878.4
Fine 0.50 902,938 1387.0 1884.2
Finest 0.40 1,713,997 1392.2 1887.9

2.4 Validation

To secure the reliability of the numerical model, this model was applied to the pintle injector hot firing test conditions of a previous study [16] using the same LOx/GCH4 propellant combination. Numerical analysis was performed for eight experimental cases from the previous study, and the results were compared with the experimental results. Fig. 5 displays the validation results based on the total momentum ratio (TMR), which is expressed by Eq. (3). The results indicate that the numerical model agrees well with the experimental data across all cases, with an average accuracy of 97.9% and a maximum error of 3.6%.

https://cdn.apub.kr/journalsite/sites/jpne/2025-005-02/N0640050202/images/jpne_2025_52_017_F5.jpg
Fig. 5.

Validation results of the numerical model.

(3)
TMR=ρrVr2ArρaVa2Aa

3. Results and Discussion

Fig. 6 shows the combustion efficiency versus IMR, with the pintle tip type and mixture ratio for each case also indicated. The combustion efficiency was represented by c* efficiency, which was expressed by Eq. (4). The following describes the effect of each factor on combustion efficiency.

(4)
ηc*=c*cideal *

https://cdn.apub.kr/journalsite/sites/jpne/2025-005-02/N0640050202/images/jpne_2025_52_017_F6.jpg
Fig. 6.

Combustion efficiency versus momentum ratio.

3.1 Effects of Momentum Ratio

First, the c* efficiency was inversely proportional to the IMR across the range of 0.3 to 4.4. Although the maximum c* efficiency was not observed within this IMR range, it can be inferred that the maximum value would occur at an IMR less than 1. This result contradicts TRW’s hot firing test results [10,11] but aligns with the findings of Kim et al. and Vasques et al. [16,17]. It was reconfirmed that the peak c* efficiency in gas-liquid pintle injectors can occur at an IMR other than 1. Therefore, further exploration of c* efficiency across a wider range of IMRs is required.

3.2 Effects of Blockage Factor

The BF varied depending on the pintle tip type, increasing from 0.174 to 0.347 as the pintle tip type changed from #1 to #3. Correspondingly, c* efficiency also increased. This result can be interpreted from two perspectives. First, it can be attributed to an increase in the interaction area between the oxidizer and fuel due to the increased BF as shown in Fig. 7[25], which subsequently enhanced mixing performance and c* efficiency. Secondly, the increased BF and consequently decreased IMR may have resulted in more efficient liquid radial flow atomization by the axial gaseous flow, thus improving the c* efficiency.

Fig. 7 shows the mole fraction contours of O2 and CH4 for different BFs at a fixed mixture ratio of 3.40. For Type #1, which has a small BF (0.174), the injected O2 impinges on the combustor inner wall, with a significant portion proceeding along the wall towards the nozzle. Simultaneously, CH4 is injected axially and flows along the central axis of the combustor, indicating insufficient mixing between O2 and CH4. On the other hand, for Type #3, which has a relatively larger BF of 0.347, O2 and CH4 are uniformly distributed in the front part of the combustor. Due to their active mixing and reaction, the mole fraction of each significantly decreases towards the rear of the combustor. From these observations, it can be concluded that an increase in BF can enhance propellant mixing performance, thereby improving combustion efficiency. However, it is difficult to distinguish whether the primary cause is the increased interaction area between the oxidizer and fuel, or the momentum ratio reaching an optimal value. Further research is necessary to investigate this by fixing either the momentum ratio of the BF while varying the other.

https://cdn.apub.kr/journalsite/sites/jpne/2025-005-02/N0640050202/images/jpne_2025_52_017_F7.jpg
Fig. 7.

Contours of O2 (top) and CH4 (bottom) for different blockage factor cases.

3.3 Effects of Mixture Ratio

The c* efficiency also varied with the mixture ratio. For all pintle tip types, higher c* efficiency was observed at lower mixture ratios. This indicates that at higher mixture ratios, the interaction between the oxidizer and fuel is insufficient. Conversely, it can be inferred that better interaction occurs at lower mixture ratios, leading to higher c* efficiency.

This can also be observed from the distribution of oxidizer and fuel in the combustor. Fig. 8 shows the mole fraction contours of O2 and CH4 for different mixture ratios when using pintle tip Type #2. When the mixture ratio is 2.72, O2 and CH4 are distributed in the front part of the combustor, and their mole fractions distinctly decrease towards the rear as the chemical reaction occurs. However, at a mixture ratio of 4.08, O2 reaches near the combustor inner wall in the front section and is partially present even close to the nozzle. CH4 also remains distributed along the combustor’s central axis, indicating that the mixing and reaction between O2 and CH4 are limited in this case. Therefore, it can be considered that the combustion efficiency of pintle injectors is influenced not only by the mixture ratio itself but also by the momentum ratio, which varies with the mixture ratio.

https://cdn.apub.kr/journalsite/sites/jpne/2025-005-02/N0640050202/images/jpne_2025_52_017_F8.jpg
Fig. 8.

Contours of O2 (left) and CH4 (right) for different mixture ratio cases.

4. Conclusion

In this study, the effects of injection conditions on combustion efficiency of a multi-slit type gas-liquid pintle injector were explored. In summary, the following conclusions were derived.

(1) The c* efficiency increased as the momentum ratio decreased. Notably, the highest c* efficiency was not observed at a momentum ratio of 1. This stands in conflict with TRW’s results, implying that the momentum ratio exhibiting maximum combustion efficiency in gas-liquid pintle injectors is not unity. While a wide range of momentum ratios was investigated, the momentum ratio with the highest combustion efficiency could not be identified, thus necessitating further research.

(2) The blockage factor significantly affected the combustion efficiency of the multi-slit type gas-liquid pintle injector. As the blockage factor increased from 0.174 to 0.347, the combustion efficiency also increased. This indicates that the mixing between O2 and CH4 improved due to the increased blockage factor. Furthermore, with a fixed propellant mass flow rate and mixture ratio, an increased blockage factor causes a decrease in the momentum ratio. This reduced momentum ratio may have improved the interaction between the oxidizer and fuel.

(3) The mixture ratio also affected the combustion efficiency of the multi-slit type gas-liquid pintle injector. Regardless of the blockage factor, the combustion efficiency increased as the mixture ratio decreased from 4.08 to 2.72. This suggests that, within the cases explored in this study, the radial flow momentum was excessively high compared to the axial flow momentum. This imbalance prevented sufficient interaction between the oxidizer and fuel at higher mixture ratios, leading to higher combustion efficiency at lower mixture ratios.

The results of this study raise questions about TRW’s design criteria for gas-liquid pintle injectors and highlight the necessity for further research into combustion efficiency under more comprehensive and broader injection conditions, such as momentum ratio, blockage factor, and mixture ratio. In addition, this study was conducted based on RANS simulation; therefore, high-fidelity numerical simulations are necessary for a more accurate prediction of combustion efficiency.

Acknowledgements

This work was supported by the Institute of Advanced Aerospace Technology at Seoul National University, and by the National Research Foundation of Korea (NRF) and the Korea institute of Energy Technology Evaluation and Planning (KETEP) grants funded by the Korea Government (MSIT, MOTIE) (RS-2023-00249707, RS-2023-00254668).

References

1

Dressler, G. A., “Summary of deep throttling rocket engines with emphasis on Apollo LMDE,” 42nd AIAA/ASME/SAE/ASEE Joint Propulsion Conference, Sacramento, California, 2006.

10.2514/6.2006-5220
2

Lei, J. P., Lan, X. H., Zhang, R. H., and Chen, W., “The development of 7500 N variable thrust engine for Chang’E-3,” Scientia Sinica Technologica, Vol. 44, No. 6, 2014, pp. 569-575.

10.1360/092014-52
3

Seedhouse, E., SpaceX: Making commercial spaceflight a reality, Springer, New York, 2013.

10.1007/978-1-4614-5514-1
4

Heister, S. D., Handbook of atomization and sprays, Springer, New York, 2011.

5

Gromski, J. M., Majamaki, A. N., Chianese, S. G., Weinstock, V. D., and Kim, T. S., “Northrop Grumman TR202 LOX/GH2 deep throttling engine project status,” 46th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Nashville, Tennessee, 2010.

10.2514/6.2010-6725
6

Suzuki, L., Nakayama, S., and Sakai, K., “ZERO by Interstellar Technologies Inc.: Lowering the cost of access to space from Japan,” Small Satellite Conference 2024, Logan, Utah, 2024.

7

Hämmerli, B., Liggieri, G., Hegre, J. S., Ruud, J. I., Maddali, P., Worsley, E., Shaw, M., and Haynes, A., “Design, manufacturing and testing of a 6kn green throttleable development engine,” 9th Space Propulsion Conference, Glasgow, Scotland, 2024.

8

Cieśliński, D., Okniński, A., Ranachowski, M., Molik, M., Bartkowiak, B., Matysek, K., Kutnik, K., Gorgeri, A., Ogórek, G., and Borys, M., “On development of green storable liquid rocket engine with thrust variation,” Acta Astronautica, Vol. 234, 2025, pp. 560-574.

10.1016/j.actaastro.2025.05.020
9

Dressler, G. A., and Bauer, J. M., “TRW Pintle Engine geritage and performance characteristics,” 36th AIAA/ASME/SAE/ASEE Joint Propulsion Conference, Las Vegas, Nevada, 2000.

10.2514/6.2000-3871
10

Mageean, J. V., “Space storage propellant performance final report,” NASA-CR-72486, 1969.

11

Carter, W. A., “Gas-liquid space storable propellant performance final report,” NASA-CR-72708, 1970.

12

Zhao, F., Zhang, H., Zhang, H., Bai, B., and Zhao, L., “Review of atomization and mixing characteristics of pintle injectors,” Acta Astronautica, Vol. 200, 2022, pp. 400-419.

10.1016/j.actaastro.2022.08.042
13

Hwang, D. K., Ryu, C., and Kwon, S., “Analysis on the research and development cases of combustion devices with liquid-Liquid pintle injector,” Journal of the Korean Society of Propulsion Engineers, Vol. 24, No. 6, 2020, pp. 126-142.

10.6108/KSPE.2020.24.6.126
14

Sakaki, K., Kakudo, H., Nakaya, S., Tsue, M., Kanai, R., Suzuki, K., Inagawa, T., and Hiraiwa, T., “Performance evaluation of rocket engine combustors using ethanol/liquid oxygen pintle injector,” 52nd AIAA/SAE/ASEE Joint Propulsion Conference, Salt Lake City, Utah, 2016.

10.2514/6.2016-5080
15

Sakaki, K., Funahashi, T., Nakaya, S., Tsue, M., Kanai, R., Suzuki, K., Inagawa, T., and Hiraiwa, T., “Longitudinal combustion instability of a pintle injector for a liquid rocket engine combustor,” Combustion and Flame, Vol. 194, 2018, pp. 115-127.

10.1016/j.combustflame.2018.04.017
16

Kim, D. H., Heo, S., Kim, I., Hwang, D., Kang, C., Lee, S., Ahn, K., and Yoon, Y., “High pressure spray and combustion characteristics of throttleable pintle injector,” Journal of the Korean Society of Propulsion Engineers, Vol. 26, No. 2, 2022, pp. 60-71.

10.6108/KSPE.2022.26.2.060
17

Vasques, B. B., and Haidn, O. J., “Effect of pintle injector element geometry on combustion in a liquid oxygen/liquid methane rocket engine,” 7th European Conference for Aeronautics and Aerospace Sciences (EUCASS), Milan, Italy, 2017.

18

Heo, S., Jeong, J., and Yoon, Y., “Spray characteristics of a multi-slit type throttleable pintle injector with different slit Heights,” Acta Astronautica, Vol. 215, 2024, pp. 475-492.

10.1016/j.actaastro.2023.12.023
19

Heo, S., Jeong, J., Kim, D. H., Yoon, Y., and Lee, B. J., “Spray angle prediction model for a multi-slit type throttleable gas-liquid pintle injector,” Acta Astronautica, Vol. 231, 2025, pp. 113-126.

10.1016/j.actaastro.2025.01.067
20

Heo, S., “Spray and combustion characteristics of a multi-slit type throttleable gas-liquid pintle injector,” Ph.D. Dissertation, Department of Aerospace Engineering, Seoul National University, Seoul, Korea, 2025.

10.1063/5.0274886
21

Magnussen, B. F., “The eddy dissipation concept: A bridge between science and technology,” ECCOMAS Thematic Conference on Computational Combustion, Lisbon, Portugal, 2005.

22

Kazakov, A., and Frenklach, M., “Reduced reaction sets based on GRI-Mech 1.2,” 1994.

23

Reitz, R. D., “Modeling atomization processes in high-pressure vaporizing sprays,” Atomisation and Spray Technology, Vol. 3, No. 4, 1987, pp. 309-337.

24

Radhakrishnan, K., Lee, K., and Koo, J., “A detailed modelling on spray atomisation and combustion of LOX/GCH4 in variable area pintle injector rocket engines,” Combustion Theory and Modelling, Vol. 25, No. 4, 2021, pp. 718-750.

10.1080/13647830.2021.1937702
25

Heo, S., Kim, I., Yoon, Y., and Lee, B. J., “Spray structure of a multi-slit type throttleable gas-liquid pintle injector,” Physics of Fluids, Vol. 37, No. 8, 2025, 083328.

10.1063/5.0274886
페이지 상단으로 이동하기