Abstract

To explore the effect of different bend pipes on the propagation characteristics of premixed methane-air explosion, the experimental explosion pipe system and numerical model were established. By adopting the comparative analysis of experiments and numerical modeling, it conducted researches on the overpressure evolution of gas explosion shock wave in pipes with different bends and obtained the expressions of attenuation coefficient of shock wave overpressure. The results showed that the change of pipe direction accelerated the attenuation of gas explosion shock wave. The propagation attenuation of gas explosion in the bend pipe was mainly affected by the bending angle and initial peak overpressure before bending. With the increase of the bending angle, the attenuation coefficient of gas explosion shock wave gradually increased. For the same bending angle, the attenuation coefficient of gas explosion shock wave increased with the increase of gas volume. The obtained coupling relationships between attenuation coefficient, bending angle, and initial peak overpressure before bending were useful for estimating the overpressure value after the bend. The results presented in this paper have important significance for the assessment of structures that have been damaged in the mine laneway of gas exploration accidents, further enriching the gas exploration spread theory.

1. Introduction

Gas explosions are one of the main hazards in coal mining, gas transportation, and related fields [1, 2]. Pipes or roadways that have different bend pipes are widely used in these industries as gas-flow channels. Gas explosions usually occur in such regions and obviously have a negative impact on industrial safety [3]. Studies on the effect of different bend pipes on premixed methane-air explosion in pipes are helpful in the design and layout of gas-explosion-suppression systems in coal mines and play an important role in preventing gas explosions or effectively mitigating the damage.

Numerous theoretical studies have examined the propagation of gas explosion shock waves in pipes with different bend pipes. Yang and Li [4, 5] established a mathematical model of bending based on the propagation laws of shock waves by means of a theoretical derivation. Zhou [6] analyzed the formation mechanism and propagation characteristics of impinging airflow with the theories of gas dynamics and fluid mechanics and obtained the relationships between the initial energy of impinging airflow in the turning and bifurcation roadway. Gao et al. [7] studied the propagation of the explosion wave and pointed out that the destruction of the gas explosion was mainly caused by the high pressure of high-speed airflow and shock waves. The results of the former are representative to some extent, but the influence of the mass force and wall friction is ignored in the derivation, and the models have certain limitations. In contrast, the latter has mainly been used to determine the head loss of liquids in agricultural machinery and petrochemical pipelines, but is not suitable for gas.

The propagation of gas explosion in different pipeline structures or with different obstacles has been studied extensively through experiments and numerical simulations [811]. Lin et al. [1215] conducted a series of studies on the influence of obstacles, bifurcated pipes, single bend, U-shaped pipe and Z-shaped pipe, and T-shaped pipes in confined spaces on the shock wave and flame wave of a gas explosion. Yu et al. [1619] studied the influences of the number, size, shape, and position of obstacles on the propagation of a blast wave by means of experiments and numerical simulations. Sun and Zheng [20, 21] qualitatively studied the effects of pipe diameters and sizes on the combustion and explosion characteristics of premixed gases using AutoGas software. Yin et al. [22] studied the propagation performance of a gasoline and air two-phase detonation wave in pipes with variable cross-sections through numerical simulations and found that the detonation wave produced an eddy current at the inflection point of the variable cross-section area and exhibited quenching. Jia et al. [23, 24] studied the gas exploration shock wave spread law via a one-way bifurcation pipeline and obtained computational formulas for the shock wave attenuation coefficients of branch pipelines and straight pipelines and the shock wave shunt coefficients of branch pipelines. Frolov et al. [25, 26] studied the deflagration to the detonation process in U-shaped structural pipes, while Zhu et al. [27, 28] found that the turning and bifurcation structures produce turbulence and studied the mechanism of flame acceleration. Obviously, the abovementioned studies have focused on the propagation characteristics of pressure waves according to built-in obstacles, bifurcated pipes, and pipes with different structures. However, there are few studies on the relationships between the attenuation rate of overpressure propagation, bending angle, and initial overpressure. The influence of the bending angle on gas explosion propagation and the peak overpressure attenuation coefficient is still unclear.

Therefore, the experimental explosion pipe system with a bend was established, and premixed methane-air explosion experiments were conducted. The overpressure propagation characteristics of three different premixed methane-air volumes and seven different bend pipes were studied experimentally. The aim was to determine the premixed methane-air explosion propagation laws in pipes with different bend pipes. The results can provide important theoretical guidelines for the prevention and control of gas explosions in different bend pipes.

