Surfactant systems for drag reduction: Physico-chemical properties and rheological behaviour

The first part of the work presents an overview of the physical chemistry of surfactants which in aqueous solutions reduce the frictional loss in turbulent pipe flow. It is shown that these surfactants form rodlike micelles above a characteristic concentraionc _t. The experimental evidence for rodlike micelles are reviewed and the prerequisites that the surfactant system must fulfill in order to form rodlike micelles are given. It is demonstrated by electrical conductivity measurements that the critical concentration for the formation of spherical micelles shows little temperature dependence, whereasc _t increases very rapidly with temperature. The length of the rodlike micelles, as determined by electric birefringence, decreases with rising temperature and increases with rising surfactant concentration. The dynamic processes in these micellar systems at rest and the influence of additives such as electrolytes and short chain alcohols are discussed.In the second part, the rheological behaviour of these surfactant solutions under laminar and turbulent flow conditions are investigated. Viscosity measurements in laminar pipe and Couette flow show the build-up of a shear induced viscoelastic state, SIS, from normal Newtonian fluid flow. A complete alignment of the rodlike micelles in the flow direction in the SIS was verified by flow birefringence. In turbulent pipe flow, drag reduction occurs in these surfactant systems as soon as rodlike micelles are present in the solution. The extent and type of drag reduction, i.e. the shape of the friction factor versus Reynolds number curve, depends directly on the size, number and surface charge of the rodlike micelles. The friction factor curve of each surfactant investigated changes in the same characteristic way as a function of temperature. For each surfactant, independent of concentration, an upper absolute temperature limit,T _L, for drag reduction exists which is caused by the micellar dynamics.T _L is influenced by the hydrophobic chain length and the counter-ion of the surfactant system. A first attempt is made to explain the drag reduction of surfactants by combining the results of these rheological measurements with the physico-chemical properties of the micellar systems.     

Polymer Induced Drag Reduction in a Turbulent Pipe Flow Subjected to a Coriolis Force

The skin friction factor f in a turbulent wall-bounded flow can be greatly reduced by using polymer solutions. In this paper we discuss experimental results on the effect of the Coriolis force on turbulent drag reduction. To study this, a horizontal smooth-walled pipe with internal diameter 25?mm is placed on a horizontal table rotating about its vertical axis. The rotation is made non-dimensional with friction velocity and pipe diameter, to form the Rotation number Ro. For a range of bulk Rotation number (Ro _b ) between 0 and 0.6 for two different Reynolds numbers (Re _b = 15 & 30 × 10³), the pressure drop is measured, from which the average friction factor f is obtained. Additionally the effect of four different polymer concentrations has been investigated. The single-phase results show that the friction factor increases monotonic but gradual with Rotation. With polymer additives a drag reduction is found that increases with concentration, but which is not affected by the rotation.     

添加剂湍流减阻流动与换热研究综述

添加剂湍流减阻是指在液体的管道湍流中添加少量的高分子聚合物或某种表面活性剂从而使湍流阻力大大降低的现象.从其被发现至今,经过近半个世纪的研究(实验研究、理论分析、数值模拟和实际系统的应用研究),尽管对这一现象及其实际应用价值已有了较为深入的认识,但仍有许多方面尚有欠缺,例如对湍流减阻的机理仍然在探索中.本文归纳评述了高分子聚合物或表面活性剂添加剂湍流减阻流动与换热现象的研究现状,从湍流减阻剂的特性、减阻剂的湍流减阻机理、湍流减阻发生时的换热机理、减阻流动速度场分布和换热控制等几个方面综述了添加剂湍流减阻流动与换热特性,并综述了湍流减阻剂在实际工业系统中的应用情况,在对添加剂湍流减阻机理、有湍流减阻发生时的对流换热机理等的理解方面进行了新的总结.      

Effect of local forcing on a turbulent boundary layer

