Zonal RANS-LES modeling for turbines in aeroengines

The cost of large-eddy simulation (LES) modeling in various zones of gas turbine aeroengines is outlined. This high cost clearly demonstrates the need to perform hybrid Reynolds-averaged Navier-Stokes-LES (RANS-LES) over the majority of engine zones because the Reynolds number is too high for pure LES. The RANS layer is used to cover over the fine streaks found in the inner part of the boundary layer. The hybrid strategy is applied to various engine zones, which is shown to typically give much greater predictive accuracy than pure RANS simulations. However, the cost estimates show that the RANS layer should be disposed within the low-pressure turbine zone. Also, the nature of the flow physics in this zone makes LES most sensible. © 2012 by Begell House, Inc.

Mixed nonlinear LES for DES suitable flows

It is shown that a new mixed nonlinear/eddy viscosity LES model reproduces profiles better than a number of competing nonlinear and mixed models for plane channel flow. The objective is an LES method that produces a fully resolved turbulent boundary layer and could be applied to a variety of aerospace problems that are currently studied with RANS, RANS-LES, or DES methods that lack a true turbulent boundary layer. There are two components to the new model. One an eddy viscosity based upon the advected subgrid scale energy and a relatively small coefficient. Second, filtered nonlinear terms based upon the Leray regularization. Coefficients for the eddy viscosity and nonlinear terms come from LES tests in decaying, isotropic turbulence. Using these coefficients, the velocity profile matches measurements data at Reτ ≈ 1000 exactly. Profiles of the components of kinetic energy have the same shape as in the experiment, but the magnitudes differ by about 25%. None of the competing LES gets the shape correct. This method does not require extra operations at the transition between the boundary layer and the interior flow.

Turbine blade tip heat transfer in low speed and high speed flows

In this paper, high and low speed tip flows are investigated for a high-pressure turbine blade. Previous experimental data are used to validate a CFD code, which is then used to study the tip heat transfer in high and low speed cascades. The results show that at engine representative Mach numbers the tip flow is predominantly transonic. Thus, compared to the low speed tip flow, the heat transfer is affected by reductions in both the heat transfer coefficient and the recovery temperature. The high Mach numbers in the tip region (M>1.5) lead to large local variations in recovery temperature. Significant changes in the heat transfer coefficient are also observed. These are due to changes in the structure of the tip flow at high speed. At high speeds, the pressure side corner separation bubble reattachment occurs through supersonic acceleration which halves the length of the bubble when the tip gap exit Mach number is increased from 0.1 to 1.0. In addition, shock/boundary-layer interactions within the tip gap lead to large changes in the tip boundary-layer thickness. These effects give rise to significant differences in the heat-transfer coefficient within the tip region compared to the low-speed tip flow. Compared to the low speed tip flow, the high speed tip flow is much less dominated by turbulent dissipation and is thus less sensitive to the choice of turbulence model. These results clearly demonstrate that blade tip heat transfer is a strong function of Mach number, an important implication when considering the use of low speed experimental testing and associated CFD validation in engine blade tip design. Copyright © 2009 by ASME.

Large eddy simulation of turbine internal cooling ducts

Abstract Large-Eddy Simulation (LES) and hybrid Reynolds-averaged Navier–Stokes–LES (RANS–LES) methods are applied to a turbine blade ribbed internal duct with a 180° bend containing 24 pairs of ribs. Flow and heat transfer predictions are compared with experimental data and found to be in agreement. The choice of LES model is found to be of minor importance as the flow is dominated by large geometric scale structures. This is in contrast to several linear and nonlinear RANS models, which display turbulence model sensitivity. For LES, the influence of inlet turbulence is also tested and has a minor impact due to the strong turbulence generated by the ribs. Large scale turbulent motions destroy any classical boundary layer reducing near wall grid requirements. The wake-type flow structure makes this and similar flows nearly Reynolds number independent, allowing a range of flows to be studied at similar cost. Hence LES is a relatively cheap method for obtaining accurate heat transfer predictions in these types of flows.

Scaling of transitionally-rough flow structures in DNSes of channels with riblets up to Reτ=550

The linear, drag-reducing effect of vanishingly small riblets breaks down once their size is in the transitionally-rough regime. We have previously reported that this breakdown is caused by the additional Reynolds stresses produced by the appearance of elongated spanwise rollers just above the riblet surface. These rollers are related with the Kelvin--Helmholtz instability of free shear layers, and to similar structures appearing over other rough and porous surfaces. However, because of the limited Reτ=180 in our previous DNSes, it could not be determined whether those structures scaled in inner or outer units. Furthermore, it is questionable if results in the transitionally-rough regime at Reτ=180 can be extrapolated to configurations of practical interest. At such small Reynolds numbers, roughness of transitional size can perturb a large portion of the boundary layer, which is not the case in most industrial and atmospheric applications. To clarify these issues we have conducted a set of DNSes at Reτ=550. Our results indicate that the spanwise rollers scale in wall units, and support the validity of the extrapolation to configurations of practical interest.

