Free-surface fluctuations in hydraulic jumps: Experimental observations

A hydraulic jump is the rapid and sudden transition from a high-velocity supercritical open channel flow to a subcritical flow. It is characterised by the dynamic interactions of the large-scale eddies with the free-surface. New series of experimental measurements were conducted in hydraulic jumps with Froude numbers between 3.1 and 8.5 to investigate these interactions. The dynamic free surface measurements were performed with a non-intrusive technique while the two-phase flow properties were recorded with a phase-detection probe. The shape of the mean free surface profile was well defined and the turbulent fluctuation profiles highlighted a distinct peak of turbulent intensity in the first part of the jump roller, with free-surface fluctuation levels increasing with increasing Froude number. The dominant free-surface fluctuation frequencies were typically between 1 and 4 Hz. A comparison between the acoustic sensor signals and conductivity probe data suggested that the air–water “free-surface” detected by the acoustic sensor corresponded to about the boundary between the turbulent shear layer and the upper free-surface layer. Simultaneous measurements of free surface and bubbly flow fluctuations for Fr = 5.1 indicated that the frequency ranges of both sensors were similar (F < 5 Hz) whatever the position downstream of the toe. The present results highlighted that the dynamic free-surface measurements can be conducted successfully using acoustic displacement meters, and the time-averaged depth measurements was a physical measure of the free-surface location in hydraulic jumps.     

Turbulent flow in a ribbed channel: Flow structures in the vicinity of a rib

PIV measurements are performed in a channel with periodic ribs on one wall. The emphasis of this study is to investigate the flow structures in the vicinity of a rib in terms of mean velocities, Reynolds stresses, probability density functions (PDF), and two-point correlations. The PDF distribution of u′ is bimodal in the separated shear layer downstream of the rib. The maximum Reynolds shear stresses occur at the leading edge of the rib. Based on quadrant analysis, it is found that ejection motions make a dominant contribution to the Reynolds shear stress in this region. Moreover, topology-based visualization is applied to the separation bubble upstream of the rib. Salient critical points and limit cycles are extracted, which gives clues to the physical processes occurring in the flow.     

Stability of a laminar rivulet liquid flow in a cylindrical duct in the approximation of one-dimensional waves

The problem of the stability of a viscous laminar liquid flow with a liquid free surface in an inclined duct is theoretically considered. Since the dependence of the flow rate on the free-surface height is not monotonic (the highest flow rate in a cylindrical duct is observed at H_*=1.7R), primary attention is given to the region H>H_*. It is proved that there is aw region of instability: for an arbitrarity low Reynolds number, there is a free-surface level above which the flow becomes unstable against one-dimensional disturbances. When the height of the liquid layer is close to the vertical dimension of the duct, the one-dimensional disturbances propagate mainly upstream (for moderate Reynolds numbers). Hence it follows that there is not steady regime of liquid flow from a fully filled duct with an open end. Kutateladze Institute of Thermal Physics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 3, pp. 90–96, May–June, 1999.     

Oscillation of shallow flow past a cavity: Resonant coupling with a gravity wave

Self-excited oscillations of flow past a cavity are generated in a shallow free-surface system. The shear layer past the cavity opening has two basic forms: a separated free-shear flow; and a shear flow along a slotted plate. Instabilities of these classes of shear flows can couple with the fundamental gravity-wave mode of the adjacent cavity. The dimensionless frequencies of both types of oscillations scale on the length of the cavity opening, rather than the gap distance between the slats, i.e., a large-scale instability is always prevalent. A technique of high-image-density particle image velocimetry allows acquisition and interpretation of global, instantaneous images of the flow pattern, including patterns of vorticity and Reynolds stress correlation. Use of a cinema approach provides representations of the timewise evolution of the global, instantaneous flow structure, and thereby definition of the amplitude peaks and phase angles of the coupled fluctuations via auto- and cross-spectral techniques. These methods, along with global, averaged representations of the fluctuating flow field, provide insight into the onset of fully coupled (phase-locked) oscillations of the shear flow past the resonator cavity. The common, as well as the distinctive, features of the resonant-coupled instability of the shear flow past the slotted plate are characterized, relative to the corresponding coupled instability of the free-shear layer. Varying degrees of resonant coupling between the unstable shear layer and the adjacent resonator are attained by variations of the inflow velocity, which yield changes of the predominant oscillation frequency, relative to the resonant frequency of the adjacent cavity. Well-defined, coherent oscillations are indeed attainable for the case of the shear flow along the slotted plate, though their amplitude is significantly mitigated relative to the case of a free-shear layer. The degree of organization of the self-excited, resonant-coupled oscillation and the manner in which it varies with open area ratio and geometry of the plate are interpreted in terms of the flow structure on either side of, and within, the slotted plate; these features are compared with the corresponding structure of the free-shear layer oscillations.     

