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.     

减阻工况下壁面周期扰动对湍流边界层多尺度的影响

通过在平板壁面施加不同频率振幅的压电陶瓷振子周期性扰动,进行了湍流边界层主动控制减阻的实验研究.在压电陶瓷振子最大减阻工况下(80 V和160Hz),使用单丝边界层探针对压电振子自由端下游2mm处进行测量,得到不同法向位置流向速度信号的时间序列.通过对比施加控制前后的多尺度分析,发现压电振子产生的扰动只对近壁区产生影响,使得近壁区大尺度脉动降低,小尺度脉动强度增大,而对边界层的外区则基本没有影响.进一步对大尺度和小尺度的脉动信号进行条件平均,发现压电振子产生的扰动对小尺度脉动的影响在时间相位上并不均匀,小尺度脉动强度在大尺度脉动为正时比在大尺度脉动为负时具有更明显的增加.这表明壁面周期扰动主要通过使大尺度高速扫掠流体破碎为小尺度结构,来影响相应的高壁面摩擦事件,从而达到减阻效果.      

壁面温度在多参数下对平板边界层的影响研究

为了得到壁面温度在不同来流速度、不同湍流强度条件下对边界层转捩与减阻的影响规律,本文采用Transitionk-kl-ω模型对低来流速度下无压力梯度的光滑平板进行了数值模拟。结果表明,随着来流速度的升高,壁温升高所起到的减阻效果更好,即高来流速度对壁面温度更为敏感。当来流处于中高湍流强度下时,壁温升高能起到推迟转捩的作用,且随着湍流强度的升高,转捩推迟的效果越好,但减阻效果正好相反;当来流处于低湍流强度下时,壁温升高会使得转捩提前发生。壁温升高抑制了边界层内流体的脉动程度,使得层流的稳态不易被破坏,流动更加稳定;同时,壁温升高使得边界层内流体的速度梯度减小,从而降低了壁面摩擦系数,故壁温升高能起到推迟边界层转捩与减阻的作用。     

水下湍流减阻途径分析

综述了壁湍流边界层分层结构及近壁区湍流猝发过程,对现有典型水下湍流减阻技术进行了较深入地总结.根据各类减阻方法的特点对水下湍流减阻方法进行了分类,并分析了各类减阻方法对壁湍流流场的影响特征,总结了维持或延长层流、干扰湍流"猝发"、壁面隔离、增大湍流阻尼和改变壁面物性等5种水下湍流减阻可行途径,为开展水下湍流减阻研究及发现新的减阻方法提供参考.     

Reynolds-averaged modeling of polymer drag reduction in turbulent flows

A novel eddy viscosity model for predicting friction drag reduction induced by polymers in turbulent wall-bounded flows is presented. The approach is based on the elliptic relaxation model modified to account for the modified Reynolds-stress equilibrium established by the presence of elastic polymer chains in the fluid. The increased wall damping of the turbulent fluctuations is obtained by modifying the pressure–strain redistribution term. Polymer solutions are represented using the Finite Extensibility Non-linear elastic FENE-P dumbbell model; only one transport equation for the elongation of the polymer chains is considered. The model reproduces the level of drag reduction observed over a wide range of rheological parameters. In addition, both the mean velocity and the turbulent fluctuations are predicted with good accuracy. The approach is computationally attractive because of its limited increase in computational cost in comparison with its Newtonian counterpart.     

Use of particle image velocimetry to study heterogeneous drag reduction

Large polymer filaments can form when drag reducing polymers are injected through wall slots. The presence of these structures enhances the performance of the drag reducing function by mechanisms which are not understood. This paper shows how particle image velocimetry (PIV) techniques can be used to study changes in the configuration of the injected polymer and in the structure of the velocity field with increasing drag reduction. The filaments are found to behave as solid bodies which break up in high shear regions close to a boundary. The breakup process provides an explanation of why the filaments are not observed close to a wall and offers the possibility of providing a heterogeneous distribution of small aggregates of polymers which could be more effective than uniformly distributed molecules as suggested by Hoyer and Gyr (J Non-Newton Fluid Mech 65:221–240, 1996; J Fluids Eng 120:818–823, 1998), Dunlop and Cox (Phys Fluids 20:203–213, 1977) and Vlachogiannis et al. (Phys Fluid 15:3786–3794, 2004). PIV measurements show dramatic qualitative changes in the velocity patterns at maximum drag reduction.     

Turbulence measurements over longitudinal micro-grooved surfaces

