Prediction of Reynolds Number Effects on Low-Pressure Turbines Using a High-Order ILES Method

2020 ◽  
Vol 142 (3) ◽  
Author(s):  
M. Bolinches-Gisbert ◽  
David Cadrecha Robles ◽  
Roque Corral ◽  
Fernando Gisbert

Abstract A comprehensive comparison between implicit large eddy simulations (ILES) and experimental results of a modern high-lift low-pressure turbine airfoil has been carried out for an array of Reynolds numbers (Re). Experimental data were obtained in a low-speed linear cascade at the Polytechnic University of Madrid using hot-wire anemometry and laser-Doppler velocimetry (LDV). The numerical code is fourth-order accurate, both in time and space. The spatial discretization of the compressible Navier–Stokes equations is based on a high-order flux reconstruction approach while a fourth-order Runge–Kutta method is used to march in time the simulations. The losses, pressure coefficient distributions, and boundary layer and wake velocity profiles have been compared for an array of realistic Reynolds numbers. Moreover, boundary layer and wake velocity fluctuations are compared for the first time with experimental results. It is concluded that the accuracy of the numerical results is comparable to that of the experiments, especially for integral quantities such as the losses or exit angle. Turbulent fluctuations in the suction side boundary layer and the wakes are well predicted too. The elapsed time of the simulations is about 140 h on 40 graphics processor units. The numerical tool is integrated within an industrial design system and reuses pre- and post-processing tools previously developed for another kind of applications. The trend of the losses with the Reynolds number has a sub-critical regime, where the losses scale with Re−1, and a supercritical regime, where the losses scale with Re−1/2. This trend can be seen both in the simulations and in the experiments.

Author(s):  
Marc Bolinches-Gisbert ◽  
David Cadrecha Robles ◽  
Roque Corral ◽  
Fernando Gisbert

Abstract A comprehensive comparison between Implicit Large Eddy Simulations (ILES) and experimental results of a modern highlift low-pressure turbine airfoil has been carried out for an array of Reynolds numbers (Re). Experimental data were obtained in a low-speed linear cascade at the Polithecnic University of Madrid using hot-wire anemometry and LDV. The numerical code is fourth order accurate, both in time and space. The spatial discretization of the compressible Navier-Stokes equations is based on a high-order Flux Reconstruction approach while a fourth order Runge-Kutta method is used to march in time the simulations. The losses, pressure coefficient distributions, and boundary layer and wake velocity profiles have been compared for an array of realistic Reynolds numbers. Moreover, boundary layer and wake velocity fluctuations are compared for the first time with experimental results. It is concluded that the accuracy of the numerical results is comparable to that of the experiments, especially for integral quantities such as the losses or exit angle. Turbulent fluctuations in the suction side boundary layer and the wakes are well predicted also. The elapsed time of the is about 140 hours on 40 Graphics Processor Units. The numerical tool is integrated within an industrial design system and reuses pre- and post-processing tools previously developed for another kind of applications. The trend of the losses with the Reynolds number has a sub-critical regime, where the losses scale with Re−1, and a supercrital regime, where the losses scale with Re−1/2. This trend can be seen both, in the simulations and the experiments.


Author(s):  
D. S. Martínez ◽  
E. Pescini ◽  
F. Marra ◽  
M. G. De Giorgi ◽  
A. Ficarella

