The Inﬂuence of Pipe Bending Curvature on H 2 -O 2 Gaseous Detonation Wave Front

: The inﬂuence of different bend curvatures on the detonation wave propagation was analyzed by an advanced numerical simulation system. The mechanism of propagation properties is revealed by cellular structure, internal and external boundary pressure distribution, propagation process of detonation wave and chemical reaction. The cellular structure and detonation wave front of bend with different curvature are very different. The simulation results show that the detonation wave with regular cell structure propagating through the curved parts induces detonation cell size increased by diffraction near the inner wall while detonation reﬂected on the bottom surface resulting in decrease of cell size. Detonation wave was affected by the rarefaction wave and compression wave in the bent pipe. The pressure distribution of the bend shows that the peak pressure in the 45 0 curvature is the largest, which should be paid more attention in industrial design. The chemical reaction could indicate the propagation characteristics of detonation wave, and different propagation characteristics have different proﬁles of chemical components.


Introduction
Gaseous detonation is very complicated. It is usually considered to be composed of induced shock waves followed by a reaction zone [1]. Gaseous detonation technology has a good application prospect in engine [2][3][4][5][6]. Most of the experimental studies and numerical simulation in the gaseous detonation wave focus on cell structure evolution and the influence of initial conditions in a straight pipe [7][8][9][10][11][12][13][14].
In daily life, combustible gas pipelines often have the connection sections of bends which is more complex, while experimental research and numerical simulation of the propagation mechanism in pipe have been reported in recent years, but there are still many deficiencies to be further studied [15][16][17][18]. The cellular structure evolution of detonation wave in bends is more complex, when the detonation wave go through the bend, different curvature leads to different cellular structure, pressure distribution and chemical reaction morphology. Study on the evolution process of detonation cellular structure in bends is helpful to understand the mechanism of interaction among the induced shock, chemical reaction zone and wall geometry. Numerical simulation is a good approach to understand the specific physical and chemical mechanisms which can put forward some constructive opinions in industrial pipeline design, prevent and restrain accidents as well as protect the property and life safety of country and citizen. There is a lot of literature following this subject. Edwards [19] studied the shock wave propagation in 90-degree bend by experiment and points out the relationship between the shock wave velocity and the curvature of the pipe. Frolov [20,21] investigated the detonation wave propagation characteristics of the U-bend tubes by experiment and numerical simulation Uchida [22] explored a comparative study on detonation wave propagation in 90 • straight tube and bend and points out two locations in pipe where pressure values are greater than those in the straight tube. One 2 of 22 location is in the outer wall surface of bend and the other appears in the downstream of the bend exit because of the effect of transverse waves. Kudo, Y. [23] studied the effects of initial pressure and radius of curvature on oblique detonation wave in a rectangular bend. Wang, L.Q. [24] showed the detonation behavior of detonation wave with different diameters and shapes. Sugiyama, Y. [25] numerically investigated the characteristics of detonation wave propagating in two-dimensional curved channel and discusses the limit of stable propagation. Two propagation modes are discovered, one is stable propagation mode, the other is a kind of repeated decay and initiation process model. Shi, X. [26] researched the propagation characteristics and different combustion modes of hydrogen/air mixture in a closed chamber. Gwak, Min-cheol [27] performed the characteristics of DDT (Detonation Deflagration Transition) in S bend by numerically and the acceleration or inhibit of the flame is caused by complex constraints of geometrical boundary in S tube. Zhenhua Pan [28] analyzed three detonation propagation modes and the effect of initial pressure in a helical detonation channel. Otsuka [29] performed numerical simulation on detonation wave propagation in the U-bend and found that detonation quenching near the entrance of U tube while-initiation in the U tube exit. It is found that the mall pipe use of the U-bend with small channel width and curvature radius can induce fast DDT. Liang, Z. [30] conducted experiments on the structural response of detonation in the 90-degree bend. S. Gallier [31] showed no difference between numerical simulation and experimental results in hydrogen-oxygen-argon mixtures. Goodwin [32] studied the deflagration to detonation transition (DDT) in rectangular tubes with different obstruction ratios. Two mechanisms are revealed, one with DDT occurring due to the creation of a hot spot in a gradient of reactivity that forms behind a Mach reflection and another in which DDT arises from energy focusing as shocks converge. Mark Short [33] examined numerically the influence of detonation front curvature on the cellular structure and the evolution dynamics from unstable detonation to stable detonation. They found that for sufficiently wide arcs, the angular speed and structure of the detonation front was influenced only by a flow region near the inner arc surface known as the detonation driving zone (DDZ). Smaller inner arc radii and lower activation energies increase the range of stable solutions. Li Jian [34] used Additive Runge-Kutta methods to research the change of detonation cellular structure and proved that the stiffly accurate additive Runge-Kutta methods can capture the discontinuity well and describe the detonation complex wave configurations exactly, especially the typical characteristics such as triple wave structure, cellular cell, diffraction and reflection. Tomotaka Niibo [35] studied the effect of mesh resolution and high order accuracy scheme on detonation cellular structure. With the decrease of mesh size and the higher order accuracy scheme, irregular cellular structure is obtained and unburned gas pockets and circulation structures caused by Kelvin-Helmholtz instability appear behind the detonation front. Xueqiang Yuan [36] analyzed the diffraction of a cellular detonation wave by numerical and theoretical method. They found that diffracted detonation can sustain coupling for a certain distance when propagating over curved wall and curved wall can promote reinitiation and even make detonation propagate stably without decoupling. Zhenjuan Xia [37] simulated the effects of curvature radius, domain width and stagnation pressure on propagation process of the detonation wave. The results indicate that the decoupling and reinitiation occur repeatedly during the unstable propagation mode of the RDW, and the angular velocities of leading shock wave vary accordingly. Xia [38] carried out a numerical study on the rotating detonation wave in plane-radial structure and compared with the experimental results. Due to the influence of curvature, the detonation wave is strengthened near the outer concave boundary and weakened near the inner convex one.
In conclusion, most cellular structures of detonation waves propagating in the bend were obtained by experiment [39][40][41], mainly using the 90 • angle. Meanwhile the different mechanisms of propagation of detonation waves and change rule in different bending angles are not revealed. The numerical simulation is an important method to explore deto-nation mechanism. Therefore, this paper solves these problems by numerical simulation of gas detonation.

