Wall shear stress modified by bubbles in a horizontal channel flow of silicone oil in the transition region

We performed laboratory experiments on bubbly channel flows using silicone oil, which has a low surface tension and clean interface to bubbles, as a test fluid to evaluate the wall shear stress modification for different regimes of bubble migration status. The channel Reynolds numbers of the flow ranged from 1000 to 5000, covering laminar, transition and turbulent flow regimes. The bubble deformation and swarms were classified as packing, film, foam, dispersed, and stretched states based on visualization of bubbles as a bulk void fraction changed. In the dispersed and film states, the wall shear stress reduced by 9% from that in the single-phase condition; by contrast, the wall shear stress increased in the stretched, packing, and foam states. We carried out statistical analysis of the time-series of the wall shear stress in the transition and turbulent-flow regimes. Variations of the PDF of the shear stress and the higher order moments in the statistic indicated that the injection of bubbles generated pseudo-turbulence in the transition regime and suppressed drag-inducing events in the turbulent regime. Bubble images and measurements of shear stress revealed a correlated wave with a time lag, for which we discuss associated to the bubble dynamics and effective viscosity of the bubble mixture in wall proximity.     

Turbulent shear stress profiles in a bubbly channel flow assessed by particle tracking velocimetry

Particle tracking velocimetry (PTV) is applied to a bubbly two-phase turbulent flow in a horizontal channel at Re = 2 × 10⁴ to investigate the turbulent shear stress profile which had been altered by the presence of bubbles. Streamwise and vertical velocity components of liquid phase are obtained using a shallow focus imaging method under backlight photography. The size of bubbles injected through a porous plate in the channel ranged from 0.3 to 1.5 mm diameter, and the bubbles show a significant backward slip velocity relative to liquid flow. After bubbles and tracer particles are identified by binarizing the image, velocity of each phase and void fraction are profiled in a downstream region. The turbulent shear stress, which consists of three components in the bubbly two-phase flow, is computed by analysis of PTV data. The result shows that the fluctuation correlation between local void fraction and vertical liquid velocity provides a negative shear stress component which promotes frictional drag reduction in the bubbly two-phase layer. The paper also deals with the source of the negative shear stress considering bubble’s relative motion to liquid.     

Skin friction reduction by large air bubbles in a horizontal channel flow

Microbubble and air film methods are believed to be applicable to skin friction reduction in ships. Small bubbles are dispersed into the turbulent boundary layer in the former case, and wide air layers cover the wall surface in the latter case. Previous studies did not specifically address the intermediate case between the microbubble and air film conditions. This study is concerned with the possibility and mechanism of drag reduction using relatively large air bubbles compared to the boundary layer thickness in a horizontal turbulent channel flow. The relationship between local skin friction and the bubble’s interfacial structure is investigated by synchronizing the measurement of wall-shear stress with the image acquisition of bubbles. The bubble sizes range from 2 to 90 mm approximately. As a result, a negative correlation between the local skin friction and the local void fraction is confirmed by the time-resolved measurement. A new observation is the fact that the local skin friction decreases drastically in the rear part of individual large bubbles, and rapidly increases after the bubble’s rear interface passes. This characteristic underlies the bubble-size dependency of the average skin friction in the intermediate bubble size condition.     

An experimental study on the effect of air bubble injection on the flow induced rotational hub

Modification of shear stress due to air bubbles injection in a rotary device was investigated experimentally. Air bubbles inject to the water flow crosses the neighbor of the hub which can rotate just by water flow shear stresses, in this device. Increasing air void fraction leads to decrease of shear stresses exerted on the hub surface until in high void fractions, the hub motion stopped as observed. Amount of skin friction decrease has been estimated by counting central hub rotations. Wall shear stress was decreased by bubble injection in all range of tested Reynolds number, changing from 50,378 to 71,238, and also by increasing air void fraction from zero to 3.06%. Skin friction reduction more than 85% was achieved in this study as maximum measured volume of air fraction injected to fluid flow while bubbles are distinct and they do not make a gas layer. Significant skin friction reduction obtained in this special case indicate that using small amount of bubble injection causes large amount of skin friction reduction in some rotary parts in the liquid phases like as water.     

Numerical investigation of turbulent channel flow controlled by spatially oscillating spanwise Lorentz force

A formulation of the skin-friction drag related to the Reynolds shear stress in a turbulent channel flow is derived. A direct numerical simulation (DNS) of the turbulent control is performed by imposing the spatially oscillating spanwise Lorentz force. Under the action of the Lorentz force with several proper control parameters, only the periodically well-organized streamwise vortices are finally observed in the near-wall region. The Reynolds shear stress decreases dramatically, especially in the near-wall area, resulting in a drag reduction.     