The present work is focused on the investigation of an alternate current driven single dielectric barrier discharge plasma actuator (AC-SDBDPA) for the control of separated flow at Reynolds numbers up to 2·104. Laminar boundary layer separation typically occurs on the suction surface of the low pressure turbines (LPT) blades when operating at high altitude cruise conditions, as the Reynolds number can drop below 2.5·104. In this context, the implementation of an active boundary layer control system able to operate in suppressing separation — only at the critical Reynolds numbers — is of great interest. The SDBDPA was manufactured by means of the photolithographic technique, which ensured a thin metal deposition with high manufacturing reliability control. Actuator operation under sinusoidal voltage at 8 kV amplitude and 2 kHz frequency was considered. Investigations were performed in a closed loop wind tunnel. A curved plate with a shape designed to reproduce the suction surface of a LPT was mounted directly over the bottom wall of the test section. The SDBDPA was inserted in a groove made at the middle of the curved plate, located at the front side of the adverse pressure gradient region. The flow pattern and velocities in absence of actuation were experimentally measured by a two-dimensional (2-D) particle image velocimetry (PIV) system and a laser Doppler velocimetry (LDV) system. PIV measurements were performed in presence of actuation. Simultaneously to the velocity measurements, the voltage applied to the AC-SDBDPA and the discharge current flowing through the circuit were acquired in order to determine the power dissipated by the device. The experimental data were supported by computational fluid dynamics (CFD) simulations based on the finite volume method. In order to deeply investigate the effect of flow separation control by the AC-SDBDPA on the LPT blade performances, the viscous and unsteady Reynolds-averaged Navier-Stokes equations were solved to predict the characteristics of the flow with and without actuation. The actuation effect was modelled as a time-constant body force calculated prior to the fluid flow simulations by using the dual potential algebraic model. The experimental data were used to calibrate and successfully validate the numerical model. An unsteady RANS (URANS) approach, using the k-ω Lam and Bremhorst Low-Reynolds turbulence model was employed, accounting with the main transient flow structures. Results showed that the mixing action of the streamwise fluid with higher momentum and the boundary layer fluid with the lower momentum -due to the AC-SDBDPA-led, depending on the tested Reynolds number, to the alleviation or suppression of the boundary layer flow separation which occurred on the suction surface of the LPT blade. The validated numerical model will allow expanding the study of the actuation effect including different locations and multiple devices, saving considerably experimental efforts.


Author(s):  
Marion Mack ◽  
Roland Brachmanski ◽  
Reinhard Niehuis

The performance of the low pressure turbine (LPT) can vary appreciably, because this component operates under a wide range of Reynolds numbers. At higher Reynolds numbers, mid and aft loaded profiles have the advantage that transition of suction side boundary layer happens further downstream than at front loaded profiles, resulting in lower profile loss. At lower Reynolds numbers, aft loading of the blade can mean that if a suction side separation exists, it may remain open up to the trailing edge. This is especially the case when blade lift is increased via increased pitch to chord ratio. There is a trend in research towards exploring the effect of coupling boundary layer control with highly loaded turbine blades, in order to maximize performance over the full relevant Reynolds number range. In an earlier work, pulsed blowing with fluidic oscillators was shown to be effective in reducing the extent of the separated flow region and to significantly decrease the profile losses caused by separation over a wide range of Reynolds numbers. These experiments were carried out in the High-Speed Cascade Wind Tunnel of the German Federal Armed Forces University Munich, Germany, which allows to capture the effects of pulsed blowing at engine relevant conditions. The assumed control mechanism was the triggering of boundary layer transition by excitation of the Tollmien-Schlichting waves. The current work aims to gain further insight into the effects of pulsed blowing. It investigates the effect of a highly efficient configuration of pulsed blowing at a frequency of 9.5 kHz on the boundary layer at a Reynolds number of 70000 and exit Mach number of 0.6. The boundary layer profiles were measured at five positions between peak Mach number and the trailing edge with hot wire anemometry and pneumatic probes. Experiments were conducted with and without actuation under steady as well as periodically unsteady inflow conditions. The results show the development of the boundary layer and its interaction with incoming wakes. It is shown that pulsed blowing accelerates transition over the separation bubble and drastically reduces the boundary layer thickness.


1985 ◽  
Vol 160 ◽  
pp. 281-295 ◽  
Author(s):  
F. A. Milinazzo ◽  
P. G. Saffman