Numerical Model
In this paper, the phenomenon of detonation reflection is described by Euler equation [42] without considering dissipation effects such as viscosity, heat conduction and diffusion. An Additive Runge-Kutta (ARK) time step method is used to integrate the governing equations to avoid rigidity of chemical sources. A 9 species, 48 reaction mechanism is adopted. The reacting species are H 2 , O 2 , H, O, OH, HO 2 , H 2 O 2 , H 2 O and Ar. A fifth-order Weighted Essentially Non-Oscillatory scheme is used to integrated on convective terms. Either with commercial software or self-developed software, the grid subdivision is crucial for discrete area in calculation. At the same time, the grid subdivision is closely related to the calculation format. Thus, the grid subdivision quality has direct influence on the accuracy of the calculation result. The object discussed in this paper is smooth bend with different curvatures and the structure is not regular, which is different from wedge surface and treated as series of micro yuan of wedge surface with continuously changing angle. This is a more complex structure. The wall boundary conditions are difficult to give out in high precision calculation with rectangular grid format and rules. Thus, the body-fitted curve coordinate system is adopted and the curve coordinate will be transformed to the isometric partitioning coordinates in calculation. This can simplify the treatment to geometric boundary conditions. In the curve coordinates, the reaction Euler equation is as follows:

∂Ũ ∂τ
where J is Jacobi determinant; ζ x , ζ y , η x , η y are transforming coefficients of curvilinear coordinates and calculating coordinates, respectively; others are similar to traditional Euler equations. The whole bent pipe is composed of a turning angle section, a straight pipe section before the turning angle and a straight pipe section after the turning angle. The turning angle among θ = 15 • ,20 • ,30 • ,40 • ,45 • , 60 • respectively. The channel width is five times of the cell width and divided into1500 × 300 grids with size of 0.1 mm. The initial conditions are set the same as [43]: the concentration is stoichiometric H 2 /O 2 mixture diluted with 70% argon. The numerical simulation we consider is the evolution of the cellular detonation wave propagating through realistic smooth pipe bends of different angle in two space dimensions at an initial pressure and temperature of 6.67 kPa and 298 K, respectively. In order to simplify the computational initialization process, a cellular detonation is located on the left. The entrance is extrapolation boundary, the exit is Neumann boundary and other boundaries are solid wall. Figure 1 shows a schematic view of the bend.

