Shallow cavity flow tone experiments: onset of locked-on states

Fully turbulent inflow past a shallow cavity is investigated for the configuration of an axisymmetric cavity mounted in a pipe. Emphasis is on conditions giving rise to coherent oscillations, which can lead to locked-on states of flow tones in the pipe–cavity system. Unsteady surface pressure measurements are interpreted using three-dimensional representations of amplitude–frequency, and velocity; these representations are constructed for a range of cavity depth. Assessment of these data involves a variety of approaches. Evaluation of pressure gradients on plan views of the three-dimensional representations allows extraction of the frequencies of the instability (Strouhal) modes of the cavity oscillation. These frequency components are correlated with traditional models originally formulated for cavities in a free-stream. In addition, they are normalized using two length scales: inflow boundary-layer thickness and pipe diameter. These scales are consistent with those employed for the hydrodynamic instability of the separated shear layer, and are linked to the large-scale mode of the shear layer oscillation, which occurs at relatively long cavity length. In fact, a simple scaling based on pipe diameter can correlate the frequencies of the dominant peaks over a range of cavity depth.The foregoing considerations provide evidence that pronounced flow tones can be generated from a fully turbulent inflow at very low Mach number, including the limiting case of fully developed turbulent flow in a pipe. These tones can arise even for the extreme case of a cavity having a length over an order of magnitude longer than its depth. Suppression of tones is generally achieved if the cavity is sufficiently shallow.     

The onset of turbulence in physiological pulsatile flow in a straight tube

 An empirical correlation for the onset of turbulence in physiological pulsatile flow is presented. We pumped three different test fluids of kinematic viscosity 0.008–0.035 cm²/s through four straight tubes 0.4–3.0 cm in diameter. A Scotch yoke mechanism provided an oscillatory sine wave flow component of known stroke volume and frequency. We adjusted the mean flow independently until we detected signal instabilities from hot film wall shear stress probes. The critical peak Reynolds number was found to correlate with the Womersley parameter and the Strouhal number as a power law function with a root-mean-square (rms) error of 15.2%. Experimental measurements of the laminar velocity profile are compared to theoretical predictions from Poiseuille’s law and Womersley’s solution. Received: 30 October 1995/Accepted: 7 April 1997     

Effects of sound on flow separation from blunt flat plates

The effect of sound on the flow around plates with semicircular or square leading edges and square trailing edges located in a low turbulence open jet has been studied. In all circumstances the length of the leading edge separation bubbles associated with square leading edge plates was found to oscillate. When sound was applied to the flow around these plates, the leading edge shear layers reattached closer to the leading edge and the oscillations in bubble length occurred at the applied sound frequency, generating patches of concentrated vorticity in the boundary layers. These vorticity patches moved downstream near the plate surface and then beyond the trailing edge to form vortex cores in a street with a Strouhal number equal to the applied sound value. Sometimes these vortex streets are unstable and break down into streets with Strouhal numbers approaching those observed without sound. These effects of sound were not observed in the flow around plates with semicircular leading edges. Without sound, square leading edge plates of intermediate length did not shed regular vortex streets.     

Numerical analysis of characteristic features of shallow and deep cavity in supersonic flow

Flow past open cavities are numerically simulated at a Mach number of 1.5, and Reynolds number, based on initial momentum thickness at the front lip of cavity, of 3333 for variable depths (D) with constant length (L). The dominant frequency of oscillation shows a sudden jump when there is a transition from shallow (L/D > 1) to deep cavity (L/D < 1). The vorticity thickness displays two different growth rates along the length of cavity: (1) initial lower spreading rate, followed by (2) higher spreading rate. The lower spreading rate of shear layer is dictated by the type of cavity (either shallow or deep), while the higher spreading rate is directly related to the amplitude of oscillations. Proper orthogonal decomposition (POD) is implemented to visualise the coherent structures based on their energy content. The first two POD spatial structures in the shallow cavity represent vortex shedding, while in the deep cavity, they comprise vortex pairing interactions as in mixing layer. The higher POD modes contain coherent structures at mixed frequencies. The behaviour of coherent structures associated with a temporal frequency is further investigated using dynamic mode decomposition (DMD). The higher DMD modes confirm the dominance of mixing layer behaviour in the deep cavity.     

Instabilités de l'écoulement produisant le bruit de fenteFlow instabilities producing the slot-tone

