Study on Pressure Characteristics and Its Evolution Law at the Ellipsoidal End Cover Pole of Cylindrical Explosion Containment Vessels

To explore the postposition of the maximum pressure at the pole of the ellipsoidal end cover of cylindrical explosion containment vessels and to reveal the mechanism of the load evolution, the experimental method was used to measure the pressure curve at the pole under different charges, and the numerical simulation method was used to analyze the evolution law of the explosion flow field within the end cover. The results show that the end cover pole was subjected to three types of pressure: the primary explosion wave, the secondary shock wave and the convergence wave. In addition, the pressure peaks increased in sequence. The evolution of the flow field in the end cover was affected by the amount of charge and the aspect ratio of the vessel. When the scaled distance due to a small charge increased, or when the aspect ratio of the vessel was reduced, the time interval between the convergence wave and the secondary shock wave at the end cover pole decreased gradually. When the scaled distance increased to 4.05 m·kg−1/3, the convergence wave at the pole superimposed on the secondary shock wave. As the aspect ratio of the vessel ranged from 1.75 to 2.50, the time interval between the two peaks was about 150 μs. However, if the aspect ratio was less than 1.40, the convergence wave and the secondary shock wave were fused through complex interaction.


Introduction
As a limiting device for blast wave and explosion products, the explosion containment vessel has wide applications in both military and civil fields. It was found that the maximum pressure of the end cover pole of the container is not the first pressure peak, but a subsequent one. W. H. Zhu et al. [1] used the TVD (Total Variation Diminishing) numerical method, combined with the analytical solution of a strong explosion, to analyze the internal load distribution of a cylindrical explosion container. They proposed that due to the interaction between waves, the pressure characteristics at the end cover pole were very complex, and multiple pressure peaks might be present. S. Kevorkian et al. [2] observed that, in the internal explosion experiment of a cuboid explosion container, the center of the cover plate at the far end of the explosion was affected by four simultaneous reflections by the chamber wall. It was suggested that the second pressure peak corresponding to the largest value appeared a few microseconds later than the first one. C. F. Zhao et al. [3] analyzed the internal explosion loading and the structural stability of the cylindrical reinforced concrete containment by numerical simulation. It was found that the explosion shock wave was repeatedly reflected in the vessel, and the dome of the container may have been subjected to the strongest impact. By the means of experiment and numerical simulation, C. Q. Wu et al. [4] studied the pressure distribution of an explosion in an irregular confined blast chamber. They found that the second pressure peak had a high intensity and even exceeded the first peak at gauges near the chamber wall or the corner. Z. Wang et al. [5] used a TNT (2,4,6-Trinitrotoluene) equivalent numerical model to analyze the explosion impact load inside the steel container, and found that due to the convergence effect of shock waves, the maximum load occurred in the pole area of the dome and all occurred in the second or third peak of the waveform. V. Pickerd et al. [6], through experiments and numerical simulation, found that the second peak was much higher than the first peak, which was on the gauges at the edge and corner of the cubic explosive vessel. W. B. Wang et al. [7] used numerical models with different zoom ratios to analyze the pressure load distribution of a single-layer spherical lattice shell structure, in which the maximum pressure at the pole always appeared at the second peak. Z. R. Li et al. [8] experimentally measured the pressure, which was loaded on the cover plate of a cylindrical vessel, with an aspect ratio of 2:1, indicating that the pressure load on it presented multiple wave peaks, and the maximum peak of the pressure load at the center of the cover plate was three times higher than the first. In addition, A. F. Demchuk [9], V. M. Kornev [10], R. R. Karpp [11], T. A. Duffey [12], Q. Dong [13][14][15], W. X. Liu [16] and other scholars have conducted fruitful studies on the distribution of the internal load, vibration response, failure analysis and the design criteria of the explosion containment vessel.
To explore the postposition phenomenon of the maximum pressure at the end cover pole, but also to reveal its load evolution mechanism, this paper introduced an experimental result on the pressure load of the pole of the ellipsoidal end cover and a numerical simulation on the evolution and distribution of explosive flow field over various charges and aspect ratios. This study was intended to provide references for the study of the distribution law of explosion, which loads in the container, and the design of the anti-detonation strength of the explosion containment vessel.