An experimental study is performed to analyze flow structures behind local suction and blowing in a flat-plate turbulent boundary layer. The local forcing is given to the boundary layer flow by means of a sinusoidally oscillating jet issuing from a thin spanwise slot at the wall. The Reynolds number based on the momentum thickness is about Re _θ =1700. The effects of local forcing are scrutinized by altering the forcing frequency (0.011 ≤ f ⁺≤ 0.044). The forcing amplitude is fixed at A ₀=0.4. It is found that a small local forcing reduces the skin friction and the skin friction reduction increases with the forcing frequency. A phase-averaging technique is employed to capture the large-scale vortex evolution. An organized spanwise vortical structure is generated by the local forcing. The cross-sectional area of vortex and the time fraction of vortex are examined by changing the forcing frequency. An investigation of the random fluctuation components reveals that turbulent energy is concentrated near the center of vortical structures. Received: 17 March 2000/Accepted: 3 April 2001     

Drag characteristics of longitudinal and transverse riblets at low dimensionless spacings

A surface grooved with microscopic riblets aligned parallel to the flow is an effective means to reduce the turbulent skin friction up to 10% compared to a smooth surface. The maximum drag reduction is found for a dimensionless rib spacing s ⁺ in the range of 15–17. For s ⁺ < 10, a linear behaviour of the drag reduction curve is predicted by viscous theory. This linear slope of the drag reduction curve is in contradiction to Schlichting’s postulation of a hydraulically smooth behaviour of small-scale roughness in a turbulent flow. This regime of evanescent dimensionless rib spacings is investigated experimentally by direct wall shear stress measurements in a fully developed channel flow. Additionally, a numerical calculation of the viscous flow over riblets was carried out to predict the drag reducing behaviour. The experimental results show a linear drag reducing behaviour down to s ⁺ = 0.3, which is in good agreement with the numerical results of the viscous simulation. The postulation of Schlichting’s hydraulically smooth regime of a rough surface was not confirmed, neither for a riblet surface nor for a surface geometry with grooves oriented perpendicular to the flow. In the latter case, the drag increases as a quadratic function of the roughness height.     

Investigation of turbulence structures in a drag-reduced turbulent channel flow with surfactant additive by stereoscopic particle image velocimetry

In the present study, we employed stereoscopic particle image velocimetry (PIV) to investigate the characteristics of turbulence structures in a drag-reduced turbulent channel flow with addition of surfactant. The tested drag-reducing fluid was a CTAC/NaSal/Water (CTAC: cetyltrimethyl ammonium chloride; NaSal: sodium salicylate) system at 25°C. The weight concentration of CTAC was 30 ppm. Stereoscopic PIV measurement was performed for a water flow (Re=1.1×10⁴) and a CTAC solution flow (Re=1.5×10⁴ with 54% drag reduction) in both the streamwise–spanwise and wall-normal-spanwise planes, respectively. The three-dimensionality of hairpin vortex structures in the near-wall region for wall-bounded turbulent flow was reproduced by conditionally averaging the stereoscopic two-dimensional-three-component velocity fields. A series of wall-normal vortex cores were found to align with the near-wall low-speed streaks with opposite vorticity signals at both sides of the streaks and with the vorticity decreased on average by about one order of magnitude in CTAC solution flow compared with water flow; the spanwise spacing between the near-wall low-speed streaks in the solution flow is increased by about 46%. The streamwise vorticity of the vortex cores appearing in the wall-normal-spanwise plane was also decreased by the use of drag-reducing surfactant additives.     

Application of the Dorodnitsyn finite element method to swirling boundary layer flow

The Dorodnitsyn finite element method for turbulent boundary layer flow with surface mass transfer is extended to include axisymmetric swirling internal boundary layer flow. Turbulence effects are represented by the two-layer eddy viscosity model of Cebeci and Smith¹ with extensions to allow for the effect of swirl. The method is applied to duct entry flow and a 10 degree included-angle conical diffuser, and produces results in close agreement with experimental measurements with only 11 grid points across the boundary layer. The introduction of swirl (w_e/u_e = 0.4) is found to have little effect on the axial skin friction in either a slightly favourable or adverse pressure gradient, but does cause an increase in the displacement area for an adverse pressure gradient. Surface mass transfer (blowing or suction) causes a substantial reduction (blowing) in axial skin friction and an increase in the displacement area. Both suction and the adverse pressure gradient have little influence on the circumferential velocity and shear stress components. Consequently in an adverse pressure gradient the flow direction adjacent to the wall is expected to approach the circumferential direction at some downstream location.     

Direct numerical simulation of turbulent flow and heat transfer in a square duct with natural convection

