International Journal of Aerospace Engineering
Volume 2020 (2020), Article ID 8825542, 18 pages
https://doi.org/10.1155/2020/8825542

Research Article

Numerical Investigation of the Wake Vortex-Related Flow Mechanisms in Transonic Turbines

Guoliang Wang, Ning Ge, and Dongdong Zhong

College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

Correspondence should be addressed to Ning Ge; gening@nuaa.edu.cn

Received 8 May 2020; Revised 3 July 2020; Accepted 11 July 2020; Published 1 August 2020

Academic Editor: Hao J. Xia

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

Abstract.

As the core equipment of the power generation system, a gas turbine is an indispensable energy-converting device in the national industry. The flow inside a high-pressure turbine (HPT) is highly unsteady, which has a great influence on the aerothermal performance and structural strength. To better clarify the flow mechanism and guide the advanced design, the basic flow characteristics of transonic turbines are investigated in the paper by a modified scale-adaptive simulation (SAS) model based on the shear stress transport (SST) turbulence model. The numerical results reveal the formation and development of the secondary flow structures such as wake vortex, pressure wave, shock wave, and the interactions among them. The length and frequency characteristics of wake are in good agreement with the large eddy simulation (LES) and the experimental data. Based on the detailed flow information, the local loss analysis is performed using the entropy generation rate. In summary, the wake vortex-related flow is the main origin of unsteadiness and entropy loss in high-pressure turbine cascade.

1. Introduction

Due to highly complex geometries and viscous effect, various secondary flows which cover a range of length scales and frequencies are formed in a gas turbine. The understanding of its flow mechanism is critical to the advanced design. The high cost of experimental equipment and the limitations of measurement technology make numerical research increasingly important.

Simulation technology is an important means to support the research and development of high-performance HPT. Over the past 30 years, significant progress has been made in computer technology and computational fluid dynamics (CFD) technology, especially in the full-scale three-dimensional CFD capability for turbomachinery. As a result, the design cycle of a gas turbine engine has been greatly reduced. In order to further improve the aerothermal performance and reduce the fuel consumption, the total pressure ratio and the inlet temperature of HPT are increasing, but its size becomes more compact. Therefore, in precise design of modern turbines, more attention should be paid to the flow details, rather than the average parameters obtained by Reynolds averaged N-S simulation. However, the existing theory for turbine aerodynamic design is mainly based on statistical modeling, ignoring the inherent unsteady flow characteristics. High-fidelity numerical simulation is an effective means to study the complex flow phenomena inside turbine components, and accurate simulation provides a solid foundation for accurate loss prediction. In recent years, the global hybrid models such as scale-adaptive simulation and delayed detached-eddy simulation have received considerable attentions, especially in engineering practices. Because they combine the advantages of RANS and LES methods, the pure LES and DNS are costly to be applied to high Reynolds number flows due to unacceptable CPU requirements [1].

The systematic research on vortex shedding started from the explanation of von Kármán, and the phenomenon is since associated to his name. The alternate vortex shedding is considered as responsible for a remarkable number of collapses of civil structures and damage to industrial equipment. Lienhard et al. analyzed a large number of experimental results of circular cylinders and summarized the law of vortex shedding at different Re [2]. An important feature of the flow past a circular cylinder is that the flow pattern depends on the Reynolds number. The knowledge of the mechanism of cylinder wake improved the understanding of turbine wake. Gehrer et al. [3] carried out experimental research by means of laser Doppler velocimetry (LDV) on the wake vortex of a linear turbine cascade and simulated it with different unsteady Reynolds averaged N-S simulation (URANS) models. As a result, the URANS method had insufficient resolution for unsteady wake vortices. Experiments of a linear turbine cascade by Sieverding et al. [4] revealed that the distances between the adjacent wake vortices are not equal, which is quite different from the flow past a circular cylinder. This is because the flow conditions are the same on both sides of a cylinder. By using the DDES method based on the S-A one-equation model [5], Lin et al. carried out a detailed numerical study on an HPT guide vane [6]. The periodic wake vortex shedding is an important origin of unsteadiness, but the unsteady effects related to the wake vortex have not been fully understood, because of the high Reynolds number, transonic Mach numbers, and high frequency, which impose great challenges for both experimental facilities and CFD solvers [6]. Their numerical results confirmed the findings of Sieverding et al., that is, the distance between an upstream pressure-side shedding vortex and a suction-side shedding vortex is always smaller than the distance between the latter and the downstream pressure-side shedding vortex. Egorov and Menter [7, 8] applied the SST-SAS model to a simplified blade geometry of MTU aeroengines. The strong mixing process of gas and cold jets was correctly predicted, and the film cooling effectiveness of the SAS was more consistent with the experiments, compared with URANS. Morata et al. [9] compared the performance of the LES approach in HPT using structured and unstructured grids, respectively. The LES results can be used as an assessment criterion for the global hybrid model.

