Generation and scattering of sound in vortical flow

We present solutions to scattering problems for unsteady disturbances to a mean swirling flow in an annular duct with a rigid 'splitter'. This situation has application to rotor-stator interaction noise in aeroengines, where the flow downstream of the fan is swirling and bifurcates into the by-pass duct and the engine core. We also consider the trailing edge extension of this problem. Inviscid mean flow in a cylindrical annulus is considered, with both axial and swirling (azimuthal) velocity components. The presence of vorticity in the mean flow couples the acoustic and vorticity modes of irrotational flow. Instead we have one combined spectrum of acoustic-vorticity waves in which the 'sonic' and 'nearly-convected' modes are fully coupled. In addition to the aeroacoustics application the results offer insight into the behaviour of these acoustic-vorticity waves, and the precise nature of the coupling between the two types of mode. Two regimes are discussed in which progress has been made, one for a specialised mean flow, uniform axial flow and rigid body swirl, and a second regime in which the frequency is assumed large, valid for any axisymmetric mean flow. The Wiener-Hopf technique is used to solve the scattering problems mathematically, and we present numerical evaluations of these solutions. Several new effects are seen to arise due to the mean vorticity, in particular the generation of sound at a trailing edge due to the scattering of a nearly convected disturbance, in contrast to the way a convected gust silently passes a trailing edge in uniform mean flow.

Normal shock boundary layer control with various vortex generator geometries

Various vortex generators which include ramp, split-ramp and a new hybrid concept "ramped-vane" are investigated under normal shock conditions with a diffuser at Mach number of 1.3. The dimensions of the computational domain were designed using Reynolds Average Navier-Stokes studies to be representative of the flow in an external-compression supersonic inlet. Using this flow geometry, various vortex generator concepts were studied with Implicit Large Eddy Simulation. In general, the ramped-vane provided increased vorticity compared to the other devices and reduced the separation length downstream of the device centerline. In addition, the size, edge gap and streamwise position respect to the shock were studied for the ramped-vane and it was found that a height of about half the boundary thickness and a large trailing edge gap yielded a fully attached flow downstream of the device. This ramped-vane also provided the largest reduction in the turbulent kinetic energy and pressure fluctuations. Additional benefits include negligible drag while the reductions in boundary layer displacement thickness and shape factor were seen compared to other devices. © 2010 by Sang Lee.

Normal shock boundary layer control with various vortex generator geometries

Various vortex generators which include ramp, split-ramp and a new hybrid concept "ramped-vane" are investigated under normal shock conditions with a diffuser at Mach number of 1.3. The dimensions of the computational domain were designed using Reynolds Average Navier-Stokes studies to be representative of the flow in an external-compression supersonic inlet. Using this flow geometry, various vortex generator concepts were studied with Implicit Large Eddy Simulation. In general, the ramped-vane provided increased vorticity compared to the other devices and reduced the separation length downstream of the device centerline. In addition, the size, edge gap and streamwise position respect to the shock were studied for the ramped-vane and it was found that a height of about half the boundary thickness and a large trailing edge gap yielded a fully attached flow downstream of the device. This ramped-vane also provided the largest reduction in the turbulent kinetic energy and pressure fluctuations. Additional benefits include negligible drag while the reductions in boundary layer displacement thickness and shape factor were seen compared to other devices. © 2011 Elsevier Ltd.

Low-order modelling of ducted flames with temporally varying equivalence ratio in realistic geometries