Drag reduction effectiveness and shear stability of polymer-polymer and polymer-fibre mixtures in recirculatory turbulent flow of water

The turbulent drag reduction caused by polymer-polymer and polymerfibre mixtures has been measured in recirculatory flow of water. Shear stability studies have also been made on a number of drag reducing polymers, asbestos fibres and their mixtures in recirculatory turbulent flow of water. Reynolds numbers ranged from 20,000 to 57,000. Both positive and negative deviations from linear additive behaviour have been observed in drag reduction caused by the polymer-polymer mixtures depending upon their compositions, flow rate and polymer species in the mixture. The drag reduction by the mixtures has been predicted by using simple mixture rule equations including an interaction parameter. This interaction parameter is believed to depend upon the polymer interaction in the polymer mixture. The random coil size and rigidity of the polymer molecules appear to be responsible for the synergism observed in the drag reduction caused by the mixture. In general, mixtures having larger solvation number seem to give positive synergism.Synergism in drag reduction by the polymer-fibre mixtures has also been observed. The simple mixture law equation with interaction parameter is also applicable in predicting the drag reduction by the mixtures as above. The random coil size of the polymer molecules and the rigidity of the polymer-fibre system appear to be responsible for the synergism observed in drag reduction. In the shearstability studies it has been observed that the decrement in drag reduction (DR) is higher than the the decrement in absolute viscosity in most cases. Carboxymethyl cellulose is found to be the most shear stable polymer followed by guar gum, xanthan gum and polyacrylamide. The mixtures exhibiting synergism in causing drag reduction are found to be more shear stable.     

Influence of curvature on drag reduction by opposition control in turbulent flow along a thin cylinder

Opposition controlled fully developed turbulent flow along a thin cylinder is analyzed by means of direct numerical simulations. The influence of cylinder curvature on the skin-friction drag reduction effect by the classical opposition control (i.e., the radial velocity control) is investigated. The curvature of the cylinder affects the uncontrolled flow statistics; for instance, skin-friction coefficient increases while Reynolds shear stress (RSS) and turbulent intensity decrease. However, the control effect in the case of a small curvature is similar to that in channel flow. When the curvature is large, the maximum drag reduction rate decreased. However, the optimal location of the detection plane is the same as that in a flat plate. Further, the drag reduction effect is achieved even on a high detection plane where the drag increases in the flat plate. Although a difference in the drag reduction effect can be observed with a change in the curvature, its mechanism considered in this analysis based on the transport of the Reynolds stress is similar to that of the flat plate.     

Air–water flow in a vertical pipe: experimental study of air bubbles in the vicinity of the wall

This study deals with the influence of bubbles on a vertical air–water pipe flow, for gas-lift applications. The effect of changing the bubble size is of particular interest as it has been shown to affect the pressure drop over the pipe. Local measurements on the bubbles characteristics in the wall region were performed, using standard techniques, such as high-speed video recording and optical fibre probe, and more specific techniques, such as two-phase hot film anemometry for the wall shear stress and conductivity measurement for the thickness of the liquid film at the wall. The injection of macroscopic air bubbles in a pipe flow was shown to increase the wall shear stress. Bubbles travelling close to the wall create a periodic perturbation. The injection of small bubbles amplifies this effect, because they tend to move in the wall region; hence, more bubbles are travelling close to the wall. A simple analysis based on a two-fluid set of equations emphasised the importance of the local gas fraction fluctuations on the wall shear stress.     

Parametric Study on a Sinusoidal Riblet for Drag Reduction by Direct Numerical Simulation

The direct numerical simulation of fully developed turbulent channel flow with a sinusoidal riblet surface has been carried out at the friction Reynolds number of 110. Lateral spacing of adjacent walls in a sinusoidal riblet is varied sinusoidally in the streamwise direction. The average lateral spacing of a sinusoidal riblet is larger than the diameter of a quasi-streamwise vortex and its wetted area is smaller than that of ordinary straight-type riblets. We investigate the effect of sinusoidal riblet design parameters on the drag reduction rate and flow statistics in this paper. The parametric study shows that the maximum total drag reduction rate is approximately 9.8% at a friction Reynolds number of 110. The riblet induces downward and upward flows in the expanded and contracted regions, respectively, which contribute to periodic Reynolds shear stress. However, the random Reynolds shear stress decreases drastically as compared with the flat surface case, resulting in the reduction of total drag owing to the sinusoidal riblet. We also performed vortex tracking to discuss the motion of the vortical structure traveling over the sinusoidal riblet surface. Vortex tracking and probability analysis for the core of the vortical structure show that the vortical structure is attenuated owing to the sinusoidal riblet and follows the characteristic flow. These results show that the high skin-friction region on the channel wall is localized at the expanded region of the riblet walls. In consequence, the wetted area of the riblet decreases, resulting in the drag-reduction effect.     