A modified technique based on the Menter SST-SAS model was proposed in the previous work [10] and improved the spatial resolution of the original model in turbomachinery. As a continuation of the above work, in this paper, an in-house CFD solver integrated with the improved SAS model is employed to delineate the unsteady phenomena and local loss related to vortex shedding, and the underlying flow mechanisms of transonic HPT vanes are analyzed in different working conditions.

2. Numerical Methods

2.1. SST-SAS Model Description

In the course of the DESIDER project, Egorov and Menter [11] proposed the SST-SAS model and implemented in the commercial ANSYS CFX solver. The governing equations of the SAS model differ from those of the SST RANS model by an additional source term in the transport equation for the turbulence eddy frequency . The governing equations can be expressed as

The key to the SAS method is the introduction of the von Kármán length scale, , into the scale-determining equation of the RANS turbulence model. The second length scale depends on the local flow topology, which has different responses in steady flow regions and unsteady flow regions. The von Kármán scale is defined as the ratio of the first velocity derivative divided by the second velocity derivative, and the definition given by Menter is as follows:

where represents the velocity vector, and the average strain rate replaces the first velocity derivative .

But when the flow unsteadiness is “insufficient,” the scale-resolving capability of the original SST-SAS model may not be fully activated [12]. In order to solve the response problem to unsteadiness, a redefined von Kármán scale is proposed by theoretical analysis. According to the Helmholtz velocity decomposition theorem [13], in a three-dimensional vector field, can be resolved into a sum of an irrotational vector and a solenoidal vector. The average strain rate represents the tensile and shear deformation, and the average vorticity represents the rotation. Thus, the mathematical formulas of should contain both effects, and the Decaying Homogenous Isotropic Turbulence (DHIT) case [14] is used to quantitatively determine its two constants, and . The complete equation for in the present SAS can be written as

where is the control volume size.

The constant parameters used in the SST-SAS model are

2.2. Research Objects

2.2.1. Case 1: VKI LS89

The VKI LS89 vane from the von Kármán Institute for Fluid Dynamics is taken as the first research object [15]. A large number of experimental and numerical studies have been carried out on the transonic turbine guide vane, and its design parameters are given in Table 1. In the present study, the computational domain and mesh details at the leading and the trailing edges are shown in Figure 1. The total mesh cells are 7 million, which are divided into 76 blocks for parallel computation. The mesh has 60 evenly spaced points in the spanwise direction with a span length equal to 14.8% of the chord length, and the first layer near wall satisfies the simulation requirement of less than 1.0.

Table 1: Characteristic dimensions of research objects [15, 16].

Parameters
VKI LS89
CNRC

Chord length (mm)
67.647
204

Axial chord length (mm)
36.985
84.7

Pitch to chord ratio (-)
0.850
0.725

Throat to chord ratio (-)
0.2207
0.1515

Flow inlet angle (degree)
0
0

Stagger angle (degree)
55.0
65.5

Trailing edge diameter (mm)
1.42
6.35

Figure 1: Computational domain and mesh details for VKI LS89.

2.2.2. Case 2: CNRC Cascade

The second research object is from the Gas Turbine Laboratory of the National Research Council of Canada (CNRC) [16], and the parameters of the CNRC cascade are given in Table 1. Since there is no corresponding experimental schlieren in Case 1, the experimental schlieren images are used for quantitative verification. The CNRC cascade consisted of a set of 6 blades, whose profile represented the midspan section of a high-pressure turbine vane originally designed for a highly-loaded, single stage turbine with a stage pressure ratio of 3.8. The experimental tests were carried out in a large-scale, continuous flow, transonic linear cascade wind tunnel, and the goal of the CNRC project was to study the important flow phenomena that affect the performance of turbines, such as transition, wake and calmed region, shock/boundary layer interaction, and vortex shedding. The computational domain is shown in Figure 2. Because of the periodicity in the span, in order to save computing resources, a span length equal to 4.9% of its chord length is selected for simulations. The total mesh cell size is 1.45 million, where the mesh has 22 evenly spaced points in the spanwise direction. In the first layer near wall, the maximum is less than one.

Figure 2: Computational domain and mesh details for the CNRC cascade.

2.3. Computational Setup

