Next Article in Journal
Lifespan Development of the Human Brain Revealed by Large-Scale Network Eigen-Entropy
Next Article in Special Issue
Second Law Analysis for Couple Stress Fluid Flow through a Porous Medium with Constant Heat Flux
Previous Article in Journal
Life on the Edge: Latching Dynamics in a Potts Neural Network
Previous Article in Special Issue
A Numerical Study on Entropy Generation in Two-Dimensional Rayleigh-Bénard Convection at Different Prandtl Number
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Second-Law Analysis of Irreversible Losses in Gas Turbines

1
Center of Applied Space Technology and Microgravity (ZARM), the University of Bremen, 28359 Bremen, Germany
2
Institute of Engineering Thermophysics, Chinese Academy of Sciences, Beijing 100190, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
*
Authors to whom correspondence should be addressed.
Entropy 2017, 19(9), 470; https://doi.org/10.3390/e19090470
Submission received: 26 July 2017 / Revised: 29 August 2017 / Accepted: 31 August 2017 / Published: 4 September 2017
(This article belongs to the Special Issue Entropy in Computational Fluid Dynamics)

Abstract

:
Several fundamental concepts with respect to the second-law analysis (SLA) of the turbulent flows in gas turbines are discussed in this study. Entropy and exergy equations for compressible/incompressible flows in a rotating/non-rotating frame have been derived. The exergy transformation efficiency of a gas turbine as well as the exergy transformation number for a single process step have been proposed. The exergy transformation number will indicate the overall performance of a single process in a gas turbine, including the local irreversible losses in it and its contribution to the exergy obtained the combustion chamber. A more general formula for calculating local entropy generation rate densities is suggested. A test case of a compressor cascade has been employed to demonstrate the application of the developed concepts.

1. Introduction

A gas turbine is a type of combustion engine, which usually comprises an upstream compressor, a downstream turbine, and a combustion chamber in between. The design of a high performance gas turbine is a key technology in modern industry due to its significant applications in aircraft, electrical generators, and ships.
Compared to reciprocating engines, gas turbines have the advantages of a higher power-weight ratio, a smaller volume, lower toxic emissions, etc. However, the efficiency of gas turbines is usually lower than that of reciprocating engines. Generally, the efficiency of a gas turbine can be improved by increasing the pressure ratio [1] and reducing the engine weight [2]. However, according to the second law of thermodynamics, the efficiency of a gas turbine is ultimately determined by the irreversibility in its flow and temperature field. For example, if a gas turbine has a pressure ratio of 20, its efficiency according to an ideal Brayton cycle is about 60%, whereas the efficiency even for a modern gas turbine is no more than 40%. A gas turbine may be combined with a steam plant to form a “combined cycle” system. In recent years, the efficiency of the gas turbine combined cycle system’s power generation has been increased to 60% [3]. However, this only increases the system’s overall efficiency, while the efficiency of the gas turbine is not really improved. The key to improving the efficiency of a gas turbine still lies in how to reduce the losses due to irreversibility.
In an ideal gas turbine, the working gases follow the Brayton cycle, which is composed of the procedures of isentropic compression, isobaric (constant pressure) combustion, isentropic expansion, and isobaric heat rejection: see the cycle 0–1–2–3–0 in Figure 1. The efficiency for the ideal Brayton cycle is:
η i d e a l = Q ˙ c b c Q ˙ 3 0 Q ˙ c b c = 1 T 0 T 1
where Q ˙ c b c and Q ˙ 3 0 are the heat generation and release rates in the processes of isobaric combustion (1–2) and isobaric heat rejection (3–0).
However, irreversibility occurs in a real gas turbine. The T-S diagram for a real Brayton cycle is indicated by 0–1’–2’–3’–0 in Figure 1, in which only the irreversibility of the compression (0–1’) and expansion (2’–3’) processes are taken into account. The efficiency of the real Brayton cycle is
η r e a l = Q ˙ c b c Q ˙ 3 0 Q ˙ 1 2 .
When the compressor and the turbine are considered to be adiabatic, the entropy will increase in processes 0–1’ and 2’–3’ due to irreversibility, leading to higher entropy at state 3’ than at state 3. As a result, the heat release in process 3’–0 is larger than in the ideal cycle (3–0), i.e., Q ˙ 3 0 > Q ˙ 3 0 . Thus, the efficiency of a real gas turbine is always lower than that of an ideal one.
The second-law analysis (SLA) is a very helpful tool for understanding these irreversible processes. In the SLA, irreversibility in both flow and temperature fields are accounted for. However, the SLA was usually used for analyzing a thermal system, while the detailed flow and temperature fields were traditionally studied within the disciplines of fluid mechanics and heat transfer, in which the SLA still has not received much attention. Entropy almost never appears in the text books of fluid dynamics and heat transfer: see [4,5] as examples. This concept is ignored perhaps due to the reason that the irreversibility effect is considered to be not important in these two disciplines. However, turbulent flow and heat transfer are typical irreversible processes due to the dissipation in the flow field and irreversibility in the temperature fields. Another possible reason is that it is very difficult to calculate or measure the local entropy generation rate accurately due to model and experimental errors.
Some progress has been made in employing the SLA in flow and heat transfer problems since the 1980s. Bejan [6,7] laid the foundation with respect to analyzing and optimizing thermal systems with the SLA approach. Kock and Herwig [8] extended this concept to a detailed analysis of turbulent flows, and identified four different mechanisms of entropy generation: dissipation in a mean and fluctuating velocity field and heat flux in a mean and fluctuating temperature field. Later, Kock and Herwig [9] developed equations for computing entropy generation rates for Reynolds Averaged Navier–Stokes Simulations (RANSs) and implemented them into computational fluid dynamics (CFD) codes. Jin and Herwig [10] indicated that model errors in RANS simulations may lead to considerable uncertainties in entropy generation results.
With the development of high performance computers in recent years, people have started to calculate the losses in blade cascades from RANS results. Orhan [11] investigated the loss mechanism of an axial turbine cascade with the SLA. Denton and Pullan [12] and Zlatinov et al. [13] analyzed the local entropy generation rate in turbines with unsteady RANSs. However, Kopriva et al. [14,15] found that RANSs have much lower accuracy than large eddy simulations (LESs). Tucker [16,17] indicated that it is particularly difficult to predict unsteady separations in turbines with RANSs. As a compromise between the computational cost and the accuracy, Lin et al. [18] studied the local entropy generation in a turbine cascade passage with a delayed detached eddy simulation (DDES) method. The losses in both flow and temperature fields were visualized according to the numerical results. With the same method, Wang et al. [19] analyzed the interaction between the corner separation and wakes in a compressor cascade. The detailed coherent structures, local losses information, and turbulence characteristics were identified according to the local entropy generation rate. Despite this progress, more systematic studies are still required for understanding the irreversible processes in gas turbines. In particular, the relationship between local entropy generation and global efficiency should be better understood.
In the present paper, we attempt to investigate several fundamental concepts with respect to the SLA of the irreversible processes in gas turbines. Entropy and exergy transport equations for compressible flows in a rotating frame will be derived in Section 2. Through the derivation, we will show the relationship between the SLA and the other laws in fluid mechanics. In Section 3, we will discuss the dimensionless coefficients for assessing irreversible processes in a gas turbine. The concepts of exergy transformation efficiency and exergy transformation number will be introduced. CFD modeling for calculating the local entropy generate rate density will be discussed in Section 4. A test case for applying the developed theories will be provided in Section 5. The conclusions are given in Section 6.