Finite Difference Schemes
The Steger warming vector flux splitting method is used to decompose the physical fluxes F and G [44], and then the semidiscrete scheme of the governing equations can be written as follows: where, 1/ 2, represent the numerical flux on grid boundary and can be obtained by the WENO reconstruction with the grid values of Fij, Gij. A addtive Runge-Kutta scheme method [45] is used to couple the Euler equations to the chemical reactions. According to the work of Araújo, Kennedy C.A., Inmaculada Higueras [46][47][48], additional Runge-Kutta explicit implicit algorithm can be used to handle the following ordinary differential equations:

Finite Difference Schemes
The Steger warming vector flux splitting method is used to decompose the physical fluxes F and G [44], and then the semidiscrete scheme of the governing equations can be written as follows: where,F ± i±1/2,j ,Ĝ ± i±1/2,j represent the numerical flux on grid boundary and can be obtained by the WENO reconstruction with the grid values of Fij, Gij. A addtive Runge-Kutta scheme method [45] is used to couple the Euler equations to the chemical reactions. According to the work of Araújo, Kennedy C.A., Inmaculada Higueras [46][47][48], additional Runge-Kutta explicit implicit algorithm can be used to handle the following ordinary differential equations: In this paper, the Euler equation is divided into the convection term and the reaction term. N in the equation is set to 2.
The reaction Euler equation is written as: It refers to both the nonrigid convection term and the rigid reaction term. The IMEX ARK2 method is adopted in which, ERK (display Runge-Kutta) method are used for the nonrigid term, and the rigid term is treated by ESDIRK (rigid diagonal explicit/implicit Runge-Kutta) method.
The above equation is written as and U i is iteratively solved.I is the unit matrix.

Detailed Chemical Kinetics Model
A detailed chemical reaction model comprised of 9 species and 48 elementary reactions was employed for describing hydrogen-oxygen-argon detonation, as shown in Table 1 where v ik and v ik are the reactant and product chemical stoichiometric coefficients, x i denotes the specie, K f ,k and K b,k are the forward and backward reaction rate constants, respectively, they can be work according to the Arrhenius law and the chemical equilibrium conditions. K f ,k and K b,k follow the Arrhenius law and satisfy the chemical equilibrium conditions: where R i = R µ /w i , R µ is the universal gas constant, P atm is the atmospheric pressure, s 0 i k b,k is the entropy at the standard state, A k is the pre-exponential constant, nk is the temperature power and E ak is the activation energy. The finite production rate ω i for the i-th specie is calculated under the assumption that all involved reactions are elementary reactions, where w i and C xi are the molecular weight and the molar concentration of i-th species (C xi = ρ i /w i ), respectively, and a ik represents the third body coefficients for the i-th species and k-th equations.