Influence of bubble size on micro-bubble drag reduction

Micro-bubble drag reduction experiments were conducted in a turbulent water channel flow. Compressed nitrogen was used to force flow through a slot injector located in the plate beneath the boundary layer of the tunnel test section. Gas and bubbly mixtures were injected into a turbulent boundary layer (TBL), and the resulting friction drag was measured downstream of the injector. Injection into tap water, a surfactant solution (Triton X-100, 20 ppm), and a salt-water solution (35 ppt) yielded bubbles of average diameter 476, 322 and 254 μm, respectively. In addition, lipid stabilized gas bubbles (44 μm) were injected into the boundary layer. Thus, bubbles with d ⁺ values of 200 to 18 were injected. The results indicate that the measured drag reduction by micro-bubbles in a TBL is related strongly to the injected gas volumetric flow rate and the static pressure in the boundary layer, but is essentially independent of the size of the micro-bubbles over the size range tested.     

Optimization of the injection process for drag-reducing additives

The primary objective of this experimental study was to determine the optimum combination of additive concentration, additive flowrate, injector angle and injector width for injecting a drag-reducing additive into a channel flow of water. The experiments were designed to keep an effective concentration of the additive in the buffer region where previous experiments have shown the additives directly affect the turbulent structures. Flow visualization of the wall-layer structures was conducted for the optimum combination of injection variables at the streamwise locations where the drag reduction peaked and where the additive became fully mixed with the channel flow of water.     

Experimental study on the effects of shear induced structure in a drag-reducing surfactant solution flow

In this paper, The drag reduction characteristics of surfactant solutions have been experimentally studied, as well as, the shear viscosities of turbulent drag-reducing surfactant solution have been measured as a function of concentration, shear rate and temperature by using AG-G2 (TA Instruments, New Castle, USA) rheometer. In comparison the rheological property with the macroscopic behavior of the solutions in turbulent channel flow, a deeper insight into the mechanisms of drag-reducing surfactant solution has been obtained. For no shear induced structure of surfactant solutions they just show features shear thinning, but the drag reduction is very significant phenomenon. Surfactant solution of the shear induced structure is not a surfactant fluid drag reduction of the necessary elements.     

Lattice Boltzmann simulations of turbulent shear flow between parallel porous walls

The effects of two parallel porous walls are investigated, consisting of the Darcy number and the porosity of a porous medium, on the behavior of turbulent shear flows as well as skin-friction drag. The turbulent channel flow with a porous surface is directly simulated by the lattice Boltzmann method （LBM）. The Darcy-Brinkman- Forcheimer （DBF） acting force term is added in the lattice Boltzmann equation to simu- late the turbulent flow bounded by porous walls. It is found that there are two opposite trends （enhancement or reduction） for the porous medium to modify the intensities of the velocity fluctuations and the Reynolds stresses in the near wall region. The parametric study shows that flow modification depends on the Darcy number and the porosity of the porous medium. The results show that, with respect to the conventional impermeable wall, the degree of turbulence modification does not depend on any simple set of param- eters obviously. Moreover, the drag in porous wall-bounded turbulent flow decreases if the Darcy number is smaller than the order of O（10-4） and the porosity of porous walls is up to 0.4.     

Numerical study of the turbulent channel flow under space-dependent electromagnetic force control at different Reynolds numbers

In this paper, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with the direct numerical simulation(DNS) methods for different Reynolds numbers. A formulation is derived to express the relation between the drag and the Reynolds shear stress. With the application of optimal electromagnetic force, the in-depth relations among characteristic structures in the flow field, mean Reynolds shear stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The regular quasi-streamwise vortex structures, which appear in the flow field, have the same period with that of the electromagnetic force.These structures suppress the random velocity fluctuations, which leads to the absolute value of mean Reynolds shear stress decreasing and the distribution of that moving away from the wall. Moreover, the wave number of optimal electromagnetic force increases,and the scale of the regular quasi-streamwise vortex structures decreases as the Reynolds number increases. Therefore, the rate of drag reduction decreases with the increase in the Reynolds number since the scale of the regular quasi-streamwise vortex structures decreases.     

