Next Article in Journal
Descriptions of Entropy with Fractal Dynamics and Their Applications to the Flow Pressure of Centrifugal Compressor
Next Article in Special Issue
Application of a Fluid–Structure Interaction Model for Analysis of the Thermodynamic Process and Performance of Boil-Off Gas Compressors
Previous Article in Journal
Limiting Uncertainty Relations in Laser-Based Measurements of Position and Velocity Due to Quantum Shot Noise
Previous Article in Special Issue
Performance Analysis of a Proton Exchange Membrane Fuel Cell Based Syngas
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Entropy
Volume 21
Issue 3
10.3390/e21030265
entropy-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Pressure Gain Characteristics and Cycle Performance in Gas Turbines Based on Interstage Bleeding Rotating Detonation Combustion

by
Lei Qi
,
Zhitao Wang
*,
Ningbo Zhao
*,
Yongqiang Dai
,
Hongtao Zheng
and
Qingyang Meng
College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China
*
Authors to whom correspondence should be addressed.
Entropy 2019, 21(3), 265; https://doi.org/10.3390/e21030265
Submission received: 28 January 2019 / Revised: 21 February 2019 / Accepted: 5 March 2019 / Published: 8 March 2019
(This article belongs to the Special Issue Advances in Applied Thermodynamics III)
Browse Figures
Versions Notes

Abstract

:
To further improve the cycle performance of gas turbines, a gas turbine cycle model based on interstage bleeding rotating detonation combustion was established using methane as fuel. Combined with a series of two-dimensional numerical simulations of a rotating detonation combustor (RDC) and calculations of cycle parameters, the pressure gain characteristics and cycle performance were investigated at different compressor pressure ratios in the study. The results showed that pressure gain characteristic of interstage bleeding RDC contributed to an obvious performance improvement in the rotating detonation gas turbine cycle compared with the conventional gas turbine cycle. The decrease of compressor pressure ratio had a positive influence on the performance improvement in the rotating detonation gas turbine cycle. With the decrease of compressor pressure ratio, the pressurization ratio of the RDC increased and finally made the power generation and cycle efficiency enhancement rates display uptrends. Under the calculated conditions, the pressurization ratios of RDC were all higher than 1.77, the decreases of turbine inlet total temperature were all more than 19 K, the power generation enhancements were all beyond 400 kW and the cycle efficiency enhancement rates were all greater than 6.72%.

1. Introduction

As one of the most important types of power plants, gas turbines are widely used in the fields of aviation, shipping, power stations, etc. Faced with a greater demand for high efficiency and low pollution, higher level gas turbines have attracted more and more attention in recent years. However, it is quite difficult to improve the cycle efficiency in conventional gas turbines due to the high entropy change during the combustion process [1,2]. In order to break through the bottleneck of efficiency improvement, various advanced combustion technologies, including detonation [3,4,5], wave rotors [6], and shockless explosion [7,8] were studied in the past few decades.
Theoretical pure detonation combustion which has a lower entropy change and self-pressure gain compared with iso-pressure combustion under similar conditions has already attracted wide attention since the beginning of the 21st century [9,10]. Due to the difficulty of the high Mach number environment in the standing oblique detonation combustion process [11] and the high frequency ignition in the pulse detonation combustion process [12], rotating detonation with its wide working range and self-sustaining propagation has been attracting increasing attention. The evolution characteristics of rotating detonation wave [13,14], the thermodynamic characteristics of rotating detonation flow fields [15,16,17,18,19] and the propulsion performance of the rotating detonator combustor (RDC) [20,21,22] have been investigated with a series of experiments and numerical studies.
In recent years, many studies on the application of rotating detonation combustion to gas turbine cycles have been gradually developed. It has been found that rotating detonation combustion could significantly improve the cycle performance of gas turbines. Wolański et al. applied RDC to a GT350 engine and found that stable continuous detonation was difficult to establish for pure Jet-A fuel. The addition of hydrogen fuel could greatly improve the stability of RDC. With the dual fuel including Jet-A and hydrogen, the thermal efficiency showed an increase of 5–7% compared with that of the base engine [23,24]. After a series of studies on the flow field characteristics of a three-dimensional RDC [25,26,27], Frolov et al. explored the feasibility of the application of RDC to gas turbines in 2016. An isolator was designed to dampen the pressure disturbances towards the upstream compressor and a total pressure gain of 15% was found during the combustion course [28]. Naples et al. improved the export structure of the RDC by adding a trailing annular ejector, which reduced the static pressure pulsation at the RDC outlet by 60–70% [29]. On this basis, a new T63 engine with a RDC using hydrogen as fuel was tested in 2017 and the cycle efficiency was improved under some conditions [30]. Combined with the design of a RDC expansion outlet and a supersonic turbine, Sousa analyzed the cycle performance of rotating detonation gas turbines using hydrogen as fuel and found the cycle efficiency could be increased by about 5% at low pressure ratios, but it showed no increase at high pressure ratios [31]. Qi et al. adopted a simplified model to analyze the cycle performance of the rotating detonation gas turbine under different working conditions. The efficiency could be increased by 5–9% compared with the base gas turbine, and equivalence ratio in RDC played an important role in design of the new cycle [32]. In addition, there were also some studies on the connection between RDCs and turbines from an experimental and numerical view, although with no analysis of the engine cycles [33,34,35].
Up to now, plenty of research on gas turbines based on rotating detonation has used a fuel with a small detonation cell size, such as hydrogen and ethylene, which cannot be applied in marine engines and power stations. The little research using industrial fuel either ignored many of the characteristics of rotating detonation flow fluid or brought in extra structures which were difficult to design. In our previous studies, the ejector was simplified as a mixing model including a series of assumptions which will cause many engineering design difficulties. In this paper, a new thermodynamic cycle model of a gas turbine based on interstage bleeding rotating detonation combustion was established using methane as fuel. The air from interstage bleeding participated in combustion and the higher pressure air at the outlet of compressor sprayed into the outlet of RDC combustor for mixing, which will avoid the need for a complex ejector design and makes better use of the air from interstage bleeding. An investigation of the pressure gain characteristics and cycle performance were carried out with a series of numerical simulations and cycle parameter calculations.

2. Methodology

2.1. Cycle Scheme