Flow around a NACA0018 airfoil with a cavity and its dynamical response to acoustic forcing

Trapping of vortices in a cavity has been explored in recent years as a drag reduction measure for thick airfoils. If, however, trapping fails, then oscillation of the cavity flow may couple with elastic vibration modes of the airfoil. To examine this scenario, the effect of small amplitude vertical motion on the oscillation of the shear layer above the cavity is studied by acoustic forcing simulating a vertical translation of a modified NACA0018 profile. At low Reynolds numbers based on the chord (O(10⁴)), natural instability modes of this shear layer are observed for Strouhal numbers based on the cavity width of order unity. Acoustic forcing sufficiently close to the natural instability frequency induces a strong non-linear response due to lock-in of the shear layer. At higher Reynolds numbers (above 10⁵) for Strouhal number 0.6 or lower, no natural instabilities of the shear layer and only a linear response to forcing were observed. The dynamical pressure difference across the airfoil is then dominated by added mass effects, as was confirmed by numerical simulations.     

Turbulence production in flow over a wavy wall

Measurements of the spatial and time variation of two components of the velocity have been made over a sinusoidal solid wavy boundary with a height to length ratio of 2a/λ = 0.10 and with a dimensionless wave number of α⁺ = (2π/λ)(v/u ^?) = 0.02. For these conditions, both intermittent and time-mean flow reversals are observed near the troughs of the waves. Statistical quantities that are determined are the mean streamwise and normal velocities, the root-meansquare of the fluctuations of the streamwise and normal velocities, and the Reynolds shear stresses. Turbulence production is calculated from these measurements. The flow is characterized by an outer flow and by an inner flow extending to a distance of about α^?1 from the mean level of the surface. Turbulence production in the inner region is fundamentally different from flow over a flat surface in that it is mainly associated with a shear layer that separates from the back of the wave. Flow close to the surface is best described by an interaction between the shear layer and the wall, which produces a retarded zone and a boundary-layer with large wall shear stresses. Measurements of the outer flow compare favorably with measurements over a flat wall if velocities are made dimensionless by a friction velocity defined with a shear stress obtained by extrapolating measurements of the Reynolds stress to the mean levels of the surface (rather than from the drag on the wall).     

Bifurcation of the stationary Ekman flow into a stable periodic flow

The two cases of stationary Ekman boundary layer flow of an incompressible fluid near i) a plane boundary and ii) a free surface with constant shear are considered. It is proven that a stable secondary flow in the form of traveling waves bifurcates from the stationary flow at a certain Reynolds number, and that the stationary flow is unstable above this number. The values of the critical Reynolds number and of the numbers that characterize the traveling wave are computed and compared with experimental values.     

Reynolds stress modelling of rectangular open‐channel flow

A Reynolds stress model for the numerical simulation of uniform 3D turbulent open‐channel flows is described. The finite volume method is used for the numerical solution of the flow equations and transport equations of the Reynolds stress components. The overall solution strategy is the SIMPLER algorithm, and the power‐law scheme is used to discretize the convection and diffusion terms in the governing equations. The developed model is applied to a flow at a Reynolds number of 77000 in a rectangular channel with a width to depth ratio of 2. The simulated mean flow and turbulence structures are compared with measured and computed data from the literature. The computed flow vectors in the plane normal to the streamwise direction show a small vortex, called inner secondary currents, located at the juncture of the sidewall and the free surface as well as the free surface and bottom vortices. This small vortex causes a significant increase in the wall shear stress in the vicinity of the free surface. A budget analysis of the streamwise vorticity is carried out. It is found that both production terms by anisotropy of Reynolds normal stress and by Reynolds shear stress contribute to the generation of secondary currents. Copyright © 2006 John Wiley & Sons, Ltd.     