In this paper, a direct numerical simulation of a fully developed turbulent flow and heat transfer are studied in a square duct with an imposed temperature difference between the vertical walls and the perfectly insulated horizontal walls. The natural convection is considered on the cross section in the duct. The numerical scheme employs a time-splitting method to integrate the three dimensional incompressible Navier-Stokes equation. The unsteady flow field was simulated at a Reynolds number of 400 based on the Mean friction velocity and the hydraulic diameter (Re _m = 6200), while the Prandtl number (Pr) is assumed 0.71. Four different Grashof numbers (Gr = 10⁴, 10⁵, 10⁶ and 10⁷) are considered. The results show that the secondary flow and turbulent characteristics are not affected obviously at lower Grashof number (Gr ≤ 10⁵) cases, while for the higher Grashof number cases, natural convection has an important effect, but the mean flow and mean temperature at the cross section are also affected strongly by Reynolds stresses. Compared with the laminar heat transfer at the same Grashof number, the intensity of the combined heat transfer is somewhat decreased.     

Numerical investigation of turbulent flow,heat transfer and entropy generation in a helical coiled tube with larger curvature ratio

Three-dimensional turbulent forced convective heat transfer and flow characteristics, and the non-dimensional entropy generation number in a helical coiled tube subjected to uniform wall temperature are simulated using the k–ε standard turbulence model. A finite volume method is employed to solve the governing equations. The effects of Reynolds number, curvature ratio, and coil pitch on the average friction factor and Nusselt number are discussed. The results presented in this paper cover a Reynolds number range of 2 × 10⁴ to 6 × 10⁴, a pitch range of 0.1–0.2 and a curvature ratio range of 0.1–0.3. The results show that the coil pitch, curvature ratio and Reynolds number have different effects on the average friction factor and Nusselt number at different cross-sections. In addition, the flow and heat transfer characteristics in a helical coiled tube with a larger curvature ratio for turbulent flow are different from that of smaller curvature ratio for laminar and turbulent flow in certain ways. Some new features that are not obtained in previous researches are revealed. Moreover, the effects of Reynolds number, curvature ratio, and coil pitch on the non-dimensional entropy generation number of turbulent forced convection in a helical coiled tube are also discussed.     

Rheological behavior and birefringence investigations on drag-reducing surfactant solutions of tallow-(tris-hydroxyethyl)-ammonium acetate/sodiumsalicylate mixtures

Rheological and flow birefringent properties of a drag-reducing mixture of tallow-(tris-hydroxiethyl)-ammonium acetate (ETHOQUAD T/13-50) and sodiumsalicylate (NaSal) have been studied as a function of the concentration and of the salt/surfactant molar ratio x. The optimum molar ratio x for drag reduction is around 2.5. It is shown that shear-induced supramicellar structures (SIS) which are believed to be responsible for friction reduction in turbulent pipe flow develop in the presence of NaSal. It was observed that SIS are also formed even if the concentration c exceeds c ^*, i.e., the concentration where the volumes of rotation of the individual rodlike micelles start to overlap. The validity of the stress optical law is discussed. A switch from a reptation-controlled stress relaxation to a kinetically controlled mechanism takes place at x 2.5 for this system.     

Drag measurements on V-grooved surfaces on a body of revolution in axial flow

In water flows with velocities of up to 9 m/s the friction drag of a body of revolution in axial flow was investigated for dependence on the body surface structure. This was done for different types of riblet film fixed on the surface with the riblet direction aligned with the flow. The lateral spacing between the triangular shaped riblets varied between 0.033 mm and 0.152 mm. In all cases the riblet spacing was equal to the riblet height. For comparison a smooth reference film was used.Depending on the Reynolds number and the non-dimensional riblet spacings ⁺, a turbulent drag reduction of up to 9% could be achieved with riblets in comparison with the flow over a smooth surface.In the region of transition to turbulent flow and with non-dimensional riblet spacings ofs ⁺10–15 drag reductions of up to 13% were obtained. It is therefore conjectured, that in addition to hampering the near wall momentum exchange, the riblets can delay the development of initial turbulent structures in time and space.     

A Reynolds stress model for turbulent flow of homogeneous polymer solutions