For Case 1, based on the SST turbulence model, the RANS, the URANS, Menter SST-SAS, Menter SST-DDES, and the present SAS were used for simulations with the same structured mesh. Except for the RANS which adopted the WENO third-order scheme, because the difference between the third-order and the fifth-order scheme is very small for RANS, all other simulations used the WENO fifth-order scheme to reconstruct the inviscid flux at interfaces. The sixth-order central difference scheme was used to calculate the viscous flux. The total temperature, total pressure, and velocity direction were given at the inlet, with the mean static pressure defined at the outlet based on the flow conditions for the VKI experiments as in Table 2. The blade wall was set as adiabatic wall with a no-slip condition. Periodic boundary condition was used in the circumferential direction and the spanwise direction. The dual-time step method was used for unsteady calculations, and the time step size is 0.25 μs.

Table 2: Boundary conditions of research objects.

Flow conditions
(K)
(Pa)
(degree)
(Pa)

VKI LS89: MUR129
409.2
184900
0
116500
0.840

VKI LS89: MUR235
413.3
182800
0
104900
0.927

CNRC: Mach 1
293.0
101300
0
67000
0.80

CNRC: Mach 2
293.0
101300
0
54602
1.00

CNRC: Mach 3
293.0
101300
0
43600
1.16

For Case 2, the numerical method and the boundary condition were given in the same way as in the first case. The isentropic exit Mach number under different operating conditions was 0.8, 1.0, and 1.16, respectively, which was obtained by adjusting the back pressure.

3. Results and Discussions

3.1. Case 1: VKI LS89

3.1.1. Numerical Method Validation

Figure 3 shows the comparisons of isentropic Mach number distributions along the blade surface of the flow condition MUR129 in Table 2 obtained by the present SAS, the experimental data, and the time-averaged result of LES. The numerical results are basically consistent with the experimental data. In addition, the fluctuation of the velocity at certain downstream locations of the trailing edge is monitored for MUR129, represented by power spectral density (PSD) obtained using a Fourier transform. The PSD curves in Figure 4 follow the -5/3 law, indicating that the inertial region of the turbulence is correctly resolved by the present SAS. Moreover, the peak frequency is 37.5 kHz, which correspond to the vortex shedding frequency, consistent with the shedding frequency of the reference [6]. The grid quantity, the turbulence model, and the numerical scheme adopted in the current case are feasible.

Figure 3: Isentropic Mach number distributions of flow condition MUR129.

Figure 4: Power spectral density of the present SAS under flow condition MUR129.

3.1.2. An Overview of Flow Field Structure

The density gradient contours can clearly delineate the flow fields, and its definition is as follows:

Figure 5 shows the dimensionless density gradient for the flow condition MUR235, reference to Table 2, obtained with the RANS, the URANS, Menter SAS, Menter DDES, the present SAS, and the time averaging of the present SAS in sequence, followed by the LES result cited from the reference. The RANS model cannot predict the existence of the pressure wave and the vortex shedding at all. The trailing edge pressure wave calculated by the URANS is weak, and the shape of wake vortex is quite different from the actual flow structure. There is a lack of resolution in the wake vortex for Menter SAS, indicating that the ability of scale resolving is not fully excited. As shown in Figures 5(d), 5(e), and 5(g), the result of the present SAS is as capable as the large eddy simulation, and the DDES model also performs well: (1) there is a weak shock wave ① at approximately 84% axial chord length of the suction surface; (2) the wake vortices are predicted in pairs; (3) with the shedding of the wake, the pressure waves are induced, and the pressure wave ④ formed on the pressure side reflects after hitting the suction surface of the adjacent vane; (4) there are interactions among the shock wave, the pressure waves, and the wake. In particular, the suction-side pressure wave ③ dissipate the wake flow, and the core of wake vortices basically disappears, which qualitatively shows the broken process of wake vortex during transport.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 5: Dimensionless density gradient contours of flow condition MUR235: (a) SST-RANS result; (b) SST-URANS result; (c) Menter SST-SAS result; (d) Menter SST-DDES result; (e) present SAS result; (f) time-averaged result of present SAS; (g) LES result [9].

Therefore, the following sections use the results of the modified SAS model for flow mechanism analysis.

The contours in Figure 6 indicate the similar flow details for the flow condition MUR129. The difference from MUR235 is that there is no shock wave in this flow condition. As the exit Mach number decreases, the dissipation of pressure waves on the wake becomes weaker, so the core of wake vortices can maintain its stability and be transported further downstream.

(a)

(b)

Figure 6: Dimensionless density gradient contours of flow condition MUR129: (a) present SAS result; (b) LES result [9].

3.1.3. Formation and Development of Wake Vortex

