Energy Loss in Pulse Detonation Engine due to Fuel Viscosity

1 School of Mechanics, Civil Engineering and Architecture, Northwestern Polytechnical University, Xi’an, Shaanxi 710072, China 2 State Key Laboratory of Structural Analysis of Industrial Equipment, Dalian University of Technology, Dalian, Liaoning 116023, China 3 Guangxi Key Laboratory of Disaster Prevention and Structural Safety, Guangxi University, Nanning, Guangxi 530004, China 4 School of Mechanical Engineering, Northwestern Polytechnical University, Xi’an, Shaanxi 710072, China


Introduction
Improving the efficiency of the pulse detonation engine (PDE) is an everlasting objective for the researchers in the associated fields.Generally, the efficiency of the PDE is determined by the energy loss due to the viscosity of the fuel, the heat dissipation of the structure, and other factors, in which the energy loss due to the viscosity of the fuel is a major factor.
In the current research findings, the effects on the efficiency of the PDE from the structure parameters of PDE, the mass fraction of the fuels, and the frequency of the detonation have been reported.For example, Bratkovich and Bussing [1] presented a PDE performance model, which considered inlet, mixer, combustor, and nozzle contributions to PDE performance attributes, referring to the ideal PDE model presented by Endo and Fujiwara [2].Radulescu and Hanson [3] presented a study to address the influence of convective heat losses on the flow field and performance of PDE, considering various working conditions; Tangirala and Dean [4] proposed a performance estimation model for the PDE with exit nozzle; Glaser [5] studied the effect on the efficiency of the PDE from the structure parameters of PDE in laboratory; and Brophy and Hanson [6] investigated the fuel distribution effects on pulse detonation engine operation and performance.Recently, Yan and colleagues [7][8][9][10] studied the fuel filled length, the frequency of the detonation, and the structure parameters of PDE effects on the performance of the PDE.
For the limitation of the measuring means, the energy loss due to the viscosity of the fuel in the PDE has not been reported in current literatures, while the energy losses due to viscosity of the fluid in some fluid dynamic models have been investigated in many literatures [11][12][13][14][15].The research results implied that the energy losses due to viscosity could not be neglected in the majority of fluid models.Thus, the energy loss due to the viscosity of the fuel in the PDE is presented by studying the energy loss due to the viscosity of the fluid in Burgers model excited by periodic impulses.Focusing on the energy loss due to the viscosity of the fluid in the PDE, a Burgers model excited by periodic impulses is introduced to describe the flow field in the PDE firstly; and then based on the generalized multisymplectic theory [16][17][18], the expression of the energy loss due to the viscosity of 2 Mathematical Problems in Engineering the fluid is derived from the first-order symmetric form of the Burgers model; finally, the numerical results of the energy loss are reported.

Energy Loss in PDE Model due to Fuel Viscosity
The practical model of the PDE in laboratory is a slender tube; thus the practical flow field in the detonation tube is threedimensional.Ignoring the effect of the boundary layer, the flow field on arbitrary cross section in the detonation tube is homogeneous approximately.With this assumption, the three-dimensional flow field in the detonation tube can be simplified to one-dimensional flow field.
In addition, the purpose of this paper is to investigate the energy loss due to the viscosity of the fuel; thus the effect of the chemical reactions among the components of the fuel can be represented by a series of periodic impulses.Then the flow fields in the detonation tube can be expressed by the following Burgers model: where ] is the kinematics viscosity of the fuel.(, ) =  ∑  =0 Ψ( + ,  −  0 −  +  0 ) ( = 0, 1, . . ., ,  = [( −  0 )/], [⋅] is the symbol of the rounding operation) are sequence of periodic impulses that describes the effect of the periodic detonations resulting from the chemical reactions among the components of the fuel. +  is the moment of the impulse and  −  0 −  +  0 is the position of the impulse that satisfies 0 ≤  −  0 −  +  0 ≤ ;  is a constant that denotes the amplitude of periodic detonation wave.Ψ is a switch function: is the length of the detonation tube. is the period of the detonation. is the propagation speed of the detonation wave. 0 is the position where the first detonation wave appears. 0 is the time lag between burning the fuel and appearing of the first detonation wave.
Equation ( 1) can be used to describe the flow filed of the transient state as well as the steady state in the detonation tube.The symmetric form of (1) can be obtained referring to the results of [17].
For simplicity, skew-symmetric matrices M, K, and symmetric matrix K are defined as Obviously, if ] = 0 and F(, ) = 0, the symmetric form (2) is a multisymplectic form [19,20], which satisfies the multisymplectic conservation law and the local energy conservation law as well as the local momentum conservation law.
Actually, the kinematic viscosity coefficient of the fuel ] in the PDE is a positive real number, which makes the symmetric form (2) generalized multisymplectic [17] with external excitation.Following the outline of [17], the local energy error of the symmetric form (2) can be obtained [17]: where the local energy error contains two parts: one is the local energy increment Δ le = (1 + ) ∑  =0 Ψ( + ,  −  0 −+ 0 ) resulting from the impulse excitation; another part is the local energy loss Δ ] = (]/2)   due to the fluid viscosity.
Superficially, the local energy loss Δ ] has nothing to do with all the detonation parameters (, , ,  0 ,  0 , and ), but actually the effect of the detonation is included in the distribution of the flow field .
From the local energy loss (3), it is easy to get the global energy error; In this paper, we consider that the components of the fuel are mixed uniformly, which means that the kinematic viscosity coefficient of the fuel ] in the PDE is independent of time and position; thus, the global energy error is similarly, the global energy error contains two parts: one is the global energy increment Δ ge =  2 (1 + ) ∑  =0 Ψ( + , − 0 −+ 0 ) resulting from the impulse excitation; another part is the global energy loss Δ ] = (]/2) ∫  0   d due to the fluid viscosity.