The interaction between unsteady heat release and acoustic pressure oscillations in gas turbines results in self-excited combustion oscillations which can potentially be strong enough to cause significant structural damage to the combustor. Correctly predicting the interaction of these processes, and anticipating the onset of these oscillations can be difficult. In recent years much research effort has focused on the response of premixed flames to velocity and equivalence ratio perturbations. In this paper, we develop a flame model based on the socalled G-Equation, which captures the kinematic evolution of the flame surfaces, under the assumptions of axisymmetry, and ignoring vorticity and compressibility. This builds on previous work by Dowling [1], Schuller et al. [2], Cho & Lieuwen [3], among many others, and extends the model to a realistic geometry, with two intersecting flame surfaces within a non-uniform velocity field. The inputs to the model are the free-stream velocity perturbations, and the associated equivalence ratio perturbations. The model also proposes a time-delay calculation wherein the time delay for the fuel convection varies both spatially and temporally. The flame response from this model was compared with experiments conducted by Balachandran [4, 5], and found to show promising agreement with experimental forced case. To address the primary industrial interest of predicting self-excited limit cycles, the model has then been linked with an acoustic network model to simulate the closed-loop interaction between the combustion and acoustic processes. This has been done both linearly and nonlinearly. The nonlinear analysis is achieved by applying a describing function analysis in the frequency domain to predict the limit cycle, and also through a time domain simulation. In the latter case, the acoustic field is assumed to remain linear, with the nonlinearity in the response of the combustion to flow and equivalence ratio perturbations. A transfer function from unsteady heat release to unsteady pressure is obtained from a linear acoustic network model, and the corresponding Green function is used to provide the input to the flame model as it evolves in the time domain. The predicted unstable frequency and limit cycle are in good agreement with experiment, demonstrating the potential of this approach to predict instabilities, and as a test bench for developing control strategies. Copyright © 2011 by ASME.

Geometry and interaction of structures in homogeneous isotropic turbulence

A strategy to extract turbulence structures from direct numerical simulation (DNS) data is described along with a systematic analysis of geometry and spatial distribution of the educed structures. A DNS dataset of decaying homogeneous isotropic turbulence at Reynolds number Reλ = 141 is considered. A bandpass filtering procedure is shown to be effective in extracting enstrophy and dissipation structures with their smallest scales matching the filter width, L. The geometry of these educed structures is characterized and classified through the use of two non-dimensional quantities, planarity' and filamentarity', obtained using the Minkowski functionals. The planarity increases gradually by a small amount as L is decreased, and its narrow variation suggests a nearly circular cross-section for the educed structures. The filamentarity increases significantly as L decreases demonstrating that the educed structures become progressively more tubular. An analysis of the preferential alignment between the filtered strain and vorticity fields reveals that vortical structures of a given scale L are most likely to align with the largest extensional strain at a scale 3-5 times larger than L. This is consistent with the classical energy cascade picture, in which vortices of a given scale are stretched by and absorb energy from structures of a somewhat larger scale. The spatial distribution of the educed structures shows that the enstrophy structures at the 5η scale (where η is the Kolmogorov scale) are more concentrated near the ones that are 3-5 times larger, which gives further support to the classical picture. Finally, it is shown by analysing the volume fraction of the educed enstrophy structures that there is a tendency for them to cluster around a larger structure or clusters of larger structures. Copyright © 2012 Cambridge University Press.

Effects of Boundary Layer Ingestion on the Aero-Acoustics of Transonic Fan Rotors

The use of boundary-layer-ingesting, embedded propulsion systems can result in inlet flow distortions where the interaction of the boundary layer vorticity and the inlet lip causes horseshoe vortex formation and the ingestion of streamwise vortices into the inlet. A previously-developed body-force-based fan modeling approach was used to assess the change in fan rotor shock noise generation and propagation in a boundary-layer-ingesting, serpentine inlet. This approach is employed here in a parametric study to assess the effects of inlet geometry parameters (offset-to-diameter ratio and downstream-to-upstream area ratio) on flow distortion and rotor shock noise. Mechanisms related to the vortical inlet structures were found to govern changes in the rotor shock noise generation and propagation. The vortex whose circulation is in the opposite direction to the fan rotation (counter-swirling vortex) increases incidence angles on the fan blades near the tip, enhancing noise generation. The vortex with circulation in the direction of fan rotation (co-swirling vortex) creates a region of subsonic relative flow near the blade tip radius which decreases the sound power propagated to the far-field. The parametric study revealed that the overall sound power level at the fan leading edge is set by the ingested streamwise circulation, and that for inlet designs in which the streamwise vortices are displaced away from the duct wall, the sound power at the upstream inlet plane increased by as much as 9 dB. By comparing the far-field noise results obtained to those for a conventional inlet, it is deduced that the changes in rotor shock noise are predominantly due to the ingestion of streamwise vorticity.

