Research Article  Open Access
JianLiang Gong, "Numerical Simulation of Coupling Thermal and Mechanical Influence on Submerged Nozzle in Solid Rocket Motor", International Journal of Aerospace Engineering, vol. 2020, Article ID 8795316, 9 pages, 2020. https://doi.org/10.1155/2020/8795316
Numerical Simulation of Coupling Thermal and Mechanical Influence on Submerged Nozzle in Solid Rocket Motor
Abstract
The coupling thermal and mechanical effect on submerged nozzles is important in the design of modern rockets upon thermal loading and aerodynamic pressure. In this paper, a simulation with the subroutine of nonuniform pressure and nonuniform heat transfer coefficient was conducted to study the thermostructural response of a submerged nozzle at the pressure 6 MPa and stagnation temperature 3200 K. Both the aerodynamic parameters and heat coefficients were obtained through analyzing the flow field. It was found that the thermal loading had an important influence on the stress of throat insert for the solid rocket motor (SRM). The hoop stress increases at first and then decreases with the increase of time for the throat insert. The ground hot firing test of SRM with a submerged nozzle was carried out. The experimental results showed that the structural integrity of the submerged nozzle is very normal during SRM operation. The present method is reasonable, which can be applied to study the thermostructural response of submerged nozzle for SRM.
1. Introduction
The SRM nozzle always acts as the energy transformation equipment where the chemical energy of the propellant turns into the kinetic energy of gas [1]. In order to shorten the length of a solid rocket motor, the submerged nozzle is adopted in high altitude solid rocket motor [2, 3]. It is classified as a submerged nozzle, as its entrance extends to within the combustion chamber of the SRM. In general, the submerged nozzle is always working under extremely harsh circumstances, such as thermal load, mechanical load, thermochemical ablation, and particleladed flow [4–7]. Hence, the structure integrity is of great importance for the normal operation of the nozzle. Hence, the simulation of the submerged nozzle is carried out to analyze the structure integrity, such as flow simulation, heat transfer simulation, and structure analysis [8].
A variety of recent works account for the growth in knowledge and techniques in the assessment of the structural behavior of SRM nozzles. Kumar et al. [9] undertook an extensive thermostructural analysis of composite structures, incorporating temperaturedependent properties, and thermal, thermochemical and mechanical loads, and subsequently performed a coupled thermostructural stress analysis of an SRM nozzle comprised of various orthotropic and isotropic materials. Morozov and Beaujardiere [10] investigated the dynamic thermostructural response of a composite rocket nozzle throat using the commercial finite element code ADINA. It is found that the dynamic response oscillates about the quasistatic response in all cases, and that, in general, the variance in stress magnitudes between the two solution techniques is significant. Tian et al. [11] studied the problems with gap design and contact stress about the SRM nozzle under thermal loading and pressure using the finite element method. Li et al. [12] had a focus on the thermal structural analysis of the nozzle throat insert; the influences of the thermoelastic stress are discussed by the computational software ABAQUS. Zheng [13] studied the thermostructural analysis and failure behaviors of carboncarbon composite throat insert. Hu et al. [14] determined the temperature and stress field of the conical trapped joint carboncarbon nozzle divergent section. Sun et al. [15] studied the thermostructural response of a typical nozzle with consideration of the structure gap based on the finite element method.YU et al. [16] considered the influence of the fit clearance of the throat insert on the deformation and stress distribution of the nozzle based on the threedimensional finite element model. The research shows that the fit clearance of throat insert does not change the overall distribution of stress, but can affect the stress on the matching surface. Cozart and Shivakumar [17] performed a finite element axisymmetric stress analysis of a 3D braided preform ablative composite rocket nozzle, incorporating a material ablation model. In the work of Yoo et al. [18], a kick motor nozzle utilizing spatially reinforced composites was subjected to finite element thermoelastic analyses following the development of a material model used to homogenize various spatial reinforcement architectures.
To sum up, a large amount of work has been done for investigating the thermostructural response. However, literature on the thermostructural response of the submerged nozzle at the initial stages of operation is scarce. This study focuses on the numerical simulation of the submerged nozzle under the condition of internal pressure and thermal loading. It addresses a method for analyzing the thermostructural response of nozzle using a finite element analysis program and discusses the development of a user subroutine which allows us to model the nonuniform pressure and nonuniform heat transfer coefficients on the wall. Furthermore, the ground hot firing test of SRM is carried out. Finally, some conclusions are drawn in the end.
2. Submerged Nozzle Model
2.1. Geometry
The threedimensional cyclic symmetry finite element analysis is performed for the submerged nozzle made of different composites. This nozzle consists of five substructures, namely, throat insert made from the punctured carboncarbon composite, the tape wound 2D silicaphenolic entrance insulator—which isolates the hot gas and metal case, the silicaphenolic liner—which isolates the hot and cooler substructures, divergent insulator made from the silicaphenolic, and metal case made from the titanium alloy, as shown in Figure 1. In Figure 1, the coordinate of axis is dimensionless by the nozzle length. Furthermore, in order to reduce the thermal stress of the throat insert, the clearance was set up in the internal interface of the throat insert. The values of interface AB, BC, DE, and EF are all 0.05 millimeters.
2.2. Assumptions
Some reasonable assumptions are considered as follows to simplify the geometry model:(1)The outer surface of the nozzle has no heat exchange(2)The pure gas steady gas is considered, the flow field is steady(3)The contact thermal resistance is totally ignored for the simulation(4)The complicated phenomena of erosion and pyrolysis behavior of the erosion and heat insulation materials are neglected(5)The radiation heat transfer is not considered
2.3. Material Properties
There are three kinds of materials for the submerged nozzle, which are titanium alloy, silicaphenolic, and carboncarbon composite. Silicaphenolic and carboncarbon materials are treated as homogeneous and orthotropic, and their properties are related to temperature. The properties of titanium alloy, silicaphenolic, and carboncarbon materials are shown in Tables 1–3, where represents temperature (), represents thermal conductivity (W.m^{1}.K^{1}), represents specific heat (W.kg^{1}.K^{1}), represents the coefficient of thermal expansion (K^{1}), represents Young’s modulus (MPa), represents shear modulus (MPa), represents Poisson’s ratio, and represents the density (kg/m^{3}).