Validation of Numerical Scheme
Diffusion reaction equation is shown in the following: ε is the stiffness factor, which shows the stiff intensity of the equation. The initial condition is set as u(0, x) = sin(2πx) and the periodic boundary condition is adopted. The analytical solution of the equation is u(t, x) = e εt sin(2π(x − t)). The fifth order WENO difference scheme is used to discretize the convection term and the explicit processing is carried out while the reaction term is treated implicitly. The results are shown in Figure 2.
analytical solution of the equation is ( , ) sin(2 ( )) t u t x e x t ε π = − . The fifth ord difference scheme is used to discretize the convection term and the explicit pro carried out while the reaction term is treated implicitly. The results are shown in  Time accuracy is a measurement of time. To test the time accuracy of the use the equation as follows: where q is the time accuracy and L is the numerical error presented by the value of the difference between the exact solution and the numerical solution. accuracy test results are shown in the following Table 2: 3.29D-5 7.9 1.63D-6 15.9 7.43D-8 where q is the time accuracy and L is the numerical error presented by the absolute value of the difference between the exact solution and the numerical solution. The time accuracy test results are shown in the following Table 2: It can be seen from the table that the fluctuation range of ARK method in different time steps is not very large, which shows that it has high time precision.

Validation of the Physical Model
In order to prove the validation of the physical model, the numerical simulation results are compared with the experimental results.
A comparison of the detonation reflection on a 30 • wedge between experiment and simulation is conducted. The numerical cellular detonation cell pattern and the experimental smoked foil record are presented together in the figure as follows. It is obvious that the wave front near the wedge is overdriven due to the shock reflection and the Mach stem In order to prove the validation of the physical model, the numerical simulation results are compared with the experimental results.
A comparison of the detonation reflection on a 30° wedge between experiment and simulation is conducted. The numerical cellular detonation cell pattern and the experimental smoked foil record are presented together in the figure as follows. It is obvious that the wave front near the wedge is overdriven due to the shock reflection and the Mach stem detonation causes smaller detonation cells behind it. The very similar characteristics exist both in the smoked foil image in Figure 3a and numerical results in Figure  3b.
(a) Experiment smoke foil photograph (b) Numerical results A comparison of the detonation reflection on a 90° bend between experiment and simulation is conducted. The numerical cellular detonation cell pattern and the experimental smoked foil record are presented together in the figure as follows. The cellular structure disappears in the detonation failure region and the cellular structure is regular in the stable detonation regular. The very similar characteristics exist both in the smoked foil image in Figure 4a and numerical results in Figure 4b.

The Effect of Resolution on the Results
We verify the influence of grid resolution on the computation. With the improvement of the grid resolution in the induction region, the triple structure can be resolved more clearly. There is no difference between 20 and 30, so we choose 20, which can be seen in the Figure 5.

The Effect of Resolution on the Results
We verify the influence of grid resolution on the computation. With the improvement of the grid resolution in the induction region, the triple structure can be resolved more clearly. There is no difference between 20 and 30, so we choose 20, which can be seen in the

The Effect of Resolution on the Results
We verify the influence of grid resolution on the computation. With the improvement of the grid resolution in the induction region, the triple structure can be resolved more clearly. There is no difference between 20 and 30, so we choose 20, which can be seen in the Figure 5.

Cell Structure
Comparing Figure 6 θ = 15°, θ = 20°, θ = 30°withθ = 40°, θ = 45°, θ = 60°, a great difference is found between the results of these two groups. In the same pipe distance, the stable cell structure of the first group is destroyed and the size of cell becomes inhomogeneous after being affected by the curvature of the pipe. However, the cellular detonation wave can spread continuously. The second group is so heavily affected by the curvature of the pipe that a wide range of cell seven disappear completely on the surface of the pipe wall.