The nature of the instability governing the self-sustained tones produced by a low Mach number plane jet impinging on a slotted plate, known as slot-tone, is identified experimentally. For a given Reynolds number, the natural shear-layer and the jet column mode frequencies of the free jet delimit the values of the measured slot-tone operating frequencies. The oscillations at lower frequencies are the result of the amplification of the jet column mode, and those at higher frequencies correspond to the shear layer instabilities. To cite this article: A. Billon et al., C. R. Mecanique 332 (2004).     

Numerical simulation of shear effect on vortex shedding behind a square cylinder

This paper is concerned with the numerical simulation of the flow structure around a square cylinder in a uniform shear flow. The calculations were conducted by solving the unsteady 2D Navier–Stokes equations with a finite difference method. The effect of the shear parameter K of the approaching flow on the vortex-shedding Strouhal number and the force coefficients acting on the square cylinder is investigated in the range K=0·0–0·25 at various Reynolds numbers from 500 to 1500. The computational results are compared with some existing experimental data and previous studies. The effect of shear rate on the Strouhal number and the force acting on the cylinder has a tendency to reduce the oscillation. The Strouhal number, mean drag and amplitude of the fluctuating force tend to decrease as the shear rate increases, but show no significant change at low shear rate. Increasing the Reynolds number decreases the Strouhal number and increases the force acting on the cylinder. At high shear rate the shedding frequencies of the fluctuating drag and lift coefficients are identical. © 1997 John Wiley & Sons, Ltd.     

Numerical Simulations of High-Speed Turbulent Cavity Flows

Detached-Eddy Simulations (DES) of flows over clean and controlled cavities with and without doors are presented in this paper. Mach and Reynolds numbers (based on cavity length) were 0.85 and one million respectively. Spectral analyses showed that the DES computations were able to correctly predict the frequencies of the Rossiter modes for both uncontrolled and controlled cases. Flow visualisations revealed that the impact of the shear layer formed along the cavity on a slanted aft wall no longer creates a large source of acoustic noise. Therefore little acoustic propagation was seen up the cavity. This was confirmed by the analysis of the cavity wall forces, which showed that the oscillations of the shear layer were reduced when the wall was slanted. This aided in reducing the overall Sound Pressure Levels throughout the cavity and far-field. Comparisons of the flow-fields suggested that the addition of the doors also aided in stabilising the shear layer, which was also shown in the analysis of the wall forces. As a result, the addition of the doors was found to affect the clean cavity configuration significantly more than the controlled one.     

Experimental investigation of cavitating structures in the near wake of a cylinder

An experimental investigation of cavitating structures in the near-wake region of a cylinder is presented. From high-speed imaging of this subcritical flow (Reynolds number of 64,000), it is found that inception of cavities occurs in the shear layer. At the developed cavitation condition, the cavities in the separated zone and the free shear layer merge. A distinct spanwise variation in cavitation activity is observed. The non-dimensionalized correlation length at inception varies from close to a non-cavitating value of about 3.5 to about 1 at developed cavitation. The non-dimensionalized length of formation, characterized by crossover of the free shear layer and the wake axis, increases from 1 to 1.8 as the cavitation number is reduced from 85% to 50% of the inception value. A frequency analysis of the cavity dynamics indicates that although the vortex shedding frequency is dominant in the shear layer, there are peaks corresponding to other frequencies in other flow regions. The presence of a sharp peak at 125 Hz, corresponding to a Strouhal number of 0.2, along with a range of frequencies, is also verified independently through measurement of fluctuating pressure at the cylinder surface.     

Revisiting URANS Computations of the Backward-facing Step Flow Using Second Moment Closures. Influence of the Numerics

Previous studies have shown that Unsteady Reynolds-Averaged Navier–Stokes (URANS) computations are able to reproduce the vortex shedding behind a backward-facing step. The aim of the present work is to investigate not only the quantitative predictions of the URANS methodology concerning the characteristic frequencies, but also the amplitude of the energy of the resolved eddies, by using the Elliptic Blending Reynolds Stress Model. This innovative low-Reynolds number second moment closure reproduces the non-viscous, non-local blocking effect of the wall on the Reynolds stresses, and it is compared to the standard k − ε and LRR models using wall-functions. Consistent with previous studies, in the 2D computations shown in the present article, the vortex shedding is captured with the correct Strouhal number, when second moment closures are used. To complete these previous analyses, we particularly focus here on the energy contained in the unsteady, resolved part and its dependency on the numerical method. This energy is less than 5% of the total energy and is strongly dependent on the mesh. Using a refined mesh, surprisingly, a steady solution is obtained. It is shown that this behaviour can be linked to the very small spatial oscillations at the step corner, produced by numerical dispersion, which act as perturbations that are sufficient to excite the natural mode of the shear layer, when the local Peclet number, comparing convection and diffusion effects, is high enough. This result suggests that URANS is not appropriate to quantitatively predict the amplitude of the large-scale structures developing in separated shear-layers, and that URANS results must be interpreted with care in terms of temporal variations of forces, temperatures, etc., in industrial applications using marginally fine meshes.     

