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.     

Air–water flow properties in step cavity down a stepped chute

For the last three decades, the research into skimming flows down stepped chutes was driven by needs for better design guidelines. The skimming flow is characterised by some momentum transfer from the main stream to the recirculation zones in the shear layer developing downstream of each step edge. In the present study some physical modelling was conducted in a relatively large facility and detailed air–water flow measurements were conducted at several locations along a triangular cavity. The data implied some self-similarity of the main flow properties in the upper flow region, at step edges as well as at all locations along the step cavity. In the developing shear layer and cavity region (i.e. y/h < 0.3), the air–water flow properties presented some specific features highlighting the development of the mixing layer downstream of the step edge and the strong interactions between cavity recirculation and mainstream skimming flows. Both void fraction and bubble count rate data showed a local maximum in the developing shear layer, although the local maximum void fraction was always located below the local maximum bubble count rate. The velocity profiles had the same shape as the classical mono-phase flow data. The air–water flow properties highlighted some intense turbulence in the mixing layer that would be associated with large shear stresses and bubble–turbulence interactions.     

Bubbly flow measurements in hydraulic jumps with small inflow Froude numbers

The transition from supercritical to subcritical open channel flow is characterised by a strong dissipative mechanism called a hydraulic jump. A hydraulic jump is turbulent and associated with the development of large-scale turbulence and air entrainment. In the present study, some new physical experiments were conducted to characterise the bubbly flow region of hydraulic jumps with relatively small Froude numbers (2.4 < Fr₁ < 5.1) and relatively large Reynolds numbers (6.6 × 10⁴ < Re < 1.3 × 10⁵). The shape of the time-averaged free-surface profiles was well defined and the longitudinal profiles were in agreement with visual observations. The turbulent free-surface fluctuation profiles exhibited a peak of maximum intensity in the first half of the hydraulic jump roller, and the fluctuations exhibited some characteristic frequencies typically below 3 Hz. The air–water flow properties showed two characteristic regions: the shear layer region in the lower part of the flow and an upper free-surface region above. The air–water shear layer region was characterised by local maxima in terms of void fraction and bubble count rate. Other air–water flow characteristics were documented including the distributions of interfacial velocity and turbulence intensity. The probability distribution functions (PDF) of bubble chord time showed that the bubble chord times exhibited a broad spectrum, with a majority of bubble chord times between 0.5 and 2 ms. An analysis of the longitudinal air–water structure highlighted a significant proportion of bubbles travelling within a cluster structure.     

Turbulence,dynamic similarity and scale effects in high-velocity free-surface flows above a stepped chute

In high-velocity free-surface flows, air entrainment is common through the interface, and intense interactions take place between turbulent structures and entrained bubbles. Two-phase flow properties were measured herein in high-velocity open channel flows above a stepped chute. Detailed turbulence measurements were conducted in a large-size facility, and a comparative analysis was applied to test the validity of the Froude and Reynolds similarities. The results showed consistently that the Froude similitude was not satisfied using a 2:1 geometric scaling ratio. Lesser number of entrained bubbles and comparatively greater bubble sizes were observed at the smaller Reynolds numbers, as well as lower turbulence levels and larger turbulent length and time scales. The results implied that small-size models did underestimate the rate of energy dissipation and the aeration efficiency of prototype stepped spillways for similar flow conditions. Similarly a Reynolds similitude was tested. The results showed also some significant scale effects. However a number of self-similar relationships remained invariant under changes of scale and confirmed the analysis of Chanson and Carosi (Exp Fluids 42:385-401, 2007). The finding is significant because self-similarity may provide a picture general enough to be used to characterise the air–water flow field in large prototype channels.     

Experimental assessment of scale effects affecting two-phase flow properties in hydraulic jumps