Finite difference solution of a Newtonian jet swell problem

A finite difference technique has been developed to study the Newtonian jet swell problem. The streamfunction and vorticity were used as dependent variables to describe the jet flow. The boundary-fitted co-ordinate transformation method was adopted to map the flow geometry into a rectangular domain. The standard finite difference method was then applied for solving the flow equations. The location of the jet free surface was updated by the kinematic boundary condition, and an adjustable parameter was included in the free-surface iteration. We could obtain numerical solutions for the Reynolds number as high as 100, and the differences between the present study and previous finite element simulations on the jet swell ratio are less than 5%.     

LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.     

Fence probe measurement of flow reversal in separating turbulent boundary layers

A small fence probe was evaluated for measurements in the time-dependent flow reversal region of the transition from boundary layer to separated flow. For moderate and high Reynolds numbers, the fence probe is demonstrated to be a usable tool for the measurement of the reverse flow associated with separation. Although the present probe pressure transducer system was limited to approximately 200?Hz, pulses of positive and negative shear stress were readily detected. At or near the location of zero surface shear stress, the measurements were limited by the signal-to-noise ratio. For the separated flow investigated, a marked reduction in the pressure gradient occurred when the fence probe indicated approximately 20?% reversal for the higher Reynolds numbers. The reversal increased to 24?% for the lower Reynolds numbers. The measurements indicate that flow reversal alone may not be adequate to identify the degree of separation. Upstream of turbulent boundary layer (intermittent) separation, the duration of the reversed shear stress was found to be very short (0.002?C0.007?s), suggesting a local, small-scale, impulse-type separation. At and beyond the location of intermittent separation, the shear stress reversal duration was an order of magnitude longer. Estimates of the maximum and minimum surface shear stress in the separation region were also obtained with the fence probe.     

Measurements of velocity,temperature and velocity fluctuation distributions in falling liquid films

The velocity, temperature and velocity fluctuation distributions within falling spindle oil films in an inclined rectangular channel were measured using hot-wire techniques and thin thermocouples. The interfacial shear was caused by cocurrent air flow.The results indicate that the liquid films are as a whole much more laminar-like than turbulent in a range of Reynolds numbers (4γ/μ) up to the experimental limit of 6000. Mixing motion occurs in the vicinity of the interface; however, the flow near the wall surface exhibits no sign of such eddy motions, as predicted by the wall law for single phase turbulent flow. Although velocity fluctuation is observed within films with interfacial shear, mean velocity profiles are approximately the same as those obtained by the laminar film prediction.     

Experimental Investigation of Separated Shear Flow under Subharmonic Perturbations over a Backward-Facing Step

Subharmonic-perturbed shear flow downstream of a two-dimensional backward-facing step was experimentally investigated. The Reynolds number was Re_h = 2.0 ×10⁴, based on free-stream velocity and step height. Planar 2D-2C particle image velocimetry was employed to measure the separating and reattaching flow in the horizontal-vertical plane in the center position. The subharmonic perturbations were generated by an oscillating flap which was implemented over the step edge and driven by periodic Ampere force. The subharmonic frequency was 55 Hz as the half of the fundamental frequency of the turbulent shear layer. As a result of the subharmonic perturbations, the size of recirculation region behind the backward-facing step is reduced and the time-averaged reattachment length is 31.0% shorter than that of the natural flow. The evolution of vortices, including vortex roll-up, growth and breakdown process, is analyzed by using phase-averaging, cross-correlation function and proper orthogonal decomposition. It is found that Reynolds shear stress is considerably increased in which the vortices roll up and then break down further downstream. In particular, rapid growth of vortices based on the “step mode” occurs at approximate half of the recirculation region, caused by in interaction between the shear layer and the recirculation region. Furthermore, the coherent structures, which are represented by a phase-correlated POD mode pair, are reconstructed in phases in order to show regular patterns of the subharmonic-perturbed coherent structures.     