Computations of two-dimensional solutions of the Navier–Stokes equations are carried out for finite-amplitude waves on steady unidirectional flow. Several cases are considered. The numerical method employs pseudospectral techniques in the streamwise direction and finite differences on a stretched grid in the transverse direction, with matching to asymptotic solutions when unbounded. Earlier results for Poiseuille flow in a channel are re-obtained, except that attention is drawn to the dependence of the minimum Reynolds number on the physical constraint of constant flux or constant pressure gradient. Attempts to calculate waves in Couette flow by continuation in the velocity of a channel wall fail. The asymptotic suction boundary layer is shown to possess finite-amplitude waves at Reynolds numbers orders of magnitude less than the critical Reynolds number for linear instability. Waves in the Blasius boundary layer and unsteady Rayleigh profile are calculated by employing the artifice of adding a body force to cancel the spatial or temporal growth. The results are verified by comparison with perturbation analysis in the vicinity of the linear-instability critical Reynolds numbers.


Author(s):  
Can Ma ◽  
Xin Yuan

This paper numerically investigates the transitional flow on a LPT (low pressure turbine) blade with fluctuating downstream potential field. A linear T106 cascade is subjected to an oscillating potential field generated by downstream moving bars. Previous experimental results in open literature showed that the unsteady downstream potential field has an obvious influence on the transitional boundary layer of LPT blade. For the numerical simulations in this paper, the unsteady Reynolds-Averaged Navier-Stokes equations are solved using the commercial software FLUENT. The transition model used in this paper is the γ-Reθ model, which has been validated against a number of transitional flows previously, including the influence of upstream wakes on the transitional boundary layer of T106 turbine blade. The simulation results are first compared to the experimental results in open literature to validate the numerical methods. Two different FSTI (free stream turbulence intensity), 1.6% and 4.0% are investigated with axial spacing between the blade and the downstream bar varying from 50% axial chord to 25% axial chord. To investigate the influence of flow compressibility, two different inlet Mach numbers, 0.02 and 0.2 are simulated. Results show that decreasing the axial spacing has an influence on the unsteady boundary layer separation and transition and the influence is enhanced at elevated inlet Mach number.


Author(s):  
Ali Mahallati ◽  
Steen A. Sjolander

The relative motion of rotor and stator blade rows causes periodically unsteady flows that influence the performance of airfoils through their effects on the boundary layer development. Part 1 of this two-part paper described the influence of Reynolds number, freestream turbulence intensity and turbulence length scales on a low-pressure (LP) high-lift turbine airfoil, PakB, under steady inlet flow conditions. The aerodynamic behaviour of the same airfoil under the influence of incoming wakes is presented in Part 2. The unsteady effects of wakes from a single upstream blade-row were measured in a low-speed linear cascade facility at Reynolds numbers of 25000, 50000 and 100000 and at two freestream turbulence intensity levels of 0.4% and 4%. In addition, eight reduced frequencies between 0.53 and 3.2, at three flow coefficients of 0.5, 0.7 and 1.0 were examined. The complex wake-induced transition, flow separation and reattachment on the suction surface boundary layer was determined from an array of closely-spaced surface hot-film sensors. The wake-induced transition caused the separated boundary layer to reattach to the suction surface at all conditions examined. The time-varying profile losses were measured downstream of the trailing edge. Profile losses increase with decreasing Reynolds number and the influence of increased freestream turbulence intensity is only evident in between wake-passing events at low reduced frequencies. At higher values of reduced frequency, the losses increase slightly and for the cases examined here, losses were slightly larger at lower flow coefficients than the higher flow coefficients. An optimum wake-passing frequency was observed at which the profile losses were a minimum.


Author(s):  
Xiao Qu ◽  
Yanfeng Zhang ◽  
Xingen Lu ◽  
Ge Han ◽  
Ziliang Li ◽  
...  