A hydraulic jump is the rapid transition from a supercritical to subcritical free-surface flow. It is characterised by strong turbulence and air bubble entrainment. New air–water flow properties were measured in hydraulic jumps with partially developed inflow conditions. The data set together with the earlier data of Chanson (Air bubble entrainment in hydraulic jumps. Similitude and scale effects, 119 p, 2006) yielded similar experiments conducted with identical inflow Froude numbers Fr ₁ = 5 and 8.5, but Reynolds numbers between 24,000 and 98,000. The comparative results showed some drastic scale effects in the smaller hydraulic jumps in terms of void fraction, bubble count rate and bubble chord time distributions. The present comparative analysis demonstrated quantitatively that dynamic similarity of two-phase flows in hydraulic jumps cannot be achieved with a Froude similitude. In experimental facilities with Reynolds numbers up to 10⁵, some viscous scale effects were observed in terms of the rate of entrained air and air–water interfacial area.     

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 time and length scale measurements in high-velocity open channel flows

In high-velocity open channel flows, the measurements of air–water flow properties are complicated by the strong interactions between the flow turbulence and the entrained air. In the present study, an advanced signal processing of traditional single- and dual-tip conductivity probe signals is developed to provide further details on the air–water turbulent level, time and length scales. The technique is applied to turbulent open channel flows on a stepped chute conducted in a large-size facility with flow Reynolds numbers ranging from 3.8E+5 to 7.1E+5. The air water flow properties presented some basic characteristics that were qualitatively and quantitatively similar to previous skimming flow studies. Some self-similar relationships were observed systematically at both macroscopic and microscopic levels. These included the distributions of void fraction, bubble count rate, interfacial velocity and turbulence level at a macroscopic scale, and the auto- and cross-correlation functions at the microscopic level. New correlation analyses yielded a characterisation of the large eddies advecting the bubbles. Basic results included the integral turbulent length and time scales. The turbulent length scales characterised some measure of the size of large vortical structures advecting air bubbles in the skimming flows, and the data were closely related to the characteristic air–water depth Y ₉₀. In the spray region, present results highlighted the existence of an upper spray region for C > 0.95–0.97 in which the distributions of droplet chord sizes and integral advection scales presented some marked differences with the rest of the flow.     

绕水翼空化流动多尺度数值研究

空化的多尺度效应是一种涉及连续介质尺度、微尺度空化泡以及不同尺度间相互转化的复杂水动力学现象, 跨尺度模型的构建是解析该多尺度现象的关键. 本文基于欧拉-拉格朗日联合算法, 通过界面捕捉法求解欧拉体系下大尺度空穴演化, 通过拉格朗日体系下离散空泡模型求解亚网格尺度离散空泡的运动及生长溃灭. 同时, 通过判断空泡与网格尺度间的关系判定不同尺度空化泡的求解模型. 基于建立的多尺度算法对绕NACA66水翼空化流动进行模拟, 将数值结果与实验进行对比, 验证了数值计算方法的准确性. 研究结果表明, 离散空泡数量与空化发展阶段密切相关, 在附着型片状空穴生长阶段, 离散空泡数量波动较小, 离散空泡主要分布在气液交界面位置; 在回射流发展阶段, 离散空泡逐渐增加并分布在回射流扰动区; 在云状空穴溃灭阶段, 离散空泡数量增多且主要分布在气液掺混剧烈的空化云团溃灭区. 在各空化发展阶段, 离散空泡直径概率密度函数均符合伽玛分布. 空化湍流流场特性对拉格朗日空泡空间分布具有重要影响, 离散空泡主要分布在强湍脉动区、旋涡及回射流发展区域.      

Structure of turbulent bubbly jets—I. Methods and centerline properties

