Comparison of DDES and URANS for Unsteady Tip Leakage Flow in an Axial Compressor Rotor

2019 ◽  
Vol 141 (12) ◽  
Author(s):  
Yangwei Liu ◽  
Luyang Zhong ◽  
Lipeng Lu

Tip leakage vortex (TLV) has a large impact on compressor performance and should be accurately predicted by computational fluid dynamics (CFD) methods. New approaches of turbulence modeling, such as delayed detached eddy simulation (DDES), have been proposed, the computational resources of which can be reduced much more than for large eddy simulation (LES). In this paper, the numerical simulations of the rotor in a low-speed large-scale axial compressor based on DDES and unsteady Reynolds-averaged Navier–Stokes (URANS) are performed, thus improving our understanding of the TLV dynamic mechanisms and discrepancy of these two methods. We compared the influence of different time steps in the URANS simulation. The widely used large time-step makes the unsteadiness extremely weak. The small time-step shows a better result close to DDES. The time-step scale is related to the URANS unsteadiness and should be carefully selected. In the time-averaged flow, the TLV in DDES dissipates faster, which has a more similar structure to the experiment. Then, the time-averaged and instantaneous results are compared to divide the TLV into three parts. URANS cannot give the loss of stability and evolution details of TLV. The fluctuation velocity spectra show that the amplitude of high frequencies becomes obvious downstream from the TLV, where it becomes unstable. Last, the anisotropy of the Reynolds stress of these two methods is analyzed through the Lumley triangle to see the distinction between the methods and obtain the Reynolds stress. The results indicate that the TLV latter part in DDES is anisotropic, while in URANS it is isotropic.

Author(s):  
Yanfei Gao ◽  
Yangwei Liu ◽  
Luyang Zhong ◽  
Jiexuan Hou ◽  
Lipeng Lu

AbstractThe standard k-ε model (SKE) and the Reynolds stress model (RSM) are employed to predict the tip leakage flow (TLF) in a low-speed large-scale axial compressor rotor. Then, a new research method is adopted to “freeze” the turbulent kinetic energy and dissipation rate of the flow field derived from the RSM, and obtain the turbulent viscosity using the Boussinesq hypothesis. The Reynolds stresses and mean flow field computed on the basis of the frozen viscosity are compared with the results of the SKE and the RSM. The flow field in the tip region based on the frozen viscosity is more similar to the results of the RSM than those of the SKE, although certain differences can be observed. This finding indicates that the non-equilibrium turbulence transport nature plays an important role in predicting the TLF, as well as the turbulence anisotropy.


Entropy ◽  
2018 ◽  
Vol 21 (1) ◽  
pp. 21 ◽  
Author(s):  
Hui Li ◽  
Xinrong Su ◽  
Xin Yuan

In unshrouded turbine rotors, the tip leakage vortices develop and interact with the passage vortices. Such complex leakage flow causes the major loss in the turbine stage. Due to the complex turbulence characteristics of the tip leakage flow, the widely used Reynolds Averaged Navier–Stokes (RANS) approach may fail to accurately predict the multi-scale turbulent flow and the related loss. In order to effectively improve the turbine efficiency, more insights into the loss mechanism are required. In this work, a Delayed Detached Eddy Simulation (DDES) study is conducted to simulate the flow inside a high pressure turbine blade, with emphasis on the tip region. DDES results are in good agreement with the experiment, and the comparison with RANS results verifies the advantages of DDES in resolving detailed flow structures of leakage flow, and also in capturing the complex turbulence characteristics. The snapshot Proper Orthogonal Decomposition (POD) method is used to extract the dominant flow features. The flow structures and the distribution of turbulent kinetic energy reveal the development of leakage flow and its interaction with the secondary flow. Meanwhile, it is found that the separation bubble (SB) is formed in tip clearance. The strong interactions between tip leakage vortex (TLV) and the up passage vortex (UPV) are the main source of unsteady effects which significantly enhance the turbulence intensity. Based on the DDES results, loss analysis of tip leakage flow is conducted based on entropy generation rates. It is found that the viscous dissipation loss is much stronger than heat transfer loss. The largest local loss occurs in the tip clearance, and the interaction between the leakage vortex and up passage vortex promotes the loss generation. The tip leakage flow vortex weakens the strength of up passage vortex, and loss of up passage flow is reduced. Comparing steady and unsteady effects to flow field, we found that unsteady effects of tip leakage flow have a large influence on flow loss distribution which cannot be ignored. To sum up, the current DDES study about the tip leakage flow provides helpful information about the loss generation mechanism and may guide the design of low-loss blade tip.


Author(s):  
Chunwei Gu ◽  
Fan Feng ◽  
Xuesong Li ◽  
Meilan Chen