In order to further analyze the characteristic length scale, the coherent structure of wake vortex is displayed by Q-criterion, and local Mach number is used for coloring. Figures 7(a) and 7(b) show the scale characteristics during vortex transportation. The vortices shedding from the pressure surface and the suction surface are denoted as VPS (pressure-side shedding vortex) and VSS (suction-side shedding vortex), respectively, and are numbered in chronological order. The vortex core is defined as the position with the lowest pressure. In each vortex pair, the intensity of VPS is stronger than that of VSS, and the streamwise distance between VPS₁ and VSS₁, , is smaller than the streamwise distance between VSS₁ and VPS₂, , which are quite different from the vortices of a circular cylinder. So when studying the effect of upstream wake vortices on the current blade row, the feasibility of using circular cylinders to generate upstream wakes is worth reexamination.

(a)

(b)

Figure 7: The length characteristics of wake vortex: (a) flow condition MUR235; (b) flow condition MUR129.

To explain the asymmetry phenomena, it is necessary to start from the generation process of wake vortex. The coherent structures displayed by Q-criterion at different moments are shown in Figure 8, and the pressure coefficient is used for coloring. First, define the period of a vortex pair as , and define the moment when the suction-side near-wall vortex starts to appear as . Then, the vortex continuously absorbs the low-momentum fluid near the trailing edge and grows over time as shown in Figure 8(a). At the same time, due to the blocking effect of VSS₁ on the main flow, a high-pressure zone appears near the trailing edge. The wake is a velocity deficit, especially just downstream of the trailing edge, and the downstream VSS₁ enters the range of the main flow on the suction side. But the main flow maintains a high speed due to inertia, thus forming a high-pressure region similar to the stagnation flow. Therefore, due to the pressure difference and the blade geometry, the near-wall vortex is divided into two parts, denoted as and , respectively. Soon at 0.53 in Figure 8(b), the vortex is divided into two parts by the same pressure difference, denoted as and , respectively. Then, the downstream VPS₁ gradually enters the main flow on the pressure side in Figure 8(c). As expected, these two pressure differences have the same squeezing effect on and , and at about 0.68 , they are divided separately, namely, , , , and . Figure 8(d) indicates that the becomes a part of VPS₁ and moves downstream following the previous vortex pair; in addition, the is related to the formation of the suction-side pressure wave. Eventually, the and vortices merge to become the main part of the new VSS, and the vortex is gradually dissipated. It is worth noting that VPS does not undergo such a complicated splitting process, so its strength is stronger than VSS. For the same reason, the asymmetric flow mechanism also causes to be smaller than . The pressure coefficient is calculated by

(a)

(b)

(c)

(d)

Figure 8: The formation of wake vortex: (a) instantaneous flow field at ; (b) ; (c) ; (d) .

As the vortices transported to the downstream, the distance gradually decreases. On the one hand, this is because the reverse pressure gradient received by VSS is larger than VPS as shown by the black dotted line in Figure 9, so that the motion of VSS is hindered and decelerated. On the other hand, the double-vortex structure formed in the vortex pair makes the two vortices closer to each other as shown in Figure 10. Ignoring the effect of the main flow, the spatial distance between the two vortices in the same vortex pair is small, so there is an interaction between them. There is an induced rotating flow field around each concentrated vortex, which means that VPS will produce an induced velocity for VSS, namely, in Figure 10, and vice versa. But the intensity of VPS is stronger, that is, the induced velocity is also greater than . Therefore, VSS moves upstream relative to VPS, and the distance between the two vortices in the main flow direction is further reduced.

Figure 9: Pressure coefficient distribution of MUR129.

Figure 10: Double-vortex structure in a vortex pair.

3.1.4. Formation and Development of Pressure Wave

Figure 11 shows the formation and development of suction-side pressure waves which are generated by the intense vortex shedding. More specifically, it should be that the vortex induces the generation of the suction-side pressure wave. As mentioned above, the vortex is divided into two parts at 0.68 , namely, and , and the three-dimensional view at this instant is displayed in Figure 11(a). As the vortex forms, it gradually moves away from the trailing edge and enters the main flow in a direction almost perpendicular to the main flow. The blockage effect on the main flow will cause pressure fluctuations to propagate in all directions, since the freestream is a subsonic flow, but the waves travel at the local speed of sound, which means that the Mach number relative to the motion of pressure waves is 1.0. It is worth noting that the main flow has absolute velocity and direction, and its Mach number near the trailing edge is around 0.9, so the pressure waves observed from Figures 11(b) and 11(c) have different moving speeds in different directions. These pressure waves move upstream over the suction surface at about 11% of the freestream velocity, affecting the aft suction surface. In addition, since the suction-side pressure wave is closely related to the vortex , its frequency is exactly equal to the frequency of the VSS. In the same way, it can be concluded that the formation of pressure-side pressure waves is closely related to VPS.