2. The Entropy and Exergy Transport Equations for Compressible Flows in a Rotating Frame

The governing equations for the flows in a cascade passage or the combustion chamber are the compressible Navier–Stokes and energy equations. The equations in a rotating frame were adopted, and thus they can be also used for the flows in rotors. The governing equations [20,21] read:
ρ t + ( ρ u i ) x i = 0
( ρ u i ) t + ( ρ u j u i ) x j = p x i + τ j i x j + ρ f i
( ρ e t ) t + ( ρ u i h t ) x i = x i ( λ θ x i ) + x j ( τ j i · u i )
The chemical reaction in the combustion chamber is not considered here for simplicity, while the combustion heat is accounted for by the heat exchange at the boundary walls. The reference frame rotates with the angular velocity of ω . The value of ω is zero when a stationary frame is under consideration. The reference frame velocity is v = ω × r , where r is the displacement from the axis to the position vector x . If x 1 is selected as the axis, we have r = ( 0 ,   x 2 , x 3 ) . f i is the component of the body acceleration rate vector F = ( 2 ω × u + ω × v ) , in which the Coriolis and centrifugal forces are taken into account.
ρ , v i , u i , θ and λ are the density, frame velocity component, relative velocity component, temperature, and thermal conductivity, respectively. The relative total energy and enthalpy are calculated by e t = e + 1 2 ( u i 2 v i 2 ) and h t = e + p ρ + 1 2 ( u i 2 v i 2 ) , respectively. τ i j = 2 ν s i j is the anisotropic part of the viscous stress tensor, where the strain rate s i j = 1 2 ( u i x j + u j x i ) .
Multiplying Equation (4) with u i and considering ρ f i u i = ( 1 2 ρ v j 2 ) / t + ( 1 2 ρ v j 2 u i ) / x i , the transport equation of the relative kinetic energy k r can be derived:
( ρ k r ) t + ( ρ u i k r ) x i = u i p x i + u i τ j i x j
where k r = 1 2 ( u i 2 v i 2 ) . Subtracting Equation (6) from Equation (5), the transport equation of the enthalpy can be obtained, i.e.,
( ρ h ) t + ( ρ u i h ) x i = d p d t + x i ( λ θ x i ) + τ j i u i x j .
Substituting Equation (7) into the fundamental equation of thermodynamics T d s = d q = d h 1 ρ d p , the entropy equation can be derived, i.e.,
ρ d s d t = ( ρ s ) t + ( ρ u i s ) x i = x i ( λ θ θ x i ) + λ θ 2 θ x i θ x i + 2 μ s i j s i j T .
The irreversibility is due to the last two terms in Equation (8), which are always positive. They are called the entropy generation rate density in the temperature field S ˙ C and in the flow field S ˙ D . Their definitions are:
S ˙ C λ θ 2 θ x i θ x i ;   S ˙ D 2 μ s i j s i j T .
For an open system, the total differential of exergy of a gas is d h E = d h T d s + d k , where T is the environmental temperature, h E is the specific exergy of a gas, and k = 1 2 ( u i + v i ) 2 is the total kinetic energy. The transport equation of k can be obtained from Equation (4), i.e.,
( ρ k ) t + ( ρ u j k ) x j = ( u i + v i ) p x i + ( u i + v i ) τ j i x j .
Substituting Equations (7), (8), and (10) into the total derivative of exergy, the balance equation of exergy can be obtained, i.e.,
ρ d h E d t = ( ρ h E ) t + ( ρ u i h E ) x i = p t + x i ( ( 1 T θ ) λ θ x i ) + x j ( τ j i u i ) + x j ( ( τ j i p δ j i ) v i ) T ( S ˙ C + S ˙ D ) .
Equations (8) and (11) were derived for compressible flows in a rotating frame, but they can also be used for incompressible flows in a non-rotating frame when ρ is a constant and v i is zero.
Integrating Equations (8) and (11) in a volume under consideration, e.g., a cascade passage or the combustion chamber, the integral forms of the entropy and exergy equations can be obtained. They are:
t V ( ρ s ) d V + A ρ u i s n i d A = S ˙ + Q ˙ S
t V ( ρ h E p ) d V + A ( ρ u i h E τ i j u j ) n i d A = T S ˙ W ˙ t + Q ˙ E
where Δ S ˙ = V ( S ˙ C + S ˙ D ) d V is the overall entropy generation rate in domain V , n i is the direction vector of a surface element, Q ˙ S = A λ θ θ x i n i d A is the exchange of entropy due to heat transfer at the boundary, W t = A ( τ j i + p δ j i ) v i n j d A is the shaft work when a rotor is accounted for, and Q ˙ E = A ( 1 T θ ) λ T x i n i d A is the increase of exergy due to heat transfer. When the compressor and turbine are assumed to be adiabatic, Q ˙ E is non-zero only in the combustion chamber.
Equation (12) shows that the entropy will always increase or remain constant in an isolated system in which Q ˙ S and A ρ u i s n i d A are zero. The corresponding destruction of exergy according to Equation (13) is V T ( S ˙ C + S ˙ D ) d V . This is accordance with the second law of thermodynamics. The derivation shows that the traditional laws in fluid mechanics and heat transfer are sufficient for satisfying the second law of thermodynamics. It is not necessary to solve the entropy or exergy equations in numerical simulations. The entropy generation rate can be calculated as post processing of the numerical results.

3. Dimensionless Coefficients for Assessing Irreversible Processes