Periodic wakes affect not only the surface boundary layer characteristics of low-pressure turbine blades and profile losses but also the vortex structures of the secondary flow and the corresponding losses. Thus, understanding the physical mechanisms of unsteady interactions and the potential to eliminate secondary losses is becoming increasingly important for improving the performance of high-lift low-pressure turbines. However, few studies have focused on the unsteady interaction mechanism between periodic wakes and endwall secondary flow in low-pressure turbines. This paper verified the accuracy of computational fluid dynamics by comparing experimental results and those of the numerical predictions by taking a high-lift low-pressure turbine cascade as the research object. Discussion was focused on the interaction mechanisms between the upstream wakes and secondary flow within the high-lift low-pressure turbine. The results indicated that upstream wakes have both positive and negative effects on the endwall flow, where the periodic wakes can decrease significantly the size of the separation bubble, prevent the formation of secondary vorticity structures at relatively high Reynolds numbers (100,000 and 150,000), and reduce the cross-passage pressure gradient of cascade. In addition, periodic wakes can improve the cascade incidence characteristic in terms of reducing the overturning and underturning of the secondary flow at downstream of the cascade all of which are beneficial for decreasing the endwall secondary losses, whereas more endwall boundary layer is involved in the main flow passage due to the wake transport, resulting in increased strength of the secondary flow at low Reynolds number of 25,000 and 50,000. Compared with the results without wakes, the total pressure loss for unsteady condition at the cascade exit decreases by 2.7% and 6.1% at high Reynolds number of 100,000 and 150,000, respectively. However, the secondary loss at unsteady flow conditions increases at low Reynolds number of 25,000 and 50,000.


Author(s):  
Muhammad Aqib Chishty ◽  
Khalid Parvez ◽  
Sijal Ahmed ◽  
Hossein Raza Hamdani ◽  
Ammar Mushtaq

The boundary layer of low-pressure turbine blades has received a great deal of attention due to advent of high lift and ultra high lift LP turbines. At cruising condition, Reynolds number is very low in engine and LP turbine performance suffers mainly from losses due to the laminar separation bubble on suction surface. In this paper, T106A low pressure turbine profile has been used to study the behavior of boundary layer and subsequently, flow is controlled using the passive technique. Unsteady Reynolds Averaged Navier Stokes equations were solved using SST Gamma-Theta transition model for turbulence closure. Hybrid mesh topology has been used to discretize the computational domain, with highly resolved structured mesh in boundary layer (Y+ < 1) and unstructured mesh in the rest of domain. Simulations were performed using commercial CFD code ANSYS FLUENT ® at Reynolds number 91000 (based on inlet velocity and chord length) and turbulence intensity of 0.4%. To study the effect of dimple on the flow separation, dimple has been positioned at different axial location on the suction side. It was found that shifting the dimple downstream results in controlled flow and reduced loss coefficient as compared to the case when no dimple is applied.


1973 ◽  
Vol 59 (2) ◽  
pp. 391-396 ◽  
Author(s):  
N. C. Freeman ◽  
S. Kumar

It is shown that, for a spherically symmetric expansion of a gas into a low pressure, the shock wave with area change region discussed earlier (Freeman & Kumar 1972) can be further divided into two parts. For the Navier–Stokes equation, these are a region in which the asymptotic zero-pressure behaviour predicted by Ladyzhenskii is achieved followed further downstream by a transition to subsonic-type flow. The distance of this final region downstream is of order (pressure)−2/3 × (Reynolds number)−1/3.


1991 ◽  
Vol 113 (4) ◽  
pp. 608-616 ◽  
Author(s):  
H. M. Jang ◽  
J. A. Ekaterinaris ◽  
M. F. Platzer ◽  
T. Cebeci

Two methods are described for calculating pressure distributions and boundary layers on blades subjected to low Reynolds numbers and ramp-type motion. The first is based on an interactive scheme in which the inviscid flow is computed by a panel method and the boundary layer flow by an inverse method that makes use of the Hilbert integral to couple the solutions of the inviscid and viscous flow equations. The second method is based on the solution of the compressible Navier–Stokes equations with an embedded grid technique that permits accurate calculation of boundary layer flows. Studies for the Eppler-387 and NACA-0012 airfoils indicate that both methods can be used to calculate the behavior of unsteady blade boundary layers at low Reynolds numbers provided that the location of transition is computed with the en method and the transitional region is modeled properly.


Sign in / Sign up

Export Citation Format

Share Document