(a)

(b)

(c)

Figure 11: The formation and development of suction-side pressure wave (density gradient contours and Q-criterion isosurfaces colored by ): (a) ; (b) ; (c) .

3.1.5. Unsteady Interactions Related to Wake Vortex

Combined with the instantaneous three-dimensional perspective of Figure 12, the evolution of wake vortex is explored. The shedding vortex just leaving the trailing edge exhibits periodicity in the spanwise direction, and the shape of the vortex core is similar to that of a cylinder with two-dimensional features. The tube around the vortex core exhibits three-dimensional features but maintains periodicity in the span. When the wake vortex moves to the area affected by the shock wave, the spatial size of the vortex pair increases, and the presence of the shock causes the vortex tube to appear three-dimensional rolled up along the span, marked with the white dotted line. As a comparison in Figure 12(b), there is almost no radial mixing process in the same region under the flow condition MUR129. Then, the vortex is transported to the influence area of suction-side pressure waves. As the pressure wave passes, accompanied by continuous violent oscillations, the large-scale tube is strongly dissipated and broken into multiscale structures, regardless of MUR129 or MUR235. It can therefore be concluded that the interactions of the shock wave, pressure wave, and wake vortex accelerate the instability and fragmentation of the vortex.

(a)

(b)

Figure 12: Three-dimensional development of wake vortex: (a) flow condition MUR235; (b) flow condition MUR129.

3.1.6. Time-Averaged Flow Analysis

In order to provide a design reference for the HPT vane in the existing RANS design system, it is necessary to quantify the impact of unsteady flow on the aerodynamic performance. First, the time-averaged flow field of the present SAS under the flow condition MUR129 is compared with the RANS result as shown in Figure 13, and the centerline direction of the trailing edge is marked with a red dot-dash line. The comparison of the flow condition MUR235 is shown in Figures 5(a) and 5(f). It can be concluded that the main difference between the unsteady time-averaged flow and the RANS flow is in the wake region. On the one hand, the width of the unsteady time-averaged wake is wider, indicating that the mixing with the main flow is stronger, and on the other hand, the unsteady time-averaged wake has obvious asymmetry about the centerline.

(a)

(b)

Figure 13: Comparison of density gradient contours of MUR129: (a) SST-RANS result; (b) time-averaged result of present SAS.

The trailing edge and the free shear layers on both sides of the trailing edge form a triangular area, which is called the base region, and the pressure in the base region is closely related to the wake vortices; increasing the base pressure coefficient can effectively reduce the trailing edge loss [17]. Therefore, in the design process of turbine blades, correct estimation of the base pressure coefficient is essential for predicting profile loss. Figure 14 shows the base pressure coefficient distribution obtained by the RANS and the time averaging, where the horizontal axis represents the arc length from the blade surface to the center point of the trailing edge (dimensionless by the trailing edge diameter). There are three main differences. Firstly, under different flow conditions, the RANS results are always higher than the time-averaged results, which can underestimate the importance of wake loss. Secondly, the RANS have a base pressure plateau, but the time-averaged results show significant nonuniformity. Thirdly, the RANS results are symmetrical about the trailing edge centerline, whereas the time-averaged results are not. The profile loss can be quantified based on the stagnation pressure loss coefficient. For the flow condition MUR235, the loss coefficients of the RANS and the present SAS are 5.86% and 6.32%, respectively, and the relative error is 7.85%. For another flow condition MUR129, the RANS and the unsteady time-averaged loss coefficients are 4.50% and 4.89%, respectively, and their relative error is 8.67%. Combined with the analysis in Section 3.1.7, it can be concluded that the error mainly comes from the shedding and development of wake vortex, because the RANS system overestimates the base pressure coefficient and ignores the nonuniformity and asymmetry of the base pressure. The stagnation pressure loss coefficient is calculated by

Figure 14: Comparison of base pressure coefficient distribution.

3.1.7. Local Loss Analysis

There are many definitions of loss coefficient in regular use for blade rows, such as stagnation pressure loss coefficient, enthalpy loss coefficient, and entropy loss coefficient [17]. But those coefficients are used for global loss analysis, which represent the total losses of the whole region. The classical breakdown of loss into profile loss, end-wall loss, and leakage loss continues to be widely used in advanced design, so it is necessary to reflect the information about local loss generation in high-fidelity simulation. Essentially, losses are the results from viscous irreversibility and thermal irreversibility, and the entropy generation rate can quantify irreversibility and reveal the specific location of local loss [18, 19]. The definition of entropy generation rate is