An attempt is made in the present paper to apply DES (Detached Eddy Simulation), which is based on S-A model of RANS, for investigating the flow field around a subsonic compressor rotor with a tip clearance of 2% blade height. Comparison of the results by DES and S-A model shows that DES model can capture more intensive vortex flow, such as tip leakage flow, double leakage flow, as well as interaction between the leakage flow and wake flow downstream of the rotor passage. DES model predicts more complicated flow at the separation region near the hub. DES simulation for different operation conditions also reveals interesting details. The shedding angle and strength of the tip leakage flow changes with the blade loading. The starting point of the leakage vortex moves towards the leading edge when the blade loading increases. Double leakage is observed only at the design and higher loading conditions, and is not at a lower loading condition. The tip leakage vortex splits into two branches downstream of the rotor blade due to interaction with the wake flow. Instantaneous results show unsteadiness of the tip leakage vortex. Alternating regions of higher and lower loss is found along the time-averaged leakage vortex trajectory. Obvious is also the unsteadiness in the separation region near the hub.


2021 ◽  
pp. 1-400
Author(s):  
Xiao He ◽  
Fanzhou Zhao ◽  
Mehdi Vahdati

Abstract Detached Eddy Simulation (DES) and its variants are emerging tools for turbomachinery simulations. In this paper, the state-of-the-art upgrades of DES are reviewed, and their capabilities in predicting compressor tip leakage flow are discussed. The upgrade with the best potential is identified as the Delayed DES (DDES) method with the grid spacing FKHΔhyb, which unlocks the physics of the Kelvin-Helmholtz instability in compressor tip leakage flow. The upgraded grid spacing FKHΔhyb is compared against the widely used default one Δmax in a backward-facing step and a low-speed axial compressor rotor. Results show that the DDES method with FKHΔhyb predicts both the main flow field and the turbulence field with reasonably good accuracy. However, the original DDES method with Δmax predicts a delayed transition to turbulence, which leads to an inaccurate prediction of the main flow field when using a coarse mesh. The findings in this paper highlight the future opportunities for using the DDES-FKHΔhyb method to predict tip-driven compressor stall and generate a turbulence database for turbulence model development.


Author(s):  
Yueqing Zhuang ◽  
Hui Liu

Since the unsteadiness of tip leakage flow has profound effects on both aerodynamic performance and stall margin of axial compressors, it is important to accurately predict the transient tip flow at affordable computational cost. Limited by the high requirement of grid resolution of wall turbulence flow, large eddy simulation (LES) method is greatly restricted in engineering application. In the present work, a Reynolds-stress-constrained large eddy simulation (CLES) method has been introduced, in which the whole domain is simulated using LES while Reynolds stress constraint is enforced on the subgrid-scale (SGS) stress model for near-wall regions aiming at reducing the near-wall grid resolution. The CLES simulations have been performed to investigate the flow behaviors of the unsteady tip leakage flow in a transonic compressor NASA Rotor 67 at near-stall conditions. Reliability assessments have been conducted through comparisons of experimental measurements and numerical results obtained by RANS, DES, CLES as well as LES, respectively. Both the total pressure ratio and isentropic efficiency calculated by CLES agree well with experiment. The turbulence statistical results show three distinct high flow fluctuation regions near the blade tip. The first one is a long and narrow strip ahead of the leading edge of the rotor caused by the movement of the passage shock wave. The second one is formed on the suction side from the leading-edge of the rotor blade due to the oscillation of the tip leakage vortex. And the third one, which occupies most of the blade passage from the middle part of the rotor blade, is generated under multiple factors. The frequency characteristic of the unsteady tip leakage flow has been analyzed. The energy spectrums of the local transient pressure signals are highly related with the local unsteady flow features. The originating mechanisms of the flow unsteadiness in the rotor tip leakage flow have also been discussed, and the results show that the flow unsteadiness is mainly caused by a combined interaction effect of the double leakage flow, the tip leakage vortex flow spilled from the adjacent blade passage, as well as the involved main flow.


Energies ◽  
2019 ◽  
Vol 12 (5) ◽  
pp. 954 ◽  
Author(s):  
Xiangyu Su ◽  
Xiaodong Ren ◽  
Xuesong Li ◽  
Chunwei Gu

Tip leakage vortex (TLV) in a transonic compressor rotor was investigated numerically using detached-eddy simulation (DES) method at different working conditions. Strong unsteadiness was found at the tip region, causing a considerable fluctuation in total pressure distribution and flow angle distribution above 80% span. The unsteadiness at near choke point and peak efficiency point is not obvious. DES method can resolve more detailed flow patterns than RANS (Reynolds-averaged Navier–Stokes) results, and detailed structures of the tip leakage flow were captured. A spiral-type breakdown structure of the TLV was successfully observed at the near stall point when the TLV passed through the bow shock. The breakdown of TLV contributed to the unsteadiness and the blockage effect at the tip region.


