1. Introduction
Liquid rocket engine suffers high temperature and pressure difference which is the reason for failure and life limitation, shown in
Figure 1. With the rapid development of reusable technology for liquid rocket engines(LREs), life prediction of the key component has to be considered after each launching procedure. There are many life prediction methods for liquid rocket engine thrust chambers, such as stress-strain methods, beam models, unified visco-plastic models, and finite element analyses.
H. J. Kasper [
1] developed an analytical method that modeled the cyclic creep failure phenomenon for the Space Shuttle Main Engine(SSME). The research provided advice for increasing the life of thrust chambers. A.K. Asraff [
2,
3] performed a stress analysis and cyclic and creep life prediction of a cryogenic rocket engine thrust chamber materials by ANSYS finite element analysis software package. The failure mode of stainless steel and copper alloys are creep and low cycle fatigue (LCF). The research also explored the different alloy plasticity models and details of cyclic stress analysis carried out for the double-walled rocket engine thrust chamber. W. Schwarz [
4] developed a model that quantitatively accounted for the crack failure mechanism of the thrust chamber’s hot-gas-side wall and verified it using 3D finite element simulation data. Tong Jun [
5] aimed at the problem of the muster and ablation on the inner wall of the thrust chamber throat section, a three-dimensional analysis of modeling finite element was conducted to calculate the steady creep and accelerated creep deformation of the structure in the high-temperature phase. Cheng Cheng [
6] investigated the thermo-mechanical response of channel wall nozzles under cyclic working loads, and the finite volume fluid-thermal coupling calculation method and the finite element thermal-structural coupling analysis technique are applied. A simplified nonlinear damage accumulation approach has been suggested to estimate the fatigue service life of the channel wall nozzle, which is obviously conservative to the Miner’s life. Zhang sheng [
7] established a structural model based on finite elements to evaluate the cooling channel of thrust chamber reliability. The parameter sensitivity was analyzed and demonstrated the dangerous point of the thrust chamber’s inner wall. M. Ferraiuolo [
8,
9] performed numerical investigations on the influence of the closeout geometry and materials on the number of cycles to failure of the thrust chamber. The impact of the creep phenomenon on the service life was also illustrated by Viscoplastic models. M. Ferraiuolo [
10] also researched a global linear model and a local nonlinear model to improve the accuracy of the numerical simulations. Citarella R. [
11] investigated the main failure phenomena that could occur during the thrust chamber’s service life using non-linear finite element analyses. The most critical areas appeared on the inner surface of the chamber, which was consistent with experimental tests from literature conducted on similar geometries. Di Liu [
12] developed a Finite Element Method (FEM) based on experimental data in order to evaluate the structural failure risk of the regenerative cooling thrust chamber inner wall. The multi-cycle thermo-structural analysis demonstrated the evolution of the stress-strain curve and damage analysis for dangerous points.
However, except for the beam model, most of the life prediction methods are quite difficult to conduct life analyses for because of complex mathematical expressions and fail to demonstrate the thinning of the inner wall [
13]. The beam model is a simplified method to reveal the life of reusable liquid rocket engines, which contains the Porowski clamped beam model, the creep-modified model, and the sandwich beam model [
14,
15].
The Porowski beam model is a simplified model for the thrust chamber inner walls, which considers the ligament as a clamped beam. This model was first proposed by J.S. Porowski and M. Badlani for the Space Shuttle Main Engine (SSME) and is still widely used in thrust chamber life prediction because of its convenience until now. Porowski [
16,
17] developed the clamped beam model to predict the life of the SSME main combustion chamber, which first offers the deformation and strain range of the coolant channel on the side of the gas wall and takes into account the plastic ratchet effect. The calculation results are verified through finite element analysis, and the results are relatively satisfactory. Badlani [
18] improved the previous clamped beam life model, increased the deformation caused by creep, and considered the impact of fatigue and ratchet damage in the life prediction. Combining the relevant parameters of SSME to calculate thinning of the coolant channel liner under one cycle, it’s almost consistent with the results of the finite element analysis, which verifies the effectiveness and accuracy of this method. The creep-modified beam model can estimate the life of the inner wall after a working cycle, and the life calculation results are more conservative than other methods.
M. Niino [
19], in order to improve the performance of high-pressure thrust chamber closeout and increase the life of the thrust chamber, developed the cold isostatic pressing (CIP) forming method. The Porowski beam model was used to predict the life of the thrust chamber, and the lifetime was 37 cycles for case A and 390 cycles for case B, respectively. The firing test was performed to verify the results. In-Kyung Sung [
20] analyzed the life of the SSME inner wall and calculated the temperature field, plastic strain and deformation by the Porowski beam model and other methods. The estimated life of the Porowski beam model is 67 cycles. Tao Chen [
21] used the creep-modified Porowski beam model to estimate the life of a reusable liquid engine thrust chamber in the throat section. The life prediction results offered a reference for the design of reusable thrust chambers. Marco Pizzarelli [
22,
23] developed an algebraic model for structural and life analysis of regeneratively cooled thrust chambers. The life of instability is calculated by the Porowski beam model. The calculation results are verified by experimental data from NASA Lewis Research Center.
Although the Porowski beam model is widely used in thrust chamber life prediction, life analysis of this model is infrequent in many articles. Life analysis is critical in reusable thrust chamber design and manufacture. The present work aims at the life analysis of LRE thrust chambers using the Porowski beam model by derivation and simulation.
2. Basic Equations and Procedures of the Porowski Beam Model
In the Porowski beam model, the ligament of the inner wall is considered a clamped beam, in which both ends are constrained. The inner wall of the thrust chamber suffers from high temperature and pressure differences between the hot-gas-side wall and coolant channel wall, which causes incremental damage and thinning after each working cycle. The simplified structure of the thrust chamber is shown in
Figure 2.
2.1. Inelastic Strain
The average temperatures of the inner wall ligament and closeout wall were denoted by T
i and T
0, respectively. The inelastic strain range of the Porowski model and creep-modified model are given by Equations (1) and (2). T
i and T
0 can be calculated from Equations (3) and (4). The coefficient of
is taken as 0.35 for the rectangular coolant channel in Equation (4).
During the ignition and steady-state working conditions, there is also a temperature drop across the ligament, which causes bending, because the inner wall ends are clamped. Conservatively assuming that all of the available elastic energy goes into plastic straining of the ligament, the correction due to thermally induced bending is given by:
Then, the total inelastic strain for the complete cycle is:
The coefficient of thermal expansion is , is the yield stress, is the creep relaxation stress, for Young’s modulus, for Poisson’s ratio. and are hot-gas-side wall temperature and coolant channel wall temperature, which can be obtained by heat transfer analysis, and is temperature drop across the inner wall ligament.
2.2. Deflection
Bending deflection , shear deflection and creep-induced deflection of the ligament are given by Equations (7) and (9) as follows.
The pressure difference between the coolant channel and thrust chamber is
, and
is the length of the coolant channel.
and
are the material constants,
is the firing time,
is the thickness of the ligament and
has been tabulated for various
-values in reference [
18].
The total deflection per cycle of the Porowski beam model and creep-modified model are
and
, respectively. This is accomplished with the assumption that the deflection during each firing cycle remains constant.
2.3. Ligament Deformation
Experimental evidence from NASA Lewis Research Center [
24,
25] shows that the distorted shape of the inner wall can be approximated by a linear variation in thickness thinning. Thus using the linear thinning model developed in the Porowski model, the thinning after N cycles,
, is then given by:
Based on the linear variation assumption, the thickness of the inner wall ligament can be described by
and
after N cycles which are given by Equations (11) and (12). The symbol
means the thickness of the end of the inner wall ligament, which increases during the firing cycle, and the symbol
means the thickness of the middle section of the inner wall ligament, which decreases during the firing cycle.
2.4. Fatigue and Creep Rupture Damage
According to the linear thinning model of the inner wall, the effective strain range in the section of minimum ligament for entering the fatigue curve is then:
The hoop strain of the inner wall ligament in the minimum section is calculated by:
The ultimate strength of inner wall material is , and is the effective strain range.
The average hoop strain
and the axial strain
in the minimum ligament section are given by:
2.5. Life Prediction
The life of the thrust chamber contains plastic instability life and fatigue-creep life according to the Porowski model based on different failure modes. Life of instability is given by the following Equation (18) and
is the ratio of
and
. The critical thickness at the life of instability is given by (18):
For some inner wall materials, the number of cycles at which ligament thinning stops
is related to the strain hardening parameter q. The following empirical criterion is obtained from reference [
17]:
Fatigue curve of thrust chamber inner wall material matched by Manson Universal Slopes (MUS) Equation.
The fatigue life of this material is
,
stands for the reduction in area.
and
for both materials are shown in
Section 3.2. The stress-to-rupture curves of OFHC copper and Narloy-Z using data matched using the MUS equation are shown in
Figure 3.
Therefore, the procedure of life prediction of the thrust chamber is shown in
Figure 4. Firstly, the parameters of the combustion chamber should be acquired.
4. Life Analysis
The life analysis of the thrust chamber contains pressure difference and temperature differences between the inner wall, structural parameters and material parameters. Except for pressure difference analysis, the rest of the chamber pressure difference is set at 18.0 MPa.
The failure of this LOX/Kerosene rocket engine is plastic instability, according to
Figure 4, so the deflection per cycle and the life of instability are analyzed.
The figure legend with (P) is the result of the Porowski beam model, and the figure legend with (C) stands for the creep-modified model. Two kinds of inner wall materials, OFHC and Narloy-Z, are analyzed as a comparison. The deflections of both materials calculated by the model are given in
Table 4.
4.1. Pressure Difference
The pressure difference range of this LOX/Kerosene liquid rocket engine thrust chamber set in this research is 10–30 MPa because of its full-flow staged combustion cycle. The relationship between pressure difference, deflection per cycle and life of instability are shown in the pictures as follows. It can be seen from
Figure 6 and
Figure 7 that the deflection and life of instability under the Porowski beam model approximately obey a linear relationship. The deflection of OFHC copper is higher than the Narloy-Z alloy. However, the life of instability seems different, and OFHC copper seems better than the Narloy-Z alloy. The failure mode of the thrust chamber ligament is instability, but the Narloy-Z alloy has better resistance against fatigue and creep failure. The following results may also be for the same reason.
4.2. Temperature Difference
The temperature difference between the ligament also influences the deflection per cycle and the life of instability. For regeneratively cooled thrust chamber ligament, changes in the temperature difference influence the total deflection and instability life. Usually, the thickness variation of the inner wall affects the heat transfer coefficient of the thrust chamber wall, and then the temperature difference of the inner wall changes. According to the heat transfer calculation results mentioned above, the temperature difference between the inner wall at the throat section is 135.41 K (243.73 °F). For different inner wall materials, the temperature difference is not the same, so the change interval of the inner wall temperature difference range is taken 100–500 K for analysis.
With the increase of the inner wall temperature difference, the total deflection decreases first and then increases because of the quadratic function of ΔT. Compared with Narloy-Z alloy, the variation in the total deflection of OFHC copper is more obvious in
Figure 8. From the results, it can be concluded that under the pressure difference of 18 MPa, the maximum instability life of OFHC and Narloy-Z thrust chamber is about 38 lifetimes in
Figure 9. The results have the same tendency for the two beam models for the life of instability because the deflectionδ
3 is quite small in the calculation. Therefore, controlling the temperature difference of the inner wall ligament is an important measure to prolong the life of the thrust chamber of reusable liquid rocket engines.
4.3. Thickness of Ligament
The inner wall suffers from high pressure and temperature difference during engine combustion. The thickness of the inner wall ligament is generally in the order of millimeters, which is extremely vulnerable to damage during each firing cycle. The thickness increase in the ligament intensifies the structural strength, but the temperature difference of the ligament increases due to heat transfer coefficient variation. During the analysis of instability life, the temperature difference caused by the variation of thickness is ignored.
Increasing the inner wall thickness significantly reduces the total deflection and increases the instability life in
Figure 10. The life of OFHC in the Porowski beam is higher than that of Narloy-Z, and the life of Narloy-Z is higher in the creep-modified model is shown in
Figure 11. The reason is that the change in the inner wall thickness has a great influence on the total deflection, and when other conditions remain invariant, the life of instability is inverse to the total deflection.
4.4. Width of Coolant Channel and Rib
The width of the coolant channel affects the life of the inner wall ligament in the thrust chamber, and the width of the coolant channel is usually in the order of millimeters. In the analysis of the total deflection and instability life of the inner wall, the interval of width is from 1 mm to 4 mm. The width of the rib only influences the instability life of the thrust chamber wall according to the Porowski beam model and the creep-modified model.
The total deflection increases with the width of the coolant channel while the instability life of the inner wall decreases. In the creep-modified model, the total deflection of the inner wall increases significantly when the groove width is 3.5 mm in
Figure 12, so the life of instability decreases rapidly in
Figure 13. The cause of this increase is the coolant width power in δ
3. With the increase of rib width, the instability life of the thrust chamber wall shows a tendency to decrease in
Figure 14, and the same reason could explain the similar results of the two models in
Figure 9. However, Porowski beam models are unable to analyze the relationship between rib width and deflection.
4.5. Rib Filling Ratio
The rib filling ratio, also called the width ratio of the rib and coolant channel which, is defined as the width of the coolant channel divided by the width of the rib,
. The rib filling ratio was referenced in many articles; however, no quantified analyses were performed. In this study, the life of instability, t
min, and t
max of the inner wall ligament were analyzed when the rib filling ratio changed. The asymptotic lines were deduced simply with the following equations.
From Equation (26), we know that only when , which means the width of the rib and coolant channel are the same, then the maximum thickness of the ligament is achieved at .
With the increase of the ratio of the rib and coolant channel, the effective strain range decreases in
Figure 15. The reason is that the ratio of t
max and t
min value is contained in the effect strain range according to Equations (13) and (14). Therefore, it can be seen that the ratio of t
max and t
min value has the same reduction with the effective strain range when the rib width ratio is increasing.
Figure 16,
Figure 17 and
Figure 18, the maximum ligament thickness of Narloy-Z is lower than OFHC, and the minimum ligament thickness is higher than OFHC under different cycles. This suggests that Narloy-Z has a slower rate of deformation than OFHC. Hence, Narloy-Z is more suitable for thrusting inner walls.
5. Conclusions
With the development of reusable liquid rocket engines, life prediction is receiving increasing attention in aerospace. This research performed a quantitative analysis of life prediction based on the Porowski beam model and the creep-modified model for a LOX/Kerosene rocket engine thrust chamber.
From the life analysis, we can draw a conclusion that the pressure and temperature difference, structural parameters and material parameters have a significant impact on the deflection per cycle and life of instability. With the increasing pressure difference, the deflection increases as well, and the life of instability decreases. The temperature difference between deflection and life seems to have a relationship with the quadratic function. During the design of reusable liquid rocket engine thrust chambers, pressure and temperature difference have to be strictly constrained to extend the lifetime. The life of instability goes up when the thickness of the inner wall ligament grows. When the width of the rib and coolant channel increase, the life of instability reduces clearly. The rib filling ratio plays an important role in ligament thinning; the asymptotic line of tmax, tmin and effective strain range appears when the rib filling ratio is augmented. This situation of widths of the same value in ribs and coolant channels should be avoided in structural design and manufacturing. From the analysis, Narloy-Z alloy is better than OFHC copper for thrust chamber inner wall material because of the lower deformation rate.
Life analysis conducted in this research offered an appropriate reference for the design and manufacturing of reusable liquid rocket engine thrust chambers.