The structure of dilute bubbly turbulent round jets, injected vertically upward in still water, was studied both theoretically and experimentally. All measurements were nonintrusive, including mean and fluctuating phase velocities, bubble number intensities, bubble-size distributions and calibration of the motion of individual bubbles. Predictions from three analyses were compared with measurements: (1) locally homogeneous flow analysis, where velocity differences between the phases were neglected; (2) deterministic separated flow analysis, where relative velocity was considered but bubble/turbulence interactions were ignored; and (3) stochasic separated flow analysis, where both relative velocity and bubble/turbulence interactions were considered using random-walk methods. This paper describes theoretical and experimental methods, flow structure near the source and mean properties along the jet axis. Effects of relative velocity were important almost everywhere in the flow; therefore, only the separated flow models yielded satisfactory predictions of bubble velocities along the axis. A companion paper treats mean and fluctuating properties in other regions of the flow.     

Air transport in chute flows

The air bubble rise velocity in still water depends mainly on the bubble size and is basically influenced by buoyancy, viscosity and surface tension. In high-speed flows the number of forces acting on air bubbles increases with turbulence, non-hydrostatic pressure gradient, shear forces, bubble clouds and free-surface entrainment. Air bubbles in these flows are used for cavitation protection of hydraulic structures such as chutes, spillways and bottom outlets. Here, air is normally added by means of aerators upstream of regions where the cavitation number falls below a critical value mainly to reduce the sonic velocity of the fluid and cushion the cavitation bubble collapse process. The distance between successive aerators depends basically on the bubble rise velocity. Until today, the bubble rise velocity in high-speed flows was not thoroughly investigated because of limited laboratory instrumentation. The present project focused on the streamwise development of air concentrations in high-speed flows along a 14 m long model chute. The bubble rise velocity was indirectly derived from the air detrainment gradient of the air concentration contour lines downstream of an aeration device. It accounts for the main hydraulic parameters chute slope, Froude number and air concentration. It is demonstrated that the bubble rise velocity in high-speed flow and stagnant water differ significantly due to fracturing processes, turbulence, and the ambient air concentration.     

Measurement of bubble size, velocity and concentration in flashing flow behind a sudden constriction

In the present paper the results of investigations in flashing flow behind a sudden constriction in vertical upflow are described. Flow visualization, laser-Doppler and phase-Doppler anemometry have been used to measure local bubble and fluid velocities, local bubble sizes and void fractions. The measurements were performed in the midplane of a two-dimensional channel with a 2:1 stepwise constriction.It was found that bubble nucleation takes place in the recirculation zone immediately behind the constriction, which is the location of the lowest static pressure. These bubbles are transported downstream by the mean flow field, while undergoing further growth. No additional nucleation was observed downstream of the recirculation zone. A periodic, cloudwise behaviour of the bubble formation was found which could be explained by the interaction between the bubble growth and the mean flow field. This interaction results in strong disturbances of the mean flow field, which show up as an increase of the fluctuating bubble velocity by a factor of 3 compared to single-phase measurements in a region of 10 step heights behind the constriction. However, these fluctuations appear more like a periodic change in the mean velocity rather than a higher turbulence level. The measured arithmetic mean bubble diameters rise from approx. 50 μm in the recirculation region to about 70–80 μm 50 step heights downstream. Maximum local bubble number density and void fraction were found to be 160001/cm³ and 0.8%, respectively.     

Transport of salts and micron-sized particles entrained from a boiling water pool

The entrainment of soluble (KI, CsI) and non-soluble (Al₂O₃) substances through droplets, which are produced by disintegrating steam bubbles at the surface of a boiling water pool, is determined in a pilot-scale facility. Integral measurements are conducted at steady-state conditions in an atmosphere of either pure steam or an air–steam mixture. The ratio of the entrained liquid mass flow and the gas mass flow through the pool, the entrainment factor, is determined for air–steam ratios between 0 and 0.47 kg/kg in the gas atmosphere and at constant total pressures between 2 and 6 bar. The influence of the vertical temperature profile in the gas atmosphere on the convective velocity field is demonstrated by phase Doppler anemometry and particle image velocimetry measurements at a location 2.1 m above the pool surface. The influences of nucleation and natural convection are demonstrated during slow de-pressurization of the facility at rates below 420 Pa/s.     