Figures 15(a) and 15(b) are the local entropy generation rates corresponding to viscous irreversibility and thermal irreversibility of the flow condition MUR235, respectively. Similarly, Figures 16(a) and 16(b) show the viscous entropy generation rate and the thermal entropy generation rate of MUR129, respectively. For linear cascades, the only loss is profile loss, and its components come from wake vortex, boundary layer, pressure wave, shock wave, etc. In particular, the periodic wake vortex is the most important origin of profile loss, and the pressure waves induced by wake vortex strongly dissipate the downstream vortices. In addition, these graphs visually compare viscous loss and thermal loss, the former is significantly larger than the latter. In order to quantitatively analyze these two losses, the Bejan number (Be) is introduced here, which represents the percentage of thermal loss in total local losses [19]. The distribution of Be in Figures 17 and 18 illustrates that in areas with high local losses, such as wake vortex and boundary layer, the Be is between 0.1% and 1%, while in other low local-loss areas, the Be is about 90%. It can be summarized that the viscous effect plays a leading role in the loss of the VKI LS89 cascade.

(a)

(b)

Figure 15: Comparison of local losses due to viscous and heat-transfer effects (MUR235): (a) viscous entropy generation rate; (b) thermal entropy generation rate.

(a)

(b)

Figure 16: Comparison of local losses due to viscous and heat-transfer effects (MUR129): (a) viscous entropy generation rate; (b) thermal entropy generation rate.

Figure 17: Local Bejan number distribution of MUR235.

Figure 18: Local Bejan number distribution of MUR129.

3.1.8. Energy Separation in Wakes

In addition to the entropy loss, vortex shedding may cause adverse effects which include high-frequency sound propagation, vibration, and local high heat transfer and energy separation into hot and cold regions [20]. The stagnation temperature on the wake centerline is found to be 30 K lower than that of the incoming fluid, while the stagnation temperature at the edges of the wake can be up to 30 K higher than incoming fluid, as shown in Figure 19(a), and Figure 19(b) is the stagnation temperature of the profile. It was demonstrated that the temperature separation is a thermoacoustic effect where the wake centerline emerges colder than the surrounding fluid, which leads to the redistribution of the stagnation temperature and the generation of local hot spots. The nonuniform total temperature distribution is a source of entropy production, and more importantly, it should be carefully considered in the turbine cooling design.

(a)

(b)

Figure 19: Time-resolved total temperature distribution in wakes: (a) overview; (b) profile.

3.2. Case 2: CNRC Cascade

3.2.1. Comparison and Analysis of Experimental Schlieren

The experimental schlieren photographs from [16, 20] are selected to quantitatively analyze the relative error of simulation results and experimental data in different flow conditions. The density gradient contours are displayed in Figures 20–22.

Figure 20: Density gradient contours of : (a) present SAS; (b) experimental schlieren [20].

(a)

(b)

Figure 21: Comparison of density gradient contours of : (a) present SAS; (b) experiment [16, 20].

(a)

(b)

Figure 22: Comparison of density gradient contours of : (a) present SAS; (b) experiment [20].

For the supersonic condition in Figure 20, the exit Mach number is equal to 1.16, the downstream disturbance cannot propagate upstream, and there is no pressure wave in the turbine nozzle passage. Instead, there are two strong oblique shock waves generated at the trailing edge, and the passage shock wave impinges on the suction surface of the adjacent blade. The incident point of passage shock is well predicted, and the relative error of its length and angle characteristics is within 2%. After experiencing the shock/boundary-layer interaction, a reflected shock wave is formed. Under the supersonic flow condition, the origin of the vortex shedding moved from the vane trailing edge to the confluence of the two trailing edge shear layers, and the predicted vortex shedding frequency is about 25.9 kHz, which is consistent with the experimental frequency in the literature [16].

For in reference [16] (or in reference [20]), the experimental flow structure is shown in Figure 21(b). There is a normal shock wave at the throat, and supersonic flow decelerates to subsonic flow through this shock wave. The pressure waves and their reflected waves are located downstream of the normal shock wave. The present SAS model reproduces the flow details well, such as the distribution of wake vortex and wave system. The distances between different vortices are marked in the figures, that is, , , and in Figure 21(a) and , , and in Figure 21(b), and the maximum error between the distances is about 2%. For in Figure 22, the flow structure similar to the flow condition MUR129 of LS89 can be observed. The wake vortices on the SS and PS appear in pairs and shed off alternately. With the shedding of the wake vortices, the pressure waves on both sides are induced and they propagate upstream, and the pressure-side pressure waves reflect on the adjacent vane.

