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.

A canonical normal SBLI flow relevant to external compression inlets

The normal shock wave / boundary layer interaction (normal SBLI) is important to the operation and performance of a supersonic inlet, and the normal SBLI is particularly prominent in external compression inlets. To improve our understanding of such interactions, it is helpful to make use of fundamental flows which capture the main elements of inlets, without resorting to the level of complexity and system integration associated with full-geometry inlets. In this paper, several fundamental fiow-fleld configurations have been considered as possible test cases to represent the normal SBLI aspects found in typical external compression inlets, and it was found that the spillage-diffuser more closely retains the basic flow features of an external compression inlet than the other configurations. In particular, this flow-fleld allows the normal shock Mach number as well as the amount and rate of subsonic diffusion to be all held approximately constant mid independent of the application of flow control. In addition, a survey of several external compression inlets was conducted to quantify the flow and geometric parameters of the spillage-diffuser relevant to actual inlets. The results indicated that such a flow may be especially relevant if the terminal Mach number is about 1.3 to 1.4, the confinement parameter is around 10%, the width around twice or three times the height, and with the area expansion just downstream of the shock on the conservative side of the stall limit for incompressible diffusers. © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

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.

Joint experimental and numerical approach to three-dimensional shock control bump research

Previous studies of transonic shock control bumps have often been either numerical or experimental. Comparisons between the two have been hampered by the limitations of either approach. The present work aims to bridge the gap between computational fluid dynamics and experiment by planning a joint approach from the outset. This enables high-quality validation data to be produced and ensures that the conclusions of either aspect of the study are directly relevant to the application. Experiments conducted with bumps mounted on the floor of a blowdown tunnel were modified to include an additional postshock adverse pressure gradient through the use of a diffuser as well as introducing boundary-layer suction ahead of the test section to enable the in-flow boundary layer to be manipulated. This has the advantage of being an inexpensive and highly repeatable method. Computations were performed on a standard airfoil model, with the flight conditions as free parameters. The experimental and computational setups were then tuned to produce baseline conditions that agree well, enabling confidence that the experimental conclusions are relevant. The methods are then applied to two different shock control bumps: a smoothly contoured bump, representative of previous studies, and a novel extended geometry featuring a continuously widening tail, which spans the wind-tunnel width at the rear of the bump. Comparison between the computational and experimental results for the contour bump showed good agreement both with respect to the flow structures and quantitative analysis of the boundary-layer parameters. It was seen that combining the experimental and numerical data could provide valuable insight into the flow physics, which would not generally be possible for a one-sided approach. The experiments and computational fluid dynamics were also seen to agree well for the extended bump geometry, providing evidence that, even though thebumpinteracts directly with the wind-tunnel walls, it was still possible to observe the key flow physics. The joint approach is thus suitable even for wider bump geometries. Copyright © 2013 by S. P. Colliss, H. Babinsky, K. Nubler, and T. Lutz. Published by the American Institute of Aeronautics and Astronautics, Inc.

CO2-gasification of a lignite coal in the presence of an iron-based oxygen carrier for chemical-looping combustion

Chemical-looping combustion (CLC) has the inherent property of separating the product CO2 from flue gases. Instead of air, it uses an oxygen carrier, usually in the form of a metal oxide, to provide oxygen for combustion. All techniques so far proposed for chemical looping with solid fuels involve initially the gasification of the solid fuel in order for the gaseous products to react with the oxygen carrier. Here, the rates of gasification of coal were compared when gasification was undertaken in a fluidised bed of either (i) an active Fe-based oxygen carrier used for chemical looping or (ii) inert sand. This enabled an examination of the ability of chemical looping materials to enhance the rate of gasification of solid fuels. Batch gasification and chemical-looping combustion experiments with a German lignite and its char are reported, using an electrically-heated fluidised bed reactor at temperatures from 1073 to 1223 K. The fluidising gas was CO2 in nitrogen. The kinetics of the gasification were found to be significantly faster in the presence of the oxygen carrier, especially at temperatures above 1123 K. A numerical model was developed to account for external and internal mass transfer and for the effect of the looping agent. The model also included the effects of the evolution of the pore structure at different conversions. The presence of Fe2O3 led to an increase in the rate of gasification because of the rapid oxidation of CO by the oxygen carrier to CO2. This resulted in the removal of CO and maintained a higher mole fraction of CO2 in the mixture of gas around the particle of char, i.e. within the mass transfer boundary layer surrounding the particle. This effect was most prominent at about 20% conversion when (i) the surface area for reaction was at its maximum and (ii) because of the accompanying increase in porosity and pore size, intraparticle resistance to gas mass transfer within the particle of char had fallen, compared with that in the initial particle. Excellent agreement was observed between the rates predicted by the numerical model and those observed experimentally. ©2013 Elsevier Ltd. All rights reserved.

Global modes with multiple saddle points

Significant progress has been made towards understanding the global stability of slowly-developing shear flows. The WKBJ theory developed by Patrick Huerre and his co-authors has proved absolutely central, with the result that both the linear and the nonlinear stability of a wide range of flows can now be understood in terms of their local absolute/convective instability properties. In many situations, the local absolute frequency possesses a single dominant saddle point in complex X-space (where X is the slow streamwise coordinate of the base flow), which then acts as a single wavemaker driving the entire global linear dynamics. In this paper we consider the more complicated case in which multiple saddles may act as the wavemaker for different values of some control parameter. We derive a frequency selection criterion in the general case, which is then validated against numerical results for the linearized third-order Ginzburg-Landau equation (which possesses two saddle points). We believe that this theory may be relevant to a number of flows, including the boundary layer on a rotating disk and the eccentric Taylor-Couette-Poiseuille flow. © 2014 Elsevier Masson SAS. All rights reserved.

On the laminar-turbulent transition of the rotating-disk flow: The role of absolute instability

© 2014 Cambridge University Press. This paper describes a detailed experimental study using hot-wire anemometry of the laminar-turbulent transition region of a rotating-disk boundary-layer flow without any imposed excitation of the boundary layer. The measured data are separated into stationary and unsteady disturbance fields in order to elaborate on the roles that the stationary and the travelling modes have in the transition process. We show the onset of nonlinearity consistently at Reynolds numbers, R, of ∼ 510, i.e. at the onset of Lingwood's (J. Fluid Mech., vol. 299, 1995, pp. 17-33) local absolute instability, and the growth of stationary vortices saturates at a Reynolds number of ∼ 550. The nonlinear saturation and subsequent turbulent breakdown of individual stationary vortices independently of their amplitudes, which vary azimuthally, seem to be determined by well-defined Reynolds numbers. We identify unstable travelling disturbances in our power spectra, which continue to grow, saturating at around R=585, whereupon turbulent breakdown of the boundary layer ensues. The nonlinear saturation amplitude of the total disturbance field is approximately constant for all considered cases, i.e. different rotation rates and edge Reynolds numbers. We also identify a travelling secondary instability. Our results suggest that it is the travelling disturbances that are fundamentally important to the transition to turbulence for a clean disk, rather than the stationary vortices. Here, the results appear to show a primary nonlinear steep-fronted (travelling) global mode at the boundary between the local convectively and absolutely unstable regions, which develops nonlinearly interacting with the stationary vortices and which saturates and is unstable to a secondary instability. This leads to a rapid transition to turbulence outward of the primary front from approximately R=565 to 590 and to a fully turbulent boundary layer above 650.