Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2014, Article ID 735926, 5 pages http://dx.doi.org/10.1155/2014/735926

Research Article

Energy Loss in Pulse Detonation Engine due to Fuel Viscosity

Weipeng Hu,1,2,3 Ziehen Deng,1,2 and Gongnan Xie4

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

Correspondence should be addressed to Weipeng Hu; wphu@nwpu.edu.cn Received 7 April 2014; Accepted 5 May 2014; Published 29 May 2014 Academic Editor: Zhijun Zhang

Copyright © 2014 Weipeng Hu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Fluid viscosity is a significant factor resulting in the energy loss in most fluid dynamical systems. To analyze the energy loss in the pulse detonation engine (PDE) due to the viscosity of the fuel, the energy loss in the Burgers model excited by periodic impulses is investigated based on the generalized multisymplectic method in this paper. Firstly, the single detonation energy is simplified as an impulse; thus the complex detonation process is simplified. And then, the symmetry of the Burgers model excited by periodic impulses is studied in the generalized multisymplectic framework and the energy loss expression is obtained. Finally, the energy loss in the Burgers model is investigated numerically. The results in this paper can be used to explain the difference between the theoretical performance and the experimental performance of the PDE partly. In addition, the analytical approach of this paper can be extended to the analysis of the energy loss in other fluid dynamic systems due to the fluid viscosity.

1. 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-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-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 [1618], the expression of the energy loss due to the viscosity of

the fluid is derived from the first-order symmetric form of the Burgers model; finally, the numerical results of the energy loss are reported.

2. 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 three-dimensional. 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:

dtu + udxu- vdxxu = f(t,x) xe[0,L], (1)

where v is the kinematics viscosity of the fuel.