Experimental Setup
The explosion containment vessel and the measures used in this experiment are shown in Figure 1a. The total length L of the vessel was 166.4 cm, the thickness 2δ of the vessel wall was 2.2 cm, the inner major axis diameter D a of the ellipsoidal end cover was 80 cm, the inner minor axis diameter D b was 40 cm and the aspect ratio (L-4δ)/(D a ) was 2.025. The internal shape of the upper end cover of the container and the ellipsoidal container bottom were the same. Therefore, a pressure-measuring point P pole was set at the bottom of the container, as shown in Figure 1b. When the spherical charge was detonated in the center of the container, its pressure curve should have been the same as that of the upper end of the container. Four kinds of compressed TNT spherical charges with 27 g, 65 g, 106 g and 155 g were set, with a density of 1.50 g·cm −3 . The TNT charge ball was hoisted in the center of the container and the central detonation method was adopted.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 2 of 10 explosion in an irregular confined blast chamber. They found that the second pressure peak had a high intensity and even exceeded the first peak at gauges near the chamber wall or the corner. Z. Wang et al. [5] used a TNT (2,4,6-Trinitrotoluene) equivalent numerical model to analyze the explosion impact load inside the steel container, and found that due to the convergence effect of shock waves, the maximum load occurred in the pole area of the dome and all occurred in the second or third peak of the waveform. V. Pickerd et al. [6], through experiments and numerical simulation, found that the second peak was much higher than the first peak, which was on the gauges at the edge and corner of the cubic explosive vessel. W. B. Wang et al. [7] used numerical models with different zoom ratios to analyze the pressure load distribution of a single-layer spherical lattice shell structure, in which the maximum pressure at the pole always appeared at the second peak. Z. R. Li et al. [8] experimentally measured the pressure, which was loaded on the cover plate of a cylindrical vessel, with an aspect ratio of 2:1, indicating that the pressure load on it presented multiple wave peaks, and the maximum peak of the pressure load at the center of the cover plate was three times higher than the first. In addition, A. F. Demchuk [9], V. M. Kornev [10], R. R. Karpp [11], T. A. Duffey [12], Q. Dong [13][14][15], W. X. Liu [16] and other scholars have conducted fruitful studies on the distribution of the internal load, vibration response, failure analysis and the design criteria of the explosion containment vessel.
To explore the postposition phenomenon of the maximum pressure at the end cover pole, but also to reveal its load evolution mechanism, this paper introduced an experimental result on the pressure load of the pole of the ellipsoidal end cover and a numerical simulation on the evolution and distribution of explosive flow field over various charges and aspect ratios. This study was intended to provide references for the study of the distribution law of explosion, which loads in the container, and the design of the anti-detonation strength of the explosion containment vessel.

Experimental Setup
The explosion containment vessel and the measures used in this experiment are shown in Figure  1(a). The total length L of the vessel was 166.4cm, the thickness 2δ of the vessel wall was 2.2cm, the inner major axis diameter Da of the ellipsoidal end cover was 80cm, the inner minor axis diameter Db was 40cm and the aspect ratio (L-4δ)/(Da) was 2.025. The internal shape of the upper end cover of the container and the ellipsoidal container bottom were the same. Therefore, a pressure-measuring point Ppole was set at the bottom of the container, as shown in Figure 1(b). When the spherical charge was detonated in the center of the container, its pressure curve should have been the same as that of the upper end of the container. Four kinds of compressed TNT spherical charges with 27g, 65g, 106g and 155g were set, with a density of 1.50g⋅cm -3 . The TNT charge ball was hoisted in the center of the container and the central detonation method was adopted.

Experimental Results
Based on previous numerical simulation studies [5,7], it was generally believed that the initial pressure stage at the end cover pole was divided into two parts: first, the explosion shock wave inside the charge was encountered at the pole and the reflected pressure was generated. Then, the convergence collision wave of the shock wave was formed at the pole, and the convergence pressure was significantly higher than the first pressure peak. Z. R. Li et al. [8] believed that the central area of the end cover was first affected by the blast wave, and then by the superposition of the wall reflection wave, the secondary shock wave and then by the cover reflection wave. In this test, the time history curve of pressure was measured at the measuring point P pole at the end cover pole, as shown in Figure 2. However, what differed from previous research results was the presence of three increasing pressure peaks in a row at the initial pressure stage of the end cover pole. For example, in Figure 2b, at the initial stage (700 µs-1500 µs), there are three pressure peaks, 0.99 MPa, 2.41 MPa and 13.89 MPa, respectively. In order to explore this phenomenon and reveal the mechanism of load evolution, numerical simulation was used for the analysis and discussion in the following article.