Characteristics of Oscillations in Supersonic Open Cavity Flows

Characteristics of oscillations in supersonic open cavity flows are investigated numerically using hybrid RANS/LES (Reynolds-Averaged Navier-Stokes/Large Eddy Simulation) method. The oscillation regimes and feedback mechanisms for the supersonic cavity flows are identified and analyzed. The calculation captures a mixed shear-layer/wake oscillation mode in the flow of Ma = 1.75, where these two modes occur alternately. The shear-layer mode and wake mode are driven by vortex convection-acoustic feedback and absolute instability, respectively. In particular, the results indicate that the feedback-acoustic-wave in the shear-layer mode is probably generated by the reflection of the downstream-traveling pressure wave, associated with the shed vortex in the shear layer, on the aft wall. The cavity flow of Ma = 2.52 is then simulated to see the influence of Mach number. It is found that the increase of Mach number may decrease the amplitude of the fluctuations in the shear layer, inhibiting the transition to wake mode. Furthermore, the influence of upstream injection is also studied, where the results show that the injection only weakens the oscillations and faintly shifts the resonant frequencies.     

Time resolved flow-field measurements of a turbulent mixing layer over a rectangular cavity

High Reynolds number, low Mach number, turbulent shear flow past a rectangular, shallow cavity has been experimentally investigated with the use of dual-camera cinematographic particle image velocimetry (CPIV). The CPIV had a 3 kHz sampling rate, which was sufficient to monitor the time evolution of large-scale vortices as they formed, evolved downstream and impinged on the downstream cavity wall. The time-averaged flow properties (velocity and vorticity fields, streamwise velocity profiles and momentum and vorticity thickness) were in agreement with previous cavity flow studies under similar operating conditions. The time-resolved results show that the separated shear layer quickly rolled-up and formed eddies immediately downstream of the separation point. The vortices convect downstream at approximately half the free-stream speed. Vorticity strength intermittency as the structures approach the downstream edge suggests an increase in the three-dimensionality of the flow. Time-resolved correlations reveal that the in-plane coherence of the vortices decays within 2–3 structure diameters, and quasi-periodic flow features are present with a vortex passage frequency of ~1 kHz. The power spectra of the vertical velocity fluctuations within the shear layer revealed a peak at a non-dimensional frequency corresponding to that predicted using linear, inviscid instability theory.     

Wall Effect on the Convective–Absolute Boundary for the Compressible Shear Layer

The linear stability of inviscid compressible shear layers is studied. When the layer develops at the vicinity of a wall, the two parallel flows can have a velocity of the same sign or of opposite signs. This situation is examined in order to obtain first hints on the stability of separated flows in the compressible regime. The shear layer is described by a hyperbolic tangent profile for the velocity component and the Crocco relation for the temperature profile. Gravity effects and the superficial tension are neglected. By examining the temporal growth rate at the saddle point in the wave-number space, the flow is characterized as being either absolutely unstable or convectively unstable. This study principally shows the effect of the wall on the convective–absolute transition in compressible shear flow. Results are presented, showing the amount of the backflow necessary to have this type of transition for a range of primary flow Mach numbers M ₁ up to 3.0. The boundary of the convective–absolute transition is defined as a function of the velocity ratio, the temperature ratio and the Mach number. Unstable solutions are calculated for both streamwise and oblique disturbances in the shear layer. Received 9 May 2001 and accepted 21 August 2001     

Aeolian tones of a plate cascade

A model for the aeroacoustic resonance effects (aeolian tones) excited around a plate cascade in a gas flow is suggested. Methods of calculating the frequencies of natural acoustic oscillations near the cascade are developed. The effect of the cascade geometry and the Mach number of the main flow on the frequencies, abundance, and modes of the natural oscillations is investigated. Anomalous acoustic oscillations near a cyclic plate cascade are shown to exist and are studied. It is shown that there always exist no less than two natural oscillation frequencies in the gas flow near any nontrivial cyclic plate cascade. It has been found that the natural oscillation frequencies can be combined in bundles such that in the case where the number of plates in a period is large the frequencies pertaining to each bundle occupy a certain interval with arbitrary density. The natural oscillations are classified with respect to the form of the eigenfunctions; the classification is based on the theory of representations of groups of locally plane symmetries of the cyclic plate cascade in the solution space. The correctness of the proposed model of the aeroacoustic resonance effects (aeolian tones) excited near a plate cascade in a gas flow is supported by a comparison with the available experimental and theoretical data. On the basis of the investigation performed, some previously unknown physical phenomena are predicted. Thus, the existence of frequency zones or main-flow Mach number ranges on which aeroacoustic resonance phenomena exist near a cyclic cascade with a large number of plates in a period is proved; it is shown that for certain frequencies of the natural oscillations near the cyclic plate cascade the resonance oscillations may be localized in the vicinity of the source; and the existence of narrow-band wave packets slowly propagating along the cascade is demonstrated. Novosibirsk, e-mail: sukhinin@hydro.nsc.ru. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 171–186, March–April, 2000.     