Figure 1 shows the cycle scheme of a gas turbine with two combustion modes. Above is the gas turbine cycle based on interstage bleeding rotating detonation combustion, and below is the conventional gas turbine cycle. Compared with the conventional gas turbine cycle, the bleeding air from a 3/4 compressor entered into the RDC and participates in the pressure gain combustion, the higher pressure air from the compressor outlet enters into a mixer and cools the gas from the RDC outlet. The design complexity of the mixer could be effectively reduced as the air from the compressor outlet had a pressure that is higher or close to that of the air from interstage bleeding. In terms of the cycle scheme, the actual T-S diagrams which respectively represent the thermodynamic cycles of the conventional gas turbine and rotating detonation gas turbine are studied in this section, as shown in Figure 2.

2.2. Physical Model and Computing Method

In this section, the model information of gas turbine cycle based on interstage bleeding rotating detonation combustion will be mainly described, including the compressor model, combustion model, mixer model, turbine model and other necessary equations. Meanwhile, to visually investigate the performance of the rotating detonation gas turbine cycle, a conventional gas turbine cycle based on approximate iso-pressure combustion (called reference gas turbine cycle, of which power parameters referred to a certain type of gas turbine in practical application) was brought in for reference. The difference of computing method only exists in the combustion progress to analyze the influence of combustion mode on the cycle performance.
  • Compressor model (1→2)
As the mass flow rate of air g a , compressor inlet temperature T 1 * , compressor inlet pressure P 1 * , compressor efficiency η c , compressor pressure ratio π c , compressor temperature ratio τ were known, compressor outlet temperature T 2 * and compressor outlet pressure P 2 * of compressor could be obtained by:
η c = H 2 s * − H 1 * H 2 * − H 1 *
T 2 * = T 1 * τ
P 2 * = P 1 * π c
where H 1 * and H 2 * are the enthalpies at state points of 1 and 2. H 2 s * is the ideal enthalpy in isentropic condition.
Compressor pressure ratio π c , compressor outlet temperature T 2 * , compressor outlet pressure P 2 * , interstage bleeding temperature T 2 m * , interstage bleeding pressure P 2 m * , mass flow rate of air g a were calculated at the 100% load (design condition) ~50% load of the reference gas turbine, as shown in Table 1. As the power generations of reference gas turbine L E , 0 at different loads were known, the state parameters of compressor outlet in different conditions could be obtained through cycle computations.
  • Combustion model (2m→3c/2→30)
In this paper, RDC(2m→3c) was simplified to a two-dimensional rectangular model, of which the length and width were respectively 800 mm and 200 mm (these sizes were preliminarily determined by referring to [19,36,37,38]). A more detailed simplified method and its reasonability can be found in [39,40], where the effectiveness of this method has been verified. The density-based solver was chose to solve the two-dimensional unstable Euler Equation. The laminar finite rate model was selected due to the fact viscosity is neglected. In the previous studies such as [16,41,42,43,44], the rotating detonation phenomenon could be effectively described by a laminar finite rate model. In addition, to effectively capture the shock waves and detonation waves, the flux term was dispersed by the advection upstream splitting method (AUSM) [45]. The convective term and time term were, respectively, dispersed in the third-order upwind scheme and using the four-step Runge–Kutta method with the second-order accuracy [46]:
∂ Φ ∂ t + ∂ U ∂ x + ∂ V ∂ y = Ω
where the dependent variable vector Φ , the convective flux vectors U and V , and the source vector Ω are defined as:
Φ = ( ρ , ρ u , ρ v , ρ e t , ρ j ) T
U = ( ρ u , ρ u 2 + P , ρ u v , ( ρ e t + P ) u , ρ j u ) T
V = ( ρ v , ρ u v , ρ v 2 + P , ( ρ e t + P ) v , ρ j v ) T
Ω = ( 0 , 0 , 0 , 0 , ω j ) T
where ρ is the density of the premixture, ρ j is the density of species j , u and v are respectively the velocity in X and Y direction, P is the pressure that can be given by the thermal equation of state for a perfect gas. The total internal energy is defined as:
e t = e + 1 2 ( u 2 + v 2 )
where a caloric equation of state for the internal energy of a reacting mixture:
e = e ( ρ j , T )
Figure 3 shows the temperature shadowgraph of the RDC at t = 6800 μs. The temperature characteristics reflected in the temperature contour and extra important pressure characteristic can both be reflected in the temperature shadowgraph.
It could be seen that there were detonation wave A, oblique shock wave B, mixing region C, discontinuity region D and premixed gas region E, which is in good agreement with the published results [47,48,49]. The independence test and model validation have already been carried out in the previous study [19,32], so the grid size and time step were respectively selected as 1 mm and 0.2 μs. In this study, the propagation speed of detonation wave was 1800.4 m/s in the design condition and the theoretical Chapman Jouguet (CJ) propagation velocity of detonation wave is 1808.0 m/s under the same conditions, so the error is only 0.8%. With the parameters calculated in the numerical simulation, the combustion thermal efficiency η r pressurization ratio π 3 c , the nonuniformities of outlet static pressure NUP and outlet static temperature NUT can be obtained using the following equations:
π 3 c = P 3 c * / P 2 m *
η r = ( H 3 c * − H 2 m * ) / Q i n
N U P = ∑ i = 1 n ( P i − P ¯ ) 2 / ( n − 1 ) / P ¯
N U T = ∑ i = 1 n ( T i − T ¯ ) 2 / ( n − 1 ) / T ¯
where P 3 c * and H 3 c * are the outlet total pressure and outlet total enthalpy at the state point of 3c, Q i n represents the heat release of fuel, and P ¯ , T ¯ are average static pressure and temperature of RDC outlet.
In the conventional approximate iso-pressure combustion progress (2→30) of the reference gas turbine cycle, there was an important parameter—the total pressure recovery coefficient σ —which was applied in the followed equation and defined as 0.934:
σ = P 3 , 0 * / P 2 *
  • Mixer model (3c→3)