Experimental Results
Based on previous numerical simulation studies [5,7], it was generally believed that the initial pressure stage at the end cover pole was divided into two parts: first, the explosion shock wave inside the charge was encountered at the pole and the reflected pressure was generated. Then, the convergence collision wave of the shock wave was formed at the pole, and the convergence pressure was significantly higher than the first pressure peak. Z. R. Li et al. [8] believed that the central area of the end cover was first affected by the blast wave, and then by the superposition of the wall reflection wave, the secondary shock wave and then by the cover reflection wave. In this test, the time history curve of pressure was measured at the measuring point Ppole at the end cover pole, as shown in Figure  2. However, what differed from previous research results was the presence of three increasing pressure peaks in a row at the initial pressure stage of the end cover pole. For example, in Figure 2(b), at the initial stage (700μs-1500μs), there are three pressure peaks, 0.99MPa, 2.41MPa and 13.89MPa, respectively. In order to explore this phenomenon and reveal the mechanism of load evolution, numerical simulation was used for the analysis and discussion in the following article.

Modeling and Parameters
As the TNT charges detonated in the center of the container, the distribution of the internal flow field was symmetrical. Therefore, a two-dimensional model could be used to describe the distribution evolution law. In order to ensure the calculation accuracy of the flow field inside the container, the mesh size of the solid was limited to less than 0.6cm, and the mesh size of the air was 0.1cm. For the convenience of calculation, the flange, bolt and other structures independent of the internal flow field were simplified in the model, and the 1/4 containment vessel model was adopted as shown in Figure  3. The mapping method was used to map the calculation results of the one-dimensional explosive shock wave to the two-dimensional model, which ensured the calculation accuracy while

