Development characteristics of drag-reducing surfactant solution flow in a duct

Development characteristics of dilute cationic surfactant solution flow have been studied through the measurements of the time characteristics of surfactant solution by birefringence experiments and of the streamwise mean velocity profiles of surfactant solution duct flow by a laser Doppler velocimetry system. For both experiments, the concentration of cationic surfactant (oleylbishydroxymethylethylammonium chloride: Ethoquad O/12) was kept constant at 1000 ppm and the molar ratio of the counter ion of sodium salicylate to the surfactants was at 1.5. From the birefringence experiments, dilute surfactant solution shows very long retardation time corresponding to micellar shear induced structure formation. This causes very slow flow development of surfactant solution in a duct. Even at the end of the test section with the distance of 112 times of hydraulic diameter form the inlet, the flow is not fully developed but still has the developing boundary layer characteristics on the duct wall. From the time characteristics and the boundary layer development, it is concluded that the entry length of 1000 to 2000 times hydraulic diameter is required for fully developed surfactant solution flow.List of abbreviations and symbols A₁, A₂ Coefficients for time constant fitting [-] - B Breadth of the test duct [m] - C₁, C₂ Coefficients for time constant fitting [-] - D Pipe diameter [m] - D_H Hydraulic diameter [m] - g Impulse response function [Pa] - H Width of the test duct [m] - n Index of Bird-Carreau model [-] - Re Reynolds number (=U_mD_H/) - Re_D Pipe Reynolds number (=U_mD/) - Re_x Streamwise distance Reynolds number (=U₀x/) - T Absolute temperature [K] - t Time [s] - t_a Retardation time [s] - t_b Build-up time [s] - t_x Relaxation time [s] - t_x1, t_x2 Relaxation time for double time constant fitting [s] - t Time constant in Bird-Carreau model [s] - U Time mean velocity [m/s] - U_m Bulk mean velocity [m/s] - U_max Maximum velocity in a pipe [m/s] - U₀ Main flow velocity [m/s] - u Friction velocity [m/s] - x, y Coordinates [m] - Shear rate [s^–1] - Mean shear rate [s^–1] - n Birefringence [-] - 99% boundary layer thickness [m] - Solution viscosity [Pa·s] - _P, _S Surfactant and solvent viscosity [Pa·s] - ₀, Zero and infinite viscosity of Bird-Carreau model [Pa·s] - Characteristic time in Maxwell model [s] - Water kinematic viscosity [m²/s] - Density [kg/m³] - Solution shear stress [Pa] - _P, _S Surfactant and solvent shear stress [Pa] - Time in convolution [s]     

On pipe diameter effects in surfactant drag-reducing pipe flows

Remarkable power saving in a fluid transport system is possible if the surfactant drag reduction technology is used. Application of surfactant drag reduction to district heating and cooling systems has been investigated in the past. The establishment of the scale-up law in drag-reducing pipe flows is one of the most important problems in this application. Main purpose of this study is aimed to develop a reliable scale-up law in surfactant drag-reducing flows. As the basic data of surfactant solutions, both non-Newtonian viscosity and viscoelasticity were experimentally determined. A turbulent eddy diffusivity model based on the Maxwell model was employed to estimate the drag reduction of surfactant solutions. The predictions by the turbulence model developed in this study with proper rheological characteristics of surfactant solutions has resulted in a reliable estimation of the pipe diameter effect in surfactant drag-reducing flows over the pipe diameter range from 11 to 150mm. Received: 30 June 1997 Accepted: 29 December 1997     

The effects of salts on the rheological characteristics of a drag-reducing cationic surfactant solution with shear-induced micellar structures