2. Experimental Study on Shock Wave of Gas Explosion via Different Bend Pipes

In the premixed methane-air area, the spread of shock wave will be coupled with, and interact with, the combustion wave. In the area containing only air, the gas combustion finishes, and the shock wave decreases. According to the principles of explosion mechanics and hydromechanics, the propagation of gas explosion pressure in a pipe varies according to the structure of the pipe. When the shock wave of gas explosion propagates through the bend points, the reflected wave will form. The purpose of this paper is to study the propagation characteristics and attenuation law of the shock wave with different bend pipes under different conditions.

2.1. Experimental System

The explosion experimental device is shown in Figure 1. The system components mainly include an explosion pipe system, air-gas distribution system, ignition system, and high-frequency data acquisition system. The explosion pipe system consists of a blasting chamber with a volume of 0.18 m3 and a pipe of diameter 180 mm. The pipe is made of 16Mn steel, which has high temperature resistance, corrosion resistance, and pressure resistance up to 20 MPa.

The experimental pipe measures 17 m in length and consists of two parts, A and B. Part A is filled with a 9.5% concentration of premixed methane-air combustible gas, and part B contains only air.

The air-gas distribution system is a key component of gas explosion experiments and mainly consists of a methane cylinder, vacuum pump, air compressor, circulating pump, vacuum meter, and regulating valve. The methane cylinder is filled with high-purity methane gas (99.99% purity). The vacuum pump (2XZ-2 rotary vane vacuum pump, Shanghai Hujing Industrial Pump Plant) is shown in Figure 2(a); this ensures the vacuum status of the blasting chamber. The air compressor (550-24 (1524) oil-free air compressor, Shanghai Jaguar Compressor Manufacturing Co., Ltd.) is shown in Figure 2(b). After each experiment, the air compressor conducts high-pressure ventilation of the explosion pipe network system for 30 min to discharge residual exhaust gas from the pipeline. The circulating pump shown in Figure 2(c) is operated for 20 min to ensure the uniform and sufficient mixing of methane. A precise vacuum manometer is used to control the methane and air inlets accurately. The principle of the gas distribution is to adopt Dalton’s law of partial pressures. Dalton’s law points out that the pressure of ideal gas mixture () is equal to the sum of the partial pressures () of each component gas, namely,

where is the total pressure of mixed gas, is the partial pressure of each component gas, and is the volume fraction of each component gas, %.

According to Dalton’s law of partial pressures, the partial pressure values of each component gas are calculated. Then, each gas component is filled into the device according to the value. After the gas preparation is completed, start the circulating pump, so as to ensure that the gas is mixed evenly and in a static state.

The ignition system is mainly composed of ignition electrodes, fuses, wires, and regulators, as shown in Figure 2(d). The ignition end remains closed when the other end is open. The ignition system is located at one end of the pipe and controls the ignition electrode in the middle of the flange plates at the front end of the pipe via the ignition device. Electric sparks generated by a 36 V direct current are used for ignition.

The high-frequency data acquisition system includes dynamic data acquisition units, pressure sensors, a controller for the flammable gas deflagration/detonation testing device, and a working host. The dynamic data acquisition unit (TST5206, see Figure 2(e)) has a sampling rate range of 0–20 MSPS and an accuracy of 0.1% FS. The pressure sensor model (CYG1401) has a measuring range of 0-3 MPa and an accuracy of 0.5% FS. The signals acquired by the sensors are transmitted to the working host by a 32-way dynamic data acquisition unit. The controller of the flammable gas deflagration/detonation test device, which was developed by the Northeast University Safety Engineering Research Center, is shown in Figure 2(f). Its main functions are to display the vacuum degree of the pipeline, start and stop the vacuum pump, start and stop the circulating pump, and control the ignition device.

2.2. Experimental Scheme

When shock wave of gas explosion propagates at the bend, its propagation attenuation law is not only related to the change of bending angle but also related to initial overpressure of shock wave before bending. In the experiment, experimental devices of various bend pipes are built to study the relationships between the overpressure of explosion shock wave, bending angle, and initial overpressure before bending.

According to the bend characteristics of underground roadways in mines, seven bending angles, 30°, 45°, 60°, 90°, 120°, 135°, and 150°, are designed in the experiment. Figure 3 shows the experimental device diagrams of 30°, 60°, 90°, and 135° bend angles. As shown in Figure 4, six pressure sensors are arranged along the central line of the pipe. The equivalent distances of the measuring points from the explosion source are listed in Table 1. The peak overpressures at measuring points P1, P2, ⋯, P5 are expressed as , , ⋯, , respectively.