Mode transition and oscillation suppression in supersonic cavity flow

Supersonic flows past two-dimensional cavities with/without control are investigated by the direct numerical simulation(DNS). For an uncontrolled cavity, as the thickness of the boundary layer declines, transition of the dominant mode from the steady mode to the Rossiter Ⅱ mode and then to the Rossiter Ⅲ mode is observed due to the change of vortex-corner interactions. Meanwhile, a low frequency mode appears. However, the wake mode observed in a subsonic cavity flow is absent in the current simulation.The oscillation frequencies obtained from a global dynamic mode decomposition(DMD)approach are consistent with the local power spectral density(PSD) analysis. The dominant mode transition is clearly shown by the dynamic modes obtained from the DMD. A passive control technique of substituting the cavity trailing edge with a quarter-circle is studied. As the effective cavity length increases, the dominant mode transition from the Rossiter Ⅱ mode to the Rossiter Ⅲ mode occurs. With the control, the pressure oscillations are reduced significantly. The interaction of the shear layer and the recirculation zone is greatly weakened, combined with weaker shear layer instability, responsible for the suppression of pressure oscillations. Moreover, active control using steady subsonic mass injection upstream of a cavity leading edge can stabilize the flow.     

Mixing Characteristics in a Supersonic Combustor with Gaseous Fuel Injection Upstream of a Cavity Flameholder

Non-reacting experiments and numerical simulations have been performed to investigate the mixing characteristics in a supersonic combustor with gaseous fuel injection upstream of a flameholding cavity in a supersonic vitiated air flow with stream Mach number 1.7. Using helium as simulated fuel, the acetone vapor is adulterated into the fuel jet. The fuel distribution in spanwise and streamwise direction is imaged by the planar laser-induced fluorescence (PLIF) measurement. According to the similarity of experimental observations with different cavities, the typical L/D = 7 cavity with aft wall angle 45° is chosen and the flowfield with the injection is calculated by Large Eddy Simulation. Experimental and numerical results have shown that most of the fuel flow away upon the open cavity with the lifting counter-rotating vortex structures induced by the transverse jet. Only a small portion of the fuel is convected into the cavity shear layer by the vortex interaction of the jet with cavity shear layer, and then transported into the cavity due to the cavity shear layer motion and the interaction of the shear layer with the cavity trailing edge.     

LES and DES Investigations of Turbulent Flow over a Sphere at Re = 10,000

Large Eddy Simulation (LES) using a dynamic Smagorinsky type subgridstress (SGS) model and Detached Eddy Simulation (DES) are applied toprediction and investigation of the flow around a sphere at a Reynoldsnumber of 10⁴ in the subcritical regime. In this regime the boundarylayers at separation are laminar, and transition to turbulence occursfarther downstream in the separated shear layers via Kelvin–Helmholtz(K–H) instabilities. The dynamic eddy viscosity model of Germano et al.(Physics of Fluids 3 (1991) 1760–1765) is used in the LES, while the current implementation of the DESemploys a formulation based on the Spalart–Allmaras (S–A) model. DES isa hybrid approach in which the closure is a modification to theproduction/destruction term of the original Reynolds-AveragedNavier–Stokes (RANS) model, reducing to RANS in the attached regions,and to LES away from the wall. In the present work where we simulate theflow over a sphere in the subcritical regime in which the boundarylayers at separation are laminar, DES can be viewed as LES with adifferent SGS model. Effects of the discretization scheme used toapproximate the convective terms are considered, along with sensitivityof predictions to changes in the additional model coefficient, C _DES, in the DES formulation. DES and LES yield similar predictions of the wakestructure, large-scale vortex shedding and the Strouhal numberassociated with the low frequency mode in the wake. Predictions ofquantities such as the drag coefficient, wake frequencies, position oflaminar separation on the sphere, and the mean pressure andskin-friction distributions along the sphere are in good agreement withthe measurements of Achenbach (Journal of Fluid Mechanics 54 (1972) 565–575). Predictions of the primaryReynolds shear stress, turbulent kinetic energy, eddy viscosity, andturbulent dissipation for the two models are also similar. In addition,both models successfully resolve the formation of the vortex tubes inthe detached shear layers along with the value of the Strouhal numberassociated with the high frequency instability mode, provided that thelevel of numerical dissipation introduced by the discretization schemeremains sufficiently low. Flow physics investigations are focused onunderstanding the wake structure in the subcritical regime.     