Author(s):  
Edward A. Miller ◽  
Michael J. Cave ◽  
David M. Williams ◽  
Khandan Thayalakhandan

Abstract Computational fluid dynamics (CFD) of industrial-scale, axial compressor geometries has traditionally been performed using steady state methods such as the mixing plane approach. With the surge in the development of large-scale, massively-parallel computing platforms, fully 3D unsteady approaches are rapidly growing in popularity. The fully 3D, unsteady approach involves building a full 3D domain for each blade row, and then coupling the stationary and rotating domains using a sliding interface. In the literature, there are various methods for solving this 3D unsteady problem, such as the Unsteady Reynolds Averaged Navier-Stokes (URANS) and the Detached Eddy Simulation (DES) methods. While these methods are well documented for a variety of real-world problems, there have been limited efforts to compare the effectiveness of these methods for fully 3D, unsteady turbomachinery problems. In this study, the first stage of an industrial-scale axial compressor was simulated using: i) the URANS approach, and ii) the DES approach. The compressor geometry consisted of an inlet housing, inlet guide vanes (IGV), a rotor, and a stator. The RANS model for both simulations was the k-epsilon model. For both of these cases, sliding mesh interfaces were located between the IGV and rotor, and between the rotor and stator. The results of the URANS and DES approaches were time-averaged and their predictions were compared. Throughout the study, our goal was to provide important insights into the performance of the URANS and DES approaches, and to highlight the essential differences.


Author(s):  
Hui Li ◽  
Xiutao Bian ◽  
Xinrong Su ◽  
Xin Yuan

Abstract The complex leakage flow structure in the tip region of unshrouded rotor is a main source of turbine aerodynamic loss. Due to the complex turbulence characteristics of the tip leakage flow, the widely used Reynolds Averaged Navier-Stokes (RANS) approach may fail to accurately predict the multi-scale turbulent flow and the related loss. In order to effectively improve the turbine efficiency, more insights into the turbulence characteristics and the loss mechanism in the tip leakage flow are required. In this work, a Delayed Detached Eddy Simulation (DDES) study is conducted to simulate the flow inside a high pressure turbine blade, with emphasis on the tip region. DDES results are in good agreement with the experiment and the comparison with RANS results verifies the advantages of DDES in resolving finer flow structures of leakage flow, also in capturing the complex turbulence characteristics. The snapshot Proper Orthogonal Decomposition (POD) method is used to extract the dominant flow features. The flow structures and the distribution of Reynolds stress help to reveal the process of leakage flow and its interaction with the secondary flow. Meanwhile, it is found that the separation vortex (SV) forms from leading edge to trailing edge, and the strong interactions between tip leakage vortex (TLV) and passage secondary vortex (PSV) significantly enhance the turbulence intensity. Based on the DDES results, loss analysis of tip leakage flow is conducted based on entropy generation rates. For the leakage flow related loss, the largest local entropy generation rate occurs at 50 % of axial chord, and the interaction between the leakage vortex and up passage vortex promotes the loss generation. To sum up, the current DDES study about the tip leakage flow provides helpful information about the loss generation mechanism and may guide the design of low-loss blade tip.


Author(s):  
Chenkai Zhang ◽  
Jun Hu ◽  
Zhiqiang Wang ◽  
Wei Yan ◽  
Chao Yin ◽  
...  

To deepen the knowledge of tip leakage flow/vortex flow structure in the tip clearance of axial compressor rotors, this paper presents steady numerical studies on a subsonic rotor. The rotor and its related low-speed large-scale repeating-stage axial compressor are used for low-speed model testing of a modern high-pressure compressor. Results were first compared with available experimental data to validate adopted numerical method. Then complex endwall flow structure and flow loss mechanism at design operating point were studied. At last, comparisons were made for tip leakage vortex structure, interface of the leakage flow/main flow, endwall blockage and loss between design and near-stall operating points. Results show that only the spilled flows below 62.5% clearance height at the leading edge will roll into tip leakage vortex for this rotor. In addition, tip leakage vortex plays a secondary important role for higher positions, where secondary leakage flow occurs and occupies broader chordwise range. Although tip leakage vortex would expand and strongly mix with the mainflow when it propagates downstream, which leads to a rapid reduction of the normalized streamwise vorticity, the value of the normalized helicity shows that concentrated vortex feature is still maintained.


Sign in / Sign up

Export Citation Format

Share Document