Effects of Boundary Layer Ingestion on the Aero-Acoustics of Transonic Fan Rotors

The use of boundary-layer-ingesting, embedded propulsion systems can result in inlet flow distortions where the interaction of the boundary layer vorticity and the inlet lip causes horseshoe vortex formation and the ingestion of streamwise vortices into the inlet. A previously-developed body-force-based fan modeling approach was used to assess the change in fan rotor shock noise generation and propagation in a boundary-layer-ingesting, serpentine inlet. This approach is employed here in a parametric study to assess the effects of inlet geometry parameters (offset-to-diameter ratio and downstream-to-upstream area ratio) on flow distortion and rotor shock noise. Mechanisms related to the vortical inlet structures were found to govern changes in the rotor shock noise generation and propagation. The vortex whose circulation is in the opposite direction to the fan rotation (counter-swirling vortex) increases incidence angles on the fan blades near the tip, enhancing noise generation. The vortex with circulation in the direction of fan rotation (co-swirling vortex) creates a region of subsonic relative flow near the blade tip radius which decreases the sound power propagated to the far-field. The parametric study revealed that the overall sound power level at the fan leading edge is set by the ingested streamwise circulation, and that for inlet designs in which the streamwise vortices are displaced away from the duct wall, the sound power at the upstream inlet plane increased by as much as 9 dB. By comparing the far-field noise results obtained to those for a conventional inlet, it is deduced that the changes in rotor shock noise are predominantly due to the ingestion of streamwise vorticity.

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.

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.

Flow physics of a normal-hole bled supersonic turbulent boundary layer

An experimental study on normal hole bleed in a supersonic turbulent boundary layer has been conducted. A combination of LDV, Schlieren imagery and oil flow visualization were used to provide a better understanding of the three-dimensional flow field surrounding a supersonic bleed array. Experiments were performed at Mach numbers of 1.8 and 2.5, while previously published results at Mach numbers of 1.3 and 1.5 were also incorporated. The bleed system was capable of removing up to approximately 10% of the incoming boundary layer through a tunnel-spanning array of discrete holes with diameters the same order of magnitude of boundary layer displacement thickness. Inspection of boundary layer profiles downstream of the bleed region indicates that vorticity generated by the discrete holes can have a substantial influence on changes to the boundary layer shape factor and skin friction coefficient, through modification of the lower 20% of the boundary layer. This vorticity was visualized through oil-flow visualization, and LDV measurements, showing the development of two vortices off each bleed hole, and corresponding upwash and downwash regions with far-reaching three dimensional effects. © 2013 by J. M. Oorebeek and H. Babinsky.

Wave propagation in a slowly-varying duct with mean vortical swirling flow

The propagation of unsteady disturbances in a slowlyvarying cylindrical duct carrying mean swirling flow is investigated using a multiple-scales technique. This is applicable to turbomachinery flow behind a rotor stage when the swirl and axial velocities are of the same order. The presence of mean vorticity couples acoustic and vorticity equations which produces an eigenvalue problem that is not self-adjoint unlike that for irrotational mean flow. In order to determine the amplitude variation along the duct, an adjoint solution for the coupled system of equations is derived. The solution breaks down where a mode changes from cut on to cut off. In this region the amplitude is governed by a form of Airy's equation, and the effect of swirl is to introduce a small shift in the origin of the Airy function away from the turning-point location. The variation of axial wavenumber and amplitude along the duct is calculated. In hard-walled ducts mean swirl is shown to produce much larger amplitude variation along the duct compared with a nonswirling flow. Mean swirl also has a large effect in ducts with finite-impedance walls which differs depending on whether modes are co-rotating with the swirl or counter rotating. © 2001 by A.J. Cooper, Published by the American Institute of Aeronautics and Astronautics, Inc.