Experiments with three-dimensional riblets as an idealized model of shark skin

The skin of fast sharks exhibits a rather intriguing three-dimensional rib pattern. Therefore, the question arises whether or not such three-dimensional riblet surfaces may produce an equivalent or even higher drag reduction than straight two-dimensional riblets. Previously, the latter have been shown to reduce turbulent wall shear stress by up to 10%. Hence, the drag reduction by three-dimensional riblet surfaces is investigated experimentally. Our idealized 3D-surface consists of sharp-edged fin-shaped elements arranged in an interlocking array. The turbulent wall shear stress on this surface is measured using direct force balances. In a first attempt, wind tunnel experiments with about 365,000 tiny fin elements per test surface have been carried out. Due to the complexity of the surface manufacturing process, a comprehensive parametric study was not possible. These initial wind tunnel data, however, hinted at an appreciable drag reduction. Subsequently, in order to have a better judgement on the potential of these 3D-surfaces, oil channel experiments are carried out. In our new oil channel, the geometrical dimensions of the fins can be magnified 10 times in size as compared to the initial wind tunnel experiments, i.e., from typically 0.5 mm to 5 mm. For these latter oil channel experiments, novel test plates with variable fin configuration have been manufactured, with 1,920–4,000 fins. This enhanced variability permits measurements with a comparatively large parameter range. As a result of our measurements, it can be concluded, that 3D-riblet surfaces do indeed produce an appreciable drag reduction. We found as much as 7.3% decreased turbulent shear stress, as compared to a smooth reference plate. However, in direct comparison with 2D riblets, the performance of 3D-riblets is still inferior by about 1.7%. On the other hand, it appears conceivable, with a careful design of the fin shape (possibly supported by theory), that this inferiority in performance might be reduced. Nevertheless, at present, it seems to be rather unlikely, that 3D-riblets can significantly outperform 2D-riblets. Finally, one interesting finding remains to be mentioned: The optimum drag reduction for short 3D-riblets occurs at a lower rib height than for longer 3D-riblets or for infinitely long 2D-riblets. The same observation had been made previously on shark scales of different species with differing rib lengths, but no explanation could be given. Received: 1 March 1999/Accepted: 16 July 1999     

Influence of wavy structured surfaces and large scale polymer structures on drag reduction

Drag reduction was studied for turbulent flow over a structured wall that contained 600 sinusoidal waves with a wavelength of 5 mm and an amplitude of 0.25 mm. A concentrated solution of a co-polymer of polyacrylamide and sodium acrylate was injected into the flow through wall slots. Laser Doppler velocimetry was used to measure turbulence. A fluorescence technique was developed that enabled us to demonstrate the existence, under certain circumstances, of large gelatinous structures in the injected polymer solution and in the flow channel.At maximum drag reduction, the Reynolds shear stress was zero and the velocity field was the same as found for a smooth surface. Larger drag reductions could be realized for a wavy wall because the initial drag was larger. The influences of polymers on the turbulent fields are similar for smooth and wavy boundaries. These results are of interest since the interaction with the wall can be quite different for water flow over smooth and wavy boundaries (which are characterized as being completely rough). An important effect of polymers is a decreasing relative importance of high frequency fluctuations with increasing drag reduction that is characterized by a cut-off frequency. This cut-off is the same for smooth and wavy walls at maximum drag reduction. The sensitivity of drag reduction to the method of preparing and delivering the polymer solution suggests that aggregation of polymers could be playing an important role for the system that was studied. For example, drag reduction was enhanced when large polymer structures are present.     

Evaluation of flush mounted hot-film sensors for skin friction reduction measurements in viscoelastic polymer solutions

The purpose of this investigation was to evaluate the performance of flush mounted hot-film sensors for mean wall shear stress measurement in turbulent flows of dilute drag reducing polymer solution. A series of pipe flow expriments were conducted over a range of Reynolds numbers and polymer solution concentrations to compare the level of skin friction drag reduction measured by hot-film sensors with values calculated from pipe pressure drop. It is shown that water calibrated hot-film sensors consistently underestimate the wall shear stress suggesting that Reynolds analogy is not valid in dilute polymer solutions. The Newtonian form of the relationship between the wall shear stress and the heat transfer remains valid in dilute polymer solutions. However, multiplicative and additive factors in the relationship are shown to increase linearly with the logarithm of the polymer concentration.     

Image analysis applied to study on frictional-drag reduction by electrolytic microbubbles in a turbulent channel flow