The power spectral density curves in Figure 23 reveal that the turbulent inertial regions are correctly resolved. All the above results illustrate that the numerical methods can reproduce the flow field of real experimental flow conditions and it is reliable for analyzing the flow mechanisms inside transonic HPT.

Figure 23: Turbulent velocity power spectral density of the CNRC cascade.

3.2.2. Local Loss Analysis

In this section, the local entropy generation under supersonic condition is analyzed. Figures 24(a) and 24(b) are the entropy generation rates due to viscous effect and thermal effect, respectively, and Figure 24(c) shows the corresponding local Bejan number distribution. It can be found that although the is greater than 1, the conclusions are consistent with that obtained from the subsonic condition. First, the entropy generation rates are on the same order of magnitude as in Case 1, and then, the Be is between 0.1% and 1% in wakes and boundary layer, while Be is close to 90% in other areas.

(a)

(b)

(c)

Figure 24: Comparison of local losses due to viscous and heat-transfer effects (): (a) viscous entropy generation rate; (b) thermal entropy generation rate; (c) local Bejan number distribution.

4. Conclusions

A modified SAS model based on the Menter SST-SAS model was applied to the numerical simulation and flow mechanism analysis of the transonic turbine. The conclusions are as follows:

(1)The wake vortex-related flow is the main origin of unsteadiness and losses in transonic HPT cascade, and local entropy generation analysis shows that the viscous loss is dominant(2)Comparing the RANS results with the time-averaged results of the present SAS, it can be found that the RANS overestimates the base pressure coefficient and cannot predict the nonuniformity and asymmetry of the base pressure, which leads to errors in the prediction of profile loss(3)The formation of suction-side shedding vortex has gone through a complex process of division and reunion, whereas the pressure-side wake vortex has not, due to the asymmetric geometries and pressure conditions. The vortices VPS and induced the generation of pressure-side and suction-side pressure waves, respectively(4)There is energy separation in the wake vortex, and the stagnation temperature difference between the hot region and the cold region is as high as 15%(5)The high-fidelity simulation based on the present SAS model can reveal the length scale and time scale characteristics of the high-pressure turbine. The CPU hours of the present SAS is about 2.5% of the structured LES in the reference

Future works include continued flow mechanisms in other turbomachinery applications, in particular for multistage turbomachinery design and aerothermal performance analysis.

Abbreviations

SAS: | Scale-adaptive simulation
RANS: | Reynolds averaged N-S simulation
DDES: | Delayed detached-eddy simulation
LES: | Large eddy simulation
DNS: | Direct numerical simulation
PSD: | Power spectral density
VPS: | Pressure-side shedding vortex
VSS: | Suction-side shedding vortex.

Nomenclature

: | Chord length
: | Trailing edge diameter
Re: | Reynolds number
Ma: | Mach number
: | Pressure
: | Temperature
: | Pressure coefficient
: | Stagnation pressure loss coefficient
: | von Kármán length scale
: | Velocity
: | The first velocity derivative
: | The second velocity derivative
: | Period of a vortex pair
: | Grid cell size
: | Entropy generation rate
Be: | Bejan number
: | Induced velocity.

Greek Letters

: | von Kármán constant
: | Density
: | Velocity direction
: | Viscosity.

Subscripts

is: | Isentropic
ax: | Axial
ex: | Blade exit
inl: | Inlet
out: | Outlet
CV: | Control volume
visc: | Viscous
therm: | Thermal.

Superscripts

:
Total.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

Dongdong Zhong developed an in-house CFD parallel solver NUAA-Turbo 2.0 under the instruction of Ning Ge. Guoliang Wang proposed a modified SAS model based on the Menter SST-SAS model and applied it to the simulations and mechanism analysis of turbomachinery. Guoliang Wang wrote the paper, and Ning Ge provided insightful opinions. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This paper is supported by the AECC Shenyang Engine Research Institute of China. The research is partly supported by the National Key Project GJXM92579.

References

P. R. Spalart, “Strategies for turbulence modelling and simulations,” International Journal of Heat and Fluid Flow, vol. 21, pp. 252–263, 2010.

J. H. Lienhard, Synopsis of Lift, Drag, and Vortex Frequency Data for Rigid Circular Cylinders, Technical Extension Service: Washington State University, Bulletin 300, USA, 1966.

A. Gehrer, H. Lang, N. Mayrhofer, and J. Woisetschläger, “Numerical and experimental investigation of trailing edge vortex shedding downstream of a linear turbine cascade,” in Proceedings of the ASME Turbo Expo 2000: Power for Land, Sea, and Air. Volume 1: Aircraft Engine; Marine; Turbomachinery; Microturbines and Small Turbomachinery, Munich, Germany, 2000.