The expectation that drag reducing longitudinal microgrooves modify coherent structures in the near wall region of a turbulent boundary layer is examined experimentally. The experiments were conducted in both water channel and wind tunnel facilities. Results of direct drag measurement and local flow detail are given. The results of the measurements are discussed in the light of various theories proposed to explain the drag reduction properties of microgrooves.     

The effect of modeled drag reduction on the wall region

The low-dimensional model derived for the wall region of a turbulent boundary layer (Aubry et al., 1988) is applied to a drag-reduced flow. In agreement with some experimental results, drag reduction is modeled by thickening the wall region, which is achieved by applying stretching transformations to the original flow. By application of a Galerkin projection, a set of ordinary differential equations (ODEs) is obtained whose structure is identical to the set corresponding to the unmodified flow. The coefficients of the ODEs are modified in a nontrivial way. The bifurcation diagrams plotted for different values of the stretching parameter are different in detail but the structure is globally the same. In particular, the intermittent behavior which Aubry et al. identified with the cyclic bursting events experimentally observed is still present. The scenario by which intermittency appears through a subcritical Hopf bifurcation in which a heteroclinic cycle is created and disappears through a bifurcation to traveling waves is identical. These results hold for values of the stretching between 1 and 2.65, the value at which the top of the buffer layer reaches the centerline of the pipe. This is in agreement with experimental results for flows whose drag is reduced but which still display intermittency. The bifurcations occur in the stretched flow at increased levels of dissipation (relative to the unstretched flow), consistent with theoretical pictures of drag reduction, in which the increase of scale is due to stabilization by an increase of dissipation in the turbulent part of the flow. Moreover, this method is a systematic way to perturb the coefficients of the ODEs of Aubry et al. (1988). Under this kind of perturbation, the behavior of the solution (in the part of the bifurcation diagram physically relevant) is found to be extremely robust.     

Numerical prediction of turbulent boundary layer noise from a sharp-edged flat plate

An efficient hybrid uncorrelated wall plane waves–boundary element method (UWPW-BEM) technique is proposed to predict the flow-induced noise from a structure in low Mach number turbulent flow. Reynolds-averaged Navier-Stokes equations are used to estimate the turbulent boundary layer parameters such as convective velocity, boundary layer thickness, and wall shear stress over the surface of the structure. The spectrum of the wall pressure fluctuations is evaluated from the turbulent boundary layer parameters and by using semi-empirical models from literature. The wall pressure field underneath the turbulent boundary layer is synthesized by realizations of uncorrelated wall plane waves (UWPW). An acoustic BEM solver is then employed to compute the acoustic pressure scattered by the structure from the synthesized wall pressure field. Finally, the acoustic response of the structure in turbulent flow is obtained as an ensemble average of the acoustic pressures due to all realizations of uncorrelated plane waves. To demonstrate the hybrid UWPW-BEM approach, the self-noise generated by a flat plate in turbulent flow with Reynolds number based on chord Re_c = 4.9 × 10⁵ is predicted. The results are compared with those obtained from a large eddy simulation (LES)-BEM technique as well as with experimental data from literature.     

A study with particle-image velocimetry of the influence of drag-reducing polymers on the structure of turbulence

Particle-image velocimetry has been used to study the effect of drag-reducing polymers on the structure of turbulence in a channel flow, under conditions of 41% and 55% drag reduction. The fluctuating velocity fields in the x-y plane and in one x-z plane were measured. The striking features of these results are the damping of small scales and the repression of fluctuations of the velocity component normal to the wall. The role of the wall in creating turbulence diminishes greatly at large drag reductions; Warholic et al. (1999) have shown that a turbulent flow with zero Reynolds stress exists at maximum drag reduction. Velocity fields presented for conditions approaching this critical behavior are of particular interest. Received: 23 April 1999/Accepted: 12 February 2001     

Tomographic PIV investigation on coherent vortex structures over shark-skin-inspired drag-reducing riblets