A mixing process happens in the mixer between the gas of the RDC outlet and the air of the compressor outlet. Two assumptions were made in the mixer model, one was that the all mechanical losses were ignored, the other was that the development of an oblique shock wave stopped in the mixer. Thus, the key parameters in this model can be calculated by:
g a   =   g a , m +   g a , c
( g a , c + g f ) P 3 c * + g a , m P 2 * = ( g a + g f ) P 3 *
where g a , c is the mass flow rate of the combustion air, g a , m is the mass flow rate of the mixing air, g f is the mass flow rate of the fuel, and P 3 * , T 3 * , π 3 respectively represent the total pressure, total temperature and composite pressure ratio at the state point of 3.
  • Turbine model (3→4/30→40)
As the loss coefficient of exhaust pressure ξ o u t was known, and expansion ratio ( π t ) and efficiency ( η t ) of turbine were:
π t = P 3 * P 2 * π c ξ o u t
η t = H 3 * − H 4 * H 3 * − H 4 s *
where H 3 * and H 4 * are the enthalpies at state points of 3 and 4. H 4 s * is the ideal enthalpy in isentropic condition (the equations were the same in the 30→40 case)
Based on the above models and the necessary iterative computations, the following thermodynamic cycle performance parameters can be obtained:
Compressor input power L C :
L C = g a , c ∫ T 1 * T 2 * C P ( T ) d T   +   g a ∫ T 1 * T 2 m * C P ( T ) d T
Turbine output power L T :
L T = ( g a + g f ) ∫ T 4 * T 3 * C P ( T ) d T
Power generation L E :
L E = L T   −   L C
Cycle efficiency η :
η = L E / Q i n
Cycle efficiency increasing rate w Δ η :
w Δ η = η − η 0 η 0 × 100 %
In addition, as the fuel of the gas turbine can be completely combusted, the gas composition mainly included N2, O2, CO2, and H2O. In view of the additivity of C P for the mixed gas, the following equation was used:
C P ( T ) = 8.3145 ( A 1 + A 2 T + A 3 T 2 + A 4 T 3 + A 5 T 4 )
where A is coefficient related to composition, obtained in [50].
The enthalpy of each point was calculated by:
H = ∫ T r e f T C P ( T ) d T
Entropy of any point S :
S = ∫ T ref T C P ( T ) T d T − R g ln P P r e f
where T r e f is the reference temperature, 288.15 K, and P r e f is the reference pressure, 0.1013 MPa.
The efficiency characteristics of the compressor and turbine were referenced from universal characteristic lines, obtained in [51,52].

3. Results and Discussion

3.1. Flow Field Characteristic Research of RDC

Figure 4 presents entropy shadowgraphs at t = 6800 μs. In the axial direction, the height of the detonation wave was about 95 mm, the height of the oblique shock wave was about 90 mm, and the height of the transition region was about 15 mm at different compressor pressure ratios. The six path lines passed through the midpoint of detonation wave (type 1 of path lines) and oblique shock wave (type 2 of path lines) in the axial direction, respectively. The path lines were obtained by transformation of coordinates, which was explained in details in [19,41]. Path lines ①, ③, ⑤ passed through detonation wave at about Y = 50 mm and a weak shock wave at about Y = 150 mm. Path lines ②, ④, ⑥ passed through the detonation wave at about Y = 7.8 mm, weak shock waves at about Y = 25 mm, 150 mm and the oblique shock wave at about Y = 155 mm.
The changes of static temperature and static pressure along path lines ①, ③, ⑤ are shown in Figure 5a. Path lines passed through the detonation wave at about Y = 47.5 mm, where the pressure and temperature both displayed jumps and then gradually decreased. At Y = 150 mm, the pressure and temperature shows a smaller jump due to the fact the path lines are passing through a weak shock wave, where the pressure increased by 8.07%, 5.22%, 8.27% and the temperature increased by 1.66%, 1.16%, 1.75%. In all, the pressure decreased by 27.34%, 26.96% 30.44% from the beginning to the end, respectively.
Figure 5b presents the changes of static temperature and static pressure along path lines ②, ④, ⑥. Pressure and temperature had the first big jumps when the path lines passed through the detonation wave at about Y = 7.8 mm. Then the two parameters dropped extremely. Path lines ②, ④, ⑥ passed through the first weak shock wave at Y = 25 mm. Although the two parameters had no significant change, the slope of downtrend slowed down obviously. Path lines ②, ④, ⑥ passed through the second weak shock wave at Y = 150 mm, where the two parameters decreased less than 1%. Pressure and temperature had their second big jumps when the path lines passed through the oblique shock wave at Y = 155 mm. In all, the pressure increased by 96.05%, 99.45%, 96.30% from the beginning to the end, respectively.
Figure 5c shows the change of the entropy along the six path lines. The same as the change trend of temperature, big jumps happened when the path lines passed through the detonation wave. At the moment, the entropy change reached about 1900 J/kg·K along path lines ①, ③, ⑤ and about 1850 J/kg·K along path lines ②, ④, ⑥. Another jump happened when the path lines ②, ④, ⑥ passed through the oblique shock wave. As for path lines ②, ④, ⑥, the major entropy change was rooted in the detonation region and the second entropy change led by the oblique shock wave was about 6.61% of the total. In all, the entropy change reached about 1900 J/kg·K and 2000 J/kg·K for the two styles of path lines from the beginning to the end, respectively. In addition, although the weak shocks caused sudden changes in the pressure and temperature diagram, it had no significant effect on the entropy change.
As there was a significant difference in the pressure changes along the two types of path lines, mass-weighted average parameters of the RDC cross-section in the axial direction were analyzed, as shown in Figure 6a. Combined with the heights of the detonation wave, oblique shock wave and transtion region as mentioned above, the total pressure P * and total temperature T * increased significantly due to the existence of the detonation wave, and the total pressure reached the highest value at the end point of the detonation wave at Y = 95 mm. Entering into the transtion region including deflagration combustion, the total pressure started to decrease and the total temperature reached its peak at the end of transtion region when the combustion was completed. In the region including the oblique shock wave, the total pressure continued to decrease but the total temperature remain unchanged. In general, there were significant pressurization characteristics in the RDC, of which the instantaneous pressurization ratios respectively reached 1.7757, 1.8005, and 1.8150 at different compressor ratios. Figure 6b presents the entropy of the RDC cross-section in the axial direction. The entropy kept increasing from the inlet to the outlet of RDC. The slope was the largest in the detonation combustion region, decreased in the transition region and reached a minimum in the oblique shock wave region.