Air–water flow characteristics in high-velocity free-surface flows with 50% void fraction

High-velocity free-surface flows are complex two-phase flows and limited information is available about the interactions between air and water for void fractions of about 50%. Herein a detailed experimental study was conducted in the intermediate flow region (C ∼ 50%) on a stepped spillway and the microscopic air–water flow characteristics were investigated. The results showed differences in water and droplet chord times with comparatively larger number of air chord times (0–2 ms), and larger number of water chord times (2–6 ms). A monotonic decrease of particle chord modes was observed with increasing bubble count rates. Several characteristic time scales were identified based upon inter-particle arrival time analyses of characteristic chord time classes as well as spectral analyses of the instantaneous void fraction signal. Chord times of 3–5 ms appeared to be characteristic time scales of the intermediate flow region having similar time scales compared to the local correlation and integral turbulent time scales and to time scales associated with bubble break-up and turbulent velocity fluctuations. A further characteristic time scale of 100 ms was identified in a frequency analysis of instantaneous void fraction. This time scale was of the same order of magnitude as free-surface auto-correlation time scales suggesting that the air–water flow structure was affected by the free-surface fluctuations.     

Structure of turbulent bubbly jets—II. Phase property profiles

This is the second part of a two-part study reporting structure measurements in bubbly turbulent round jets in a still environment. Measurements are compared with three theoretical approaches: (1) locally homogeneous flow analysis, where velocity differences between the phases were neglected; (2) deterministic separated flow analysis, where relative velocity was considered but bubble/turbulence interactions were ignored; and (3) stochastic separated flow analysis, where both relative velocity and bubble/turbulence interactions were considered using random-walk methods. This part of the study considers measurements and predictions of mean and fluctuating phase velocities and mean bubble number intensities at several axial stations. Locally homogeneous flow analysis was not very satisfactory since effects of relative velocity were important for present test conditions. Deterministic separated flow analysis was also ineffective, since neglecting turbulent dispersion caused the width of the bubble-containing region to be underestimated. In contrast, the stochastic separated flow analysis yielded reasonably good predictions.     

Free-surface fluctuations and turbulence in hydraulic jumps

A hydraulic jump is the highly turbulent transition between a high-velocity impinging flow and a turbulent roller. The jump flow is characterised by some substantial air bubble entrainment, spray and splashing. In the present study, the free-surface fluctuations and air-water properties of the hydraulic jump roller were investigated physically for relatively small Froude numbers (2.4 < Fr₁ < 5.1) and relatively large Reynolds numbers (6.6 × 10⁴ < Re < 1.3 × 10⁵). The shape of the mean free surface profile was well defined, and the time-averaged free-surface elevation corresponded to the upper free-surface, with the quantitative values being close to the equivalent clear-water depth. The turbulent fluctuation profiles exhibited a maximum in the first part of the hydraulic jump roller. The free-surface fluctuations presented some characteristic frequencies between 1.4 and 4 Hz. Some simultaneous free-surface measurements at a series of two closely located points yielded the free-surface length and time scales of free-surface fluctuations in terms of both longitudinal and transverse directions. The length scale data seemed to depend upon the inflow Froude number, while the time scale data showed no definite trend. Some simultaneous measurements of instantaneous void fraction and free-surface fluctuations exhibited different features depending upon the phase-detection probe sensor location in the different regions of the roller.     

Water entrainment due to spillway surface jets

Strong flow entrainment has been observed downstream of spillways constructed with flow deflectors. This water entrainment has important environmental and ecological impacts because it improves the mixing of powerhouse and spillway flows, but may negatively impact fish migration or create adverse flow conditions.