The main steps in the gas explosion experiments are the installation and debugging of the pipeline system, gas tightness test, gas distribution and charging, circulation mixing, ignition, and data acquisition. In the experiments, the temperature is 298 K, and atmospheric pressure is 101.325 kPa.

This paper derives the propagation and attenuation law of shock wave after gas explosion accident in general air districts, so the concept of the overpressure is that the absolute pressure of shock wave removes the standard atmospheric pressure, which is a parameter characterizing shock wave propagation. The shock wave overpressure is measured in Pascals. The experiments mainly focus on shock wave attenuation laws under different initial overpressures and distinct bend pipes. Therefore, 0.23m3, 0.28m3, and 0.31m3 of combustible gas are added to the pipe at a 9.5% concentration under the same bent configuration to change initial peak overpressure of shock wave. Three tests are conducted for each experimental condition, and the mean value is taken as the result for that configuration, and a total of 21 successful experiments are required.

2.3. Experimental Results and Analysis
2.3.1. Characteristics of Shock Wave Propagation via Bend Pipes

Figure 5 shows the overpressure-time variation curves of measuring points P1, P2, P3, P4, and P5 under 60° bend angle and 0.28m3 gas volume. P3 has a peak overpressure at 0.047 s before bending. After the shock wave propagates for a certain distance, another peak overpressure appears. It can be seen that explosion shock wave reflects at the corner, and the reflected shock wave propagates backwards, which shows that the bend structure has a certain influence on the reflection of explosion shock wave.

According to the propagation process of explosion shock wave, the peak overpressure of measuring point P1 is 0.4749 MPa. With the increase of propagation distance, the peak overpressure at P2 increases to 0.9833 MPa; then, the peak overpressure at P5 decreases to 0.3732 MPa. The peak overpressure of explosion shock wave presents a rising and then descending trend, and other operating conditions show the same change trend. In addition, negative pressure occurs at the outlet of the pipe (measuring point P5).

Peak overpressure at P3 is 0.7162 MPa, and that at P4 is 0.5552 MPa. After the 60° bend angle, the peak overpressure attenuates by 22.48%. Compared with 30° and 150° bend angles, the peak overpressure attenuates by 12.73% and 51.29% under the same gas volume, respectively, which means that attenuation amplitude of peak overpressure is different with different bend angles. In addition, the initial overpressure before bending also has certain influence on the attenuation of shock wave. For the 60° bend angle, when the initial pressure before bending is 0.6026 MPa and 0.8245 MPa, respectively, the peak overpressure attenuates by 18.86% and 21.86%, respectively. It can be seen that attenuation amplitude of peak overpressure varies with initial overpressure before bending.

Due to the influence of local resistance at the bend, the total energy decreases, and the change of the pipe direction accelerates the attenuation of shock wave. It can be seen that propagation attenuation of gas explosion in bend pipes is mainly affected by the bend angle and initial overpressure before bending. Table 2 shows experimental peak overpressures of gas explosion at each measuring point when gas volume is 0.23m3, 0.28m3, and 0.31m3, respectively, in different bend angles.

2.3.2. Attenuation Characteristics of Explosion Shock Wave in Different Bend Pipes

The attenuation degree of explosion shock wave at the bend is described by the attenuation coefficient, which is equal to the ratio of peak overpressure before bending () and peak overpressure after bending (). The attenuation coefficient of peak overpressure at the bend is determined by the following formula, namely,

According to equation (2), the attenuation coefficient of explosion shock wave at each bend is calculated as shown in Table 3. It can be seen from the table that peak overpressure of gas explosion has an obvious downward trend at the bend, and the downward range increases obviously with the increase of angle. The peak overpressure before bending also has certain influence on the attenuation coefficient; that is, it increases with the increase of initial peak overpressure.

As can be seen from Table 3 and Figure 6, due to the different reflection effects of bend pipes, peak overpressures of gas explosion are different before and after bending. For the pipes with the same bend angles, initial peak overpressure before bending and peak overpressure after bending both increases with the increase of gas volume. With the increase of bending angle, the attenuation coefficient of explosion shock wave gradually increases. For the same bend angle, the attenuation coefficient of explosion shock wave increases with the increase of gas volume, but the increase is small. From 30° to 90°, the average value of the attenuation coefficient increases from 1.130 to 1.354, and the increase amplitude is relatively gentle. From 90° to 150°, the average value of the attenuation coefficient increases from 1.354 to 2.046, with a steep increase.