Cell Structure
Comparing Figure 6 θ = 15 • , θ = 20 • , θ = 30 • with θ = 40 • , θ = 45 • , θ = 60 • , a great difference is found between the results of these two groups. In the same pipe distance, the stable cell structure of the first group is destroyed and the size of cell becomes inhomogeneous after being affected by the curvature of the pipe. However, the cellular detonation wave can spread continuously. The second group is so heavily affected by the curvature of the pipe that a wide range of cell seven disappear completely on the surface of the pipe wall. great difference is found between the results of these two groups. In the same pipe distance, the stable cell structure of the first group is destroyed and the size of cell becomes inhomogeneous after being affected by the curvature of the pipe. However, the cellular detonation wave can spread continuously. The second group is so heavily affected by the curvature of the pipe that a wide range of cell seven disappear completely on the surface of the pipe wall.  Figure 6. Cellular structure of different bend angles. Figure 6 shows different propagation characteristics of detonation wave in the bend pipe, which can be explained by the change of cellular structure. The cellular structure plays an important role in the investigation of detonation propagation mechanism. The  Figure 6 shows different propagation characteristics of detonation wave in the bend pipe, which can be explained by the change of cellular structure. The cellular structure plays an important role in the investigation of detonation propagation mechanism. The detonation wave propagation is able to form a regular cellular structure and this is mainly driven by the chemical reaction. The disappearance of cellular structure indicates that the energy supplied by the chemical reaction is insufficient to maintain the structure and the chemical reaction surface is gradually decoupled from the detonation wave front. The increase of cellular structure shows that the driving energy supplied by chemical reaction is less than the requirement of the stable cell and the detonation wave is under driven. The decrease of cellular structure implies that the energy supplied by the chemical reaction is greater than the needs of stable cell and the detonation wave is over driven. The decrease of the chemical energy is generated by the rarefaction wave on the up wall, which causes the pressure and temperature drop behind the detonation, while the increase is due to the compression wave on the lower wall, which causes the pressure and temperature increases behind the detonation. The disappearance of the cellular structure also indicates the detonation is extinguished. The stable structure must have periodic transverse wave structure. However, the rarefaction wave generated in the upper wall makes the transverse wave decay rapidly near the edge of the wall. Thus, the quenching area is generated.
A comparison between Figure 6 θ = 15 • , θ = 20 • and θ = 30 • shows that the whole pipe is divided into A, B and C regions. The structure in A region is regular and the propagation of detonation is stable. Cellular structure in B region becomes irregular, where the cell size increases in the upper half while it decreases in the lower half. The reason is that the detonation wave is diffracted on the upper wall, and the pressure of the detonation wave reduced by the rarefaction makes it instable. This leads to the increasing of cellular structure. However, the Maher reflection happens between the pipe and detonation wave in the lower wall and induces increase of the pressure and temperature correspondingly after the Maher. With the increase of the detonation wave intensity, the cell structure is reduced. In contrast to B region, the cell size in C region decreases in the upper part while it increases in the lower part. The result is that along with the propagation of detonation waves spreading to downstream of the bend, the three wave structure interacts among the others at the same time, the Maher reflection between the detonation and the pipe is stopped. Unlike θ = 15 • and θ = 20 • , the results of θ = 30 • show that there is a critical detonation in the B region (the transition curvature between the presence of quenching and absence of quenching).
Comparison between Figure 6 θ = 40 • , θ = 45 • and θ = 60 • shows that a quenching area exists in B region. The cellular structure disappears in this area. With the increase of the curvature of the bend in θ = 40 • , θ = 45 • and θ = 60 • , the no cell area moves to the upstream bend. θ = 40 • has the smallest area in the farthest distance, while θ = 60 • has the largest area in the nearest distance. Figure 7 reveals the effect of the curvature of the bend on the quenching position. Then, as the detonation wave leaves the elbow section, θ = 40 • , θ = 45 • and θ = 60 • in the downstream position close to the upper wall has a range of pressure concentration, in which the detonation is reinitiated. The reason is when the detonation wave comes out gradually from the elbow section, the pressure experiences significant reinitiation because of the collision of three wave structure. Similar to DDT, a strong detonation wave is formed and then the steady detonation wave is restored slowly. pressure experiences significant reinitiation because of the collision ture. Similar to DDT, a strong detonation wave is formed and then t wave is restored slowly.

Internal and External Boundary Pressure Distribution
In order to own the evolution process of the detonation wave distribution in different curvature pipes is shown in Figure 8 from θ the 6 pictures, θ = 15° and θ = 20° is in a similar situation, the sa

