A Comparison of Unsteady RANS and DES for Simulating an Axial Compressor Stage

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.

2021 ◽  
Author(s):  
David Vanpouille ◽  
Dimitrios Papadogiannis ◽  
Stéphane Hiernaux

Abstract Surge margin is critical for the safety of aeronautical compressors, hence predicting it early in the design process using CFD is mandatory. However, close to surge, steady-state Reynolds Averaged Navier-Stokes (RANS) simulations are proven inadequate. Unsteady techniques such as Unsteady RANS (URANS) and Large Eddy Simulation (LES) can provide more reliable predictions. Nevertheless, the accuracy of such methods are dependent on the method used to handle the rotor/stator interfaces. The most precise method, the sliding mesh, requires simulating the full annulus or a periodic sector, which can be very costly. Other techniques to reduce the domain exist, such as the phase-lagged approach or geometric blade scaling, but introduce restrictive assumptions on the flow at near-stall conditions. The objective of this paper is to investigate the near-stall flow of a low-pressure compressor using unsteady methods of varying fidelity: URANS with the phase lag assumption, URANS on a periodic sector and a high-fidelity LES on a smaller periodic sector achieved using geometric blade scaling. Results are compared to experimental measurements. An overall good agreement is found. Results show that the tip leakage vortex is not the origin of the stall on the studied configuration and a hub corner separation is initiated. LES further validates the (U)RANS flow predictions and brings additional insight on unsteady flow separations.


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.


2005 ◽  
Vol 109 (1102) ◽  
pp. 631-638 ◽  
Author(s):  
R. M. Ashworth

AbstractIt is argued that acoustic resonance phenomena in open cavities such as weapons bays cannot be adequately predicted through numerical solution of Reynolds averaged Navier-Stokes (RANS) equations. The requirement to resolve the growth of the shear layer instability from the lip of the cavity inevitably implies that turbulence further downstream is resolved while also being modelled thus making RANS over dissipative. Large eddy simulation (LES) models only unresolved scales and a hybrid method combining RANS near walls with LES in the cavity appears a practical alternative to pure RANS. This paper compares computations of the M219 cavity configuration made with unsteady RANS and with the hybrid method known as detached eddy simulation (DES). It is shown that whilst unsteady RANS and DES give very similar predictions for the 1stand 3rdmodes of the acoustic resonance the 2ndmode (which is dominant near the centre of the cavity) is absent in the RANS results but well predicted by DES. The 2ndmode is thought to arise from an interaction with vortical structures in the shear layer which are suppressed in the highly dissipative RANS method. The 4thmode, which is much weaker than the other three modes, is over-predicted by DES and under-predicted by a smaller amount in RANS.


2019 ◽  
Vol 33 (20) ◽  
pp. 1950228
Author(s):  
Chunlai Tian ◽  
Tairan Chen ◽  
Tian Zou

Unsteady cavitating flow with high Reynolds number and significant instability commonly exists in fluid machinery and engineering system. The high-resolution approaches, such as direct numerical simulation and large eddy simulation, are not practical for engineering issues due to the significant cost in the computational resource. The objective of this paper is to provide the approach with Detached-Eddy Simulation (DES) model based on the Reynolds-averaged Navier–Stokes (RANS) equations for predicting unsteady cavitating flows. The credibility of the approach is validated by a set of numerical examples of its application: the unsteady cavitating flows around the two-dimensional (2D) Clark-Y hydrofoil and the three-dimensional (3D) blunt body. It is found that the calculated cavity shapes, cavity lengths and unsteady characteristics by DES model agree well with the experimental measurements and observations. Further analysis indicates that the turbulent eddy viscosity around the cavity and wake region is well predicted by the DES model, which results in the development of large-scale vortexes, and further cavitation instability. The DES model, which exhibits a significantly unsteady 3D behavior, is a more comprehensive turbulence model for unsteady cavitating flows.


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.


2019 ◽  
Author(s):  
Κωνσταντίνος Διακάκης