Modeling and Parameters
As the TNT charges detonated in the center of the container, the distribution of the internal flow field was symmetrical. Therefore, a two-dimensional model could be used to describe the distribution evolution law. In order to ensure the calculation accuracy of the flow field inside the container, the mesh size of the solid was limited to less than 0.6 cm, and the mesh size of the air was 0.1 cm. For the convenience of calculation, the flange, bolt and other structures independent of the internal flow field were simplified in the model, and the 1/4 containment vessel model was adopted as shown in Figure 3. The mapping method was used to map the calculation results of the one-dimensional explosive shock wave to the two-dimensional model, which ensured the calculation accuracy while maintaining simpler calculations [17][18][19]. The symmetrical boundaries of the two planes, both xoz and yoz, were set, respectively. The initial volume fraction method, *INITIAL_VOLUME_FRACTION_GEOMETRY, was used to fill the explosive at the origin of the one-dimensional air grid.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 4 of 10 maintaining simpler calculations [17][18][19]. The symmetrical boundaries of the two planes, both xoz and yoz, were set, respectively. The initial volume fraction method, *INITIAL_VOLUME_FRACTION_GEOMETRY, was used to fill the explosive at the origin of the one-dimensional air grid.
. The body and end cover of the explosion containment vessel were made of Q345R steel and the vessel wall was 2.2cm thick. The deformation of these solids under pressure loads had little effect on the flow field [20,21]. Therefore, the solid deformation was ignored, and *MAT_RIGID could be selected for the material model. The material density ρ0 was 7.85g⋅cm -3 , the Young's modulus E was 209GPa and the Poisson's ratio ν was 0.28.
A high energy explosive combustion model was used to describe the explosive, with the keyword *MAT_HIGH_EXPLOSIVE_BURN. The combined JWL (Jones-Wilkins-Lee) equation represents the pressure of the explosion product. The form is as follows [22]: The units of A, B, and E are the same as the pressure, and R1, R2, ω, and V are dimensionless values, where E is the energy density and V is the relative volume. When entering the parameters, the initial energy density E0 and the initial relative volume V0 should be determined, and the specific parameters are shown in Table 1.
Regarding air as an ideal gas, the parameters were set to C0 = C1 = C2 = C3 = C6 = 0. The polytropic exponent γ is usually set to 1.4, thus C4 = C5 = γ -1 = 0.4 [22]. The initial density ρ0, the initial energy density E0 and the initial relative volume V0 of air were set to 1.29g⋅cm -3 , 0.25MPa and 1.0, respectively. The parameters are shown in Table 2. The body and end cover of the explosion containment vessel were made of Q345R steel and the vessel wall was 2.2 cm thick. The deformation of these solids under pressure loads had little effect on the flow field [20,21]. Therefore, the solid deformation was ignored, and *MAT_RIGID could be selected for the material model. The material density ρ 0 was 7.85 g·cm −3 , the Young's modulus E was 209 GPa and the Poisson's ratio ν was 0.28.
A high energy explosive combustion model was used to describe the explosive, with the keyword *MAT_HIGH_EXPLOSIVE_BURN. The combined JWL (Jones-Wilkins-Lee) equation represents the pressure of the explosion product. The form is as follows [22]: The units of A, B, and E are the same as the pressure, and R 1 , R 2 , ω, and V are dimensionless values, where E is the energy density and V is the relative volume. When entering the parameters, the initial energy density E 0 and the initial relative volume V 0 should be determined, and the specific parameters are shown in Table 1.  The material model of air was *MAT_NULL, and its state equation was *EOS_LINEAR_POLYNOMIAL. The equation can be written as follows: Regarding air as an ideal gas, the parameters were set to C0 = C1 = C2 = C3 = C6 = 0. The polytropic exponent γ is usually set to 1.4, thus C4 = C5 = γ −1 = 0.4 [22]. The initial density ρ 0 , the initial energy density E 0 and the initial relative volume V 0 of air were set to 1.29 g·cm −3 , 0.25 MPa and 1.0, respectively. The parameters are shown in Table 2.  [23,24].  Figure 4 shows the explosion pressure curve of a 65 g TNT charge in the containment vessel with an aspect ratio of 2.025. The scaled distance between this charge and the pole was 2.01 m·kg −1/3 . The numerical simulation and experimental results were well matched and the maximum pressure was reached at the third peak.   Figure 4 shows the explosion pressure curve of a 65g TNT charge in the containment vessel with an aspect ratio of 2.025. The scaled distance between this charge and the pole was 2.01m⋅kg -1/3 . The numerical simulation and experimental results were well matched and the maximum pressure was reached at the third peak.  Figure 5 illustrates the nephogram of the internal explosion in the container under this working condition. After the charge explosion at 0μs, the explosion shock wave reached the vicinity of the vessel wall at about 210μs, and then because of the blocking effect of the wall being moved along the vessel wall towards the bottom, the shock wave formed the reflection shock wave. At about 495μs, a clear Mach reflection and a three-wave intersection phenomenon appeared. It is worth noting that the reflection waves at the detonation center torus were gradually converging toward the center of the container. At 530μs, the reflection shock waves from the torus successfully converged at the center of the container and generated the secondary shock wave that propagated to all directions of the container. At the same time, the Mach wave continued to propagate along the end cover inner wall after passing the corner between the barrel and the end cover, as shown in figure 765μs. At around 840μs, the first explosion shock wave had already acted on the end cover, and the first pressure peak appeared at the pole of the end cover, as shown in Figure 4 and Figure 5, respectively. The secondary shock wave propagated faster in the air after the primary explosion wave was compressed, while the Mach wave propagated slower due to the wall retardation. Therefore, before the Mach wave at 1135μs, the secondary shock wave reached the end cover pole, forming the second pressure peak. At the following moment 1295μs, Mach reflection waves from all directions converged and collided at the pole of the end cover, creating a third pressure jump. Over time, the new shock wave formed at the pole and travelled around.