Internal and External Boundary Pressure Distribution
In order to own the evolution process of the detonation wave in tube, the pressure distribution in different curvature pipes is shown in Figure 8 from θ = 15 • to θ = 60 • . In the 6 pictures, θ = 15 • and θ = 20 • is in a similar situation, the same trend of the pressure is shown on the upper wall and lower wall. The pressure of the outer wall increases and then decreases, while the inner wall increases slowly, which explains why the small curvature bending has hardly any effect on the propagation of detonation. At the same time, the pressure of the inner and outer walls is not changed greatly. The pressure on the inner wall of diagram θ = 30 • decreases gradually, and with the increase of curvature, the rarefaction effect on the inner wall is strengthened gradually. However, there is a critical detonation point different to θ = 15 • and θ = 20 • where the pressure of inner wall suddenly increased in the late stage of detonation propagation and then decreased. All θ = 40 • , θ = 45 • and θ = 60 • are similar and have a big jump in the downstream of the inner wall where the pressure change increased more than 7.3 atm. Owing to the rarefaction of the inner wall, the detonation wave is extinguished, while the detonations are overdriven on the outer wall. A larger pressure gradient is generated between inner and outer walls and moves to the upper wall along with the propagation of detonation wave. Then a strong transverse wave is produced and compacts with the inner wall leading to a reflection which reinitiates the unburned gas rapidly. the inner wall where the pressure change increased more than 7.3 atm. Owing to the rarefaction of the inner wall, the detonation wave is extinguished, while the detonations are overdriven on the outer wall. A larger pressure gradient is generated between inner and outer walls and moves to the upper wall along with the propagation of detonation wave. Then a strong transverse wave is produced and compacts with the inner wall leading to a reflection which reinitiates the unburned gas rapidly.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 14 of 23 From the pressure in Figure 8 θ = 40°, θ = 45° and θ = 60°, it can be found that a peak pressure on outer wall exists when the detonation wave begins entering the pipe, which indicates that the detonation has a Maher reflection increasing the pressure. As the detonation wave propagates through the downstream of the pipe, the pressure on the outer wall decreases gradually, while a peak pressure appears on the inner wall of the From the pressure in Figure 8 θ = 40 • , θ = 45 • and θ = 60 • , it can be found that a peak pressure on outer wall exists when the detonation wave begins entering the pipe, which indicates that the detonation has a Maher reflection increasing the pressure. As the detonation wave propagates through the downstream of the pipe, the pressure on the outer wall decreases gradually, while a peak pressure appears on the inner wall of the downstream pipe. This peak pressure is much higher than that in the upstream pipe. It shows the beginning of the detonation. At the same time, it is shown that when the detonation wave propagates in the bend with a certain curvature, the range of pressure increase appears both in the upstream and downstream of the bend, but the pressure on the inner wall of downstream is much larger than the outer wall. Therefore, considering the pressure of the material surface in the industrial design, these two positions of the bend pipe should be taken into account and the inner wall of downstream especially needs an extra reinforcement.
From the maximum pressure curve in Figure 9, it is obvious that the 30 • curvature is critical value. This also verifies that reinitiation will appear when the degree of pipe is more than 30 • , which is also shown in Figure 6. The peak of the pressure in the 45 • bend pipe is the largest one and this fits to the result in Figure 8. The results show that the pressure reinitiation at this ratio is larger than that other curvature on the inner wall. Based on this conclusion, the curvature of bend surface considered in actual engineering damage should try to avoid 45 • .  Figure 10 is the detonation wave propagation process in differen ing, in which high brightness white represents high pressure region tween θ = 15°, θ = 30° and θ = 60° with three kinds of curvature is p found the propagation processes are different. The small curvature be effect on the forward propagation of detonation wave and the detona  Figure 10 is the detonation wave propagation process in different curvature bending, in which high brightness white represents high pressure region. Comparison between θ = 15 • , θ = 30 • and θ = 60 • with three kinds of curvature is performed and it is found the propagation processes are different. The small curvature bend has hardly any effect on the forward propagation of detonation wave and the detonation wave front is tidy. The position of the detonation wave front near the upper wall in θ = 30 • is relatively large, and the wave front becomes a bit flat, while the influence of the upper and lower walls of the detonation wave in θ = 60 • are relatively large, making the whole wave front distorted and generating a high pressure transverse wave near the upper wall of the bend downstream.