Στην παρούσα διατριβή μελετήθηκε η μετάβαση της ροής από στρωτή σε τυρβώδη καθώς και η συμπεριφορά ροών μεγάλων αριθμών Reynolds στα πλαίσια προσομοίωσης τους με μεθόδους υψηλής πιστότητας.Για την προσομοίωση ροών με μετάβαση εξετάστηκαν μέθοδοι με υπολογισμό οριακού στρώματος και μέθοδοι με εξισώσεις μεταφοράς. Αυτές περιλαμβάνουν την μέθοδο e N καθώς και τα μοντέλα γ-Re θ , γ και AFT. Όλες οι μέθοδοι δοκιμάστηκαν σε αεροτομές, πτέρυγες και άτρακτο γενικής μορφής, σε εφαρμογές οι οποίες προέρχονταν από τους τομείς της αεροναυτικής και της αιολικής ενέργειας. Οι συγκρίσεις αφορούσαν κατά κύριο λόγο σε αεροδυναμικά φορτία και θέσεις μετάβασης. Στα πλαίσια διδιάστατων προσομοιώσεων, η μέθοδος e N με υπολογισμό οριακού στρώματος και το μοντέλο AFT έδωσαν πιο ακριβή αποτελέσματα από τις υπόλοιπες μεθόδους. Το μοντέλο γ-Re θ είναι μια καλή εναλλακτική, αρκεί ο αριθμός Reynolds να μην υπερβαίνει τα 6 εκατομμύρια. Πέραν αυτού του ορίου, η ακρίβεια των αποτελεσμάτων του μοντέλου μειώνεται σημαντικά. Ωστόσο, η μέθοδος e N και το μοντέλο AFT δεν δύνανται να χρησιμοποιηθούν για την μοντελοποίηση τρισδιάστατης μετάβασης στο πλαίσιο τρισδιάστατων προσομοιώσεων. Σε αυτές τις περιπτώσεις, το μοντέλο γ-Re θ εμπλουτισμένο με όρους εγκάρσιας ροής μπορεί να δώσει καλά αποτελέσματα, αρκεί ο αριθμός Reynolds να είναι στα αποδεκτά για το μοντέλο όρια. Όσον αφορά στις μεθόδους προσομοίωσης τύρβης υψηλής πιστότητας, εξετάστηκαν οι μέθοδοι Large Eddy Simulation (LES) και Detached Eddy Simulation (DES). Για τις προσομοιώσεις LES χρησιμοποιήθηκε το μοντέλο μικρών κλιμάκων του Smagorinsky. Η εφαρμογή του DES περιελάμβανε τις μεθόδους Delayed DES (DDES) και Improved Delayed DES (IDDES). Το ενδιαφέρον εστιάστηκε στην μοντελοποίηση ροών με μεγάλη αποκόλληση. Τόσο το LES όσο και το DES ήταν σε θέση να δώσουν πιο ακριβή αποτελέσματα από τους απλούς, μη-μόνιμους Reynolds Averaged Navier Stokes υπολογισμούς (Unsteady RANS) σε σύγκριση με πειράματα και υπολογιστικά αποτελέσματα από τη βιβλιογραφία. Το μοντέλο DES θεωρείται λιγότερο απαιτητικό σε υπολογιστικούς πόρους λόγω της μοντελοποίησης του οριακού στρώματος η οποία οδηγεί σε μικρότερες απαιτήσεις πλέγματος κοντά στην στερεή επιφάνεια. Ωστόσο, το DES δεν αναμένεται να μπορεί να δώσει αξιόπιστα αποτελέσματα σε ροές όπου η παρουσία και η εξέλιξη μικρών κλιμάκων τύρβης στο οριακό στρώμα είναι σημαντική, και που το μοντέλο LES πλεονεκτεί εκ κατασκευής. Σχετικά σημειώνεται ότι οι LES προσομοιώσεις δεν έφτασαν στα υπολογιστικά τους όρια όσον αφορά στο πλέγμα. Για να παραχθούν αξιόπιστα αποτελέσματα σε αυτές τις περιπτώσεις πρέπει να χρησιμοποιηθεί LES με πυκνό υπολογιστικό πλέγμα.


Author(s):  
Julien Marty ◽  
Cédric Uribe

The present study focuses on the impact of the underlying RANS turbulence model in the Zonal Detached Eddy Simulation (ZDES) method when used for secondary flow prediction. This is carried out in light of three issues commonly investigated for hybrid RANS/LES methods (detection and protection of attached boundary layer, emergence, and growth of resolved turbulent fluctuations and accurate prediction of separation front due to progressive adverse pressure gradient). The studied configuration is the first rotor of a high pressure compressor. Three different turbulence modelings (Spalart and Allmaras model (SA), Menter model with (SST) and without (BSL) shear stress correction) are assessed as ZDES underlying turbulence model and also as turbulence model of unsteady RANS simulations. Whatever the underlying turbulence model, the ZDES behaves well with respect to the first two issues as the boundary layers appear effectively shielded and the RANS-to-LES switch is close downstream of trailing edges and separation fronts leading to a quick LES treatment of wakes and shear layers. Both tip leakage and corner flows are strongly influenced by the Navier–Stokes resolution approach (unsteady RANS vs. ZDES) but the underlying turbulence modelling (SA vs. SST vs. BSL) impacts mainly the junction flow near the hub for both approaches. ZDES underlying turbulence model choice appear essential since it leads to quite different corner flow separation topologies and so to inversion of the downstream stagnation pressure radial gradient.


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.


Author(s):  
Wei Ma ◽  
Feng Gao ◽  
Xavier Ottavy ◽  
Lipeng Lu ◽  
A. J. Wang

Recently bimodal phenomenon in corner separation has been found by Ma et al. (Experiments in Fluids, 2013, doi:10.1007/s00348-013-1546-y). Through detailed and accurate experimental results of the velocity flow field in a linear compressor cascade, they discovered two aperiodic modes exist in the corner separation of the compressor cascade. This phenomenon reflects the flow in corner separation is high intermittent, and large-scale coherent structures corresponding to two modes exist in the flow field of corner separation. However the generation mechanism of the bimodal phenomenon in corner separation is still unclear and thus needs to be studied further. In order to obtain instantaneous flow field with different unsteadiness and thus to analyse the mechanisms of bimodal phenomenon in corner separation, in this paper detached-eddy simulation (DES) is used to simulate the flow field in the linear compressor cascade where bimodal phenomenon has been found in previous experiment. DES in this paper successfully captures the bimodal phenomenon in the linear compressor cascade found in experiment, including the locations of bimodal points and the development of bimodal points along a line that normal to the blade suction side. We infer that the bimodal phenomenon in the corner separation is induced by the strong interaction between the following two facts. The first is the unsteady upstream flow nearby the leading edge whose angle and magnitude fluctuate simultaneously and significantly. The second is the high unsteady separation in the corner region.


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.


Sign in / Sign up

Export Citation Format

Share Document