2.3.3. Coupling Relationships between Attenuation Coefficient, Initial Peak Overpressure, and Bending Angle

In order to further explore the relationships between the attenuation coefficient of peak overpressure, bending angle, and initial peak overpressure at the bend, the coupling relationships between them are established based on a large number of experimental data.

Assuming that the functional relation between attenuation coefficient , initial peak overpressure , and bending angle can be obtained:

where , , and are the unsolved coefficients and exponentials of the nonlinear model.

Taking the natural logarithm of both sides of equation (2), we can get the following result:

According to the experimental parameters and result data in Table 3, equation (4) is used as the mathematical model between the attenuation coefficient of peak overpressure, initial peak overpressure, and bending angle. Using the MATLAB software for multiple nonlinear regression, the fitting coefficients are obtained: , , and . By substituting above coefficients and indexes into equation (4), the attenuation coefficient can be described as

Figure 7 shows the comparison curves between the experiment value and model fitting value of the attenuation coefficient. It can be seen from the curve that the relative error between the experimental value and model fitting value is within ±10%. Therefore, equation (5) can be used as a formula to solve the ventilation network in the later catastrophic process.

3. Numerical Simulations of Gas Explosion Shock Waves via Different Bend Pipes

Using the FLUENT software, the shock wave overpressure of gas explosion spread law in pipes with different bend pipes was investigated numerically. The experimental results were incorporated into the numerical simulations, and mathematical models were established under the experimental conditions. The attenuation degree of shock wave overpressure at different angles was analyzed, and the attenuation formula obtained from the experiment in Section 2 was verified.

3.1. Numerical Models and Simulation Methods

In a mathematical simulation, the gas explosion process can generally be treated as the thermal expansion process of an ideal gas, and the gas dynamics can be captured by the equations of mass conservation, momentum conservation, energy conservation, and component equilibrium equations [29]. Generally, the Reynolds-averaged Navier-Stokes equations are used as the governing equations [30, 31]. The main equations are as follows: where is the density, is the static pressure, is the particle velocity, is the dynamic viscous coefficient, is the enthalpy, is the fuel mass fraction, is the rate for volume combustion, is the turbulent kinetic energy, and is the heat source term.

Turbulence is an important factor in explosive gas combustion. Using the turbulence model, the turbulent kinetic energy equation and the turbulent kinetic dissipation rate equation are as follows [32]: where is the effective viscous coefficient; is the laminar viscous coefficient; is the turbulent viscous coefficient; is the Kronecker operator; , , , , and are the turbulent model constants, respectively; and is obtained as follows:

The explosion of combustible gases is a turbulent combustion process accompanied by complex chemical reactions. It is necessary to consider both turbulent effects and chemical reaction dynamics. Therefore, the EBU-Arrhenius combustion model is adopted [33]. The turbulent flame velocity is usually calculated by the flame front model, and the turbulent flame velocity is given by

The turbulence intensity and turbulence characteristic scale are calculated as follows:

(14)

where is the model constant, usually 0.52, is the laminar flame velocity [m/s], is the heat transfer coefficient of temperature [m2/s], is the turbulence time scale, is the time scale of the chemical reaction, is the flame stretching factor (when there is no stretching, ), and is a constant, usually 0.37.

In the numerical simulations, the finite volume method is adopted. The SIMPLEC algorithm [34] for pressure–velocity coupling is selected, and the subrelaxation method is adopted to control the convergence rate.

3.2. Analysis of Shock Wave Overpressure Propagation in Pipes

In order to verify the experimental results, seven geometric models of bend pipes are established in total, including 30°, 45°, 60°, 90°, 120°, 135°, and 150°. The size of the pipe is exactly the same as the experimental system. The pipe is closed at one end and open at the other end. Within the range of m (, 10 m, and 11 m) at the closed end, it is filled with premixed methane-air combustible gas with a concentration of 9.5%; that is, the gas volume is 0.23m3, 0.28m3, and 0.31m3, respectively. The ignition source is at the closed end, and initial ignition temperature is 1400 K. Two monitoring sections are set up at monitoring point A before bending and monitoring point B after bending to monitor the overpressure variation law of explosion shock wave and airflow velocity with time, to analyze the attenuation degree of peak overpressure of shock wave at different angles, and to verify the attenuation formula obtained in Section 2.3 experiment. Both monitoring points A and B are 1 m away from the bend.