Evolution Process of Detonation Wave in Bend
found the propagation processes are different. The small curvature bend has hardly any effect on the forward propagation of detonation wave and the detonation wave front is tidy. The position of the detonation wave front near the upper wall in θ = 30° is relatively large, and the wave front becomes a bit flat, while the influence of the upper and lower walls of the detonation wave in θ = 60° are relatively large, making the whole wave front distorted and generating a high pressure transverse wave near the upper wall of the bend downstream.

Details of Chemical Reactions
The variation of the curvature of the bend has a huge influence on the evolution of detonation wave. The propagation mechanism is shown in different forms when propagating indifferent curvature bends. The performance of the chemical reaction in the process of detonation wave propagation in different time and different propagation mechanism is shown in Figure 11.

Details of Chemical Reactions
The variation of the curvature of the bend has a huge influence on the evolution of detonation wave. The propagation mechanism is shown in different forms when propagating indifferent curvature bends. The performance of the chemical reaction in the process of detonation wave propagation in different time and different propagation mechanism is shown in Figure 11.

Details of Chemical Reactions
The variation of the curvature of the bend has a huge influence on the evolution of detonation wave. The propagation mechanism is shown in different forms when propagating indifferent curvature bends. The performance of the chemical reaction in the process of detonation wave propagation in different time and different propagation mechanism is shown in Figure 11 In Figure 11 θ = 15°, we can find that when detonation wave propagates in a small curvature pipe, the OH composition is mainly distributed in the wave front and the lower wall, indicating that this local chemical reaction is much more intense and the detonation wave continues to propagate forward without any distortion in the whole propagation process. This is in good agreement with the analyzed results of pressure profile. Figure θ = 15° is similar to Figure θ = 30°, except for some disturbances appearing to the surface wave front on the inner wall in the process of detonation wave propagation and making the wave front unsmooth. The detonation wave still propagates forward. Figure 11 θ = 60° is very different from θ = 15° and θ = 30°; it can be found the detonation is greatly affected when it is entering the bend. An air mass surrounded by detonation wave appears near the upper wall. As the detonation wave continues to move forward, the wave front becomes very in homogeneous after a large disturbance. Moreover, the air mass is increased and the wave front becomes more inhomogeneous. At the same time, the surface is gradually depressed, indicating that there is no chemical reaction and the detonation wave and chemical reaction is decoupled, which prompts the detonation quenching near the inner wall. By analyzing the OH content map, we find that the detonation wave moving toward the inner wall is gradually formed near the lower wall and ignites the unreacted gas on the inner surface. This drives the whole detonation wave front to move forward continuously. Some local high concentration points of OH group is found near the inner wall, indicating the formation of hot spots which evidently shows that the chemical reaction is more intense and the detonation wave is formed. The content of OH in the inner wall still stays at a high level in the later stage of the detonation, which shows that the chemical reaction on the surface of the in- Figure 11. Profiles of OH components with different propagation characteristics.
In Figure 11 θ = 15 • , we can find that when detonation wave propagates in a small curvature pipe, the OH composition is mainly distributed in the wave front and the lower wall, indicating that this local chemical reaction is much more intense and the detonation wave continues to propagate forward without any distortion in the whole propagation process. This is in good agreement with the analyzed results of pressure profile. Figure θ = 15 • is similar to Figure θ = 30 • , except for some disturbances appearing to the surface wave front on the inner wall in the process of detonation wave propagation and making the wave front unsmooth. The detonation wave still propagates forward. Figure 11 θ = 60 • is very different from θ = 15 • and θ = 30 • ; it can be found the detonation is greatly affected when it is entering the bend. An air mass surrounded by detonation wave appears near the upper wall. As the detonation wave continues to move forward, the wave front becomes very in homogeneous after a large disturbance. Moreover, the air mass is increased and the wave front becomes more inhomogeneous. At the same time, the surface is gradually depressed, indicating that there is no chemical reaction and the detonation wave and chemical reaction is decoupled, which prompts the detonation quenching near the inner wall. By analyzing the OH content map, we find that the detonation wave moving toward the inner wall is gradually formed near the lower wall and ignites the unreacted gas on the inner surface. This drives the whole detonation wave front to move forward continuously. Some local high concentration points of OH group is found near the inner wall, indicating the formation of hot spots which evidently shows that the chemical reaction is more intense and the detonation wave is formed. The content of OH in the inner wall still stays at a high level in the later stage of the detonation, which shows that the chemical reaction on the surface of the inner wall is strongly affected by the shock and is consistent with the result of the pressure curve. Figure 12 is the concentration of chemical reactions in the bend pipe. Hydrogen and oxygen in reactants (a), (c) and products (H) change fastest. (b), (d) and (E) represent the intermediate products of O, H and OH and indicate the chemical reaction. The concentration of these intermediate products has the maximum value at the reaction front and then gradually decreased. The pressure profile of (E) shows that a high pressure area is formed near the lower wall affected by the compression wave, while the unreacted gas near the upper wall is formed by rarefaction wave. The contents of (f) and (g) reach the maximum value at the reaction front and then decay very fast. Comparing (f), (g), (b), (d) and (E), it is found that H 2 O 2 and HO 2 in the chemical reaction flow field are compositions with the least content.