Turbulent separated reattached flow in a two-dimensional curved-wall diffuser

A turbulent separation-reattachment flow in a two-dimensional asymmetrical curved-wall diffuser is studied by a two-dimensional laser doppler velocimeter. The turbulent boundary layer separates on the lower curved wall under strong pressure gradient and then reattaches on a parallel channel. At the inlet of the diffuser, Reynolds number based on the diffuser height is 1.2×10⁵ and the velocity is 25.2m/s. The results of experiments are presented and analyzed in new defined streamline-aligned coordinates. The experiment shows that after Transitory Detachment Reynolds shear stress is negative in the near-wall backflow region. Their characteristics are approximately the same as in simple turbulent shear layers near the maximum Reynolds shear stress. A scale is formed using the maximum Reynolds shear stresses. It is found that a Reynolds shear stress similarity exists from separation to reattachment and the Schofield-Perry velocity law exists in the forward shear flow. Both profiles are used in the experimental work that leads to the design of a new eddy-viscosity model. The length scale is taken from that developed by Schofield and Perry. The composite velocity scale is formed by the maximum Reynolds shear stress and the Schofield-Perry velocity scale as well as the edge velocity of the boundary layer. The results of these experiments are presented in this paper.     

Self-aeration and turbulence in a stepped channel: Influence of cavity surface roughness

The strong interactions between free-surface flows and atmospheric surroundings may lead to substantial air–water mixing with void fractions ranging from zero in clear-water to 100%. In this study, the air–water flow properties were studied in a large stepped water channel operating at large Reynolds numbers. Interactions between free-surface and cavity recirculation were systematically investigated in the skimming flow regime. Some surface roughness was introduced on the cavity walls and identical experiments were performed with several configurations. Basic results demonstrated some influence of step surface roughness on the flow properties leading to some counter-intuitive finding. The presence of cavity roughness was associated with higher flow velocities and comparatively lower turbulence levels. Distributions of bubble/droplet chords spanned over several orders of magnitude without significant influence of the cavity roughness. The distributions of turbulence levels and bubble count rates showed some correlation and highlighted strong interactions between entrained particles (bubbles, drops) and the flow turbulence.     

Experimental Investigation of Three-Dimensional Vortex Structures Downstream of Vortex Generators Over a Backward-Facing Step

An experimental investigation of vortex generators has been carried out in turbulent backward-facing step (BFS) flow. The Reynolds number, based on a freestream velocity U₀ = 10 m/s and a step height h = 30 mm, was Re_h = 2.0 × 10⁴. Low-profile wedge-type vortex generators (VGs) were implemented on the horizontal surface upstream of the step. High-resolution planar particle image velocimetry (2D-2C PIV) was used to measure the separated shear layer, recirculation region and reattachment area downstream of the BFS in a single field of view. Besides, time-resolved tomographic particle image velocimetry (TR-Tomo-PIV) was also employed to measure the flow flied of the turbulent shear layer downstream of the BFS within a three-dimensional volume of 50 × 50 × 10 mm³ at a sampling frequency of 1 kHz. The flow control result shows that time-averaged reattachment length downstream of the BFS is reduced by 29.1 % due to the application of the VGs. Meanwhile, the Reynolds shear stress downstream of the VGs is considerably increased. Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) have been applied to the 3D velocity vector fields to analyze the complex vortex structures in the spatial and temporal approaches, respectively. A coherent bandwidth of Strouhal number 0.3 < St_h < 0.6 is found in the VG-induced vortices, and moreover, Λ-shaped three-dimensional vortex structures at St_h = 0.37 are revealed in the energy and dynamic approaches complementarily.     

Parametric study of separation and transition characteristics over an airfoil at low Reynolds numbers