3.2. Variations of RDC Characteristic Parameters at Different Compressor Pressure Ratios

The mass flow rate of air was divided into two parts: the mass flow rate of bleeding air from the 3/4 compressor which will enter into the RDC was defined as g a , c , and the mass flow rate of the higher pressure air from the compressor outlet was defined as g a , m . Table 2 presents the mass flow rates and distribution proportions of the two parts of air when the equivalent ratio of the RDC is 1.
Figure 7 shows the variations of outlet total pressure P 3 c * and outlet total temperature T 3 c * of RDC at different compressor pressure ratios. All the following data of the RDC outlet were calculated using the mass-weighted and time average method [32]. The two parameters both had a downtrend with the decrease of compressor pressure ratio. When the compressor pressure ratio decreased from 13.66 to 10.02, P 3 c * and T 3 c * decreased from 1.2561 MPa, 2513.3 K to 1.0383 MPa, 2482.95 K, respectively.
In order to further investigate the performance of RDC, the pressurization ratio and combustion thermal efficiency of the RDC were calculated under different working conditions according to the Equations (11) and (12), as shown in Figure 8. It can be seen that as the compressor pressure ratio decreased, η r of the RDC did not change significantly, all being beyond 0.9971. Additionally, π 3 c gradually increased with the decrease of compressor pressure ratio. At π c = 13.66, π 3 c reached 1.7791. When π c decreased to 10.02, π 3 c increased by 4.3%, reaching 1.8545. From the view of variation of flow field parameters, higher temperature jumps and pressure jumps in the detonation wave would promote a higher pressure gain in the RDC [53]. Therefore, variations of temperature and pressure jumps were calculated along each of the two types of path lines at different compressor pressure ratios, as shown in Figure 9. Choices of the typical path lines were same as those in Section 3.1, respectively, passing through the midpoint of the detonation wave (type 1 of path lines in Figure 9) and the midpoint of the oblique shock wave (type 2 of path lines in Figure 9). With the decrease of compressor pressure ratio, pressure and temperature jumps at the detonation wave both had an uptrend, which was in accordance with the variation trend seen in Figure 8.
Furthermore, the entropy change was analyzed to explain the reason for the variation of pressure gain characteristics. As the pressure and temperature both had obvious influences on the entropy change, the entropy change of iso-pressure combustion at the same compressor pressure ratio was used for reference to eliminate the influence of temperature, so the entropy change will only be reflected in the change of pressure. Figure 10 shows the related parameters of the entropy change. As the compressor pressure ratio decreased, both the entropy change of RDC S3c and the entropy change of iso-pressure combustion S3i increased. The difference between them ΔS3c (ΔS3c = S3i − S3c) also had the same uptrend, which meant that the available energy loss savings in the rotating detonation combustion gradually increased with the decrease of compressor pressure ratio, and was finally reflected in the higher and higher total pressure of the RDC outlet. The change of entropy change increased from 176.14 J/kg·K to 188.83 J/kg·K when the compressor pressure ratio decreased from 13.66 to 10.02.
In addition, the non-uniformities of the outlet static pressure and outlet static temperature were calculated at different compressor pressure ratios. As shown in Figure 11, the nonuniformities of outlet static pressure NUP and outlet static temperature NUT both had uptrends with the decrease of compressor pressure ratio, but with a small variation range. However, the non-uniformity of outlet static pressure has already reached to 0.46, which sets a quite high demand for the design of the mixer and first stage blade of the turbine.

3.3. Investigation of Cycle Performance in Gas Turbine Based on Interstage Bleeding Rotating Detonation Combustion

Figure 12 shows the parameters related to the turbine inlet total temperature and turbine inlet total pressure of the rotating detonation gas turbine and reference gas turbine at different pressure ratios.
It can be seen in Figure 12a that as the compressor pressure ratio decreased from 13.66 to 10.22, the turbine inlet total temperature of the rotating detonation gas turbine decreased from 1382.90 K to 1165.69 K. Moreover, there was a decreasing difference compared with those of the reference gas turbine, but all beyond 19 K. As shown in Figure 12b, with the decrease of compressor pressure ratio, the difference of turbine inlet total pressure between the reference gas turbine and rotating detonation gas turbine gradually increased, but with a decreasing slope. The reason for this phenomenon was that when the compressor pressure ratio was too low, the pressurization ratio of the RDC would be greater than the pressure ratio of the last 1/4 compressor. At that moment, the pressure of air from the compressor outlet is less than the pressure of gas from the RDC outlet, which could have a negative effect on the pressure in the mixer. The occurrence of this phenomenon led to the uncertainty about the change of cycle efficiency, which will be analyzed in details in the last part of this paper. Figure 13 shows the difference of entropy change between the reference gas turbine and the rotating detonation gas turbine, including the differences in the whole working progress ΔSW, compressor ΔS2, combustor with mixer ΔS3, and turbine ΔS4. With the decrease of compressor pressure ratio, ΔS2 gradually decreased, ΔS4 increased, ΔS3 first increased and then decreased with a turning point at π c = 11.44. The turbine played a decisive role in the entropy change of whole working progress, as the uptrend of ΔSW was similar to ΔS4. Although a turning point appeared in the variation trend of the combustor with mixer, there was no indication that the slope of ΔSW was slowing down.
Figure 14 shows compressor input power and turbine output power of the rotating detonation gas turbine and the reference gas turbine at different compressor pressure ratios. It can be clearly seen there is no significant difference of turbine output power between the two gas turbines, but the compressor input power of the rotating detonation gas turbine was less than that of the reference gas turbine, especially at high compressor pressure ratio. Figure 15 presents the variation of power generation at different compressor pressure ratios. With the decrease of compressor pressure ratio, the difference ΔLE between the power generation of the rotating detonation gas turbine LE and the power generation of the reference gas turbine LE,0 gradually decreased, but the power generation enhancement rate w Δ L E had the opposite trend. ΔLE reached 498.95 kW at π c = 13.66, and the power generation of the rotating detonation gas turbine reached 7915.20 kW at that moment. When the compressor pressure ratio decreased to 10.02, a difference of power generation of about 441.18 kW still can be obtained. Meanwhile, the power generation enhancement rate increased from 6.72% to 11.89%. Figure 16 shows the related parameters of cycle efficiency at different compressor pressure ratios. At π c = 13.66, the cycle efficiency of the rotating detonation gas turbine η and reference gas turbine η 0 were respectively 0.3493 and 0.3273, with an efficiency enhancement rate ( w Δ η ) of about 6.72%. The two cycle efficiencies decreased to 0.3073 and 0.2731 at π c = 10.02, but with an increasing efficiency enhancement rate of about 11.89%. In all the calculated cases, the efficiency enhancement rates of the rotating detonation gas turbine were all greater than 6.72% and increased with the decrease of compressor pressure ratio. In addition, the slope of the efficiency increase rate was approximately stable and had no downward trend with the decrease of compressor pressure ratio. In the previous discussion, the air pressure from the compressor outlet was less than the pressure of gas from the RDC outlet when the compressor pressure ratio was too low, which could have a negative effect on the pressure in the mixer. However, the changes of efficiency enhancement rate showed that the high pressure gain characteristic at low compressor pressure ratio still played a leading role in the variation of cycle efficiency.
The new style gas turbine cycle based on interstage bleeding rotating detonation combustion finally brought an increase about 0.03 in the cycle efficiency of the gas turbine. There was still a gap compared with that described in [32], but this does not mean the research was meaningless. Firstly, from the viewpoint of structural design, enormous convenience was brought in the design of mixer as the pressure of the mixing air (pressurizing in the last 1/4 compressor) was greater than or close to the RDC outlet pressure.
Moreover, the air of interstage bleeding for anti-surge could be better utilized, reducing the loss in the cycle. In summary, it was believed that interstage bleeding rotating detonation combustion could play an important role as an auxiliary combustion system in the gas turbine cycle. Therefore, a gas turbine cycle with combined combustion type (including conventional approximate iso-pressure combustion and rotating detonation combustion) should be further investigated in the future.