Most studies found in the literature attempt to explain this entrainment with turbulent mixing. Both reduced-scale hydraulic models and single-phase, isotropic RANS models grossly under-predict the degree of entrainment observed in prototypes. In this paper, an anisotropic model that accounts for the bubble volume fraction and attenuation of the normal velocity fluctuations at the free surface is presented. The model adequately predicts the main mechanisms causing water entrainment and compares well against experimental data for a round surface jet and for Brownlee Dam at model scale. It is shown that appropriate entrainment can only be captured if the turbulence anisotropy and the two-phase nature of the flow are modelled.     

不同发射深度下导弹水下点火气水流体动力计算

从流体动力角度研究了不同发射深度下，导弹水下点火这一非定常非线性过程。整个系统分为外部水流场、喷管流场和燃气泡流场三个区域加以考虑。水流场采用不可压势流模型，用边界元方法求解；喷管内流场采用非定常一元流动模型，用特征线差分法求解，并设置了激波检测功能；燃气泡采用基于质量和能量守恒的零维计算模型。在时间域中用步进方法实现了三个流场的耦合求解。给出了四种发射深度下的数值计算结果，展示了导弹水下点火的一     

An Experimental Study of the Modulation of the Bubble Motion by Gas–Liquid-Phase Interaction in Oscillating-Grid Decaying Turbulence

The purpose of the present research is to understand dynamic bubble–liquid interaction in a bubbly flow based on the experimental results of the modulation of the bubble motion in oscillating-grid decaying turbulence. By comparing the experimental results obtained from stagnant water and those from oscillating-grid decaying turbulence, we discussed and described detailed process of the modulation of the bubble motion in a water vessel. We discussed the enhancement of the transition of the bubble motion from 2D to 3D by combining the liquid-phase motion obtained through particle imaging velocimetry/laser-induced fluorescence (PIV/LIF) measurement and the bubble wake motion captured through the LIF/HPTS (8-hydroxypyrene-1, 3, 6-trisulfonic acid) method, under both conditions (in the stagnant water and in the oscillating-grid decaying turbulence) in which the initial bubble formation and the bubble motion (gravity-center motion and surface oscillation) were considered to be the same. In addition, by using PIV/LIF measurement along with an infrared shadow technique, we simultaneously obtained the bubble motion (2D zigzagging motion in stagnant water, and 3D motion in the decaying turbulence) and the standard deviation of the liquid-phase motion (the bubble Reynolds number: 775; the turbulent Reynolds number: 62.2). Taking all of the results together, the modulation of the bubble motion in the decaying turbulence, and the dynamic interaction between the bubble and the liquid-phase motion were experimentally and carefully investigated. Consequently, the enhancement and the modulation of the bubble wake motion were considered to be triggered by the collapse of the symmetric property of the bubble–liquid (i.e. ambient liquid-phase turbulence) interaction.     

An experimental study of a turbulent shear layer at a clean and contaminated free-surface

An experimental study was performed to measure the flow properties of a vertically-orientated shear layer in the vicinity of a free-surface. The effect of surface contamination on the near surface flow field was also determined. Digital Particle Image Velocimetry was used to measure instantaneous and averaged velocity, vorticity, and Reynolds stresses. Results show that the presence of surfactants can cause directional shifts of the shear layer, as well as an overall damping of the turbulence in the vicinity of the free-surface, except in the vicinity of a Reynolds ridge where an increase in Reynolds stress was observed.     

Interfacial aeration and bubble count rate distributions in a supercritical flow past a backward-facing step

An example of high-velocity open channel flows is a supercritical flow past an abrupt drop. In such a geometry, the basic air–water flow properties were measured, including distributions of void fraction and bubble count rate, and local air and water chord size distributions, at and downstream of the backward-facing step. The bubble count rate distributions were compared with a conceptual model of streamwise distribution of air and water chords which yields a quasi-parabolic relationship between bubble count rate and void fraction. The proposed model was an attempt to explain the experimental relationship between bubble count rate and void fraction, rather a meticulous breakdown of the complex air–water structure.