We investigate frictional-drag reduction with electrolytic microbubbles based on image measurement of a turbulent flow in a water channel at Re = 4800 (based on the half channel height). Microbubbles with a diameter ranging 30–200 μm can reduce frictional drag by as much as 30% relative to single-phase flow even at low void fractions (α ≈ 3 × 10⁻⁴); however, drag reduction is only effective within a limited downstream distance from an electrode array. Arrangement of the optical system allows us to measure the bubble-production rate by water electrolysis from images near the wall and to trace the motion of bubbles. We also measure velocity fields using particle-tracking velocimetry based on a shallow depth-of-field approach by segregating tracer particles from microbubbles. Vertically oscillating microbubbles likely represent interaction with vortical structures near the wall, and bubbles approaching the wall appear to induce negative streamwise velocity relative to the surrounding fluid. We relate the wall friction with the double integral of the Reynolds-stress profile and show that its variation due to microbubbles decreases the drag on the wall. Microbubbles tend to coalesce downstream resulting in a fewer bubbles but with greater size; accordingly, the oscillatory motion diminishes, and the frictional drag rather increases.     

The development of the structure of water – air bubble grid turbulence

The development of the turbulent flow field generated due to the interaction of grid turbulence with a swarm of bubbles is investigated experimentally, in a vertical channel of rectangular cross section. Void fraction and streamwise mean and rms velocity distributions have been measured at several distances from the grid, with an optical probe and Laser Doppler Velocimetry, in relation to the air flow rate ratio. The obtained results indicate that close to the grid the void fraction and velocity distributions are dictated by the bubble injectors’ location on the grid. Downstream the void fraction distribution changes to a double peak pattern. The velocity distribution is characterized by a shear layer between the wall area and the central area of the channel. The extend of this shear layer is increasing as the distance from the grid and the gas flow rate ratio are increasing, and is associated with a corresponding increase of the turbulence fluctuations. Autocorrelation and spectra measurements at the centre of the channel show a reduction of the flow scales for low void fraction. Consistently, power spectra distributions indicate that bubbles cause a redistribution of energy manifested by the relative enhancement of the intermediate scales’ energy content and a consequent reduction in the larger scales. These trends are gradually alleviated and reversed at large distances from the grid, as the air flow rate is increased.     

Numerical and experimental investigation of turbulent characteristics in a drag-reducing flow with surfactant additives

Cetyltrimethyl ammonium chloride (CTAC) surfactant additives, because of their long-life characteristics, can be used as promising drag-reducers in district heating and cooling systems. In the present study we performed both numerical and experimental tests for a 75 ppm CTAC surfactant drag-reducing channel flow. A two-component PIV system was used to measure the instantaneous streamwise and wall-normal velocity components. A Giesekus constitutive equation was adopted to model the extra stress due to the surfactant additives, with the constitutive parameters being determined by well-fitting apparent shear viscosities, as measured by an Advanced Rheometric Expansion System (ARES) rheometer. In the numerical study, we connected the realistic rheological properties with the drag-reduction rate. This is different from previous numerical studies in which the model parameters were set artificially. By performing consistent comparisons between numerical and experimental results, we have obtained an insight into the mechanism of the additive-induced drag-reduction phenomena.

Our simulation showed that the addition of surfactant additives introduces several changes in turbulent flow characteristics: (1) In the viscous sublayer, the mean velocity gradient becomes gentler due to the viscoelastic forces introduced by the additives. The buffer layer becomes expanded and the slope of the velocity profile in the logarithmic layer increases. (2) The locations where the streamwise velocity fluctuation and Reynolds shear stress attain their maximum value shifted from the wall region to the bulk flow region. (3) The root-mean-square velocity fluctuations in the wall-normal direction decrease for the drag-reducing flow. (4) The Reynolds shear stress decreases dramatically and the deficit of the Reynolds shear stress is mainly compensated by the viscoelastic shear stress. (5) The turbulent production becomes much smaller and its peak-value position moves toward the bulk flow region. All of these findings agree qualitatively with experimental measurements.

Regarding flow visualization, the violent streamwise vortices in the near wall region become dramatically suppressed, indicating that the additives weaken the ejection and sweeping motion, and thereby inhibit the generation of turbulence. The reduction in turbulence is accomplished by additive-introduced viscoelastic stress. Surfactant additives have dual effects on frictional drag: (1) introduce viscoelastic shear stress, which increases frictional drag; and (2) dampen the turbulent vortical structures, decrease the turbulent shear stress, and then decrease the frictional drag. Since the second effect is greater than the first one, drag-reduction occurs.