Kelvin-Helmholtz-like instabilities in turbulent flows over complex surfaces

Riblets are small surface protrusions aligned with the flow direction, which confer an anisotropic roughness to the surface [6]. We have recently reported that the transitional-roughness effect in riblets, which limits their performance, is due to a Kelvin–Helmholtz-like instability of the overlying mean flow [7]. According to our DNSs, the instability sets on as the Reynolds number based on the roughness size of the riblets increases, and coherent, elongated spanwise vortices begin to develop immediately above the riblet tips, causing the degradation of the drag-reduction effect. This is a very novel concept, since prior studies had proposed that the degradation was due to the interaction of riblets with the flow as independent units, either to the lodging of quasi-streamwise vortices in the surface grooves [2] or to the shedding of secondary streamwise vorticity at the riblet peaks [9]. We have proposed an approximate inviscid analysis for the instability, in which the presence of riblets is modelled through an average boundary condition for an overlying, spanwise-independent mean flow. This simplification lacks the accuracy of an exact analysis [4], but in turn applies to riblet surfaces in general. Our analysis succeeds in predicting the riblet size for the onset of the instability, while qualitatively reproducing the wavelengths and shapes of the spanwise structures observed in the DNSs. The analysis also connects the observations with the Kelvin–Helmholtz instability of mixing layers. The fundamental riblet length scale for the onset of the instability is a ‘penetration length,’ which reflects how easily the perturbation flow moves through the riblet grooves. This result is in excellent agreement with the available experimental evidence, and has enabled the identification of the key geometric parameters to delay the breakdown. Although the appearance of elongated spanwise vortices was unexpected in the case of riblets, similar phenomena had already been observed over other rough [3], porous [1] and permeable [11] surfaces, as well as over plant [5,14] and urban [12] canopies, both in the transitional and in the fully-rough regimes. However, the theoretical analyses that support the connection of these observations with the Kelvin–Helmholtz instability are somewhat scarce [7, 11, 13]. It has been recently proposed that Kelvin–Helmholtz-like instabilities are a dominant feature common to “obstructed” shear flows [8]. It is interesting that the instability does not require an inflection point to develop, as is often claimed in the literature. The Kelvin-Helmholtz rollers are rather triggered by the apparent wall-normal-transpiration ability of the flow at the plane immediately above the obstructing elements [7,11]. Although both conditions are generally complementary, if wall-normal transpiration is not present the spanwise vortices may not develop, even if an inflection point exists within the roughness [10]. REFERENCES [1] Breugem, W. P., Boersma, B. J. & Uittenbogaard, R. E. 2006 J. Fluid Mech. 562, 35–72. [2] Choi, H., Moin, P. & Kim, J. 1993 J. Fluid Mech. 255, 503–539. [3] Coceal, O., Dobre, A., Thomas, T. G. & Belcher, S. E. 2007 J. Fluid Mech. 589, 375–409. [4] Ehrenstein, U. 2009 Phys. Fluids 8, 3194–3196. [5] Finnigan, J. 2000 Ann. Rev. Fluid Mech. 32, 519–571. [6] Garcia-Mayoral, R. & Jimenez, J. 2011 Phil. Trans. R. Soc. A 369, 1412–1427. [7] Garcia-Mayoral, R. & Jimenez, J. 2011 J. Fluid Mech. doi: 10.1017/jfm.2011.114. [8] Ghisalberti, M. 2009 J. Fluid Mech. 641, 51–61. [9] Goldstein, D. B. & Tuan, T. C. 1998 J. Fluid Mech. 363, 115–151. [10] Hahn, S., Je, J. & Choi, H. 2002 J. Fluid Mech. 450, 259–285. [11] Jimenez, J., Uhlman, M., Pinelli, A. & G., K. 2001 J. Fluid Mech. 442, 89–117. [12] Letzel, M. O., Krane, M. & Raasch, S. 2008 Atmos. Environ. 42, 8770–8784. [13] Py, C., de Langre, E. & Moulia, B. 2006 J. Fluid Mech. 568, 425–449. [14] Raupach, M. R., Finnigan, J. & Brunet, Y. 1996 Boundary-Layer Meteorol. 78, 351–382.

An experimental study of loss sources in a fan operating with continuous inlet stagnation pressure distortion