4. Conclusions

In this paper, a gas turbine cycle model based on interstage bleeding rotating detonation combustion was established using methane as fuel. Combined with a series of two-dimensional numerical simulations and calculations of cycle parameters, variations of pressure gain characteristic and cycle performance were investigated at different compressor pressure ratios. The main conclusions are as follows:
(1)
There were significant differences of pressure variation along different types of path lines, but the pressurization characteristics of the RDC were obvious in general.
(2)
Compared with the conventional gas turbine cycle, the non-negligible pressure gain characteristics acquired in the interstage bleeding rotating detonation combustor contributed to an obvious performance improvement in the rotating detonation gas turbine cycle. The power generation enhancements were all beyond 400 kW and the cycle efficiency enhancement rates were all greater than 6.72%. Additionally, the turbine inlet total temperature showed a decrease of over 19 K.
(3)
Due to the differences of entropy change between iso-pressure combustion and rotating detonation combustion, the pressurization ratio of the RDC increased with the decrease of compressor ratio, which would have a leading influence on the performance improvement in the rotating detonation gas turbine cycle. When the compressor pressure ratio decreased from 13.66 to 10.02, the pressurization ratio of the RDC increased from 1.7791 to 1.8545, then the power generation enhancement rate and cycle efficiency enhancement rate increased from 6.72% to 11.89%.
(4)
Compared with the reference gas turbine cycle, the difference of turbine entropy change in the rotating detonation gas turbine cycle played a decisive role in the entropy change of whole working process.

Author Contributions

L.Q. performed calculations of cycle parameters and result analysis. Y.D. and Q.M. performed the numerical simulations. H.Z. was responsible for the numerical approach. Z.W. and N.Z. wrote the manuscript. All authors have read and approved the final manuscript.

Funding

This work was supported in part by National Nature Science Foundation of China (Grant No. 51709059), the Province Nature Science Foundation of Heilongjiang (Grant No. QC2017045) and the Fundamental Research Funds for the Central Universities (Grant No. HEUCFJ170304, Grant No. HEUCFP201719, Grant No. HEUCFM180302).

Conflicts of Interest

We declare that we have no conflict of interest.