a. 210μs
b. 495μs c. 530μs d. 765μs  It is worth noting that the reflection waves at the detonation center torus were gradually converging toward the center of the container. At 530 µs, the reflection shock waves from the torus successfully converged at the center of the container and generated the secondary shock wave that propagated to all directions of the container. At the same time, the Mach wave continued to propagate along the end cover inner wall after passing the corner between the barrel and the end cover, as shown at 765 µs (Figure 5d). At around 840 µs, the first explosion shock wave had already acted on the end cover, and the first pressure peak appeared at the pole of the end cover, as shown in Figures 4 and 5, respectively. The secondary shock wave propagated faster in the air after the primary explosion wave was compressed, while the Mach wave propagated slower due to the wall retardation. Therefore, before the Mach wave at 1135 µs, the secondary shock wave reached the end cover pole, forming the second pressure peak. At the following moment 1295 µs, Mach reflection waves from all directions converged and collided at the pole of the end cover, creating a third pressure jump. Over time, the new shock wave formed at the pole and travelled around.
appeared at the pole of the end cover, as shown in Figure 4 and Figure 5, respectively. The secondary shock wave propagated faster in the air after the primary explosion wave was compressed, while the Mach wave propagated slower due to the wall retardation. Therefore, before the Mach wave at 1135μs, the secondary shock wave reached the end cover pole, forming the second pressure peak. At the following moment 1295μs, Mach reflection waves from all directions converged and collided at the pole of the end cover, creating a third pressure jump. Over time, the new shock wave formed at the pole and travelled around.

Influence of Scaled Distance on Internal Explosion Flow Field
Observing Figure 2, it can be found that with a small amount of charge, the interval time between the convergence wave peak and the secondary shock wave peak is smaller, as shown in Figure 2(a). To investigate this phenomenon, an explosion containment vessel model with a same aspect ratio of 2.025 and a smaller charge of 8g was used for numerical simulation. In conversion, the scaled distance between the charge and the pole was 4.05m⋅kg -1/3 . Figure 6 shows the pressure nephogram at the critical moment.
At 1400μs after the TNT explosion, the primary explosion wave reached near the end cover pole and caused a small pressure jump, while the secondary shock wave formed by the torus reflection wave convergence was already propagating to the pole. During the propagation process, the secondary wave and the shock wave propagating along the end cover converged toward the pole at the same time. For example, at 2000μs the two shock waves are almost equal to the wave at pole. At 2105μs, the impact of the secondary shock wave is encountered at the pole, and the shock wave propagating along the inner wall of the end cover is also about to reach the pole. Subsequently, only at 2130μs was the convergence wave at the pole compressed inward to the extreme, at which point the pressure reached the maximum. Of particular note was that, compared with the typical pressure flow field evolution, which was described in the previous section, it merely took about 25μs from the rising edge of the secondary shock wave to reach the peak of the convergence wave, instead of about 160μs. It was obvious that during the charge of 8g TNT, the secondary wave and the convergence wave arrived at the pole almost at the same time, so they almost merged into one peak in the pressure curve. In other words, it can be called a superposition wave of the convergence wave and the secondary shock wave. From the pressure curve of the amplified superimposed wave, it can also be seen that, in fact, the two waves were not completely merged, but when the pressure was not unloaded after the secondary shock wave, the convergence wave reached the pole as shown in Figure  7  In order to explore the changing law of the interval time between the convergence wave and the secondary shock wave at different scaled distances, the charge was set between 8g to 155g, and the corresponding scaled distance ranged from 1.51-4.05m⋅kg -1/3 . The simulation results are shown in Figure 7(b). When the scaled distance was small, the interval time between the convergence wave and the secondary shock wave was large, therefore three dominating peaks at the initial stage could be observed in the pressure curve, and the maximum value appeared in the third peak-the convergence wave peak. As the scaled distance grew, the interval between the convergence wave and

Influence of Scaled Distance on Internal Explosion Flow Field
Observing Figure 2, it can be found that with a small amount of charge, the interval time between the convergence wave peak and the secondary shock wave peak is smaller, as shown in Figure 2a. To investigate this phenomenon, an explosion containment vessel model with a same aspect ratio of 2.025 and a smaller charge of 8 g was used for numerical simulation. In conversion, the scaled distance between the charge and the pole was 4.05 m·kg −1/3 . Figure 6 shows the pressure nephogram at the critical moment.

