Delayed Detached Eddy Simulation of Rotating Stall for a Full Annulus Transonic Axial Compressor Stage

Author(s):  
Jiaye Gan ◽  
Hong-Sik Im ◽  
Ge-Cheng Zha

This paper solves the filtered Navier-Stokes equations to simulate stall inception of NASA compressor transonic Stage 35 with delayed detached eddy simulation (DDES). A low diffusion E-CUSP Riemann solver with a 3rd order MUSCL scheme for the inviscid fluxes and a 2nd order central differencing for the viscous terms are employed. A full annulus of the rotor-stator stage is simulated with an interpolation sliding boundary condition (BC) to resolve the rotor-stator interaction. The tip clearance is fully gridded to accurately resolve tip vortices and their effect on stall inception. The DDES results show that the stall inception of Stage 35 is initialized by a weak harmonic disturbance with the length scales of the full annulus and grows rapidly with two emerging spike like disturbance. The two spike disturbances propagate in counter rotational direction with about 42% of rotor speed. The spike stall cells cover about 6 blades. They lead to two stall cells grown circumferentially and inwardly.

Author(s):  
Hongsik Im ◽  
Xiangying Chen ◽  
Gecheng Zha

Detached eddy simulation of an aeroelastic self-excited instability, flutter in NASA Rotor 67 is conducted using a fully coupled fluid/structre interaction. Time accurate compressible 3D Navier-Stokes equations are solved with a system of 5 decoupled modal equations in a fully coupled manner. The 5th order WENO scheme for the inviscid flux and the 4th order central differencing for the viscous flux are used to accurately capture interactions between the flow and vibrating blades with the DES (detached eddy simulation) of turbulence. A moving mesh concept that can improve mesh quality over the rotor tip clearance was implemented. Flutter simulations were first conducted from choke to stall using 4 blade passages. Stall flutter initiated at rotating stall onset, grows dramatically with resonance. The frequency analysis shows that resonance occurs at the first mode of the rotor blade. Before stall, the predicted responses of rotor blades decayed with time, resulting in no flutter. Full annulus simulation at peak point verifies that one can use the multi-passage approach with periodic boundary for the flutter prediction.


Author(s):  
Purvic Patel ◽  
Yunchao Yang ◽  
Gecheng Zha

Abstract This paper utilizes the Improved Delayed Detached Eddy Simulation (IDDES) to investigate the non-synchronous vibration (NSV) mechanism of a 1.5 stage high-speed axial compressor. The NSV occurs at a part speed in the rig test. A low diffusion E-CUSP approximate Riemann solver with a third order Weighted Essentially Non-Oscillating (WENO) scheme for the inviscid flux and a second order central differencing scheme for the viscous flux are employed to solve the 3D time accurate Navier-Stokes equations. The fully conservative sliding boundary condition is used to preserve the wake-propagation. The aerodynamic instability in the tip region induces two alternating low pressure regions near the leading and the trailing edge on the suction side of the rotor blade. It is observed that the circumferential tip vortex motion in the rotor passage above 75 % span and its coupling forces cause NSV at the operating speed. This instability moves in the counter-rotating direction in the rotational frame. The NSV results using URANS simulation is also presented for comparison. The predicted frequency with the IDDES and URANS using rigid blades agrees well with the measured frequency in the rig test. In addition to the NSV, the IDDES solver also captures the dominant engine order frequencies. The tip flow structures show the vortex filament with one end on the suction side of the rotor blade and other side terminating on the casing or the pressure side of the rotor blade.


Author(s):  
Chunill Hah ◽  
Jo¨rg Bergner ◽  
Heinz-Peter Schiffer

The current paper reports on investigations aimed at advancing the understanding of the flow mechanism that leads to the onset of short-length scale rotating stall in a transonic axial compressor. Experimental data show large oscillation of the tip clearance vortex as the rotor operates near the stall condition. Inception of spike-type rotating stall is also measured in the current transonic compressor with high response pressure transducers. Computational studies of a single passage and the full annulus were carried out to identify flow mechanisms behind the spike-type stall inception in the current transonic compressor rotor. Steady and unsteady single passage flow simulations were performed, first to get insight into the interaction between the tip clearance vortex and the passage shock. The conventional Reynolds-averaged Navier-Stokes method with a standard turbulence closure scheme does not accurately reproduce tip clearance vortex oscillation and the measured unsteady pressure field. Consequently, a Large Eddy Simulation (LES) was carried out to capture more relevant physics in the computational simulation of the rotating stall inception. The unsteady random behavior of the tip clearance vortex and it’s interaction with the passage shock seem to be critical ingredients in the development of spike-type rotating stall in a transonic compressor. The Large Eddy Simulation was further extended to the full annulus to identify flow mechanisms behind the measured spike-type rotating stall inception. The current study shows that the spike-type rotating stall develops after the passage shock is fully detached from the blade passages. Interaction between the tip clearance vortex and the passage shock creates a low momentum area near the pressure side of the blade. As the mass flow rate decreases, this low momentum area moves further upstream and reversed tip clearance flow is initiated at the trailing edge plane. Eventually, the low momentum area near the pressure side reaches the leading edge and forward spillage of the tip clearance flow occurs. The flows in the affected blade passage or passages then stall. As the stalled blade passages are formed behind the passage shock, the stalled area rotates counter to the blade rotation just like the classical Emmon’s type rotating stall. Both the measurements and the computations show that the rotating stall cell covers one to two blade passage lengths and rotates at roughly 50% of the rotor speed.