The viability of Boundary Layer Ingesting (BLI) engines for future aircraft propulsion is dependent on the ability to design robust, efficient engine fan systems for operation with continuously distorted inlet flow. A key step in this process is to develop an understanding of the specific mechanisms by which an inlet distortion affects the performance of a fan stage. In this paper, detailed full-annulus experimental measurements of the flow field within a low-speed fan stage operating with a continuous 60-degree inlet stagnation pressure distortion are presented. These results are used to describe the three-dimensional fluid mechanics governing the interaction between the fan and the distortion and to make a quantitative assessment of the impact on loss generation within the fan. A 5.3 percentage point reduction in stage total-to-total efficiency is observed as a result of the inlet distortion. The reduction in performance is shown to be dominated by increased loss generation in the rotor due to off-design incidence values at its leading edge, an effect which occurs throughout the annulus despite the localised nature of the inlet distortion. Increased loss generation in the stator row is also observed due to flow separations that are shown to be caused by whirl angle distortion at rotor exit. By addressing these losses, it should be possible to achieve improved efficiency in BLI fan systems. Copyright © 2012 by ASME.

Vortical structures on three-dimensional shock control bumps

The vortical wake structure produced by a three-dimensional shock control bump (SCB) is thought to be useful for controlling transonic buffet on airfoils. However, at present the vorticity produced is relatively weak and the production mechanism is not well understood. Using a combined experimental and computational approach, a preliminary investigation on the wake vorticity for different bump geometries has been carried out. The structure of the wake for on and off-design conditions are considered, and the effects on the downstream boundary layer demonstrated. Three main vortical structures are observed: a primary vortex pair, weak inter-bump vortices and shear flow in the lambda-shock region. The effect of pressure gradients on vortex strength is examined and it is found that spanwise pressure gradients on the front section of the bump are the most significant parameter influencing vortex strength. © 2013 by S.P. Colliss et al.

Corner effects in reflecting oblique shock-wave/boundarylayer interactions

A separated oblique shock reflection on the floor of a rectangular cross-section wind tunnel has been investigated at M=2.5. The study aims to determine if and how separations occurring in the corners influence the main interaction as observed around the centreline of the floor. By changing the size of the corner separations through localised suction and small corner obstructions it was shown that the shape of the separated region in the centre was altered considerably. The separation length along the floor centreline was also modified by changes to the corner separation. A simple physical model has been proposed to explain the coupling between these separated regions based on the existence of compression or shock waves caused by the displacement effect of corner separation. These corner shocks alter the adverse pressure gradient imposed on the boundary-layer elsewhere which can lead to local reductions or increases of separation length. It is suggested that a typical oblique shock wave/boundary-layer interaction in rectangular channels features several zones depending on the relative position of the corner shocks and the main incident shock wave. Based on these findings the dependence of centre-line separation length on effective wind tunnel width is hypothesised. This requires further verification through experiments or computation. © 2013 by H. Babinsky.

Experimental study of surface roughness noise

A turbulent boundary-layer flow over a rough wall generates a dipole sound field as the near-field hydrodynamic disturbances in the turbulent boundary-layer scatter into radiated sound at small surface irregularities. In this paper, phased microphone arrays are applied to the experimental study of surface roughness noise. The radiated sound from two rough plates and one smooth plate in an open jet is measured at three streamwise locations, and the beamforming source maps demonstrate the dipole directivity. Higher source strengths can be observed in the rough plates than the smooth plate, and the rough plates also enhance the trailing-edge noise. A prediction scheme in previous theoretical work is used to describe the strength of a distribution of incoherent dipoles over the rigid plate and to simulate the sound detected by the microphone array. Source maps of measurement and simulation exhibit encouraging similarities in both source pattern and source strength, which confirms the dipole nature and the predicted magnitude of roughness noise. The simulations underestimate the streamwise gradient of the source strengths and overestimate the source strengths at the highest frequency. © 2007 by Yu Liu and Ann P. Dowling.

On high-frequency noise scattering by aerofoils in flow

Abstract A theoretical model is developed for the sound scattered when a sound wave is incident on a cambered aerofoil at non-zero angle of attack. The model is based on the linearization of the Euler equations about a steady subsonic flow, and is an adaptation of previous work which considered incident vortical disturbances. Only high-frequency sound waves are considered. The aerofoil thickness, camber and angle of attack are restricted such that the steady flow past the aerofoil is a small perturbation to a uniform flow. The singular perturbation analysis identifies asymptotic regions around the aerofoil; local 'inner' regions, which scale on the incident wavelength, at the leading and trailing edges of the aerofoil; Fresnel regions emanating from the leading and trailing edges of the aerofoil due to the coalescence of singularities and points of stationary phase; a wake transition region downstream of the aerofoil leading and trailing edge; and an outer region far from the aerofoil and wake. An acoustic boundary layer on the aerofoil surface and within the transition region accounts for the effects of curvature. The final result is a uniformly-valid solution for the far-field sound; the effects of angle of attack, camber and thickness are investigated. © 2013 Cambridge University Press.