Influence of Scaled Distance on Internal Explosion Flow Field
Observing Figure 2, it can be found that with a small amount of charge, the interval time between the convergence wave peak and the secondary shock wave peak is smaller, as shown in Figure 2(a). To investigate this phenomenon, an explosion containment vessel model with a same aspect ratio of 2.025 and a smaller charge of 8g was used for numerical simulation. In conversion, the scaled distance between the charge and the pole was 4.05m⋅kg -1/3 . Figure 6 shows the pressure nephogram at the critical moment.
At 1400μs after the TNT explosion, the primary explosion wave reached near the end cover pole and caused a small pressure jump, while the secondary shock wave formed by the torus reflection wave convergence was already propagating to the pole. During the propagation process, the secondary wave and the shock wave propagating along the end cover converged toward the pole at the same time. For example, at 2000μs the two shock waves are almost equal to the wave at pole. At 2105μs, the impact of the secondary shock wave is encountered at the pole, and the shock wave propagating along the inner wall of the end cover is also about to reach the pole. Subsequently, only at 2130μs was the convergence wave at the pole compressed inward to the extreme, at which point the pressure reached the maximum. Of particular note was that, compared with the typical pressure flow field evolution, which was described in the previous section, it merely took about 25μs from the rising edge of the secondary shock wave to reach the peak of the convergence wave, instead of about 160μs. It was obvious that during the charge of 8g TNT, the secondary wave and the convergence wave arrived at the pole almost at the same time, so they almost merged into one peak in the pressure curve. In other words, it can be called a superposition wave of the convergence wave and the secondary shock wave. From the pressure curve of the amplified superimposed wave, it can also be seen that, in fact, the two waves were not completely merged, but when the pressure was not unloaded after the secondary shock wave, the convergence wave reached the pole as shown in Figure  7  In order to explore the changing law of the interval time between the convergence wave and the secondary shock wave at different scaled distances, the charge was set between 8g to 155g, and the corresponding scaled distance ranged from 1.51-4.05m⋅kg -1/3 . The simulation results are shown in Figure 7(b). When the scaled distance was small, the interval time between the convergence wave and the secondary shock wave was large, therefore three dominating peaks at the initial stage could be observed in the pressure curve, and the maximum value appeared in the third peak-the convergence wave peak. As the scaled distance grew, the interval between the convergence wave and At 1400 µs after the TNT explosion, the primary explosion wave reached near the end cover pole and caused a small pressure jump, while the secondary shock wave formed by the torus reflection wave convergence was already propagating to the pole. During the propagation process, the secondary wave and the shock wave propagating along the end cover converged toward the pole at the same time. For example, at 2000 µs the two shock waves are almost equal to the wave at pole. At 2105 µs, the impact of the secondary shock wave is encountered at the pole, and the shock wave propagating along the inner wall of the end cover is also about to reach the pole. Subsequently, only at 2130 µs was the convergence wave at the pole compressed inward to the extreme, at which point the pressure reached the maximum. Of particular note was that, compared with the typical pressure flow field evolution, which was described in the previous section, it merely took about 25 µs from the rising edge of the secondary shock wave to reach the peak of the convergence wave, instead of about 160 µs. It was obvious that during the charge of 8 g TNT, the secondary wave and the convergence wave arrived at the pole almost at the same time, so they almost merged into one peak in the pressure curve. In other words, it can be called a superposition wave of the convergence wave and the secondary shock wave. From the pressure curve of the amplified superimposed wave, it can also be seen that, in fact, the two waves were not completely merged, but when the pressure was not unloaded after the secondary shock wave, the convergence wave reached the pole as shown in Figure 7a.
In order to explore the changing law of the interval time between the convergence wave and the secondary shock wave at different scaled distances, the charge was set between 8 g to 155 g, and the corresponding scaled distance ranged from 1.51-4.05 m·kg −1/3 . The simulation results are shown in Figure 7b. When the scaled distance was small, the interval time between the convergence wave and the secondary shock wave was large, therefore three dominating peaks at the initial stage could be observed in the pressure curve, and the maximum value appeared in the third peak-the convergence wave peak. As the scaled distance grew, the interval between the convergence wave and the secondary shock wave decreased gradually, until the convergence wave was superimposed on the secondary shock wave. When the scaled distance reached 4.05 m·kg −1/3 , there were only two obvious wave peaks in the pressure curve and the maximum pressure appeared on the second peak-the superimposed wave peak.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 10 the secondary shock wave decreased gradually, until the convergence wave was superimposed on the secondary shock wave. When the scaled distance reached 4.05m⋅kg -1/3 , there were only two obvious wave peaks in the pressure curve and the maximum pressure appeared on the second peakthe superimposed wave peak.