References

  1. Durmusoglu, Y.; Ust, Y. Thermodynamic optimization of an irreversible regenerative closed Brayton cycle based on thermoeconomic performance criterion. Appl. Math. Model. 2014, 38, 5174–5186. [Google Scholar] [CrossRef]
  2. Liu, Y.Z.; Sun, X.X.; Sethi, V.; Nalianda, D.; Li, Y.G.; Wang, L. Review of modern low emissions combustion technologies for aero gas turbine engines. Prog. Aerosp. Sci. 2017, 94, 12–45. [Google Scholar] [CrossRef]
  3. Lu, F.K.; Braun, E.M. Rotating detonation wave propulsion: Experimental challenges, modeling, and engine concepts. J. Propuls. Power 2014, 30, 1125–1142. [Google Scholar] [CrossRef]
  4. Paxson, D.E.; Perkins, H.D. Thermal load considerations for detonative combustion-based gas turbine engines. In Proceedings of the 40th AIAA Aerospace Sciences Meeting, AIAA 2014-3396, Fort Lauderdale, FL, USA, 11–14 July 2004. [Google Scholar]
  5. Kailasanath, K. Recent developments in the research on rotating-detonation wave engines. In Proceedings of the 55th AIAA Aerospace Sciences Meeting, AIAA 2017-0784, Grapevine, TX, USA, 9–13 January 2017. [Google Scholar]
  6. Chan, S.; Liu, H. Mass-based design and optimization of wave rotors for gas turbine engine enhancement. Shock Waves 2016, 2, 1–12. [Google Scholar] [CrossRef]
  7. Bobusch, B.C.; Berndt, P.; Paschereit, C.O.; Klein, R. Shockless explosion combustion: An innovative way of efficient constant volume combustion in gas turbines. Combust. Sci. Technol. 2014, 186, 1680–1689. [Google Scholar] [CrossRef]
  8. Berndt, P.; Klein, R. Modeling the kinetics of the shockless explosion combustion. Combust. Flame 2017, 175, 16–26. [Google Scholar] [CrossRef]
  9. Wintenberger, E.; Sheperd, J.E. Thermodynamic cycle of analysis for propagating detonations. Thermodynamic Cycle Analysis for Propagating Detonations. J. Propuls. Power 2006, 22, 694–698. [Google Scholar] [CrossRef]
  10. Wintenberger, E.; Sheperd, J.E. A model for the performance of air-breathing pulse detonation engines. In Proceedings of the 39th AIAA Aerospace Sciences Meeting, AIAA 2003-4511, Huntsville, AL, USA, 20–23 July 2003. [Google Scholar]
  11. Wolański, P. Detonative propulsion. Proc. Combust. Inst. 2013, 34, 125–158. [Google Scholar] [CrossRef]
  12. Yao, S.B.; Ma, Z.; Zhang, S.J.; Luan, M.Y. Reinitiation phenomenon in hydrogen-air rotating detonation engine. Int. J. Hydrog. Energy 2017, 42, 28588–28598. [Google Scholar] [CrossRef]
  13. Meng, Q.Y.; Zhao, N.B.; Zheng, H.T.; Yang, J.L.; Qi, L. Numerical investigation of the effect of inlet mass flow rates on H2/air non-premixed rotating detonation wave. Int. J. Hydrog. Energy 2018, 43, 13618–13631. [Google Scholar] [CrossRef]
  14. Lin, W.; Zhou, J.; Liu, S.J.; Lin, Z.Y. An experimental study on CH4/O2 continuously rotating detonation wave in a hollow combustion chamber. Exp. Therm. Fluid Sci. 2015, 62, 122–130. [Google Scholar] [CrossRef]
  15. Bluemner, R.; Bohon, M.D.; Paschereit, C.O. Dynamics of counter-rotating wave modes in an RDC. In Proceedings of the 56th AIAA Aerospace Sciences Meeting, AIAA 2018-4572, Cincinnati, OH, USA, 9–11 July 2018. [Google Scholar]
  16. Zhou, R.; Wang, J.P. Numerical investigation of flow particle paths and thermodynamic performance of continuously rotating detonation engines. Combust. Flame 2012, 159, 3632–3645. [Google Scholar] [CrossRef]
  17. Paxson, D.E.; Fotia, M.; Hoke, J.; Schauer, F. Comparison of numerically simulated and experimentally measured performance of a rotating detonation engine. In Proceedings of the 53rd AIAA Aerospace Sciences Meeting, AIAA 2015-1101, Kissimmee, FL, USA, 5–9 January 2015. [Google Scholar]
  18. Paxson, D.E. Impact of an exhaust throat on semi-idealized rotating detonation engine performance. In Proceedings of the 54th AIAA Aerospace Sciences Meeting, AIAA 2016-1647, San Diego, CA, USA, 4–8 January 2016. [Google Scholar]
  19. Zheng, H.T.; Qi, L.; Zhao, N.B.; Li, Z.M.; Liu, X. A thermodynamic analysis of the pressure gain of continuously rotating detonation combustor for gas turbine. Appl. Sci. Basel 2018, 8, 535. [Google Scholar] [CrossRef]
  20. Keisuke, G.; Kato, Y.; Ishihara, K.; Matsuoka, K.; Kasahara, J.; Matsuo, A.; Funaki, I. Experimental study of effects of injector configurations on rotating detonation engine performance. In Proceedings of the 54th AIAA Aerospace Sciences Meeting, AIAA 2016-5100, San Diego, CA, USA, 4–8 January 2016. [Google Scholar]
  21. Braun, J.; Saracoglu, B.H.; Paniagua, G. Unsteady performance of rotating detonation engines with different exhaust nozzles. J. Propuls. Power 2016, 33, 1–10. [Google Scholar] [CrossRef]
  22. Fotia, M.L.; Schauer, F.; Kaemming, T.; Hoke, J. Experimental study of the performance of a rotating detonation engine with nozzle. J. Propuls. Power 2015, 31, 674–681. [Google Scholar] [CrossRef]
  23. Wolański, P. Application of the continuous rotating detonation to gas turbine. Appl. Mech. Mater. 2015, 782, 3–12. [Google Scholar] [CrossRef]
  24. Li, J.M.; Teo, C.J.; Khoo, B.C.; Wang, J.P.; Wang, C. Detonation Control for Propulsion: Pulse Detonation and Rotating Detonation Engines (Shock Wave and High Pressure Phenomena), 1st ed.; Springer International Publishing AG: Cham, Switzerland, 2018; pp. 30–35. ISBN 978-3-319-68905-0. [Google Scholar]
  25. Frolov, S.M.; Aksenov, V.S.; Ivanov, V.S.; Shamshin, I.O. Large-scale hydrogen–air continuous detonation combustor. Int. J. Hydrog. Energy 2015, 40, 1616–1623. [Google Scholar] [CrossRef]
  26. Frolov, S.M.; Aksenov, V.S.; Ivanov, V.S. Three-dimensional numerical simulation of the operation of a rotating-detonation chamber with separate supply of fuel and oxidizer. Russ. J. Phys. Chem. B 2013, 7, 35–43. [Google Scholar] [CrossRef]
  27. Frolov, S.M.; Aksenov, V.S.; Ivanov, V.S. Experimental proof of Zel'dovich cycle efficiency gain over cycle with constant pressure combustion for hydrogen–oxygen fuel mixture. Int. J. Hydrog. Energy 2015, 40, 6970–6975. [Google Scholar] [CrossRef]
  28. Frolov, S.M.; Dubrovskii, A.V.; Ivanov, V.S. Three-dimensional numerical simulation of a continuously rotating detonation in the annular combustion chamber with a wide gap and separate delivery of fuel and oxidizer. Prog. Propuls. Phys. 2016, 8, 375–388. [Google Scholar] [CrossRef]
  29. Naples, A.; Hoke, J.; Schauer, F.R. Rotating detonation engine interaction with an annular ejector. In Proceedings of the 52nd AIAA Aerospace Sciences Meeting, AIAA 2014-0287, National Harbor, MD, USA, 13–17 January 2014. [Google Scholar]
  30. Naples, A.; Hoke, J.; Battelle, R.T.; Wagner, M.; Schauer, F.R. RDE Implementation into an Open-Loop T63 Gas Turbine Engine. In Proceedings of the 55th AIAA Aerospace Sciences Meeting, AIAA 2017-1747, Grapevine, TX, USA, 9–13 January 2017. [Google Scholar]
  31. Sousa, J.; Paniagua, G.; Morata, E.C. Thermodynamic analysis of a gas turbine engine with a rotating detonation combustor. Appl. Energy 2017, 195, 247–256. [Google Scholar] [CrossRef]
  32. Qi, L.; Zhao, N.B.; Wang, Z.T.; Yang, J.L.; Zheng, H.T. Pressure gain characteristic of continuously rotating detonation combustion and its influence on gas turbine cycle performance. IEEE Access 2018, 6, 70236–70247. [Google Scholar] [CrossRef]
  33. Zhou, S.B.; Ma, H.; Li, S.; Liu, D.K.; Yan, Y.; Zhou, C.S. Effects of a turbine guide vane on hydrogen-air rotating detonation wave propagation characteristics. Int. J. Hydrog. Energy 2017, 42, 20297–20305. [Google Scholar] [CrossRef]
  34. Zhou, S.B.; Ma, H.; Ma, Y.; Zhou, C.S.; Liu, D.K.; Li, S. Experimental study on a rotating detonation combustor with an axial-flow turbine. Acta Astronaut. 2018, 151, 7–14. [Google Scholar] [CrossRef]
  35. Liu, Q. Research on Fuel Reforming and Rotating Detonation Pressurized Combustion Technique of Combustor. Ph.D. Thesis, Harbin Engineering University, Harbin, China, 2013. [Google Scholar]
  36. Gavrikov, A.I.; Efimenko, A.A.; Dorofeev, S.B. A model for detonation cell size prediction from chemical kinetics. Combust. Flame 2000, 120, 19–33. [Google Scholar] [CrossRef]
  37. Bykovskii, F.A.; Zhdan, S.A.; Vedernikov, E.F. Continuous spin detonations. J. Propuls. Power 2006, 22, 1204–1216. [Google Scholar] [CrossRef]
  38. Zhe, L.; James, B.; Guillermo, P. Characterization of a supersonic turbine downstream of a rotating detonation combustor. J. Eng. Gas Turbines Power 2018, 141, 1–13. [Google Scholar] [CrossRef]
  39. Wu, D.; Liu, Y.; Liu, Y.S.; Wang, J.P. Numerical investigations of the restabilization of hydrogen-air rotating detonation engines. Int. J. Hydrog. Energy 2014, 39, 15803–15809. [Google Scholar] [CrossRef]
  40. Tsuboi, N.; Watanabe, Y.; Kojima, T.; Hayashi, A.K. Numerical estimation of the thrust performance on a rotating detonation engine for a hydrogen–oxygen mixture. Proc. Combust. Inst. 2015, 35, 2005–2013. [Google Scholar] [CrossRef]
  41. Kaemming, T.; Fotia, M.L.; Hoke, J.; Schauer, F. Thermodynamic modeling of a rotating detonation engine through a reduced-order approach. J. Propuls. Power 2017, 33, 1–9. [Google Scholar] [CrossRef]
  42. Yao, S.B.; Tang, X.M.; Wang, J.P.; Shao, Y.T.; Zhou, R. Three-dimensional numerical study of flow particle paths in rotating detonation engine with a hollow combustor. Combust. Sci. Technol. 2017, 189, 965–979. [Google Scholar] [CrossRef]
  43. Oran, E.S.; Weber, J.W., Jr.; Stefaniw, E.I. A numerical study of a two-dimensional H2-O2-Ar detonation using a detailed chemical reaction model. Combust. Flame 1998, 113, 147–163. [Google Scholar] [CrossRef]
  44. Nordeen, C.; Schwer, D.; Schauer, F.; Hoke, J.; Cetegen, B.; Barber, T. Thermodynamic modeling of a rotating detonation engine. In Proceedings of the 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, AIAA 2011-803, Orlando, FL, USA, 4–7 January 2011. [Google Scholar]
  45. Uemura, Y.; Hayashi, A.K.; Asahara, M.; Tsuboi, N.; Yamada, E. Transverse wave generation mechanism in rotating detonation. Proc. Combust. Inst. 2013, 34, 1981–1989. [Google Scholar] [CrossRef]
  46. Graf, A.P.; Riedel, U. Numerical simulation of supersonic reactive flows using explicit Runge–Kutta methods. In Proceedings of the 38th AIAA Aerospace Sciences Meeting, AIAA 2000-438, Reno, NV, USA, 10–13 January 2000. [Google Scholar]
  47. Bykovskii, F.A.; Vedernikov, E.F. Continuous detonation of a subsonic flow of a propellant. Combust. Explos. Shock Waves 2003, 39, 323–334. [Google Scholar] [CrossRef]
  48. Schwer, D.; Kailasanath, K. Numerical study of the effects of engine size on rotating detonation engines. In Proceedings of the 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, AIAA 2011-581, Orlando, FL, USA, 4–7 January 2011. [Google Scholar]
  49. Schwer, D.; Kailasanath, K. Numerical investigation of the physics of rotating-detonation-engines. Proc. Combust. Inst. 2011, 33, 2195–2202. [Google Scholar] [CrossRef]
  50. Yao, Q.; Li, S.Q.; Wang, Y. An Introduction to Combustion: Concepts and Applications, 2nd ed.; Tsinghua University Press: Beijing, China, 2009; pp. 490–516. ISBN 978-7-302-18762-2. [Google Scholar]
  51. Wang, Z.T.; Li, J.; Fan, K.; Ma, W.Q.; Lei, H.D. Prediction method for low speed characteristics of compressor based on modified similarity theory with genetic algorithm. IEEE Access 2018, 6, 36834–36839. [Google Scholar] [CrossRef]
  52. Wang, Z.T.; Li, J.; Fan, K.; Li, S.Y. The off-design performance simulation of marine gas turbine based on optimum scheduling of variable stator vanes. Math. Probl. Eng. 2017, 2017, 2671251. [Google Scholar] [CrossRef]
  53. Sousa, J.; Braun, J.; Paniagua, G. Development of a fast evaluation tool for rotating detonation combustor. Appl. Math. Model. 2017, 52, 42–52. [Google Scholar] [CrossRef]