Nature has shown us that the microstructure of the skin of fast-swimming sharks in the ocean can reduce the skin friction drag due to the well-known shark-skin effect.In the present study,the effect of shark-skin-inspired riblets on coherent vortex structures in a turbulent boundary layer(TBL) is investigated.This is done by means of tomographic particle image velocimetry(TPIV) measurements in channel fl ws over an acrylic plate of drag-reducing riblets at a friction Reynolds number of 190.The turbulent fl ws over drag-reducing riblets are verifie by a planar time-resolved particle image velocimetry(TRPIV) system initially,and then the TPIV measurements are performed.Two-dimensional(2D) experimental results with a dragreduction rate of around 4.81% are clearly visible over triangle riblets with a peak-to-peak spacing s+of 14,indicating from the drag-reducing performance that the buffer layer within the TBL has thickened;the logarithmic law region has shifted upward and the Reynolds shear stress decreased.A comparison of the spatial topological distributions of the spanwise vorticity of coherent vortex structures extracted at different wall-normal heights through the improved quadrant splitting method shows that riblets weaken the amplitudesof the spanwise vorticity when ejection(Q2) and sweep(Q4) events occur at the near wall,having the greatest effect on Q4 events in particular.The so-called quadrupole statistical model for coherent structures in the whole TBL is verified Meanwhile,their spatial conditional-averaged topological shapes and the spatial scales of quadrupole coherent vortex structures as a whole in the overlying turbulent fl w over riblets are changed,suggesting that the riblets dampen the momentum and energy exchange between the regions of near-wall and outer portion of the TBL by depressing the bursting events(Q2 and Q4),thereby reducing the skin friction drag.     

On the skin friction coefficient in viscoelastic wall-bounded flows

Analysis of the skin friction coefficient for wall bounded viscoelastic flows is performed by utilizing available direct numerical simulation (DNS) results for viscoelastic turbulent channel flow. The Oldroyd-B, FENE-P and Giesekus constitutive models are used. First, we analyze the friction coefficient in viscous, viscoelastic and inertial stress contributions, as these arise from suitable momentum balances, for the flow in channels and pipes. Following Fukagata et al. (Phys. Fluids, 14, p. L73, 2002) and Yu et al. (Int. J. Heat. Fluid Flow, 25, p. 961, 2004) these three contributions are evaluated averaging available numerical results, and presented for selected values of flow and rheological parameters. Second, based on DNS results, we develop a universal function for the relative drag reduction as a function of the friction Weissenberg number. This leads to a closed-form approximate expression for the inverse of the square root of the skin friction coefficient for viscoelastic turbulent pipe flow as a function of the friction Reynolds number involving two primary material parameters, and a secondary one which also depends on the flow. The primary parameters are the zero shear-rate elasticity number, El₀, and the limiting value for the drag reduction at high Weissenberg number, LDR, while the secondary one is the relative wall viscosity, μ_w. The predictions reproduce both types A and B of drag reduction, as first introduced by Virk (Nature, 253, p. 109, 1975), corresponding to partially and fully extended polymer molecules, respectively. Comparison of the results for the skin friction coefficient against experimental data shows good agreement for low and moderate drag reduction which is the region covered by the simulations.     

沟槽面湍流边界层结构实验研究

应用激光测速技术和氢气泡流动显示技术对沟槽面湍流边界层特性及近壁区拟序结构特征进行了精细的测量和观察。实验结果表明：与光滑面湍流边界层相比,沟槽面端流边界层的黏性底层厚度、过渡层厚度及流速分布对数公式中的积分常数C均有所增大,说明采用的沟槽面具有减阻特性。此外,无量纲低速带条间距明显减小,最多减小20%,说明无量钢低速带条平均间距的缩短与湍流减阻密切联系。     

Die Wirkung von Polymerzusätzen auf die Turbulenzstruktur in der ebenen Mischungsschicht zweier Ströme

The plane mixing layer formed between two parallel streams moving with different velocities is one of the simplest types of free turbulent boundary layers and has frequently been studied for Newtonian fluids. As a result of this and because of its good experimental accessibility this type of flow provides a good opportunity for obtaining information about the influence of drag reducing additives on the structure of free turbulence. This is all the more so because of the presence of a characteristic vortex structure which can be clearly distinguished from the overlying statistical fine turbulence. The turbulence field was investigated using an existent laser Doppler anemometer system that had been designed for space-time correlation measurements. This enabled measurements to be made of the mainstream velocity as well as of the longitudinal and transversal turbulent fluctuations and, after a simple modification, also of the Reynolds shear stresses and the cross correlation coefficients. The main result of the addition of 50 ppm of the polymer used (Separan AP30) was found to be an intensification of the Reynolds shear stresses. The resulting substantially more rapid increase (than in water) in the thickness of the shear layer can be explained theoretically; such behaviour has also been observed in free jets. On the other hand, the reduced thickness of the mixing layer in the initial region and the associated enhancement of the longitudinal fluctuations and damping of the transversal fluctuations indicate that the main shear flow induces a flow anisotropy by uncoiling and aligning the polymer molecules. The increase in the spreading angle suggests that the entrainment process at the edges of the mixing layer is intensified. This can be explained by the enhancement of the large energy carrying vortices in the turbulence spectrum. This is probably also the reason for the general increase in the correlation coefficients observed at all positions along the centreline of the flow field. However, a complete discussion of the energy transfer mechanism present here, in particular with inclusion of the fine turbulence responsible for dissipation, is only possible with the help of a detailed analysis of the vortex structure in the mixing layer. This is presented in a following paper. The relation between the degree of drag reduction and the intensity of the Reynolds shear stresses enables the direct influence of the rheological properties of the fluid on the turbulent momentum transfer to be clearly recognized.     