The effect of the counterion salt sodium salicylate (Nasal) on the transient rheological properties of a drag-reducing surfactant system tris (2-hydroxyethyl) tallowalkyl ammonium acetate (TTAA) has been studied with both rheometric and rheo-optical methods. Three types of transient behavior for N₁ and viscosity were identified in 5 mM TTAA solutions depending on the counterion concentration: induction and growth (below equimolar concentration); overshoot and growth (above equimolar concentration); and overshoot then plateau (at high concentrations of Nasal). The transient flow birefringence and orientation angle show trends similar to those of the viscosity and N₁. The second type of transient behavior suggests a two-stage alignment and shear thickening process. The SIS buildup time from the quiescent state, the rebuilding time after a strong preshear, and the relaxation time were also obtained from N₁ measurements, and show a maximum around equimolar conditions. The initial N₁ and viscosity immediately after the flow startup, on the other hand, show a maximum around a ratio of 2.5 to 3 Nasal/TTAA. For solutions with a Nasal concentration in the ratio 1.5 to 3, the steady state values of N₁ and viscosity do not show much variation with Nasal concentration over the shear rate range covered, however. The effect of an addition of sodium chloride (NaCI) to an equimolar Nasal/TTAA solution on the characteristic times and steady state values was also quantified. These rheological results provide us with tools to determine the optimal concentration ratio for practical drag reduction applications.     

Breakdown of the Reynolds analogy over drag-reducing riblets surface

An experiment was carried out in a low-speed wind tunnel to study the near-wall turbulence structure over a heated riblets surface. The results confirmed the main conclusion of previous study that riblets can enhance the convective heat-transfer rate by as much as 35 percent even within the drag- reducing region. The logarithmic profile of the mean temperature was shifted downwards, indicating that viscous-sublayer thickness of the thermal boundary layer was reduced over a riblets surface. The reason for the apparent breakdown of the Reynolds analogy seems to be related to the difference in the turbulence length scale between the thermal and momentum boundary layers in the near-wall region.     

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.     

Elastic properties of a dilute surfactant solution in the shear induced state

The first normal stress differences N ₁ of a highly dilute cationic surfactant solution are investigated in a cone-and-plate rheometer. In continuation of a previous paper (Nowak 1998), where the buildup of a shear induced structure in such a solution was attained after a reduced deformation, the N ₁ turned out to be in proportion to the square of the shear rate γ˙ reduced by a critical value γ˙ _c in a first range above γ˙ _c . At higher shear rates the N ₁ tend to lower values than predicted by this relation. Relaxation experiments were performed in the same geometry to determine the characteristic time scales of the shear induced state's decay. In the lower range above &γdot; _c the stress decay is a monoexponential process, while a second time constant has to be introduced to describe the relaxation in that range, where the N ₁ deviate from the parabolic dependence of the reduced shear rate. Received: 10 May 1999 Accepted: 15 November 2000     

Karhunen–Loeve representations of turbulent channel flows using the method of snapshots

For three‐dimensional flows with one inhomogeneous spatial coordinate and two periodic directions, the Karhunen–Loeve procedure is typically formulated as a spatial eigenvalue problem. This is normally referred to as the direct method (DM). Here we derive an equivalent formulation in which the eigenvalue problem is formulated in the temporal coordinate. It is shown that this so‐called method of snapshots (MOS) has some numerical advantages when compared to the DM. In particular, the MOS can be formulated purely as a matrix composed of scalars, thus avoiding the need to construct a matrix of matrices as in the DM. In addition, the MOS avoids the need for so‐called weight functions, which emerge in the DM as a result of the non‐uniform grid typically employed in the inhomogeneous direction. The avoidance of such weight functions, which may exhibit singular behaviour, guarantees satisfaction of the boundary conditions. The MOS is applied to data sets recently obtained from the direct simulation of turbulence in a channel in which viscoelasticity is imparted to the fluid using a Giesekus model. The analysis reveals a steep drop in the dimensionality of the turbulence as viscoelasticity is increased. This is consistent with the results that have been obtained with other viscoelastic models, thus revealing an essential generic feature of polymer‐induced drag reduced turbulent flows. Published in 2006 by John Wiley & Sons, Ltd.     

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.     

Shear-induced phase separation in cationic surfactant solutions around a rotating sphere

The results of an experimental investigation of the flow of a highly dilute cationic surfactant solution (tetradecyl trimethyl ammonium bromide with added sodium salicylate as counterion, equimolar 2.4mM) around a rotating sphere are presented. The flow and the shear-induced phase transition are visualized by means of a Toepler Schlieren optics. The buildup of the shear-induced structures occurs only above a critical shear rate. Once this critical value is exceeded the shear-induced phase separation starts after a characteristic deformation with the shear rate reduced by the critical one. Further analysis of the obtained data is performed on basis of an analytical calculation of the flow around a rotating sphere in a second order fluid (Thomas and Walters, 1964; Giesekus, 1965). From some characteristic features of the shear-induced structures as induction time and position the parameter describing the elastic properties of the fluid is estimated. Received: 6 January 1998 Accepted: 1 May 1998     