C. H. Sieverding, D. Ottolia, C. Bagnera, A. Comadoro, J. Brouckaert, and J. Desse, “Unsteady turbine blade wake characteristics,” Journal of Turbomachinery, vol. 126, no. 4, pp. 551–559, 2004.

P. R. Spalart, S. Deck, M. L. Shur, and K. D. Squires, “A new version of detached-eddy simulation, resistant to ambiguous grid densities,” Theoretical and Computational Fluid Dynamics, vol. 20, no. 3, pp. 181–195, 2006.

D. Lin, X. Yuan, and X. Su, “Local entropy generation in compressible flow through a high pressure turbine with delayed detached eddy simulation,” Entropy, vol. 19, no. 1, p. 29, 2017.

F. R. Menter and Y. Egorov, “The scale-adaptive simulation method for unsteady turbulent flow predictions. Part 1: theory and model description,” Flow Turbulence Combust., vol. 85, no. 1, pp. 113–138, 2010.

Y. Egorov, F. R. Menter, R. Lechner, and D. Cokljat, “The scale-adaptive simulation method for unsteady turbulent flow predictions. Part 2: application to complex flows,” Flow, Turbulence and Combustion, vol. 85, no. 1, pp. 139–165, 2010.

E. C. Morata, N. Gourdain, F. Duchaine, and L. Y. M. Gicquel, “Effects of free-stream turbulence on high pressure turbine blade heat transfer predicted by structured and unstructured LES,” International Journal of Heat and Mass Transfer, vol. 55, no. 21-22, pp. 5754–5768, 2012.

G. Wang, D. Zhong, N. Ge, and R. Yang, “Scale-adaptive simulations of high-pressure turbine guide vane,” in Proceedings of the ASME 2019 Gas Turbine India Conference. Volume 1: Compressors, Fans, and Pumps; Turbines; Heat Transfer; Structures and Dynamics, Chennai, Tamil Nadu, India, 2019.

Y. Egorov and F. R. Menter, “Development and application of SST-SAS turbulence model in the DESIDER project,” in Advances in Hybrid RANS-LES Modelling, S. H. Peng and W. Haase, Eds., vol. 97 of Notes on Numerical Fluid Mechanics and Multidisciplinary Design, pp. 261–270, Springer, Berlin, Heidelberg, 2008.

L. Du and F. Ning, “Scale adaptive simulation of flows past an airfoil after stall,” in Proceedings of the ASME 2012 Fluids Engineering Division Summer Meeting collocated with the ASME 2012 Heat Transfer Summer Conference and the ASME 2012 10th International Conference on Nanochannels, Microchannels, and Minichannels. Volume 1: Symposia, Parts A and B, pp. 1551–1560, Rio Grande, Puerto Rico, USA, 2012.

C. Q. Liu, Y. Q. Wang, Y. Yang, and Z. W. Duan, “New omega vortex identification method,” Science China Physics, Mechanics & Astronomy, vol. 59, no. 8, pp. 684–711, 2016.

G. Comte-Bellot and S. Corrsin, “Simple Eulerian time correlation of full-and narrow-band velocity signals in grid-generated, ‘isotropic’ turbulence,” Journal of Fluid Mechanics, vol. 48, no. 2, pp. 273–337, 1971.

T. Arts, R. M. Lambert, and A. W. Rutherford, “Aero-thermal investigation of a highly loaded transonic linear turbine guide vane cascade,” Tech. Rep., von Kármán Institute for Fluid Dynamics: Tech. Rep. VKI-174, Belgium, 1990.

T. C. Currie and W. E. Carscallen, “Simulation of trailing edge vortex shedding in a transonic turbine cascade,” Journal of Turbomachinery, vol. 120, no. 1, pp. 10–19, 1998.

J. D. Denton, “The 1993 IGTI scholar lecture: loss mechanisms in turbomachines,” Journal of Turbomachinery, vol. 115, no. 4, pp. 621–656, 1993.

H. Herwig and B. Schmandt, “How to determine losses in a flow field: a paradigm shift towards the second law analysis,” Entropy, vol. 16, no. 6, pp. 2959–2989, 2014.

A. Bejan, Entropy Generation through Heat and Fluid Flow, Wiley, New York, NY, USA, 1982.

J. P. Gostelow, W. E. Carscallen, M. Kurosaka, and A. Mahallati, “The relationship between energy separation and base drag in turbine blade wakes,” in Proceedings of the ASME Turbo Expo 2013: Turbine Technical Conference and Exposition. Volume 6A: Turbomachinery, San Antonio, Texas, USA, 2013.