Measurement and analysis of the turbulent length scales in jet flows

The characteristic length and time scales of turbulence are reported in some detail for jet flows. The objective of the work is to determine the frequency dependence of these two-point turbulent properties, which are used to model the sources necessary for noise prediction using the acoustic analogy approach. A range of jet flow conditions for single and co-axial configurations are considered so that the effect of Mach number, temperature ratio and nozzle geometry is examined. The frequency dependence of both the fixed and moving frame length scales and the convection velocities for both the turbulence and the Reynolds stress are derived using a two-point complex coherence function. At higher frequencies, the integral scales are found to be strongly isotropic and inversely proportional to the Strouhal number. A frequency-dependent Taylor scale is derived and shown to agree well with the experimental results at the higher frequencies.     

A DPIV study of a starting flow downstream of a backward-facing step

 In this paper an experimental investigation of a starting water flow downstream of a backward-facing step is described. The Reynolds number of the asymptotic steady flow is Re≈4300 based on the step height of s=2 cm and the free stream velocity of U=21.4 cm/s. Velocity measurements were performed with video-based DPIV (Digital Particle Image Velocimetry) at a sampling frequency of 25 Hz. The main purpose of this study is to reveal the temporal development of global structures which could not be analyzed with single-point probes. It was found that at initialization of the flow a regular vorticity street is formed, which collapses at a normalized time t ^* =U t/s≈17 due to vorticity interactions. After this the flow is dominated by complicated vorticity roll-up and shedding dynamics in the recirculation region. The starting phase is terminated for t ^* >40. Prior to the collapse of the vorticity street values of 9 times the steady state asymptotic wall normal stress and of twice the steady state negative wall shear stress were observed. The early increasing slope of the reattachment length is constant over a time of approximately t ^* =8. The collapse of the vorticity street and the vorticity interactions thereafter cause fluctuations both in the velocity in the free shear layer and in the reattachment length. The fully developed flow has a dominant frequency corresponding to a Strouhal number St=fs/U≈0.04. Received: 20 September 1996/Accepted: 1 April 1997     

An inclined wall jet: Mean flow characteristics and effects of acoustic excitation

 The mean velocity field of a 30° inclined wall jet has been investigated using both hot-wire and laser Doppler anemometry (LDA). Provided that the nozzle aspect ratio is greater than 30 and the inclined wall angle (β) is less than 50°, LDA measurements for various β show that the reattachment length is independent of the nozzle aspect ratio and the nozzle exit Reynolds number (in the range 6670–13,340). There is general agreement between the reattachment lengths determined by LDA and those determined using wall surface oil film visualisation technique. The role of coherent structures arising from initial instabilities of a 30° wall jet has been explored by hot-wire spectra measurements. Results indicate that the fundamental vortex roll-up frequency in both the inner and outer shear layer corresponds to a Strouhal number (based on nozzle exit momentum thickness and velocity) of 0.012. The spatial development of instabilities in the jet has been studied by introducing acoustic excitation at a frequency corresponding to the shear layer mode. The formation of the fundamental and its first subharmonic has been identified in the outer shear layer. However, the development of the first subharmonic in the inner shear layer has been severely suppressed. Distributions of mean velocities, turbulence intensities and Reynolds shear stress indicate that controlled acoustic excitation enhances the development of instabilities and promotes jet reattachment to the wall, resulting in a substantially reduced recirculation flow region. Received: 24 November 1998/Accepted: 24 August 1999     

Numerical analysis of the flow pattern and vortex breakdown over a pitching delta wing at supersonic speeds

A supersonic compressible flow over a 60° swept delta wing with a sharp leading edge undergoing pitching oscillations is computationally studied. Numerical simulations are performed by the finite volume method with the use of the k?ω turbulence model for various Mach numbers and angles of attack. Variations of flow patterns in a crossflow plane, hysteresis loops associated with the vortex core location, and vortex breakdown positions during a pitching cycle are investigated. Trends for various Mach numbers, mean angles of attack, pitching amplitudes, and pitching frequencies are illustrated.