Formation of the shear-induced state in dilute cationic surfactant solutions

In cationic surfactant solutions a change of state occurs due to mechanical stresses. In the dilute regime of rodlike micelles the formation of a so-called Shear-Induced State (SIS) occurs above a critical shear rate. In this context dilute means that there is no sterical interaction between rodlike micelles, the solution is below the overlap concentration. Employing a mathematical model, it is shown that aggregation forces are weak compared to hydrodynamic forces. The mathematical formulation is based on a model of Israelachvili which describes the chemical potential of micelles. Hydrodynamic forces are calculated with a rigid-dumbbell model. SIS formation can be explained by the destruction of rodlike micelles.     

Steady and dynamic shear properties of non-aqueous drag-reducing polymer solutions

The steady and dynamic shear properties of two non-aqueous drag-reducers (a medium molecular weight polyisobutylene and a commercial organic drag-reducer) in kerosene solutions over a wide range of temperature and concentration were presented. The intrinsic and zero-shear viscosity results were used to identify the concentrate regimes of these solutions. A characteristic time constant λ₀, which was based on the spring-bead model for dilute solutions, was employed as the scaling parameter for both steady-shear and dynamic data over a wide range of concentration and temperature. The inadequacy of the Graessley reduced-variable method in the dilute region was illustrated. The shear-thinning behaviour of these polymer solutions could be described by the Carreau model. The dynamic data followed the Zimm and Rouse-like behaviour in the low and high frequency limits. The Cox-Merz rule was obeyed in the low shear rate and frequency regions. The Carreau and the zero-frequency Maxwell time constants appeared to be related to λ₀ by a constant factor over a wide range of polymer concentrations. The finding provides a method for extrapolating viscoelastic information into the drag reduction regime, and could be useful for interpretation of drag reduction results.     

Development characteristics of fluctuating velocity field of drag-reducing surfactant solution flow in a duct

Development behavior of the fluctuating velocity of surfactant solution in a duct has been studied experimentally. The concentration of surfactants was kept constant at 1,000 ppm, mean velocity at 0.78 m/s and fluid temperature at 15 °C. Using laser Doppler velocimetry, the fluctuating streamwise velocity distributions at six cross sections, which ranged from 14 to 112 times of hydraulic diameter of the duct, were measured. From the results, the fluctuating structures of surfactant solution flow are observed to have structures different from that of turbulent water flow in the developing field. The wavelet analysis reveals that the high-level fluctuation of surfactant solution flow is characterized by periodicity rather than irregularity around the position where the fluctuation intensity takes a peak value and that the period and the scale of periodic flow structures are related to the relaxation times of the fluid. This indicates that the high-level fluctuation is deeply related to the elastic instability and has a different generation mechanism from that of turbulence observed in a Newtonian turbulent flow.     

Friction and heat transfer in drag-reducing surfactant solution flow through curved pipes and elbows

A study of drag-reducing flow in curved pipes was conducted. In contrast to earlier studies we show that if we use a modified definition of drag reduction that includes only the turbulence effects, we observe indeed the same level of drag reduction in both coiled and straight pipes. More complex results showing reduced drag reduction compared to curved pipes were achieved with elbows. Two elbows of different size and type were tested in turbulent flow of both water and drag-reducing surfactant solution. A more elaborate analysis was conducted for a half-inch threaded elbow with a ratio of curvature radius to diameter of 1.2. The pressure drop and heat transfer were measured in a section downstream from the elbow over a distance of x/D = 130 in order to investigate the hydrodynamic and thermal developments of the flow. The pressure drop coefficient of the elbow was calculated for water and a surfactant solution, based on the total increase in pressure drop in the system due to the presence of the elbow. For a larger welded elbow of 6″ diameter some drag reduction was measured for the surfactant solution.     

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.     

New Answers on the Interaction Between Polymers and Vortices in Turbulent Flows