Influence of Aspect Ratio on Internal Explosion Flow Field
In addition to the scaled distance of charge, factors such as the aspect ratio of the container also affect the internal flow field distribution [25]. Therefore, explosion containment vessel models with different aspect ratios were designed for numerical simulation and the charge was uniformly set to 65g.
The evolution of the internal explosion flow field with a container aspect ratio of 1.00 was taken as an example for analysis. As shown in Figure 8 and Figure 9(a), at about 215μs after the charge explosion, the primary explosion wave and its reflection occurred simultaneously at the pole and at the central section torus of the containment vessel. Then, the reflection point moved from the torus and the pole to the corner of the container, and finally the reflection points met near the corner, at 335μs. At the same time, all the reflection shock waves converged toward the center of the container and compressed to the extreme at about 535μs, then outward to form a series of secondary shock waves. The ellipsoidal end cover had radians, which helped the reflection shock waves to converge toward the focal point. Therefore, the reflection shock waves did not converge at the center of the container, but formed two upper and lower convergence points. Due to the influence of the location of the convergence point, at about 760μs the end cover pole first received the secondary shock wave and the reflection was formed at the central section torus. It is noteworthy that, along the inner wall of the end cover converging towards the pole, there was still a shock wave propagating. As a result of the pressure of this secondary shock wave, its reflected wave at the pole and the convergence wave were in the same order of magnitude, and they formed a complex interaction at the pole, resulting in the formation of three peaks between 760μs and 860μs, as shown in Figure 9(a). After that, the reflection shock waves from the vessel wall converged again at about 1505μs, and the formed tertiary shock wave reached the end cover pole at 1595μs. Moreover, through the complex interaction of the tertiary wave itself, its reflected wave and the convergence wave, a maximum pressure peak was finally formed at around 1645μs.

Influence of Aspect Ratio on Internal Explosion Flow Field
In addition to the scaled distance of charge, factors such as the aspect ratio of the container also affect the internal flow field distribution [25]. Therefore, explosion containment vessel models with different aspect ratios were designed for numerical simulation and the charge was uniformly set to 65 g.
The evolution of the internal explosion flow field with a container aspect ratio of 1.00 was taken as an example for analysis. As shown in Figures 8 and 9a, at about 215 µs after the charge explosion, the primary explosion wave and its reflection occurred simultaneously at the pole and at the central section torus of the containment vessel. Then, the reflection point moved from the torus and the pole to the corner of the container, and finally the reflection points met near the corner, at 335 µs. At the same time, all the reflection shock waves converged toward the center of the container and compressed to the extreme at about 535 µs, then outward to form a series of secondary shock waves. The ellipsoidal end cover had radians, which helped the reflection shock waves to converge toward the focal point. Therefore, the reflection shock waves did not converge at the center of the container, but formed two upper and lower convergence points. Due to the influence of the location of the convergence point, at about 760 µs the end cover pole first received the secondary shock wave and the reflection was formed at the central section torus. It is noteworthy that, along the inner wall of the end cover converging towards the pole, there was still a shock wave propagating. As a result of the pressure of this secondary shock wave, its reflected wave at the pole and the convergence wave were in the same order of magnitude, and they formed a complex interaction at the pole, resulting in the formation of three peaks between 760 µs and 860 µs, as shown in Figure 9a. After that, the reflection shock waves from the vessel wall converged again at about 1505 µs, and the formed tertiary shock wave reached the end cover pole at 1595 µs. Moreover, through the complex interaction of the tertiary wave itself, its reflected wave and the convergence wave, a maximum pressure peak was finally formed at around 1645 µs.  Through setting up containment vessels with different aspect ratios in the range of 1.00 to 2.50 for numerical simulation, the change of the interval time between the convergence wave and the secondary shock wave was investigated as shown in Figure 9(b). When the aspect ratio of the container was greater than 1.75, the interval between the two peaks was practically the same, about 150μs. When the aspect ratio was lower than 1.60, the interval between the two peaks decreased gradually. Until the aspect ratio came to 1.40, due to the same intensity of the secondary shock wave and its reflected wave at the pole as the convergence wave, it was impossible to distinguish the convergence wave peak from the secondary wave peak after the complex interaction. When the aspect ratio was less than 1.40, the convergence wave and the secondary shock were fused to form the complex fusion wave.