Time-resolved surface pressure measurements are used to experimentally investigate characteristics of separation and transition over a NACA 0018 airfoil for the relatively wide range of chord Reynolds numbers from 50,000 to 250,000 and angles of attack from 0° to 21°. The results provide a comprehensive data set of characteristic parameters for separated shear layer development and reveal important dependencies of these quantities on flow conditions. Mean surface pressure measurements are used to explore the variation in separation bubble position, edge velocity in the separated shear layer, and lift coefficients with angle of attack and Reynolds number. Consistent with previous studies, the separation bubble is found to move upstream and decrease in length as the Reynolds number and angle of attack increase. Above a certain angle of attack, the proximity of the separation bubble to the location of the suction peak results in a reduced lift slope compared to that observed at lower angles. Simultaneous measurements of the time-varying component of surface pressure at various spatial locations on the model are used to estimate the frequency of shear layer instability, maximum root-mean-square (RMS) surface pressure, spatial amplification rates of RMS surface pressure, and convection speeds of the pressure fluctuations in the separation bubble. A power-law correlation between the shear layer instability frequency and Reynolds number is shown to provide an order of magnitude estimate of the central frequency of disturbance amplification for various airfoil geometries at low Reynolds numbers. Maximum RMS surface pressures are found to agree with values measured in separation bubbles over geometries other than airfoils, when normalized by the dynamic pressure based on edge velocity. Spatial amplification rates in the separation bubble increase with both Reynolds number and angle of attack, causing the accompanying decrease in separation bubble length. Values of the convection speed of pressure fluctuations in the separated shear layer are measured to be between 35 and 50% of the edge velocity, consistent with predictions of linear stability theory for separated shear layers.     

Three-dimensional vortex wake structure of a flapping-wing micro aerial vehicle in forward flight configuration

The unsteady flow field past a backward-facing step in a rectangular duct is investigated by adopting time-resolved particle image velocimetry (PIV) in the Reynolds number range of 2,640–9,880 based on step height and the inlet average velocity. The PIV realizations are subjected to post-processing techniques, namely, proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD). At low Reynolds numbers, the second spatial POD modes indicate the presence of the shear layer mode, whereas this feature shifts to higher modes at higher Reynolds numbers. The corresponding temporal modes are Fourier-transformed to obtain the dominant frequency, whose Strouhal number corroborates the above observation. Short-time windows in the transverse velocity component along the shear layer are selected to investigate the temporal stability of the flow field by DMD to quantify the growth rate of the shear layer mode. The higher harmonics of this mode are also observed to grow, albeit at lesser rate. By relating to POD analysis, the most energetic structures were found to correspond to the unstable modes. The correlation between these unstable DMD modes and the Fourier-filtered flow fields for the same frequencies indicate better match for the lower operating Reynolds number case as compared to higher ones. The spatial stability analysis demonstrates the growth of the shear layer vortices, which is combined with the temporal stability analysis to evaluate the phase velocity of the identified shear layer structures. The calculated phase velocity magnitude of the shear layer is found to be reasonably below the local velocity as expected.     

Experimental investigation on the interaction between polymer solution jet and free surface

Turbulence modification with polymer additives, i.e., Tom's effect, is a well known phenomenon. In this study, Tom's effect at a free surface was experimentally investigated. The turbulence at a free surface was generated by a horizontal liquid jet. A new specklegram technique was proposed in order to quantitatively measure the free-surface shapes caused by a turbulent jet. The specklegram method was very simple and was able to evaluate the free-surface waves accurately. The measurement confirmed that the surface of a polymer solution was less wavy than that of water. The jet beneath the free surface was measured by the LDV. The turbulence diffusion of the polymer jet was much smaller than that of the water jet. The surface turbulence was found to be modified by the polymer solution indirectly. The modification was a consequence of the Tom's effect at the shear layer around the jet.     

Turbulence manipulation in air–water flows on a stepped chute: An experimental study

In a stepped channel operating with large flow rates, the flow skims over the pseudo-bottom formed by the step edges as a coherent stream. Intense three-dimensional recirculation is maintained by shear stress transmission from the mainstream to the step cavities, while significant free-surface aeration takes place. The interactions between free-surface aeration and cavity recirculation are investigated herein with seven step cavity configurations. The experiments were conducted in a large stepped channel operating at large Reynolds numbers. For some experiments, triangular vanes, or longitudinal ribs, were placed across the step cavities to manipulate the flow turbulence to enhance the interactions between the mainstream flow and the cavity recirculation region. The results showed a strong influence of the vanes on the air–water flow properties in both free-stream and cavity flows. The findings demonstrate some passive turbulence manipulation in highly turbulent air–water flows.