2.4. Mesh
A properly sized mesh can generate more accurate results and reduce the computing resources for the thermostructural simulation. A 1/12th 3D symmetric model is used in this paper for the simulation. The cylindrical coordinate system is adopted to impose symmetric constraints on the symmetric surface. Figure 2 shows the mesh generation result of the thermostructural simulation. This model is also used to simulate the nonlinear and orthotropic material properties. The finite element model is generated by a structured grid. There are 9453 elements in the finite element mesh. Figure 2 shows the mesh generation result.
2.5. Boundary Conditions
The simulation for the submerged nozzle is subjected to both thermal and mechanical excitations during the period of analysis. To solve the problem of thermalstructure, both the fluid software and structure software are employed. The first stage entails the aerodynamics of flow on the surface of the nozzle. The second stage supplies the temperature and stress distribution for structure. The flow diagram of the simulation is presented in Figure 3.
On the one hand, in order to obtain the steady flow filed, the axisymmetric numerical simulation was carried out by a finite volume method, based on a pressure solver. A standard model was used to make the gas phase equation as a closed system. On the other hand, there exist two kinds of thermal boundary conditions for structure simulation, i.e., the adiabatic condition between the outermost structural steel and the air, and the forced convection from the inner flow. For thermalstructural simulation, the heat transfer coefficient along the wall should be provided by means of fluid software. The result is shown in Figures 4–6, respectively. Figure 4 shows the temperature distribution, gradually decreasing along the wall. Figure 5 shows the force convective coefficient between hot gas and wall. As can be seen, the maximum value of the forced convection coefficient is about 8266 W/(m^{2}∙K), and the minimum value is about 324 W/(m^{2}∙K). It is noticed that the variation curve of the heat transfer coefficient presents a peak at the upstream of the throat insert. Figure 6 shows the pressure distribution, gradually decreasing along the wall, because of the gas exhaust flow.
3. Results and Discussion
3.1. Thermal Loading
For the submerged nozzle, the distribution of temperature was obtained at times 28 s based on the finite element method on the conditions of thermal loading, as shown in Figure 7. The stagnation temperature is 3200 K. The transient thermal analysis is run to get the temperature profile of the nozzle at the end of 28 seconds. It indicates that the convection heat transfer between the gas and the throat insert is very obvious, and the temperature of this material increases greatly. Because the conductivity of the carboncarbon composite is high, and the heat transfer coefficient is large as shown in Figure 5. Secondly, the outer surface of the case is low, which is 300 K environmental temperature, because the conductivity of the liner is much lower than that of the throat insert. This plays an important role in the insulation of the case. Finally, the heat transfer depth of the entrance insulator and the divergent insulator is shallow, because these convective coefficients are low as shown in Figure 5.
3.2. Mechanical Loading
For the submerged nozzle, the distribution of Von Mises stress was obtained at 28 s based on the same method on the condition of pressure, as shown in Figure 8. The stagnation of temperature and pressure is 3200 K and 6 MPa, respectively. It can be seen that the position of the maximum Von Mises stress is the outer surface of the case, and the value is about 176 MPa. Secondly, the stress of the protective materials for the nozzle is small under the pressure, and the maximum of Von Mises stress for the throat insert is only 12 MPa.
3.3. Both Thermal and Mechanical Loading
3.3.1. Discussion about Strength
On the condition of the thermal loading and pressure, the distribution of Von Mises stress was obtained at time 28 s for the submerged nozzle, as shown in Figure 9. The stagnation of temperature and pressure is 3200 K and 6 MPa, respectively. It is interesting to note that the position of the maximum Von Mises stress is the front end of the cylindrical section of the case, and the maximum value is 176 MPa. Compared with the Figure 8, the maximum Von Mises stress increased obviously, the thermal expansion of the thermal protective materials is caused by the rise of the internal temperature. The contact interface of these materials is compressed from each other due to the inconsistency of the expansion deformation. It is shown that the thermal expansion of the protective materials has an important influence on the stress of the case.
The throat insert is the one subjected to the most severe thermal and mechanical loading, which provides an appropriate case for study. Furthermore, the hoop stress is the most important stress component. Figure 10 displays the hoop stress of the throat insert at different times. It is clear that the position of the maximum hoop compressive stress is upstream at the convergent section of the throat insert, because of the maximum heat transfer coefficient in the upstream of throat insert, the temperature difference is large, and the thermal expansion deformation of throat insert is obvious. Secondly, the position of the maximum hoop tensile stress is at the inner fillet of the throat insert. In order to avoid the damage of throat insert during the test, the radius of the inner fillet can be properly increased. Thirdly, the hoop stress increases at first and then decreases with time for the throat insert. Finally, during the working operation, the maximum hoop compressive stress of throat insert is 6.1 s, and its value is 103.9 MPa, which is less than 180 MPa of compress strength. The maximum hoop stress 34.8 MPa of throat insert is 20.1 s, and its value is 34.8 MPa, which is less than 115 MPa of tensile strength. Both maximum hoop compressive stress and maximum tensile stress yield the material strength.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
The contact stress of forward and backward interface in the throat insert has an important influence on the stress distribution and the safety of nozzle during the operation. In particular, both forward gap AB and backward gap EF are very important in Figure 1. Therefore, the value of forward gap AB in Figure 1 is 0.05 millimeters, and the one of backward value EF is 0.05 millimeters. Figure 11 shows the distribution of contact stress for the throat insert at time 28 seconds. It indicates that the value of contact stress is larger at the front face AB, and the value of contact stress is smaller at the back face EF. This is the reason that the gas temperature at upstream is much higher than the one at downstream for throat insert. In order to ensure the reliability of the throat insert, the interface gap of the front face of the throat insert can be properly increased, and the interface gap of the back face can be reduced.
3.4. Experimental Verification of Structural Integrity for the Throat Insert
On the condition of combustion chamber pressure 6 MPa and stagnation temperature 3200 K, the ground hot firing test of SRM with the submerged nozzle was carried out, and the working time was 28 s. In addition, the design values of interface AB, BC, DE, and EF, as shown in Figure 1, are all 0.05 millimeters. The structural integrity of the submerged nozzle is normal. Figure 12 shows the structure of the throat insert after the test of the submerged nozzle. During the SRM operation, the performance of the submerged nozzle yields the requirement of SRM. Meanwhile, the test value of hoop strain for a metal case is at time 28 s, and the computational value is , as shown in Figure 13, Point A. It can be seen that the computational value is equivalent to the test value. Furthermore, it can be shown that the method of thermostructure response of submerged nozzle developed in this paper is reasonable. The structure design and interface clearance of throat insert is valid. The distribution of stress field and temperature field for the throat insert are sound, which yields meet the material strength.
4. Conclusion
In this paper, the thermostructure response of threedimensional submerged nozzle under the condition of thermal loading and pressure was investigated. The ground firing test of SRM has been completed. Some conclusions can be drawn as follows:(1)Thermal loading is found to have the most dominating influence on the thermal stress of nozzle(2)The hoop stress of throat insert increases at first and then decreases with the increase of time(3)In order to ensure the reliability of the submerged nozzle, the design of throat insert clearance and structure is extremely important
Data Availability
As most of the data in this manuscript were related to trade secrets, I cannot provide them completely. In the future, if necessary, I can share some data with the reviewers or readers.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This article was funded by Xi’an Modern Chemistry Research Institute, China.
References
 P. F. Melia, “Flow and ablation patterns in Titan IV SRM aft closures,” in 31st Joint Propulsion Conference and Exhibit, p. 2878, San Diego, CA, USA, 1995. View at: Publisher Site  Google Scholar
 E. R. Joyce and M. P. Snyder, “Solid rocket motor design for a high altitude composite rocket,” in 45th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, p. 5083, Denver, CO, USA, 2009. View at: Publisher Site  Google Scholar
 M. Angelone, F. Mascanzoni, and C. Milana, “Technological and programmatic development of Zefiro 40 solid rocket motor,” in 48th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, p. 4213, Atlanta, GA, USA, 2012. View at: Publisher Site  Google Scholar
 B. C. Yang, A Theoretical Study of ThermoMechanical Erosion of HighTemperature Ablatives, The Pennsylvania State University, University Park, State College, PA, USA, 1992.
 Y. W. Guan, J. Li, and Y. Liu, “Ablation characteristics and reaction mechanism of insulation materials under slag deposition condition,” Acta Astronautica, vol. 136, no. 4, pp. 80–89, 2017. View at: Publisher Site  Google Scholar
 J. Li, M. F. Guo, X. Lv, Y. Liu, K. Xi, and Y. W. Guan, “Erosion characteristics of ethylene propylene diene monomer composite insulation by hightemperature dense particles,” Acta Astronautica, vol. 145, no. 2, pp. 293–303, 2018. View at: Publisher Site  Google Scholar
 Y. Liu, Y. W. Guan, and J. Li, “Insulator ablation modes in different impact conditions of alumina droplets onto wall surfaces,” Acta Astronautica, vol. 153, no. 6, pp. 138–145, 2018. View at: Publisher Site  Google Scholar
 Q. Li, P. J. Liu, and G. Q. He, “Fluidsolid coupled simulation of the ignition transient of solid rocket motor,” Acta Astronautica, vol. 110, no. 3, pp. 180–190, 2015. View at: Publisher Site  Google Scholar
 R. R. Kumar, G. Vinod, S. Renjith, G. Rajeev, M. K. Jana, and R. Harikrishnan, “Thermostructural analysis of composite structures,” Materials Science and Engineering A, vol. 412, no. 12, pp. 66–70, 2005. View at: Publisher Site  Google Scholar
 E. V. Morozov and J. Beaujardiere, “Numerical simulation of the dynamic thermostructural response of a composite rocket nozzle throat,” Composite Structures, vol. 91, no. 4, pp. 412–420, 2009. View at: Publisher Site  Google Scholar
 S. P. Tian, G. J. Tang, D. K. Li, and D. C. Yuan, “Gap design of solid rocket motor nozzle structures,” Journal of Propulsion Technology, vol. 26, no. 5, pp. 448–451, 2005. View at: Google Scholar
 S. L. Li, F. Zhang, and B. Xiong, “Influences factor analysis of thermal structure on motor throat,” Structure and Environment Engineering, vol. 40, no. 2, pp. 56–63, 2013. View at: Google Scholar
 Q. Zheng, ThermoStructural Analysis and Failure Behavior of C/C Composite Throat, Harbin Institute of Technology, Harbin, China, 2011.
 J. H. Hu, S. H. Meng, X. L. Chang, C. H. Xu, and H. Chen, “Temperature and stress field analysis of conical trapped joint C/C nozzle divergent section,” Journal of Solid Rocket Technology, vol. 35, no. 1, pp. 64–68, 2012. View at: Google Scholar
 L. Sun, F. T. Bao, N. Zhang et al., “ThermoStructural Response Caused by Structure Gap and Gap Design for Solid Rocket Motor Nozzles,” Energies, vol. 9, no. 6, p. 430, 2016. View at: Publisher Site  Google Scholar
 X. J. Yu, S. Yu, Y. Wang, and Y. Yang, “Thermostructure coupled computation for the influence of fit clearance on the stress distribution of composite nozzle,” in 2018 Joint Propulsion Conference, p. 4486, Cincinnati, OH, USA, 2018. View at: Publisher Site  Google Scholar
 A. B. Cozart and K. N. Shivakumar, “Stress analyses of a 3D braided composite ablative nozzle,” in 40th Structures, Structural Dynamics, and Materials Conference and Exhibit, p. 1277, St. Louis, MO, USA, 1999. View at: Publisher Site  Google Scholar
 J. S. Yoo, I. H. Cho, and C. G. Kim, “Thermoelastic analysis of a kick motor nozzle incorporating spatially reinforced composites,” Journal of Spacecraft and Rockets, vol. 40, no. 1, pp. 83–91, 2003. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 JianLiang Gong. 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.