Conclusions
In this paper, the pressure load time history curve at the end cover of the ellipsoid of the cylindrical explosion containment vessel was experimentally measured under the charges of 27g, 65g, 106g and 155g TNT. Different from the previous results, it was found that there were three pressure peaks successively rising at the initial stage of the pressure at the end cover pole. In order to verify the experimental results and explore the mechanism, the method of fluid-solid coupling was used to simulate and analyze the flow field changes in the container with a two-dimensional numerical model. A series of conclusions were obtained through the numerical simulation and analysis. The study in this paper shows that the pressure peak at the end of the container increases continuously, which must be controlled in the actual container. This context expounds important practical Through setting up containment vessels with different aspect ratios in the range of 1.00 to 2.50 for numerical simulation, the change of the interval time between the convergence wave and the secondary shock wave was investigated as shown in Figure 9b. When the aspect ratio of the container was greater than 1.75, the interval between the two peaks was practically the same, about 150 µs. When the aspect ratio was lower than 1.60, the interval between the two peaks decreased gradually. Until the aspect ratio came to 1.40, due to the same intensity of the secondary shock wave and its reflected wave at the pole as the convergence wave, it was impossible to distinguish the convergence wave peak from the secondary wave peak after the complex interaction. When the aspect ratio was less than 1.40, the convergence wave and the secondary shock were fused to form the complex fusion wave.  Through setting up containment vessels with different aspect ratios in the range of 1.00 to 2.50 for numerical simulation, the change of the interval time between the convergence wave and the secondary shock wave was investigated as shown in Figure 9(b). When the aspect ratio of the container was greater than 1.75, the interval between the two peaks was practically the same, about 150μs. When the aspect ratio was lower than 1.60, the interval between the two peaks decreased gradually. Until the aspect ratio came to 1.40, due to the same intensity of the secondary shock wave and its reflected wave at the pole as the convergence wave, it was impossible to distinguish the convergence wave peak from the secondary wave peak after the complex interaction. When the aspect ratio was less than 1.40, the convergence wave and the secondary shock were fused to form the complex fusion wave.

Conclusions
In this paper, the pressure load time history curve at the end cover of the ellipsoid of the cylindrical explosion containment vessel was experimentally measured under the charges of 27g, 65g, 106g and 155g TNT. Different from the previous results, it was found that there were three pressure peaks successively rising at the initial stage of the pressure at the end cover pole. In order to verify the experimental results and explore the mechanism, the method of fluid-solid coupling was used to simulate and analyze the flow field changes in the container with a two-dimensional numerical model. A series of conclusions were obtained through the numerical simulation and analysis. The study in this paper shows that the pressure peak at the end of the container increases continuously, which must be controlled in the actual container. This context expounds important practical

Conclusions
In this paper, the pressure load time history curve at the end cover of the ellipsoid of the cylindrical explosion containment vessel was experimentally measured under the charges of 27 g, 65 g, 106 g and 155 g TNT. Different from the previous results, it was found that there were three pressure peaks successively rising at the initial stage of the pressure at the end cover pole. In order to verify the experimental results and explore the mechanism, the method of fluid-solid coupling was used to simulate and analyze the flow field changes in the container with a two-dimensional numerical model. A series of conclusions were obtained through the numerical simulation and analysis. The study in this paper shows that the pressure peak at the end of the container increases continuously, which must be controlled in the actual container. This context expounds important practical significance for the safety evaluation and protection design of explosive vessels. The main conclusions of this paper are summarized as follows: a. It was found by the explosion experiment of different charges that, at the initial stage, the ellipsoidal end cover pole of the cylindrical explosion vessel was subjected to three types of shock waves, namely the primary explosion wave, the secondary shock wave and the convergence wave, and the pressure peak values increased in sequence. The results were verified by numerical simulation.
b. In the numerical simulation, the scaled distance between the detonation center and the end cover pole was changed by adapting the charge amount. Numerical simulation was carried out under different charges and it was found that, when the scaled distance increased due to the reduction of the charge, the interval between the convergence wave and the secondary shock wave at the end cover pole decreased gradually. When the scaled distance increased to 4.05 m·kg −1/3 , the convergence wave at the pole was superimposed on the secondary shock wave.
c. By changing the length of the vessel wall, with different aspect ratios at a 65 g charge, the numerical simulation of models showed that the decrease of the aspect ratio of the container also led to the decrease of the interval between the convergence wave and the secondary shock wave at the end cover pole. When the aspect ratio was reduced to below 1.40, the convergence wave and secondary shock wave fused through complex interaction.