Profiles of Different Chemical Components
ner wall is strongly affected by the shock and is consistent with the result of the pressure curve. Figure 12 is the concentration of chemical reactions in the bend pipe. Hydrogen and oxygen in reactants (a), (c) and products (H) change fastest. (b), (d) and (E) represent the intermediate products of O, H and OH and indicate the chemical reaction. The concentration of these intermediate products has the maximum value at the reaction front and then gradually decreased. The pressure profile of (E) shows that a high pressure area is formed near the lower wall affected by the compression wave, while the unreacted gas near the upper wall is formed by rarefaction wave. The contents of (f) and (g) reach the maximum value at the reaction front and then decay very fast. Comparing (f), (g), (b), (d) and (E), it is found that H2O2 and HO2 in the chemical reaction flow field are compositions with the least content.

Conclusions
In this paper, the influence of different curvature bends on the propag teristics of the cellular detonation wave is systematically studied. The cellula wave propagated in the bends is interacted by both compression wave and tion wave in the inner and outer wall. The propagation mechanism is conclu parts: bend curvature with angle less than 30° has rarely effect on propagatio detonation wave, which can spread themselves stably and has relatively sm onation wave front and almost unchanged cell size. Curvature bend with than 30° has a great influence on the propagation of the detonation wave w decoupling, quenching and reinitiation. The whole propagation process reac value after reinitiation. The 45° bending curvature has the largest pressure in bends. The detonation wave cellular disappeared after quenching and beca ter the reinitiation. In the details of the chemical reaction, the difference o wave propagation mechanism in different bending curvature was shown in figure of OH. Meanwhile, the way the detonation wave moved forward was is obvious that the wave front in the quenching bends is largely different fro nonquenching bends.

Conclusions
In this paper, the influence of different curvature bends on the propagation characteristics of the cellular detonation wave is systematically studied. The cellular detonation wave propagated in the bends is interacted by both compression wave and the rarefaction wave in the inner and outer wall. The propagation mechanism is concluded in three parts: bend curvature with angle less than 30 • has rarely effect on propagation of cellular detonation wave, which can spread themselves stably and has relatively smoothed detonation wave front and almost unchanged cell size. Curvature bend with angle more than 30 • has a great influence on the propagation of the detonation wave which leads to decoupling, quenching and reinitiation. The whole propagation process reached the peak value after reinitiation. The 45 • bending curvature has the largest pressure in all kinds of bends. The detonation wave cellular disappeared after quenching and became larger after the reinitiation. In the details of the chemical reaction, the difference of detonation wave propagation mechanism in different bending curvature was shown in composition figure of OH. Meanwhile, the way the detonation wave moved forward was indicated. It is obvious that the wave front in the quenching bends is largely different from that in the nonquenching bends.