Using a priori analyses of direct numerical simulation (DNS) data, a Reynolds stress model (RSM) is developed to account for the influence of polymer additives on turbulent flow over a wide range of flow conditions. The Finitely Extensible Nonlinear Elastic-Peterlin (FENE-P) rheological constitutive model is utilized to evaluate the polymer contribution to the stress tensor. Thirteen DNS data sets are used to analyze the budgets of elastic stress–velocity gradient correlations as well as Reynolds stress and dissipation transport. Closures are developed in the framework of the RSM model for all the required unknown and non-linear terms. The polymer stresses, velocity profiles, turbulent flow statistics and the percentage of friction drag reduction predicted by the RSM model are in good agreement with present and those obtained from independent DNS data over a wide range of rheological and flow parameters.     

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.     

Flow field analysis of a turbulent boundary layer over a riblet surface

The near-wall flow structures of a turbulent boundary layer over a riblet surface with semi-circular grooves were investigated experimentally for the cases of drag decreasing (s ⁺=25.2) and drag increasing (s ⁺=40.6). One thousand instantaneous velocity fields over riblets were measured using the velocity field measurement technique and compared with those above a smooth flat plate. The field of view was 6.75 × 6.75 mm² in physical dimension, containing two grooves. Those instantaneous velocity fields were ensemble averaged to get turbulent statistics including turbulent intensities and turbulent kinetic energy. To see the global flow structure qualitatively, flow visualization was also carried out using the synchronized smoke-wire technique under the same experimental conditions. For the case of drag decreasing (s ⁺=25.2), most of the streamwise vortices stay above the riblets, interacting with the riblet tips frequently. The riblet tips impede the spanwise movement of the streamwise vortices and induce secondary vortices. The normalized rms velocity fluctuations and turbulent kinetic energy are small near the riblet surface, compared with those over a smooth flat plate. Inside the riblet valleys, these are sufficiently small that the increased wetted surface area of the riblets can be compensated. In addition, in the outer region (y ⁺ > 30), these values are almost equal to or slightly smaller than those for the smooth plate. For the case of drag increasing (s ⁺=40.6), however, most of the streamwise vortices stay inside the riblet valleys and contact directly with the riblet surface. The high-speed down-wash flow penetrating into the riblet valley interacts actively with the wetted riblet surface and increases the skin friction. The rms velocity fluctuations and turbulent kinetic energy have larger values compared with those over a smooth flat plate. Received: 24 March 1999/Accepted: 10 March 2000     

The effect of polymer additives on the components of the energy balance of a turbulent flow