Author(s):  
Hong-Sik Im ◽  
Xiangying Chen ◽  
Ge-Cheng Zha

This paper uses the advanced Delayed-Detached Eddy Simulation (DDES) of turbulence to simulate rotating stall inception of NASA Rotor 67. The rotor is a low-aspect-ratio transonic axial-flow fan with a tip speed of 429 m/s and a pressure ratio of 1.63. A full annulus simulation was employed with the time accurate compressible Navier-Stokes code in order to accurately capture the the formation of long-length disturbance and a short-length inception (spike). The validation for all numerical methods used in this study was accomplished by the comparisons of the CFD solutions with the test data in advance of unsteady simulations. Self-induced rotating stall development is simulated holding the same back pressure at the near stall experiment without any throttling. Spike type rotating stall occurs and rotates at roughly 50% of rotor speed counter to the rotation. After spike onset, rotating stall fully develops approximately within 2 rotor revolutions. Two distinct characteristics that can advance the mechanism of spike type rotating stall are observed. First, the passage shock is fully detached from rotor and decays during the spike inception. Consequently the shifted sonic line at the upstream of rotor allows stalling flow to propagate to the neighboring passage. Second, the trailing edge back flow contributes to the build up of a fully developed stall cell by pushing tip clearance flow toward blade leading edge and inducing tip spillage flow. Tip vortex originated from the leading edge dies out during spike inception as the swirl angle of incoming tip flow decreases, while in the unstalled passages it develops without breakdown. DDES challenge for the complete blade row reflects well the sequence of rotating stall and its unsteady behavior.


Author(s):  
Tong Li ◽  
Yibin Wang ◽  
Ning Zhao

The simple frigate shape (SFS) as defined by The Technical Co-operative Program (TTCP), is a simplified model of the frigate, which helps to investigate the basic flow fields of a frigate. In this paper, the flow fields of the different modified SFS models, consisting of a bluff body superstructure and the deck, were numerically studied. A parametric study was conducted by varying both the superstructure length L and width B to investigate the recirculation zone behind the hangar. The size and the position of the recirculation zones were compared between different models. The numerical simulation results show that the size and the location of the recirculation zone are significantly affected by the superstructure length and width. The results obtained by Reynolds-averaged Navier-Stokes method were also compared well with both the time averaged Improved Delayed Detached-Eddy Simulation results and the experimental data. In addition, by varying the model size and inflow velocity, various flow fields were numerically studied, which indicated that the changing of Reynolds number has tiny effect on the variation of the dimensionless size of the recirculation zone. The results in this study have certain reference value for the design of the frigate superstructure.


2002 ◽  
Vol 124 (2) ◽  
pp. 413-423 ◽  
Author(s):  
L. S. Hedges ◽  
A. K. Travin ◽  
P. R. Spalart

The flow around a generic airliner landing-gear truck is calculated using the methods of Detached-Eddy Simulation, and of Unsteady Reynolds-Averaged Navier-Stokes Equations, with the Spalart-Allmaras one-equation model. The two simulations have identical numerics, using a multi-block structured grid with about 2.5 million points. The Reynolds number is 6×105. Comparison to the experiment of Lazos shows that the simulations predict the pressure on the wheels accurately for such a massively separated flow with strong interference. DES performs somewhat better than URANS. Drag and lift are not predicted as well. The time-averaged and instantaneous flow fields are studied, particularly to determine their suitability for the physics-based prediction of noise. The two time-averaged flow fields are similar, though the DES shows more turbulence intensity overall. The instantaneous flow fields are very dissimilar. DES develops a much wider range of unsteady scales of motion and appears promising for noise prediction, up to some frequency limit.


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):  
M Künzelmann ◽  
R Urban ◽  
R Mailach ◽  
K Vogeler

The stable operating range of axial compressors is limited by the onset of rotating stall and surge. Mass injection upstream of the tip of an axial compressor rotor is a stability enhancement approach which can be effective in suppressing stall in tip-critical rotors, and thus increasing the operating range of compressors. In this article, investigations on active flow control related to the rotor tip gap sensitivity are discussed. The experiments were performed in a 1.5-stage low-speed research compressor. Measurements at part speed (80 per cent) and full speed (100 per cent) with varying injection rates are discussed. These tests were performed for two rotor tip clearances of 1.3 per cent and 4.3 per cent of rotor blade tip chord. Results on the compressor map, the flow field as well as transient measurements to identify the stall inception are discussed. Supplementary, the numerical results are compared to the experiments based on the configuration with the greatest benefit in operating range enhancement.


AIAA Journal ◽  
2016 ◽  
Vol 54 (8) ◽  
pp. 2377-2391 ◽  
Author(s):  
W. Riéra ◽  
J. Marty ◽  
L. Castillon ◽  
S. Deck

Sign in / Sign up

Export Citation Format

Share Document