Numerical data of polymer drag reduced flows is interpreted in terms of modification of near-wall coherent structures. The originality of the method is based on numerical experiments in which boundary conditions or the governing equations are modified in a controlled manner to isolate certain features of the interaction between polymers and turbulence. As a result, polymers are shown to reduce drag by damping near-wall vortices and sustain turbulence by injecting energy onto the streamwise velocity component in the very near-wall region.     

Characterization of micellar structure dynamics for a drag-reducing surfactant solution under shear: normal stress studies and flow geometry effects

Some surfactant solutions have been observed to exhibit a strong drag reduction behavior in turbulent flow. This effect is generally believed to result from the formation of large cylindrical micelles or micellar structures. To characterize and understand better these fluids, we have studied the transient rheological properties of an efficient drag-reducing aqueous solution: tris (2-hydroxyethyl) tallowalkyl ammonium acetate (TTAA) with added sodium salicylate (NaSal) as counter ion. For a 5/5 mM equimolar TTAA/NaSal solution, there is no measurable first normal stress difference (N ₁) immediately after the inception of shear, but N ₁ begins to increase after a well-defined induction time — presumably as shear-induced structures (SIS) are formed — and it finally reaches a fluctuating plateau region where its average value is two orders of magnitude larger than that of the shear stress. The SIS buildup times obtained by first normal stress measurements were approximately inversely proportional to the shear rate, which is consistent with a kinetic process during which individual micelles are incorporated through shear into large micellar structures. The SIS buildup after a strong preshear and the relaxation processes after flow cessation were also studied and quantified with first normal stress difference measurements. The SIS buildup times and final state were also found to be highly dependent on flow geometry. With an increase in gap between parallel plates, for example, the SIS buildup times decreased, whereas the plateau viscosity increased.     

Effects of polymer additives on the large scale structure of a two-dimensional submerged jet

The large structures on the boundary of a two-dimensional submerged jet were examined. In the first experiment, a jet of a 100 ppm aqueous polyacrylamid solution was fed into pure water (P/W). In the second experiment, the jet as well as the ambient fluid were pure water (W/W), whereas in the third experiment, polymer solution (W/P) was in the ambient only.The large structures were found to be of the same size for all three experiments. The spread-angle was lowest for P/W and highest for W/P. These effects are attributed to a strong flow which stretches the polymer molecules in the mixing layer.     

The influence of the injection system on drag reduction

The influence of the injection system for centerline injected polymer solutions (threads) on drag reduction in a turbulent pipe flow was studied using injectors of different length and grids. Compared with a short injector, the long injector showed a different behavior: the drag reduction was lower and its onset point was shifted to higher Reynolds numbers.The velocity profiles for the polymer-phase and the water-phase were measured simultaneously with a combination of laser-Doppler-velocimetry LDV and laser-induced fluorescence LIE It was found that the analysis of the LDV measurements with respect to the difference in velocity between the polymer-phase and the water-phase can give information about the mixing between both phases. For a Reynolds number of 30000 the difference between the phases is comparatively large for low drag reduction and very small for high drag reduction. The results indicate that the drag reduction achieved by injecting a concentrated polymer solution is mainly caused by a mixing process between polymer and water.     

On the boundary layer/riblets interaction mechanisms and the prediction of turbulent drag reduction

Turbulent drag reduction experienced by ribletted surfaces is the result of both (1) the interaction between riblet peaks and the coherent structures that characterize turbulent near-wall flows, and (2) the laminar sublayer flow modifications caused by the riblet shape, which can balance, under appropriate conditions, the drag penalty due to the increased wetted surface. The latter “viscous” mechanism is investigated by means of an analytical model of the laminar sublayer, which removes geometrical restrictions and allows us to take into account “real” shapes of riblet contours, affected by manufacturing inaccuracies, and to compute even for such cases a parameter, called protrusion height, related to the longitudinal mean flow. By considering real geometries, riblet effectiveness is clearly shown to be related to the difference between the longitudinal and the transversal protrusion heights. A simple method for the prediction of the performances of ribletted surfaces is then devised. The predicted and measured drag reduction data, for different riblet geometries and flow characteristics, are in close agreement with each other. The soundness of the physical interpretation underlying this prediction method is consequently confirmed.     

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.