Generalized Multisymplectic Scheme
To investigate the energy loss in PDE by generalized multisymplectic method [16,17], the generalized multisymplectic scheme for the symmetric form (2) must be constructed firstly.
Referring to [17], introduce a uniform grid {(  ,   )} in  2 with time step size Δ in the -direction and spatial step size Δ in the -direction.The approximate value of z(, ) at the point (  ,   ) is denoted by z   .
Using the Preissman box discrete method, the discrete form of the symmetric form (2) can be obtained: where  +  is the forward differences to   z and  +  is the forward differences to It has been proven that the discrete form ( 8) is generalized multisymplectic with structure-preserving property if the kinematic viscosity coefficient of the fuel ] is small enough [17].
The discrete form of the local energy error of the symmetric form (2) at the point (  ,   ) is where the local energy increment resulting from the impulse excitation is (Δ le )   = (1 + ) ∑  =0 Ψ( + ,  −  0 −  +  0 ) and the local energy loss due to the fluid viscosity is (Δ ] )   = (]/2) +1/2 +1/2  +   +   +1/2 +1/2 .Similarly, the discrete form of the global energy error of the symmetric form (2) at the   is where the global energy loss due to the fluid viscosity is (/Δ)+1/2 ] (we select the trapezoid quadrature method to calculate the numerical integration ∫  0   d) and the global energy increment resulting from the impulse excitation is

Numerical Results
In this section, some numerical results about the energy loss in PDE will be presented.In the numerical experiments, we let the step lengths be Δ = 5 × 10 We simulate the flow field in the PDE by the structurepreserving scheme (8) and calculate the energy error on each grid.The numerical results of the local energy error and the local energy loss due to the viscosity of the fuel at different time are shown in Figures 1 and 2.
The local energy error is equal to the negative local energy loss due to the viscosity of the fuel on each grid at  = 5×10 −8 s ((Δ  )   = −(Δ ] )   ), and the local energy error on each grid is negative, see Figure 1, which implies that the local energy increment resulting from the impulse excitation on each grid ((Δ le )   ) is zero (Ψ = 0) at  = 5 × 10 −8 s.The local energy error, the local energy loss due to the viscosity of the fuel, and the local energy increment resulting from the impulse excitation increase quickly on the grids near the detonation wave at  = 4 × 10 −4 s, see Figure 2, which implies that the local energy loss due to the viscosity of the fuel increases when the speed of the fuel increases.
Based on the numerical results of the local energy error, the global energy error, the global energy increment resulting from the impulse excitation, and the global energy loss due to the fluid viscosity at each time step can be easily obtained.Figure 3 implies that the global energy loss due to the fluid viscosity can not be neglected in PDE.Neglecting the global energy loss due to the fluid viscosity will exaggerate the thrust of the PDE (here, the thrust is represented by the global energy increment) by ((5.436924761525581×10 4 )/(8×10 5 ))× 100% ≈ 6.7962% at  = 5 × 10 −4 s.

Conclusion
Energy loss in PDE due to the viscosity of the fuel is an important factor that affects the efficiency of PDE.In this paper, the energy loss due to the fuel viscosity is investigated numerically.A Burgers model excited by periodic impulses is established to describe the detonation process of PDE firstly.Then, the energy loss expressions are derived from the generalized multisymplectic form of the Burgers model.Finally, the numerical results of the energy loss in PDE are reported.
From the results, it can be concluded that the energy loss

Figure 1 :
Figure 1: Local energy error and local energy loss due to the viscosity of the fuel on each grid at  = 5 × 10 −8 s.

5 Figure 2 : 2 Figure 3 :
Figure 2: Local energy error, local energy increment resulting from the impulse excitation, and local energy loss due to the viscosity of the fuel on each grid at  = 4 × 10 −4 s.