On the basis of the experimental data obtained, an analysis is made of the effect of polymer additives directly on the generation of turbulent energy, on the dissipation of the energy of the averaged motion, and on the density of the flux of the kinetic energy of the turbulence. The presence of polymer additives in the turbulent flow significantly changes the relationship between the generation of the turbulent energy and the dissipation of the energy of the averaged motion. Under the action of polymer additives, the density of the flux of kinetic energy decreases over the depth of the channel, which, in turn, brings about a decrease in the influx of energy from the averaged motion to the pulsed motion. The following definitions are adopted below: the x₁ axis of a Cartesian system of coordinates coincides with the horizontal axis of symmetry of the channel and with the direction of the averaged motion of the liquid; the x₂ axis is directed upward; the x₃ axis is perpendicular to the lateral wall of the channel; the origin of coordinates, O, coincides with the lower plane (bottom) of the channel. Further, U₀ is the mean velocity of the flow of liquid in the channel; U_i is the local component of the averaged velocity (i=1, 2, 3); H is the height of the channel; z=2x₂/H; Re is the Reynolds number of the averaged flow;v is the coefficient of kinematic viscosity; u_i is the pulsation component of the velocity (i=1, 2, 3); u_* is rate of dynamic friction; A=(₀–_p0 ^–1 is the coefficient of the lowering of the friction resistance with the flow of polymer solutions; ₀ and _p are the coefficients of the friction resistance with the motion of water and polymer solutions in a channel, respectively; c is the weight concentration of the polymer solution (%); ₀ is the friction stress at the wall; U₊ is the velocity of the flow at the axis of the channel.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 53–58, January–February, 1973.     

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.     

Structure of near wall turbulence downstream of a wall mounted protrusion: an interesting Reynolds stress suppression phenomena

A local suppression in the generation of near wall Reynolds stress is achieved by modifying the buffer region and sublayer (y ⁺ <30) of a turbulent pipe flow with a 16.4 wall unit high wall mounted protrusion. Multi-component, multi-point, time resolved laser Doppler velocimetry measurements are made in the undisturbed and modified ARL/PSU glycerin tunnel pipe flow at a Reynolds number of approximately 10000. A downstream converging flow field is produced by the divergence of the approaching mean flow around the protrusion. A pair of counter-rotating vortices, 15 wall units in diameter with common flow down, are generated by the protrusion and also contribute to the wall directed flow convergence. The convergence region is 15 wall units high and more than 100 wall units long and appears to decouple the near wall region from the outer turbulent wall layer. Locally, turbulent velocity fluctuations in the form of Reynolds stress producing events, sweeps and ejections, are retarded within this region. This results in a reduction in near wall uv Reynolds stress and local wall shear. Interestingly, the counter-rotating vortices act to increase turbulent diffusion in a manner which is uncorrelated with Reynolds stress generation.     

An experimental contribution to near-wall measurements by means of a special laser Doppler anemometry technique

 A new experimental technique for the investigation of near-wall turbulence using laser Doppler anemometry is presented, which allows an accurate measurement of the flow field very close to the wall, with good resolution and a high data rate. Such a technique is tested in a fully developed turbulent flow (with Reynolds numbers between 4,300 and 67,000) by carrying out a careful statistical analysis of the streamwise and wall-normal velocity components within the near-wall region, at distances from the wall ranging from approximately y ⁺ = 1 to y ⁺ = 100. The velocity profiles, Reynolds stresses and higher-order moments of the two-dimensional boundary layer are presented. The results, which are in agreement with the most recent data in the literature, testify the validity of the proposed experimental solution. Moreover, the accuracy of the results allows the friction velocity to be calculated as the intercept at the wall of the best linear fit of the total stress profile; in this way, an unambiguous examination of the normalized statistics is possible. Received: 17 April 2001 / Accepted: 15 August 2001     

An investigation of possible mechanisms of heterogeneous drag reduction in pipe and channel flows

When concentrated polymer solutions are injected into the core-region of a turbulent pipe or channel flow, the injected polymer solution forms a thread which preserves its identity far beyond the injection point. The resulting drag reduction is called heterogeneous drag reduction.This study presents experimental results on the mechanism of this type of drag reduction. The experiments were carried out to find out whether this drag reduction is caused by small amounts of polymer removed from the thread and dissolved in the near-wall region of the flow or by an interaction of the polymer thread with the turbulence. The friction behavior of this type of drag reduction was measured for different concentrations in pipes of different cross-sections, but of identical hydraulic diameter. The parameters of the injection, i.e. injector geometry as well as the ratio of the injection to the bulk velocity, were varied. In one set of experiments the polymer thread was sucked out through an orifice and the friction behavior in the pipe was determined downstream of the orifice. In another experiment, near-wall fluid was led into a bypass in order to measure its drag reducing properties. Furthermore, the influence of a water injection into the near-wall region on the drag reduction was studied.The results provide a strong evidence that heterogeneous drag reduction is in part caused by small amount of dissolved polymer in the near-wall region as well as by an interaction of the polymer thread with the turbulence.Nomenclature a channel height - b channel width - c _p concentration of the injected polymer solution - c _R effective polymer concentration averaged over the cross-section - d pipe or hydraulic diameter - d _i injector diameter - DR drag reduction - f friction factor - l downstream distance from injector - L length of a pipe segment - P polymer type - p differential pressure - Re Reynolds number - U bulk velocity - u ^* ratio of injection to bulk velocity - y ⁺ dimensionless wall distance - v kinematic viscosity - density of the fluid - _w wall shear stress     

A note on velocity-profile similarity of turbulent boundary layers with negligible wall stress

Summary The velocity profiles of a turbulent boundary layer with zero skin friction throughout its region of pressure rise, measured by Stratford in 1959, are analyzed in terms of a law of the wall and a velocity-defect law with a common velocity scale and a logarithmic velocity profile in the region of overlap. The analysis deviates from earlier work by Stratford and Townsend. It is shown that the flow in Stratford's boundary layer, even at the largest value of x ₁ at which measurements were taken, is not yet in a state of equilibrium. The velocity scale for turbulent boundary layers with zero skin friction is proportional to the cube root of the pressure gradient.