Finite Difference Method for the Acoustic Radiation of an Elastic Plate Excited by a Turbulent Boundary Layer: A Spectral Domain Solution

A finite difference method is developed to study, on a two-dimensional model, the acoustic pressure radiated when a thin elastic plate, clamped at its boundaries, is excited by a turbulent boundary layer. Consider a homogeneous thin elastic plate clamped at its boundaries and extended to infinity by a plane, perfectly rigid, baffle. This plate closes a rectangular cavity. Both the cavity and the outside domain contain a perfect fluid. The fluid in the cavity is at rest. The fluid in the outside domain moves in the direction parallel to the system plate/baffle with a constant speed. A turbulent boundary layer develops at the interface baffle/plate. The wall pressure fluctuations in this boundary layer generates a vibration of the plate and an acoustic radiation in the two fluid domains. Modeling the wall pressure fluctuations spectrum in a turbulent boundary layer developed over a vibrating surface is a very complex and unresolved task. Ducan and Sirkis [1] proposed a model for the two-way interactions between a membrane and a turbulent flow of fluid. The excitation of the membrane is modeled by a potential flow randomly perturbed. This potential flow is modified by the displacement of the membrane. Howe [2] proposed a model for the turbulent wall pressure fluctuations power spectrum over an elastomeric material. The model presented in this article is based on a hypothesis of one-way interaction between the flow and the structure: the flow generates wall pressure fluctuations which are at the origin of the vibration of the plate, but the vibration of the plate does not modify the characteristics of the flow. A finite difference scheme that incorporates the vibration of the plate and the acoustic pressure inside the fluid cavity has been developed and coupled with a boundary element method that ensures the outside domain coupling. In this paper, we focus on the resolution of the coupled vibration/interior acoustic problem. We compare the results obtained with three numerical methods: (a) a finite difference representation for both the plate displacement and the acoustic pressure inside the cavity; (b) a coupled method involving a finite difference representation for the displacement of the plate and a boundary element method for the interior acoustic pressure; (c) a boundary element method for both the vibration of the plate and the interior acoustic pressure. A comparison of the numerical results obtained with two models of turbulent wall pressure fluctuations spectrums - the Corcos model [3] and the Chase model [4] - is proposed. A difference of 20 dB is found in the vibro-acoustic response of the structure. In [3], this difference is explained by calculating a wavenumber transfer function of the plate. In [6], coupled beam-cavity modes for similar geometry are calculated by the finite difference method. This revised version was published online in July 2006 with corrections to the Cover Date.     

Turbulent Drag Reduction by Uniform Blowing Over a Two-dimensional Roughness

Direct numerical simulation (DNS) of turbulent channel flow over a two-dimensional irregular rough wall with uniform blowing (UB) was performed. The main objective is to investigate the drag reduction effectiveness of UB on a rough-wall turbulent boundary layer toward its practical application. The DNS was performed under a constant flow rate at the bulk Reynolds number values of 5600 and 14000, which correspond to the friction Reynolds numbers of about 180 and 400 in the smooth-wall case, respectively. Based upon the decomposition of drag into the friction and pressure contributions, the present flow is considered to belong to the transitionally-rough regime. Unlike recent experimental results, it turns out that the drag reduction effect of UB on the present two-dimensional rough wall is similar to that for a smooth wall. The friction drag is reduced similarly to the smooth-wall case by the displacement of the mean velocity profile. Besides, the pressure drag, which does not exist in the smooth-wall case, is also reduced; namely, UB makes the rough wall aerodynamically smoother. Examination of turbulence statistics suggests that the effects of roughness and UB are relatively independent to each other in the outer layer, which suggests that Stevenson’s formula can be modified so as to account for the roughness effect by simply adding the roughness function term.     

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.     

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.