Figure 8 is the simulation diagram of explosion pressure and airflow velocity propagation at , 2.0 ms, 5.4 ms, 9.8 ms, 12 ms, and 15.1 ms under the 90° bend angle and 0.23m3 gas volume. At the monitoring point A, peak overpressure is 0.613 MPa, and peak airflow velocity is 605.43 m/s. The explosion shock wave and impact airflow propagate along the straight pipe before bending, and the explosion shock wave is the positive shock wave. When , the explosion shock wave goes through a 90° right bend angle. Under the influence of the bend, the overpressure of shock wave is reflected and refracted to different degrees, forming oblique shock wave. At the bend, the overpressure of shock wave rises sharply, reaching a peak of 1.31 MPa, then reflection occurs, and the overpressure is superimposed. Part of the shock wave enters the pipe after bending and attenuates. The impact airflow collides with the wall surfaces, the airflow velocity changes, eddy currents appear, and the kinetic energy decreases. As the explosive airflow constantly rushes into the pipe after bending, the airflow velocity constantly changes, and a local eddy zone is formed at the bend.

As shown in Figures 8(e) and 8(f), with the continuous forward propagation of explosion shock wave, the front shock wave, and explosion airflow resume parallel motion, explosion shock wave changes from oblique shock wave to positive shock wave. At the monitoring point B, peak overpressure at the wave front is 0.553 MPa, and peak airflow velocity is 545.84 m/s.

Figure 9 shows the relationships between explosion overpressure and airflow velocity at different moments at the monitoring points A and B before and after the bend. There are two peak overpressures at the monitoring point A, the first peak overpressure is 0.613 MPa, and the second peak is 0.578 MPa, and it is clear that the second peak is due to the reflection of overpressure at the bend. The impact airflow velocity generated by the explosion creates the local eddy zone at the bend, which consumes a large amount of kinetic energy, and the airflow velocity decreases sharply. The airflow velocity increases slightly after reflecting back from the wall.

3.2.1. Comparison of Numerical Simulation and Experimental Results

Peak overpressure and attenuation coefficient of explosion shock wave with different gas volumes are obtained by simulating seven bend pipes with different angles of 30°, 45°, 60°, 90°, 120°, 135°, and 150°. Figure 10 shows the comparison of numerical simulation, experimental results, and model fitting values (calculated with equation (5)). It can be seen that the three results are in good agreement. Some scholars have concluded that as long as the relative error between simulation and experimental results is between -47% and+ 47%, the results obtained by numerical simulation can be applied in the engineering field [35]. Table 4 shows the relative error between the two. It can be seen that the maximum relative error is 7.13%, within ±47%, which indicates that the numerical model and mesh generation are reasonable. This further verifies the correctness of the attenuation coefficient formula obtained in the Section 2.3.

4. Conclusions

(1)The explosion shock wave not only reflects but also refracts at the bend. The local resistance produced by the bend has a certain impact on the total energy reduction of gas explosion, and the change of the pipe direction accelerates the attenuation of shock wave(2)Propagation attenuation of gas explosion in the bend pipe is mainly affected by the bend angle and initial peak overpressure before bending. Peak overpressure of gas explosion in the pipe has an obvious downward trend, and the downward trend increases with the increase of angle. Initial peak overpressure also has a certain influence on the attenuation coefficient of gas explosion; that is, the larger the initial peak overpressure is, the more obvious the attenuation coefficient of gas explosion is(3)The attenuation degree of explosion shock wave can be characterized by the attenuation coefficient in different bend pipes. With the increase of the bend angle, the attenuation coefficient of explosion shock wave gradually increases. For the same bend angle, the attenuation coefficient increases with the increase of gas volume. The coupling relationships between the attenuation coefficient of peak overpressure, bending angle, and initial peak overpressure at the bend are obtained. This allows the peak overpressure value after a bend to be estimated, which is useful for determining the source location and strength of an explosion in practical applications(4)The propagation and attenuation law of explosion shock wave are affected by the bent structure of the pipe, which presents certain complexity in the developing process. It is of great significance to find the law of gas explosion shock wave in the pipe for determining the influence scope of gas explosion disaster and disaster prevention and mitigation in coal mines

Data Availability

The data used to support the findings of the study are included in the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Nos. 51874009, 52074012, 51774014, 51874007, and 51704012), the Opening Project of State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology) (No.KFJJ20-11M), the Young Elite Scientists Sponsorship Program by China Association for Science and Technology (No.2018QNRC001), the Postdoctoral Fund of Anhui Province (No. 2017B212), and the Natural Science Foundation of the Anhui Higher Education Institutions of China (Nos. KJ2019A0133 and KJ2018A0099).