f(t,x) = A^=0 ¥S(t + nT,Vt - Vcp0 - nVT + x0) (n = 0,1,...,N,N = [(t - (p0)/T], [•] 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. t + nT is the moment of the impulse and Vt - V(p0 -nVT + x0 is the position of the impulse that satisfies 0 < Vt - Vcp0 - nVT + X0 < L; A is a constant that denotes the amplitude of periodic detonation wave. ¥ is a switch function: ¥ = {0 . L is the length of the detonation tube. T is the period of the detonation. V is the propagation speed of the detonation wave. x0 is the position where the first detonation wave appears. <p0 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].

Introducing canonical momenta dxp = u, dxq = -(1/2) dtu, and defining state variable z = (p, u, q)T,(1) becomes the first-order partial differential equations

1 2 — 0 0 2

dt z +

0 0 1 0 0 0

0 0 0 + [0 v 0

[-1 0 0 [0 0 0

where S(z) = -uq + (1/6)u3 is the Hamiltonian function and the external excitation is

F (t,x)

= (0,^f(t,x) d%,0)

0, I A^WS (t + nT, Vt - Vcp0 - nVT + x0) dx, 0

= ( 0, AL^We (t + nT, Vt - Vcp0 - nVT + x0 ),0 ) ;

here e(t + nT, Vt - V(p0 - nVT + x0) is the step function.

For simplicity, skew-symmetric matrices M, K, and symmetric matrix K are defined as

0 1 ~2

0 0 1 0 0 0 -10 0

0 0 0 0 v 0 0 0 0

Obviously, if v = 0 and F(i, x) = 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 v 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]:

A, = df

AL^We (t + nT, Vt - Vcp0 - nVT + x0)

+S(z)--zT (K + k^

+ 2-dxzr (K + K)dtz

= AL(1 + V) ^WS (t + nT, Vt - Vcp0 - nVT + x0)

+ 2dxudtu - 2dt (udxu)

= AL(1 + V) (t + nT, Vt - Vcp0 - nVT + x0)

= VzS(z) + F (t,x),

-2Ud*x U=Ale -A V

where the local energy error contains two parts: one is the local energy increment A ¡e = AL( 1 + V) WS(t + nT, Vt -V(p0-nVT+x0) resulting from the impulse excitation; another part is the local energy loss A iv = (v/2)udtxu due to the fluid viscosity.

Superficially, the local energy loss A h has nothing to do with all the detonation parameters (L, V, T, x0, %, and A), but actually the effect of the detonation is included in the distribution of the flow field u.

From the local energy loss (3), it is easy to get the global energy error;

AL(1+V) (t + nT, Vt - Vcp0 - nVT + xQ )

- 2udt* u

In this paper, we consider that the components of the fuel are mixed uniformly, which means that the kinematic viscosity coefficient of the fuel v in the PDE is independent of time and position; thus, the global energy error is

A g =1 ALYW (l + V)S(t + nT,Vt- V<p0 - nVT + x0) dx

udtx udx

= AL2 (1 + V) Y^z (t + nT, Vt - Vcp0 - nVT + x0) - - f

udtx udx

Using the Preissman box discrete method, the discrete form of the symmetric form (2) can be obtained:

1 2 — 0 0 2

St Zj+1/2 +

0 1 -0

X-+ t+l/2 s X zj

= + F (и),

where S+ is the forward differences to dtz and is the forward differences to dxz, z1-+l/2 = (z'j+1 + z'^/2, z'+1/2 =

(z+ + z!l )/2, and z^ = (Z+1 + z)+i + Z+1 + Zj)/4.

It has been proven that the discrete form (8) is generalized multisymplectic with structure-preserving property if the kinematic viscosity coefficient of the fuel v is small enough [17].

The discrete form of the local energy error of the symmetric form (2) at the point (tt, Xj) is

(Arfj = AL(l+V)Y^S(i + nT, Vi - V<p0 - nVT + x0)

-lu^S+XuM = (A k)\ -(A h)\,

where the local energy increment resulting from the impulse excitation is (A h)) = AL(1 + V) ZZ0 V8(i + nT, Vi - Vcp0 -nVT + x0) and the local energy loss due to the fluid viscosity is (A 4. = (v/2)u%%8+8lu%%.

Similarly, the discrete form of the global energy error of the symmetric form (2) at the tt is

(A g) = AL2 (1 + V) Y^e (i + nT, Vi - Vcp0 - nVT + x0)

= A ge -A g-''

similarly, the global energy error contains two parts: one is the global energy increment Age = AL2(1 + V) We(t + nT, Vt-VcpQ-nVT+xQ) resulting from the impulse excitation;

another part is the global energy loss A gv = (v/2) jQ udtxudx due to the fluid viscosity.

3. Generalized Multisymplectic Scheme

To investigate the energy loss in PDE by generalized multi-symplectic method [16,17], the generalized multisymplectic scheme for the symmetric form (2) must be constructed firstly.

Referring to [17], introduce a uniform grid {(t{, Xj)] in R2 with time step size A in the -direction and spatial step size Ax in the x-direction. The approximate value of z(t, x) at the point (t¡, xj) is denoted by z'j.

(L/Ax)-1

i+1/2o.+¡¡.+ i+1/2 t V >+1/2 C-+C-+ i+1/2

U1/2 StSXU1/2 +2 L Uj+1/2StSXUj+1/2

¡+1/2

S-+ s>+ i+1/2 ~OtOU,

(L/Ax) + 1/2 t x (L/Ax)+1/2

= (Ael -(A„I

where the global energy loss due to the fluid viscosity is (A gJ = (vAx/4)[u[+/12/28+8+xu[+/12/2 +

2y{L/hx)-l ui+ll2s+(s+ui+H2) + ix1/2 S+S+Uix1/2 1 2Lj=1 ujx1/2°t (0xUjx1/2) + U(L/Ax)x1/20t0xU(L/Ax)x1/21

(we select the trapezoid quadrature method to calculate the numerical integration J0 udtxudx) and the global energy increment resulting from the impulse excitation is (A )' = AL2(1 + V) We(i + nT, Vi - V% - nVT +

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

4. 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 Af = 5 x 10-9 s and Ax = 0.01 m, structure parameter of the PDE be L = 2 m, the detonation parameters be A = 100 m • s-2, V = 2000 m • s-1, x0 = 0.2 m, <p0 = 1 x 10-7 s, and T = 0.02 s (which implies that the detonation frequency is 50 Hz) and the kinematic viscosity coefficient of the fuel v = 0.1. The inset boundary condition considered in this paper is w(i, 0) = 100 m • s-1.

We simulate the flow field in the PDE by the structure-preserving 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 t = 5x 10-8 s ((A= -(A;v)p, 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 ((A fe)p is zero (¥ = 0)at t = 5 x 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 t = 4 x 10- 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.43692476152 5 581x104)/(8x105 ))x 100% = 6.7962% at t = 5 x 10-4 s.

0 0.5 1 1.5 2 x (m) x105

0 0.5 1 1.5 2 x (m) x104

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 t = 4 x 10-4 s.

« 2 0 -2

8 6 4 2 0

0.2 0.3

t (ms)

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

t (ms)

0 0.1 0.2 0.3

t (ms)

Figure 3: Global energy error, global energy increment resulting from the impulse excitation, and global energy loss due to the viscosity of the fuel on each grid at i = 4x 10-4 s.

5. 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

due to the fuel viscosity accounts for 6.7962% of the energy increment resulting from the impulse excitation.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The research is supported by the National Natural Science Foundation of China (11372252 and 11372253), the Science Foundation of Aviation of China (2013ZB53020), the Fundamental Research Funds for the Central Universities (3102014JCQ01035), the Open Project of Guangxi Key Laboratory of Disaster Prevention and Structural Safety (2013ZDK04), and the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment (GZ1312).

References

[1] T. E. Bratkovich and T. R. A. Bussing, "A pulse detonation engine performance model," in Proceedings of the 31st AIAA/ ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, AIAA 95-3155, pp. 1-15, San Diego, Calif, USA.

[2] T. Endo and T. Fujiwara, "Analytical estimation of performance parameters of an ideal pulse detonation engine," Transactions of the Japan Society for Aeronautical and Space Sciences, vol. 45, no. 150, pp. 249-254, 2003.

[3] M. I. Radulescu and R. K. Hanson, "Effect of heat loss on pulse-detonation-engine flow fields and performance," Journal of Propulsion and Power, vol. 21, no. 2, pp. 274-285, 2005.

[4] V. E. Tangirala and A. J. Dean, "Performance estimations of a pulse detonation engine with exit nozzle," in Proceedings of the AIAA/ASME/SAE/ASEE 42nd Joint Propulsion Conference, pp. 4724-4739, Sacramento, Calif, USA, July 2006.

[5] A. J. Glaser, Performance and Environmental Impact Assessment of Pulse Detonation Based Engine Systems, University of Cincinnati, Cincinnati, Ohio, USA, 2007.

[6] C. M. Brophy and R. K. Hanson, "Fuel distribution effects on pulse detonation engine operation and performance," Journal of Propulsion and Power, vol. 22, no. 6, pp. 1155-1161, 2006.

[7] Y. Yan, W. Fan, K. Wang, and Y. Mu, "Experimental investigation of the effect of bell-shaped nozzles on the two-phase pulse detonation rocket engine performance," Combustion, Explosion and Shock Waves, vol. 47, no. 3, pp. 335-342, 2011.

[8] W. Fan, J. Li, Q. Li, and C. Yan, "Numerical investigation on multi-cycle operation of pulse detonation rocket engine," International Journal of Turbo and Jet Engines, vol. 25, no. 3, pp. 189195, 2008.

[9] Q. Zhang, C. Yan, W. Fan, and Q. Li, "Experimental investigation of effect of partial filling on the impulse of pulse detonation engine," Chinese Science Bulletin, vol. 52, no. 20, pp. 2859-2865, 2007.

[10] Q. Li, W. Fan, C.-J. Yan, C.-Q. Hu, and B. Ye, "Experimental investigation on performance of pulse detonation rocket engine model," Chinese Journal of Aeronautics, vol. 20, no. 1, pp. 9-14, 2007.

[11] C.-H. Yap, L. P. Dasi, and A. P. Yoganathan, "Dynamic hemodynamic energy loss in normal and stenosed aortic valves," Journal of Biomechanical Engineering, vol. 132, no. 2, 2010.

[12] K. Dusling, G. D. Moore, and D. Teaney, "Radiative energy loss and v2 spectra for viscous hydrodynamics," Physical Review C, vol. 81, no. 3, Article ID 034907, 2010.

[13] V. Fester, B. Mbiya, and P. Slatter, "Energy losses of non-Newtonian fluids in sudden pipe contractions," Chemical Engineering Journal, vol. 145, no. 1, pp. 57-63, 2008.

[14] T. Mai and T. Liebner, "Loss of energy dissipation capacity from the deadzone in linear and nonlinear viscous damping devices," Earthquake Engineering and Engineering Vibration, vol. 6, no. 1, pp. 11-20, 2007.

[15] S. Norris and G. Mallinson, "Volumetric methods for evaluating energy loss and heat transfer in cavity flows," International Journal for Numerical Methods in Fluids, vol. 54, no. 12, pp. 14071423, 2007.

[16] W. P. Hu, S. M. Han, and Z. C. Deng, "Analyzing dynamic response of non-homogeneous string fixed at both ends," International Journal of Non-Linear Mechanics, vol. 47, no. 10, pp. 1111-1115, 2012.

[17] W. Hu, Z. Deng, and S. Han, "An implicit difference scheme focusing on the local conservation properties for burgers equation," International Journal of Computational Methods, vol. 9, no. 2, Article ID 1240028, 2012.

[18] W. Hu, Z. Deng, S. Han, and W. Zhang, "Generalized multi-symplectic integrators for a class of Hamiltonian nonlinear wave PDEs," Journal of Computational Physics, vol. 235, pp. 394-406, 2013.

[19] T. J. Bridges, "Multi-symplectic structures and wave propagation," Mathematical Proceedings of the Cambridge Philosophical Society, vol. 121, no. 1, pp. 147-190,1997.

[20] T. J. Bridges and S. Reich, "Multi-symplectic integrators: numerical schemes for Hamiltonian PDEs that conserve symplectic-ity," Physics Letters, Section A: General, Atomic and Solid State Physics, vol. 284, no. 4-5, pp. 184-193, 2001.

Copyright of Mathematical Problems in Engineering is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.