Large eddy simulations in low-pressure turbines: Effect of wakes at elevated free-stream turbulence

The transition of a separated shear layer over a flat plate, in the presence of periodic wakes and elevated free-stream turbulence (FST), is numerically investigated using Large Eddy Simulation (LES). The upper wall of the test section is inviscid and specifically contoured to impose a streamwise pressure distribution over the flat plate to simulate the suction surface of a low-pressure turbine (LPT) blade. Two different distributions representative of a 'high-lift' and an 'ultra high-lift' turbine blade are examined. Results obtained from the current LES compare favourably with the extensive experimental data previously obtained for these configurations. The LES results are then used to further investigate the flow physics involved in the transition process.In line with experimental experience, the benefit of wakes and FST obtained by suppressing the separation bubble, is more pronounced in 'ultra high-lift' design when compared to the 'high-lift' design. Stronger 'Klebanoff streaks' are formed in the presence of wakes when compared to the streaks due to FST alone. These streaks promoted much early transition. The weak Klebanoff streaks due to FST continued to trigger transition in between the wake passing cycles.The experimental inference regarding the origin of Klebanoff streaks at the leading edge has been confirmed by the current simulations. While the wake convects at local free-stream velocity, its impression in the boundary layer in the form of streaks convects much slowly. The 'part-span' Kelvin-Helmholtz structures, which were observed in the experiments when the wake passes over the separation bubble, are also captured. The non-phase averaged space-time plots manifest that reattachment is a localized process across the span unlike the impression of global reattachment portrayed by phase averaging. © 2013 Elsevier Inc.

Average corotation of line segments near a point and vortex identification

An easy-to-interpret kinematic quantity measuring the average corotation of material line segments near a point is introduced and applied to vortex identification. At a given point, the vector of average corotation of line segments is defined as the average of the instantaneous local rigid-body rotation over "all planar cross sections" passing through the examined point. The vortex-identification method based on average corotation is a one-parameter, region-type local method sensitive to the axial stretching rate as well as to the inner configuration of the velocity gradient tensor. The method is derived from a well-defined interpretation of the local flow kinematics to determine the "plane of swirling" and is also applicable to compressible and variable-density flows. Practical application to direct numerical simulation datasets includes a hairpin vortex of boundary-layer transition, the reconnection process of two Burgers vortices, a flow around an inclined flat plate, and a flow around a revolving insect wing. The results agree well with some popular local methods and perform better in regions of strong shearing. Copyright © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Improving intermediate pressure turbine performance by using a nonorthogonal stator

Intermediate pressure (IP) turbines in high bypass ratio civil aeroengines are characterized by a significant increase in radius and a low aspect ratio stator. Conventional aerodynamic designs for the IP turbine stator have had leading and trailing edges orthogonal to the hub and casing end walls. The IP turbine rotor, however, is stacked radially due to stress limits. These choices inevitably lead to a substantial gap between the IP stator and rotor at the outer diameter in a duct that is generally diffusing the flow due to the increasing radius. In this low Mach number study, the IP stator is redesigned, incorporating compound sweep so that the leading and trailing edges are no longer orthogonal to the end walls. Computational investigations showed that the nonorthogonal stator reduces the flow diffusion between the stator and rotor, which yields two benefits: The stator trailing edge velocity was reduced, as was the boundary layer growth on the casing end wall within the gap. Experimental measurements confirm that the turbine with the nonorthogonal stator has an increased efficiency by 0.49%, while also increasing the work output by 4.6%, at the design point. © 2014 by ASME.

An experimental study of loss sources in a fan operating with continuous inlet stagnation pressure distortion

The viability of boundary layer ingesting (BLI) engines for future aircraft propulsion is dependent on the ability to design robust, efficient engine fan systems for operation with continuously distorted inlet flow. A key step in this process is to develop an understanding of the specific mechanisms by which an inlet distortion affects the performance of a fan stage. In this paper, detailed full-annulus experimental measurements of the flow field within a low-speed fan stage operating with a continuous 60 deg inlet stagnation pressure distortion are presented. These results are used to describe the three-dimensional fluid mechanics governing the interaction between the fan and the distortion and to make a quantitative assessment of the impact on loss generation within the fan. A 5.3 percentage point reduction in stage total-to-total efficiency is observed as a result of the inlet distortion. The reduction in performance is shown to be dominated by increased loss generation in the rotor due to off-design incidence values at its leading edge, an effect that occurs throughout the annulus despite the localized nature of the inlet distortion. Increased loss in the stator row is also observed due to flow separations that are shown to be caused by whirl angle distortion at rotor exit. By addressing these losses, it should be possible to achieve improved efficiency in BLI fan systems. © 2013 by ASME.