Entropy 21 00265 g001 550
Figure 1. Cycle scheme of a gas turbine with two combustion modes.
Figure 1. Cycle scheme of a gas turbine with two combustion modes.
Entropy 21 00265 g001
Entropy 21 00265 g002 550
Figure 2. T-S diagrams of a marine gas turbine cycle with different combustion processes.
Figure 2. T-S diagrams of a marine gas turbine cycle with different combustion processes.
Entropy 21 00265 g002
Entropy 21 00265 g003 550
Figure 3. Temperature shadowgraphs of RDC.
Figure 3. Temperature shadowgraphs of RDC.
Entropy 21 00265 g003
Entropy 21 00265 g004 550
Figure 4. Entropy shadowgraphs of the RDC at 6800 μs.
Figure 4. Entropy shadowgraphs of the RDC at 6800 μs.
Entropy 21 00265 g004
Entropy 21 00265 g005 550
Figure 5. Variations of pressure, temperature and entropy along different path lines.
Figure 5. Variations of pressure, temperature and entropy along different path lines.
Entropy 21 00265 g005
Entropy 21 00265 g006 550
Figure 6. Mass-weighted average parameters of the RDC cross-section.
Figure 6. Mass-weighted average parameters of the RDC cross-section.
Entropy 21 00265 g006
Entropy 21 00265 g007 550
Figure 7. Variations of T 3 c * and P 3 c * at different compressor pressure ratios.
Figure 7. Variations of T 3 c * and P 3 c * at different compressor pressure ratios.
Entropy 21 00265 g007
Entropy 21 00265 g008 550
Figure 8. Variations of π 3 c and η r at different compressor pressure ratios.
Figure 8. Variations of π 3 c and η r at different compressor pressure ratios.
Entropy 21 00265 g008
Entropy 21 00265 g009 550
Figure 9. Variations of temperature and pressure jumps at different compressor pressure ratios.
Figure 9. Variations of temperature and pressure jumps at different compressor pressure ratios.
Entropy 21 00265 g009
Entropy 21 00265 g010 550
Figure 10. Variations of related parameters of the entropy change at different compressor pressure ratios (ΔS3c = S3i − S3c).
Figure 10. Variations of related parameters of the entropy change at different compressor pressure ratios (ΔS3c = S3i − S3c).
Entropy 21 00265 g010
Entropy 21 00265 g011 550
Figure 11. Non-uniformities of the outlet static pressure and outlet static temperature.
Figure 11. Non-uniformities of the outlet static pressure and outlet static temperature.
Entropy 21 00265 g011
Entropy 21 00265 g012 550
Figure 12. Variations of turbine inlet total temperature and turbine inlet total pressure at different compressor pressure ratios.
Figure 12. Variations of turbine inlet total temperature and turbine inlet total pressure at different compressor pressure ratios.
Entropy 21 00265 g012
Entropy 21 00265 g013 550
Figure 13. Difference of entropy change between the reference gas turbine and the rotating detonation gas turbine.
Figure 13. Difference of entropy change between the reference gas turbine and the rotating detonation gas turbine.
Entropy 21 00265 g013
Entropy 21 00265 g014 550
Figure 14. Variations of compressor input power and turbine output power at different compressor pressure ratios.
Figure 14. Variations of compressor input power and turbine output power at different compressor pressure ratios.
Entropy 21 00265 g014
Entropy 21 00265 g015 550
Figure 15. Variation of power generation at different compressor pressure ratios.
Figure 15. Variation of power generation at different compressor pressure ratios.
Entropy 21 00265 g015
Entropy 21 00265 g016 550
Figure 16. Variation of cycle efficiency at different compressor pressure ratios.
Figure 16. Variation of cycle efficiency at different compressor pressure ratios.
Entropy 21 00265 g016
Table
Table 1. Relative parameters of the compressor.
Table 1. Relative parameters of the compressor.
Load100%90%80%70%60%50%
L E , 0 [kW]7416.256674.635933.005191.384449.753708.13
π c 13.6612.9012.1711.4410.7110.02
P 2 * [Mpa]1.35731.28181.20911.13651.06410.9962
T 2 * [K]688.32677.44665.22652.89639.65627.18
P 2 m * [Mpa]0.70600.67630.64740.61800.58820.5598
T 2 m * [K]562.64555.61547.71539.71531.39522.96
g a [kg/s]26.818525.677324.614623.621822.688621.8212
Table
Table 2. Related parameters of the mass flow rate at different compressor pressure ratios.
Table 2. Related parameters of the mass flow rate at different compressor pressure ratios.
π c 13.6612.9012.1711.4410.7110.03
g f [kg/s]0.48630.44560.40550.36640.32840.2914
g a , c [kg/s]8.34837.64926.95986.28975.63715.0013
g a , m [kg/s]18.470218.028117.654917.332117.051516.8199
g a , m / g a , c 2.21242.23682.53672.75563.02493.3631