The bulk quantitates S , H , T , and P can be calculated by the statistics of the local flow quantities s , h , θ , and p , respectively. For example, the bulk temperature is calculated by T = A ρ u i θ n i d A / m ˙ , where A is the cross section under consideration. The mass flow rate is m ˙ = A ρ u i n i d A . The losses in a certain process step, such as in a cascade passage, can be evaluated with these bulk quantities.
The loss in a cascade passage was traditionally evaluated by the enthalpy loss coefficient [22], which is defined by:
ζ h = { H ˙ o u t H ˙ o u t , i s H ˙ i n l * H ˙ i n l   for   a   compressor   cascade ; H ˙ o u t H ˙ o u t , i s H ˙ o u t * H ˙ o u t for   a   turbine   cascade
where H ˙ is the rate of specific enthalpy, and the subscripts “ o u t ” and “ i s ” denote the outlet and isentropic process, respectively. The superscript * denotes the stagnation value.
In a real application, instead of Equation (14), it is more convenient to use the stagnation pressure loss coefficient [22] instead of ζ h to indicate the loss, i.e.,
ζ p = { P i n l * P o u t * P i n l * P i n l   for   a   compressor   cascade ; P i n l * P o u t * P o u t * P o u t for   a   turbine   cascade
where the subscript “ i n l ” denotes the inlet. Obviously, the irreversibility due to heat transfer is not taken into account in Equations (14) and (15). However, as shown in Figure 1, the increase of entropy due to heat transfer will also reduce the efficiency of a gas turbine. The losses in a cascade can be more accurately assessed by the SLA. Denton [23] suggested the use of an entropy loss coefficient to indicate the loss of efficiency in a cascade, i.e.,
ζ s = { T o u t Δ S ˙ H ˙ i n l * H ˙ i n l   for   a   compressor   cascade ; T o u t Δ S ˙ H ˙ o u t * H ˙ o u t for   a   turbine   cascade
where Δ s ˙ is the entropy generation rate in a cascade under consideration. The losses of exergy due to irreversibility are more accurately assessed by Equation (16) than by Equation (15). However, the contribution of the cascade to the output work is not taken into account by Equation (16). According to Equation (16), the optimized blade cascades should have low loading and work at a small incidence angle, thus the entropy generation can be minimized. This conflicts with some real applications in which higher loading blade cascades which work at higher incidence angles are preferred.
In order to assess the irreversible processes in gas turbines more comprehensively, we adopted and further developed the concepts of entropic potential and energy devaluation number, which have been proposed by Herwig and his colleagues [24,25,26] in recent years. According to these studies, the entropic potential is defined by the entropy generation rate by which the entropy of the ambient is increased when the primary energy rate E ˙ becomes part of its internal energy, i.e.,
S ˙ i r r = E ˙ T
where T is the environmental temperature. The amount of the entropic potential rate of E ˙ that is consumed by the process step i under consideration can be determined by the energy devaluation number, which is defined by
N i Δ S ˙ i S ˙ i r r = T Δ S ˙ i E ˙ .
More details can be found in [24,25,26].
In a gas turbine, the primary energy is the exergy Q ˙ E which is obtained in the combustion chamber in process 1’–2’. Q ˙ E was derived and defined in Equation (13). In order to simplify its calculation, we approximate the local wall temperature θ w in the combustion chamber with the bulk temperature T in the cross section which is enclosed by the wall surface. Thus, the primary energy can be calculated by
E ˙ = Q ˙ E = 1 2 ( 1 T T ) d Q ˙ = Q ˙ c b c T ( S 2 S 1 ) m ˙
where Q ˙ E is the exergy obtained in the combustion chamber, and m ˙ is the mass flow rate.
Time averaging and summing up the integral exergy equations (Equation (13)) of all the processes in a gas turbine and dividing it with Q ˙ E , we have
i ( T Δ S ˙ i ) Q ˙ E + i W ˙ t i Q ˙ E + Q ˙ 3 0 E Q ˙ E = 1
as the chain of energy devaluation and transportation. In this chain, T Δ S ˙ i Q ˙ E = N i is the energy devaluation number of process i defined by Equation (18). The exergy obtained in the combustion chamber Q ˙ E is devaluated by i ( T Δ S ˙ i ) , transformed to shaft work by i W ˙ t i , or transported to the environment by Q ˙ 3 0 E . Q ˙ 3 0 E occurs in process 3 0 (see Figure 1) due to the release of heat, which is not a real loss of exergy and this part of exergy can be (totally by an ideal process and partly by a real process) regained through a gas turbine combined cycle (GTCC) power generation system. Thus, the exergy transformed by the gas turbine is i W ˙ t i + E ˙ 3 0 E . Multiplying Equation (20) with Q ˙ E / Q ˙ c b c , we may define the exergy transformation efficiency of a gas turbine by:
η E ( i W ˙ t i + Q ˙ 3 0 E ) / Q ˙ c b c = Q ˙ E / Q ˙ c b c T i Δ S ˙ i / Q ˙ c b c .
η E indicates the fraction of the heat Q ˙ c b c which can be transformed to exergy. η E is higher than η r e a l defined by Equation (2) by Q ˙ 3 0 E / Q ˙ c b c .
η E is influenced by both the entropy generation rate Δ S ˙ i in each process step i and its contribution to Q ˙ E . Thus, we may define the exergy transformation number of process step i by:
η i E Q ˙ i E / Q ˙ c b c T Δ S ˙ i / Q ˙ c b c
where Q ˙ i E is the contribution of process i to Q ˙ E . That is to say, although all of the exergy Q ˙ E is obtained in the combustion chamber, the combustion chamber is not the only contributor to Q ˙ E . The compressor upstream also has important effects on Q ˙ E .
When η i E of each process is known, we have η E = i η i E as the chain of exergy transformation. In a gas turbine, we use η t u r , i E , η c o m , i E , and η c b c E to denote the exergy transformation numbers of turbine cascade i, compressor cascade i, and the combustion chamber.
η t u r , i E of turbine cascade i is determined only by the entropy generation rate in it, i.e.,
η t u r , i E = T Δ S ˙ t u r , i / Q ˙ c b c .
Compressor cascade i may have two opposite effects on exergy transformation. On the one hand, similar to a turbine cascade, the exergy is destroyed due to entropy generation. On the other hand, the exergy obtained in the combustion chamber is increased due to the increase of the temperature and the pressure through the compressor cascade. Under this consideration, η c o m , i E is calculated by
η c o m , i E = Q ˙ c o m , i E / Q ˙ c b c T Δ S ˙ c o m , i / Q ˙ c b c
where Δ S ˙ c o m , i is the entropy generation rate in compressor cascade i, and Q ˙ c o m , i E is the contribution of compressor cascade i to the exergy obtained in the combustion chamber Q ˙ E . Figure 2 shows the influence of compressor cascade i on the exergy obtained in the combustion chamber schematically: the static temperature is increased from T 1 , i 1 to T 1 , i through compressor cascade i. The potential exergy that can be obtained in the combustion chamber through a reversible process is increased from Q ˙ c b c , i 1 E = Q ˙ c b c m ˙ T T 1 , i 1 T 2 , i 1 c p T d T to Q ˙ c b c , i E = Q ˙ c b c m ˙ T T 1 , i T 2 , i c p T d T . T 2 , i 1 and T 2 , i are the temperature at the exit of the combustion chamber when the temperature at the inlet of the combustion chamber is T 1 , i 1 and T 1 , i , respectively. Thus, the contribution of compressor cascade i to Q ˙ E is calculated by
Q ˙ c o m , i E = Q ˙ c b c , i E Q ˙ c b c , i 1 E = m ˙ T T 1 , i 1 T 2 , i 1 c p T d T m ˙ T T 1 , i T 2 , i c p T d T
with T 2 , i 1 and T 2 , i being determined according to the first law of thermodynamics, by
T 1 , i 1 T 2 , i 1 c p d T = T 1 , i T 2 , i c p d T = Q ˙ c b c m ˙ .
When an ideal gas with a constant capacity c p is taken into account, Equations (25) and (26) can be simplified to
Q ˙ c o m , i E = m ˙ c p T l n ( T 1 , i 1 + Q ˙ c b c / ( c p m ˙ ) T 1 , i 1 ) m ˙ c p T l n ( T 1 , i + Q ˙ c b c / ( c p m ˙ ) T 1 , i ) .
η c b c E of the combustion chamber is calculated by
η c b c E = Q ˙ c b c E / Q ˙ c b c T Δ S ˙ c b c / Q ˙ c b c
where Δ S ˙ c b t is the entropy generation rate in the combustion chamber, and Q ˙ c b c E is the increase of exergy in the combustion chamber without the upstream compressor through a reversible process. It is calculated by
Q ˙ c b c E = Q ˙ c b c T T 1 , 0 T 2 , 0 c p T d T ;  with  T 1 , 0 T 2 , 0 c p d T = Q ˙ c b c m ˙
where T 1 , 0 = T 0 is the compressor inlet temperature (state 0 is indicated in Figure 1), and T 2 , 0 is temperature at the exit of the combustion chamber when the inlet temperature of the combustion chamber is T 0 . Obviously, we have Q ˙ c b c E + i Q ˙ c o m , i E = Q ˙ E .
When the efficiency coefficients of all the processes are known, we have
η c b c E + i η c o m , i E + i η t u r , i E = η E
as a chain of exergy transformation. The exergy transformation number η i E can be used to assess an isolated process step, since only local flow and temperature fields are needed to calculate its value.

4. Modeling of Local Entropy Generation Rate Densities

The entropy generation rate densities S ˙ D and S ˙ C can be directly calculated from direct numerical simulation (DNS) results. However, DNS requires very high computational costs and is thus not suitable for engineering applications. When RANS or other models’ results are used in the study, part of the losses cannot be calculated directly and thus must be modeled. Kock and Herwig [8,9] proposed models for calculating entropy generation rates with RANS results. In these models, S ˙ D and S ˙ C were decomposed by:
S ˙ D = 2 μ s ¯ i j s ¯ i j θ ¯ + ε θ ¯
S ˙ C = λ θ ¯ 2 θ ¯ x i θ ¯ x i + ε θ θ ¯ 2
where ε = 2 μ s i j s i j ¯ is the turbulent dissipation rate, and ε θ = λ θ x i θ x i ¯ is the temperature fluctuation dissipation rate. However, ε and ε θ are not always calculated explicitly in RANS models, and even when they are determined, their accuracies are low since they often are only used as intermediate quantities for calculating the eddy viscosity or Reynolds stresses: see [27].
Instead, we assume that, when the flow domain under consideration is sufficiently large, the produced turbulence is in balance with the dissipation, i.e.,
V ε d V V P r o d V = V ρ u i u j ¯ s ¯ i j d V
V ε θ d V V P r o θ d V = V ρ θ u i ¯ θ ¯ x i V .
With this assumption, ε and ε θ can be replaced with the turbulence production rate P r o and the temperature fluctuation production rate P r o θ . According to the eddy viscosity hypothesis, the effect of turbulence on momentum and heat transfer can be approximated in a similar way as the molecular diffusion: see [27]. Thus, S ˙ D and S ˙ C can be calculated by
S ˙ D 2 ( μ + μ T ) s ¯ i j s ¯ i j θ ¯
S ˙ C ( λ + λ T ) θ ¯ 2 θ ¯ x i θ ¯ x i
where μ T is the eddy viscosity, and λ T = μ T Pr T is the turbulent thermal conductivity. Pr T is the turbulent Prandtl number. Equations (35) and (36) can be calculated directly in all RANS models which are based on the eddy viscosity hypothesis.

5. A Test Case of Application

As an example for applying the concepts developed in Section 2, Section 3 and Section 4, we analyzed the turbulent flow in an isolated compressor cascade which was taken from the experimental database of [28] with a RANS method. Through this low cost test case, we will show how to analyze the numerical results with the concepts developed in Section 2, Section 3 and Section 4.
The computational domain is half of the cascade passage. The geometric parameters of the cascade are shown in Table 1. In order to reduce the boundary effects, the inlet and outlet regions were extended by 1.59C and 2C, respectively, where C is the length of the chord. The velocity profile at the inlet was given according to the experimental data in [28]. The time and surface averaged inlet velocity u i n l and turbulent intensity are 40 ms−1 and 0.8%, respectively. The turbulent intensity was calculated by ( 2 3 k T , i n l ) 1 / 2 / u i n l , where k T , i n l is the inlet turbulent kinetic energy. The inlet-specific dissipation rate ω T , i n l is 3.1 × 10 3 s−1. Four incidence angles were accounted for in the present study. They are 0°, 2°, 4°, and 7°. Since the flow is at a small velocity (the Mach number is smaller than 0.3), the fluid in the cascade is assumed to be incompressible with the constant density of 1.217 kg m−3.
The flow was assumed to be quasi-steady. The following steady Reynolds averaged Navier–Stokes equations were solved during the simulation:
u ¯ i x i = 0
ρ ( u ¯ j u ¯ i ) x j = p ¯ x i + x j ( ( μ + μ T ) u ¯ i x j ) .
The eddy viscosity μ T was calculated with the k-ω Shear-Stress Transport (SST) turbulence model [29].
An open source CFD software, OpenFoam v16.06+, was used to carry out the simulation. PimpleFoam was selected as the computational solver. This solver is based on a pressure correction method for incompressible flows. The second-order upwind scheme was used for spatial discretization. Body-fitted mesh, which concentrates near the wall, was adopted in the study. The dimensionless mesh spacing y + of the first grid point near the wall is ensured to be smaller than 1 to resolve the turbulent boundary layer. The mesh in the region close to the cascade’s trailing edge was refined to capture the corner separation. The mesh in a cross section is shown in Figure 3. A typical mesh has about 5.8 million grid points. The mesh independence study was performed to ensure the results are mesh resolution independent. More computational details can be found in [30].
The entropy generation in the temperature field was neglected in the present test case, since it is much smaller than the one in the flow field. The global static temperature at the inlet and outlet were approximated according to the ideal gas law, i.e., T = P / ( ρ R ) . The specific gas constant R is 287.1 J kg−1 K−1. Since the density ρ is a constant, the global temperature is proportional to the global static pressure.
The entropy generation rate density in the flow field S ˙ D was calculated with Equations (31) and (35). When Equation (31) is adopted, the turbulent dissipation rate ε must be determined corresponding to the specific turbulent model. For the k-ω SST turbulence model in use, ε is calculated by
ε = ρ β * k ω
where β * is the constant used in the k-ω SST turbulence model. It is determined empirically according to experimental and DNS data. Reference [28] suggests its value to be 0.09.
The static pressure coefficient C p at two sections ( x 3 / C = 50 % and 29.7%) are shown in Figure 4. The gap between the current CFD results and the experimental data in [28] is due to the model error. The accuracy can be further improved by using more accurate CFD methods, e.g., the large eddy simulation (LES) method. However, the current low cost CFD method is sufficient for the purpose of this study, i.e., demonstrating how to analyze cascade flows with the developed concepts.
Both model results at the incidence angle of 4° are shown in Figure 5, which indicates that the two models predict similar patterns of entropy generation. S ˙ D close to the hub is stronger according to Equation (35) (see Figure 5f) than Equation (31) (see Figure 5e), since the turbulence production is stronger than the turbulence dissipation in this region. In other words, not all of the produced turbulence is dissipated locally. Similar phenomena can be found at other incidence angles: see Figure 6. However, according to our assumption, the turbulence production is in balance with the turbulence dissipation when the domain size is sufficiently large. This assumption was validated by our numerical results: the volume-integrated entropy generation rates by Equations (31) and (35) can be found in Table 2, which are almost identical. Compared with Equation (31), Equation (35) is more general and is independent of turbulence models.
In the current test case, we are only able to compare the entropy generation rate in the flow field since the energy equation was not solved. More systematic studies, e.g., cascade flows at high Mach numbers, are still required to validate the equivalence between Equation (32) and Equation (36). Besides the entropy generation rates, the other integral quantities for calculating the dimensionless coefficients are also provided in Table 2.
Fluid properties in the combustion chamber, including the specific combustion heat Q ˙ 1 / m ˙ , the heat capacity c p , and the environmental temperature T , are required to determine the local exergy transformation number η c o m , i E , which is defined by Equation (24). The values of these parameters are shown in Table 3.
The global coefficients ζ p , ζ s , and η c o m , i according to the RANS results are shown in Table 4. Both the total pressure loss coefficient ζ p and the entropy loss coefficient ζ s indicate that losses in the cascade increase with the incidence angle. However, the exergy transformation coefficient η c o m , i E suggests that the optimal incidence angle is 2 o , at which cascade works with the best overall performance: although more losses due to irreversibility occur at 2 o than at 0 o , a larger pressure ratio is obtained through the cascade, which increases the potential exergy obtained in the combustion chamber Q c o m , i E .
The exergy transformation number η c o m , i is linearly related with the environmental temperature T , thus T does not influence the optimal results. Figure 7 shows the influence of the other reference parameters Q ˙ c b c / m ˙ and c p on η c o m , i : η c o m , i decreases with Q ˙ c b c / m ˙ and increases with c p . However, the optimal results are not influenced by Q ˙ c b c / m ˙ and c p when they are only mildly changed.

6. Conclusions

Several fundamental aspects with respect to the SLA of the turbulent flows in gas turbines were discussed in this study.
Entropy and exergy equations (Equations (8) and (11)) for compressible/incompressible flows in a rotating/non-rotating frame were derived. The derivation shows that the Navier–Stokes equations and the energy equation are sufficient to satisfy the second law of thermodynamics, thus it is not necessary to solve the entropy and exergy equations to evaluate their quantities. The entropy and exergy can be determined by the post processing of CFD simulations. However, the entropy and exergy equations and their budgets are helpful tools for analyzing the irreversible processes in gas turbines.
The exergy transformation efficiency η E of a gas turbine as well as the exergy transformation number η i E of a single process step i were proposed in this study. η i E in a turbine cascade, a compressor cascade, or the combustion chamber are suggested to be calculated by Equations (23), (24), or (28). The value of η i E indicates the overall effects of an irreversible process, including its destruction of exergy and its contribution to the potential exergy obtained in the combustion chamber. η i E can be used to assess the performance of an isolated process in a gas turbine, since only local flow and temperature fields are required to calculate its value.
The methods for calculating the local entropy generation rate densities were discussed. It was suggested to use turbulence production rates (Equations (35) and (36)) instead of the turbulence dissipation rates (Equations (31) and (32)) to calculate the local entropy generation rate densities. The assumption behind this approximation is that the turbulence production rate is in balance with the turbulence dissipation rate when the domain is sufficiently large. An advantage of Equations (35) and (36) is that they are independent from the choices of turbulence models. However, more systematic studies, e.g., LESs of cascade flows at high Mach numbers, are still required to further validate these equations.
A test case with respect to a compressor cascade has been employed for applying the concepts developed in the study. The numerical results show that the entropy generation rates calculated by Equations (35) and (31) are almost identical. The exergy transformation number suggests an optimal incidence angle at which the compressor cascade works with the best overall performance.

Acknowledgments

The authors gratefully acknowledge the support of this study by the DFG-Heisenberg program (JI 253/1) and the grants of National Natural Science Foundation of China (No. 51676184 and No. 51506195). The acknowledgement is also given to X. Ottavy of Ecole Centrale de Lyon for the data support of our test case and H. Herwig of Hamburg University of Technology for the helpful discussion with him.

Author Contributions

Yan Jin proposed the models and derived the mathematical equations. Juan Du, Zhiyuan Li, and Hongwu Zhang performed the CFD simulation. Juan Du and Yan Jin wrote the paper together. Juan Du, Hongwu Zhang, and Yan Jin had deep discussion about the model and the numerical results. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

e local specific internal energy, J kg−1
e E local specific exergy of the fluid, J kg−1
e t local relative total energy, J kg−1
f i body acceleration rate component, m s−2
F body acceleration rate vector, m s−2
E ˙ bulk specific primary energy rate, J kg−1 s−1
h local specific enthalpy, J kg−1
h E local specific exergy of a gas in an open system, J kg−1
k local kinetic energy, m2 s−2
k T , i n l inlet specific turbulent kinetic energy, m2 s−2
m ˙ mass flow rate, kg s−1
n i direction vector of a surface element
l b l d length of the blade span, m
N i energy devaluation number of process step i
N b l d blade number in a cascade
p local pressure, N m−2
P bulk pressure, N m−2
P r o local turbulence production rate, kg m−1 s−3
P r o θ local temperature fluctuation production rate, J K m−3 s−1
Q ˙ c b c heat generation rate in the combustion chamber, J s−1
Q ˙ 3 0 heat release rate in a reversible Brayton cycle, J s−1
Q ˙ 3 0 heat release rate in an irreversible Brayton cycle, J s−1
Q ˙ 3 0 E exergy of the heat release rate Q ˙ 3 0 , J s−1
Q ˙ E exergy of the heat generation rate Q ˙ c b c , J s−1
Q ˙ i E contribution of process i to Q ˙ E , J s−1
R specific gas constant, J kg−1 K−1
s local specific entropy, J kg−1 K−1
s i j strain rate tensor, s−1
S bulk specific entropy, J kg−1 K−1
S ˙ C entropy generation rate density in the flow field, J K−1 m−3 s−1
S ˙ D entropy generation rate density in the temperature field, J K−1 m−3 s−1
S ˙ i r r bulk specific entropy potential rate, J kg−1 K−1 s−1
T bulk temperature, K
u local velocity vector, m s−1
u i local velocity component, m s−1
v local frame velocity vector, m s−1
v i local frame velocity component, m s−1
W ˙ t shaft power of a rotor, J s−1
x displacement vector, m
x i coordinate vector component, m
Greek symbols
S ˙ i entropy generational rate in process i, J s−1 K−1
ε local turbulence dissipation rate, kg m−1 s−3
ε θ local temperature fluctuation dissipation rate, kg K m−1 s−3
λ thermal conductivity, J s−1 m−1 K−1
μ dynamic viscosity, kg m−1 s−1
η i d e a l efficiency of an ideal Brayton cycle
η r e a l efficiency of a real Brayton cycle
η E exergy transformation efficiency
η i E exergy transformation number of process step i
θ local temperature
τ i j viscous stress tensor, k gm−1 s−2
ω angular velocity vector, s−1
ω T ,   i n l inlet specific dissipation rate in the k- ω SST turbulence model, s−1
ζ h enthalpy loss coefficient
ζ p stagnation pressure loss coefficient
ζ s entropy loss coefficient
Subscripts
1 , 2 , 3 , 4 four states of an ideal Brayton cycle
A cross section
b l d blade number in a cascade
c b c combustion chamber
c o m compressor
C temperature field
D flow field
i process step index or vector component index
i n l inlet
i s isentropic value
j vector component index
o u t outlet
r relative value
t total value
T turbulence
t u r turbine
w wall surface
environmental value
Superscripts
E exergy
S entropy
*stagnation value
irreversible process
  ¯ Reynolds averaging

References

  1. Wang, H. Experimental and Numerical Research of Highly Loaded Axial-Flow Compressor Stages. Ph.D. Thesis, Beihang University, Beijing, China, 2009. [Google Scholar]
  2. Gao, F. Advanced Numerical Simulation of Corner Separation in a Linear Compressor Cascade. Ph.D. Thesis, Ecole Centrale de Lyon, Écully, France, 2014. [Google Scholar]
  3. Satoshi, H.; Tsukagoshi, K.; Masada, J.; Ito, E. Test Results of the World’s First 1,600C J-series Gas Turbine. Mitsubishi Heavy Ind. Tech. Rev. 2012, 49, 18–23. [Google Scholar]
  4. Batchelor, G.K. An Introduction to Fluid Dynamics; Cambridge Mathematical Library: Cambridge, UK, 2000. [Google Scholar]
  5. Incropera, F.; DeWitt, D.; Bergmann, T.; Lavine, A. Fundamentals of Heat and Mass Transfer, 6th ed.; John Wiley & Sons: New York, NY, USA, 2006. [Google Scholar]
  6. Bejan, A. Entropy Generation through Heat and Fluid Flow; John Wiley & Sons: New York, NY, USA, 1982. [Google Scholar]
  7. Bejan, A. Entropy Generation Minimization; CRC Press: Boca Raton, FL, USA, 1996. [Google Scholar]
  8. Kock, F.; Herwig, H. Local entropy production in turbulent shear flows: A high-Reynolds number model with wall functions. Int. J. Heat Mass Trans. 2004, 47, 2205–2215. [Google Scholar] [CrossRef]
  9. Kock, F.; Herwig, H. Entropy production calculation for turbulent shear flows and their implementation in CFD codes. Int. J. Heat Fluid Flow 2005, 26, 672–680. [Google Scholar] [CrossRef]
  10. Jin, Y.; Herwig, H. Turbulent flow in rough wall channels: Validation of RANS models. Comp. Fluids 2015, 122, 34–46. [Google Scholar] [CrossRef]
  11. Orhan, O.E. Investigation of the Effect of Turbulence on Entropy Generation in Turbomachinery. Ph.D. Thesis, Middle East Technical University, Ankara, Turkey, 2014. [Google Scholar]
  12. Denton, J.; Pullan, G. A numerical investigation into the sources of endwall loss in axial flow turbines. In Proceedings of the ASME Turbo Expo 2012: Turbine Technical Conference and Exposition, Copenhagen, Denmark, 11–15 June 2012; pp. 1417–1430. [Google Scholar]
  13. Zlatinov, M.B.; Tan, C.S.; Montgomery, M.; Islam, T.; Harris, M. Turbine hub and shroud sealing flow loss mechanisms. J. Turbomach. 2012, 134, 061027. [Google Scholar] [CrossRef]
  14. Kopriva, J.E.; Laskowski, G.M.; Sheikhi, M.R.H. Computational assessment of inlet turbulence on boundary layer development and momentum/thermal wakes for high pressure turbine nozzle and blade. In Proceedings of the ASME 2014 International Mechanical Engineering Congress and Exposition, Montreal, QC, Canada, 14–20 November 2014. [Google Scholar]
  15. Laskowski, G.M.; Kopriva, J.; Michelassi, V.; Shankaran, S.; Paliath, U.; Bhaskaran, R.; Wang, Q.; Talnikar, C.; Wang, Z.; Jia, F. Future directions of high-fidelity CFD for aero-thermal turbomachinery research, analysis and design. In Proceedings of the 46th AIAA Fluid Dynamics Conference, Washington, DC, USA, 13–17 June 2016. [Google Scholar]
  16. Tucker, P.G. Computation of unsteady turbomachinery flows: Part 1—Progress and challenges. Prog. Aerosp. Sci. 2011, 47, 522–545. [Google Scholar] [CrossRef]
  17. Tucker, P.G. Computation of unsteady turbomachinery flows: Part 2—LES and hybrids. Prog. Aerosp. Sci. 2011, 47, 546–569. [Google Scholar] [CrossRef]
  18. Lin, D.; Yuan, X.; Su, X.R. Local Entropy generation in compressible flow through a high pressure turbine with delayed detached eddy simulation. Entropy 2017, 19, 29. [Google Scholar] [CrossRef]
  19. Wang, H.; Lin, D.; Su, X.R.; Yuan, X. Entropy analysis of the interaction between the corner separation and wakes in a compressor cascade. Entropy 2017, 19, 324. [Google Scholar] [CrossRef]
  20. ANSYS FLUENT 13.0 User’s Guide, ANSYS Inc.: Canonsburg, PA, USA, 2010.
  21. Cariglino, F. External Aerodynamic Simulations in a rotating Frame of Reference. Ph.D. Thesis, Politecnico di Torino, Turin, Italy, 2013. [Google Scholar]
  22. Brown, L. Axial flow compressor and turbine loss coefficients: A comparison of several parameters. J. Eng. Power 1972, 94, 193–201. [Google Scholar] [CrossRef]
  23. Denton, J.D. Loss mechanisms in turbomachines. In Proceedings of the ASME 1993 International GasTurbine and Aeroengine Congress and Exposition, Cincinnati, OH, USA, 24–27 May 1993. [Google Scholar]
  24. Redecker, C.; Herwig, H. Calculating and assessing complex convective heat transfer problems: The CFD-SLA approach. In Proceedings of the 15th International Heat Transfer Conference, Kyoto, Japan, 10–15 August 2014. [Google Scholar]
  25. Wenterodt, T.; Redecker, C.; Herwig, H. Second law analysis for sustainable heat and energy transfer: The entropic potential concept. Appl. Energy 2015, 139, 376–383. [Google Scholar] [CrossRef]
  26. Herwig, H.; Redecker, C. Heat transfer and entropy. In Heat Transfer Studies and Applications; Kazi, S.N., Ed.; InTech: Rijeka, Croatia, 2015; pp. 143–161. [Google Scholar]
  27. Pope, S.B. Turbulent Flows; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  28. Ma, W. Experimental Investigation of Corner Stall in a Linear Compressor Cascade. Ph.D. Thesis, Ecole Centrale de Lyon, Écully, France, 2012. [Google Scholar]
  29. Menter, F.R. Two equation Eddy viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef]
  30. Li, Z.Y.; Du, J.; Jemcov, A.; Ottavy, X.; Lin, F. A study of loss mechanism in a linear compressor cascade at the corner stall condition. In Proceedings of the ASME Turbo Expo 2017: Turbomachinery Technical Conference and Exposition (GT2017), Charlotte, NC, USA, 26–30 June 2017. [Google Scholar]
Figure 1. The schematic diagram of ideal (0–1–2–3–0) and real (0–1’–2’–3’–0) Brayton cycles in a gas turbine. (a) The temperature–entropy diagram; (b) The schematic diagram of a gas turbine. The hot air at 3(3’) cools down in the atmosphere to become cold air at 0. Only the irreversible processes in the compressor (0–1’) and the turbine (2’–3’) are under consideration.
Figure 1. The schematic diagram of ideal (0–1–2–3–0) and real (0–1’–2’–3’–0) Brayton cycles in a gas turbine. (a) The temperature–entropy diagram; (b) The schematic diagram of a gas turbine. The hot air at 3(3’) cools down in the atmosphere to become cold air at 0. Only the irreversible processes in the compressor (0–1’) and the turbine (2’–3’) are under consideration.
Entropy 19 00470 g001
Figure 2. Influence of compressor cascade i on the exergy obtained in the combustion chamber. (a) A gas turbine with i − 1 compressor cascades: the inlet and outlet temperature of the combustion chamber is T 1 , i 1 and T 2 , i 1 ; (b) A gas turbine with i compressor cascades: the inlet and outlet temperature of the combustion chamber is T 1 , i and T 2 , i . The static temperature is increased from T 1 , i 1 to T 1 , i through compressor cascade i.
Figure 2. Influence of compressor cascade i on the exergy obtained in the combustion chamber. (a) A gas turbine with i − 1 compressor cascades: the inlet and outlet temperature of the combustion chamber is T 1 , i 1 and T 2 , i 1 ; (b) A gas turbine with i compressor cascades: the inlet and outlet temperature of the combustion chamber is T 1 , i and T 2 , i . The static temperature is increased from T 1 , i 1 to T 1 , i through compressor cascade i.
Entropy 19 00470 g002
Figure 3. The computational domain and the mesh resolution in a cross section.
Figure 3. The computational domain and the mesh resolution in a cross section.
Entropy 19 00470 g003
Figure 4. The static pressure coefficient (Cp) on the blade wall surface; comparison between the current computational fluid dynamics (CFD) results and the experimental data in [28]. (a) x 3 / C = 50 % ; (b) x 3 / C = 29.7 % .
Figure 4. The static pressure coefficient (Cp) on the blade wall surface; comparison between the current computational fluid dynamics (CFD) results and the experimental data in [28]. (a) x 3 / C = 50 % ; (b) x 3 / C = 29.7 % .
Entropy 19 00470 g004
Figure 5. The distribution of the entropy production rate density at different cross sections. The incidence angle α is 4°. (a) Equation (31), x 2 / l b l d = 50% (mid-plane); (b) Equation (35), x 2 / l b l d = 50% (mid-plane); (c) Equation (31), x 2 / l b l d = 25%; (d) Equation (35), x 2 / l b l d = 25%; (e) Equation (31), x 2 / l b l d = 5% (close to the hub); (f) Equation (35), x 2 / l b l d = 5% (close to the hub).
Figure 5. The distribution of the entropy production rate density at different cross sections. The incidence angle α is 4°. (a) Equation (31), x 2 / l b l d = 50% (mid-plane); (b) Equation (35), x 2 / l b l d = 50% (mid-plane); (c) Equation (31), x 2 / l b l d = 25%; (d) Equation (35), x 2 / l b l d = 25%; (e) Equation (31), x 2 / l b l d = 5% (close to the hub); (f) Equation (35), x 2 / l b l d = 5% (close to the hub).
Entropy 19 00470 g005
Figure 6. The distribution of the entropy production rate in a plane close to the hub ( x 2 / l b l d = 5%) for different incidence angles. (a) α = 2°, Equation (31); (b) α = 2°, Equation (35); (c) α = 4°, Equation (31); (d) α = 4°, Equation (35); (e) α = 7°, Equation (31); (f) α = 7°, Equation (35).
Figure 6. The distribution of the entropy production rate in a plane close to the hub ( x 2 / l b l d = 5%) for different incidence angles. (a) α = 2°, Equation (31); (b) α = 2°, Equation (35); (c) α = 4°, Equation (31); (d) α = 4°, Equation (35); (e) α = 7°, Equation (31); (f) α = 7°, Equation (35).
Entropy 19 00470 g006
Figure 7. Influence of the reference parameters on the exergy transformation number η c o m , i . (a) Q ˙ c b c / m ˙ varies from 3.8 × 10 5 J kg−1 to 4.2 × 10 5 J kg−1; (b) c p varies from 950 J kg−1 K−1 to 1050 J kg−1 K−1. The arrows indicate the direction of increasing Q ˙ cbc / m ˙ and c p .
Figure 7. Influence of the reference parameters on the exergy transformation number η c o m , i . (a) Q ˙ c b c / m ˙ varies from 3.8 × 10 5 J kg−1 to 4.2 × 10 5 J kg−1; (b) c p varies from 950 J kg−1 K−1 to 1050 J kg−1 K−1. The arrows indicate the direction of increasing Q ˙ cbc / m ˙ and c p .
Entropy 19 00470 g007
Table 1. Geometric parameters of the cascade.
Table 1. Geometric parameters of the cascade.
ParameterValue
Chord (m)0.15
Camber angle (°)23.22
Stagger angle (°)42.7
Pitch spacing (m)0.134
Solidity1.12
Blade span (m)0.37
Aspect ratio2.47
Design upstream flow angle (°)54.31
Design downstream flow angle (°)31.09
Table 2. Time- and volume-averaged dimensional Reynolds Averaged Navier–Stokes (RANS) results. The local entropy generation rate S ˙ D was calculated with Equation (35). S ˙ i , c o m / ( 2 N b l a ) is the entropy generation rate in a half-cascade passage, and S ˙ i , c o m is the entropy generation rate in a cascade.
Table 2. Time- and volume-averaged dimensional Reynolds Averaged Navier–Stokes (RANS) results. The local entropy generation rate S ˙ D was calculated with Equation (35). S ˙ i , c o m / ( 2 N b l a ) is the entropy generation rate in a half-cascade passage, and S ˙ i , c o m is the entropy generation rate in a cascade.
Bulk Variables α = 0 ° α = 2 ° α = 4 ° α = 7 °
S ˙ i , c o m / ( 2 N b l a ) (W K−1)
(Equation (31))
0.07050.07950.09760.108
S ˙ i , c o m / ( 2 N b l a ) (W K−1)
(Equation (35))
0.06950.07820.1010.112
m ˙ / ( 2 N b l a ) (kg s−1)0.68970.6560.6210.568
P i n l (Pa)101,009100,996101,004101,015
P i n l * (Pa)101,955101,941101,949101,960
T i n l (K)296.94296.90296.93296.96
T i n l * (K)299.72299.68299.71299.74
P o u t (Pa)101,325101,325101,325101,325
P o u t * (Pa)101,922101,896101,886101,888
T o u t (K)297.87297.87297.87298.87
T o u t * (K)299.53299.55299.52299.53
Table 3. The reference parameters for calculating the local exergy transformation number η c o m , i E . Q ˙ 1 / m ˙ and c p are the properties of the gas in the combustion chamber.
Table 3. The reference parameters for calculating the local exergy transformation number η c o m , i E . Q ˙ 1 / m ˙ and c p are the properties of the gas in the combustion chamber.
ParametersValues
T (K)288.15
Q ˙ c b c / m ˙ (J kg−1) 4 × 10 5
c p (J kg−1 K−1)1005
Table 4. Global coefficients according to the RANS results. The local entropy generation rate S ˙ D was calculated with Equation (35). The optimal exergy transformation number is shown with grey background.
Table 4. Global coefficients according to the RANS results. The local entropy generation rate S ˙ D was calculated with Equation (35). The optimal exergy transformation number is shown with grey background.
Dimensionless Coefficients α = 0 ° α = 2 ° α = 4 ° α = 7 °
ζ p (Equation (15))0.03490.04760.06670.0762
ζ s (Equation (16))0.01070.01270.01730.0210
η c o m , i (Equation (24))0.0004900.0005000.0004550.000411

Share and Cite

MDPI and ACS Style

Jin, Y.; Du, J.; Li, Z.; Zhang, H. Second-Law Analysis of Irreversible Losses in Gas Turbines. Entropy 2017, 19, 470. https://doi.org/10.3390/e19090470

AMA Style

Jin Y, Du J, Li Z, Zhang H. Second-Law Analysis of Irreversible Losses in Gas Turbines. Entropy. 2017; 19(9):470. https://doi.org/10.3390/e19090470

Chicago/Turabian Style

Jin, Yan, Juan Du, Zhiyuan Li, and Hongwu Zhang. 2017. "Second-Law Analysis of Irreversible Losses in Gas Turbines" Entropy 19, no. 9: 470. https://doi.org/10.3390/e19090470

APA Style

Jin, Y., Du, J., Li, Z., & Zhang, H. (2017). Second-Law Analysis of Irreversible Losses in Gas Turbines. Entropy, 19(9), 470. https://doi.org/10.3390/e19090470

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