Share and Cite

MDPI and ACS Style

Qi, L.; Wang, Z.; Zhao, N.; Dai, Y.; Zheng, H.; Meng, Q. Investigation of the Pressure Gain Characteristics and Cycle Performance in Gas Turbines Based on Interstage Bleeding Rotating Detonation Combustion. Entropy 2019, 21, 265. https://doi.org/10.3390/e21030265

AMA Style

Qi L, Wang Z, Zhao N, Dai Y, Zheng H, Meng Q. Investigation of the Pressure Gain Characteristics and Cycle Performance in Gas Turbines Based on Interstage Bleeding Rotating Detonation Combustion. Entropy. 2019; 21(3):265. https://doi.org/10.3390/e21030265

Chicago/Turabian Style

Qi, Lei, Zhitao Wang, Ningbo Zhao, Yongqiang Dai, Hongtao Zheng, and Qingyang Meng. 2019. "Investigation of the Pressure Gain Characteristics and Cycle Performance in Gas Turbines Based on Interstage Bleeding Rotating Detonation Combustion" Entropy 21, no. 3: 265. https://doi.org/10.3390/e21030265

APA Style

Qi, L., Wang, Z., Zhao, N., Dai, Y., Zheng, H., & Meng, Q. (2019). Investigation of the Pressure Gain Characteristics and Cycle Performance in Gas Turbines Based on Interstage Bleeding Rotating Detonation Combustion. Entropy, 21(3), 265. https://doi.org/10.3390/e21030265

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop