Stokes Mechanism of Drag Reduction

2005 ◽  
Vol 73 (3) ◽  
pp. 483-489 ◽  
Author(s):  
Promode R. Bandyopadhyay

The mechanism of drag reduction due to spanwise wall oscillation in a turbulent boundary layer is considered. Published measurements and simulation data are analyzed in light of Stokes’ second problem. A kinematic vorticity reorientation hypothesis of drag reduction is first developed. It is shown that spanwise oscillation seeds the near-wall region with oblique and skewed Stokes vorticity waves. They are attached to the wall and gradually align to the freestream direction away from it. The resulting Stokes layer has an attenuated nature compared to its laminar counterpart. The attenuation factor increases in the buffer and viscous sublayer as the wall is approached. The mean velocity profile at the condition of maximum drag reduction is similar to that due to polymer. The final mean state of maximum drag reduction due to turbulence suppression appears to be universal in nature. Finally, it is shown that the proposed kinematic drag reduction hypothesis describes the measurements significantly better than what current direct numerical simulation does.

Author(s):  
Promode R. Bandyopadhyay

The mechanism of drag reduction due to spanwise wall oscillation in a turbulent boundary layer is considered. Published measurements and simulation data are analyzed in light of Stokes’ second problem. A kinematic vorticity reorientation hypothesis of drag reduction is first developed. It is shown that spanwise oscillation seeds the near-wall region with oblique and skewed Stokes vorticity waves. They are attached to the wall and gradually align to the freestream direction away from it. The resulting Stokes’ layer has an attenuated nature compared to its laminar counterpart. The attenuation factor increases in the buffer and viscous sublayer as the wall is approached. The mean velocity profile at the condition of maximum drag reduction is similar to that due to polymer. The final mean state of maximum drag reduction due to turbulence suppression appears to be universal in nature. Finally, it is shown that the proposed kinematic drag reduction hypothesis describes the measurements significantly better than what current Direct Numerical Simulation does.


1982 ◽  
Vol 119 ◽  
pp. 423-441 ◽  
Author(s):  
M. A. Goldshtik ◽  
V. V. Zametalin ◽  
V. N. Shtern

We propose a simplified theory of a viscous layer in near-wall turbulent flow that determines the mean-velocity profile and integral characteristics of velocity fluctuations. The theory is based on the concepts resulting from the experimental data implying a relatively simple almost-ordered structure of fluctuations in close proximity to the wall. On the basis of data on the greatest contribution to transfer processes made by the part of the spectrum associated with the main size of the observed structures, the turbulent fluctuations are simulated by a three-dimensional running wave whose parameters are found from the problem solution. Mathematically the problem reduces to the solution of linearized Navier-Stokes equations. The no-slip condition is satisfied on the wall, whereas on the outer boundary of a viscous layer the conditions of smooth conjunction with the asymptotic shape of velocity and fluctuation-energy profiles resulting from the dimensional analysis are satisfied. The formulation of the problem is completed by the requirement of maximum curvature of the mean-velocity profile on the outer boundary applied from stability considerations.The solution of the problem does not require any quantitative empirical data, although the conditions of conjunction were formulated according to the well-known concepts obtained experimentally. As a result, the near-wall law for the averaged velocity has been calculated theoretically and is in good agreement with experiment, and the characteristic scales for fluctuations have also been determined. The developed theory is applied to turbulent-flow calculations in Maxwell and Oldroyd media. The elastic properties of fluids are shown to lead to near-wall region reconstruction and its associated drag reduction, as is the case in turbulent flows of dilute polymer solutions. This theory accounts for several features typical of the Toms effect, such as the threshold character of the effect and the decrease in the normal fluctuating velocity. The analysis of the near-wall Oldroyd fluid flow permits us to elucidate several new aspects of the drag-reduction effect. It has been established that the Toms effect does not always result in thickening of the viscous sublayer; on the contrary, the most intense drag reduction takes place without thickening in the viscous sublayer.


1971 ◽  
Vol 45 (3) ◽  
pp. 417-440 ◽  
Author(s):  
P. S. Virk

Further evidence of a universal maximum drag reduction asymptote is presented. In the elastic sublayer model, inferred therefrom, the mean velocity profile during drag reduction is approximated by three zones: the usual viscous sublayer, an elastic sublayer where the mixing-length constant is derived from the maximum drag reduction asymptote, and an outermost region with Newtonian mixinglength constant. Upon integration the model yields a friction factor relation, parametric in elastic sublayer thickness, which properly reproduces the known features of turbulent dilute polymer solution flow. The dependence of elastic sublayer thickness upon flow and polymeric parameters is inferred from experimental data revealing two hitherto unknown relationships: namely that on Prandtl co-ordinates, 1/f½vs. log Re f½ the difference in slope between a polymer solution and solvent is proportional to the square root of molar concentration and to the three-halves power of backbone chain links in the macromolecule. The proportionality constant in the preceding relationship is approximately the same for several different polymer species of carbon-carbon or similar skeletal structure in various thin solvents; there is an indication that this constant further depends upon the product of solvent viscosity times the cube of the effective bond length per chain link of the polymer species. Some recent results regarding the onset of drag reduction are also summarized.


1970 ◽  
Vol 37 (2) ◽  
pp. 488-493 ◽  
Author(s):  
P. S. Virk ◽  
H. S. Mickley ◽  
K. A. Smith

The maximum drag reduction in turbulent pipe flow of dilute polymer solutions is ultimately limited by a unique asymptote described by the experimental correlation: f−1/2=19.0log10(NRef1/2)−32.4 The semilogarithmic mean velocity profile corresponding to and inferred from this ultimate asymptote has a mixing-length constant of 0.085 and shares a trisection (at y+ ∼ 12) with the Newtonian viscous sublayer and law of the wall. Experimental mean velocity profiles taken during drag reduction lie in the region bounded by the inferred ultimate profile and the Newtonian law of the wall. At low drag reductions the experimental profiles are well correlated by an “effective slip” model but this fails progressively with increasing drag reduction. Based on the foregoing a three-zone scheme is proposed to model the mean flow structure during drag reduction. In this the mean velocity profile segments are (a) a viscous sublayer, akin to Newtonian, (b) an interactive zone, characteristic of drag reduction, in which the ultimate profile is followed, and (c) a turbulent core in which the Newtonian mixing-length constant applies. The proposed model is consistent with experimental observations and reduces satisfactorily to the Taylor-Prandtl scheme and the ultimate profile, respectively, at the limits of zero and maximum drag reductions.


1974 ◽  
Vol 65 (3) ◽  
pp. 439-459 ◽  
Author(s):  
Helmut Eckelmann

Hot-film anemometer measurements have been carried out in a fully developed turbulent channel flow. An oil channel with a thick viscous sublayer was used, which permitted measurements very close to the wall. In the viscous sublayer between y+ ≃ 0·1 and y+ = 5, the streamwise velocity fluctuations decreased at a higher rate than the mean velocity; in the region y+ [lsim ] 0·1, these fluctuations vanished at the same rate as the mean velocity.The streamwise velocity fluctuations u observed in the viscous sublayer and the fluctuations (∂u/∂y)0 of the gradient at the wall were almost identical in form, but the fluctuations of the gradient at the wall were found to lag behind the velocity fluctuations with a lag time proportional to the distance from the wall. Probability density distributions of the streamwise velocity fluctuations were measured. Furthermore, measurements of the skewness and flatness factors made by Kreplin (1973) in the same flow channel are discussed. Measurements of the normal velocity fluctuations v at the wall and of the instantaneous Reynolds stress −ρuv were also made. Periods of quiescence in the − ρuv signal were observed in the viscous sublayer as well as very active periods where ratios of peak to mean values as high as 30:1 occurred.


2005 ◽  
Vol 15 (6) ◽  
pp. 370-389 ◽  
Author(s):  
Anshuman Roy ◽  
Ronald G. Larson

Abstract We present a one-parameter model that fits quantitatively the mean velocity profiles from experiments and numerical simulations of drag-reduced wall-bounded flows of dilute solutions of polymers and non-Brownian fibers in the low and modest drag reduction regime. The model is based on a viscous mechanism of drag reduction, in which either extended polymers or non-Brownian fibers increase the extensional viscosity of the fluid and thereby suppress both small and large turbulent eddies and reduce momentum transfer to the wall, resulting in drag reduction. Our model provides a rheological interpretation of the upward parallel shift S+ in the mean velocity profile upon addition of polymer, observed by Virk. We show that Virk’s correlations for the dependence on polymer molecular weight and concentration of the onset wall shear stress and slope increment on the Prandtl-Karman plot can be translated to two dimensionless numbers, namely an onset Weissenberg number and an asymptotic Trouton ratio of maximum extensional viscosity to zero-shear viscosity. We believe that our model, while simple, captures the essential features of drag reduction that are universal to flexible polymers and fibers, and, unlike the Virk phenomenology, can easily be extended to flows with inhomogeneous polymer or fiber concentration fields.


Author(s):  
Dongmei Zhou ◽  
Kenneth S. Ball

This paper has two objectives, (1) to examine the effects of spatial resolution, (2) to examine the effects of computational box size, upon turbulence statistics and the amount of drag reduction with and without the control scheme of wall oscillation. Direct numerical simulation (DNS) of the fully developed turbulent channel flow was performed at Reynolds number of 200 based on the wall-shear velocity and the channel half-width by using spectral methods. For the first objective, four different grids were applied to the same computational domain and the biggest impact was observed on the logarithmic law of mean velocity profiles and on the amount of drag reduction with 28.3% for the coarsest mesh and 35.4% for the finest mesh. Other turbulence features such as RMS velocity fluctuations, RMS vorticity fluctuations, and bursting events were either overpredicted or underpredicted through coarse grids. For the second objective, two different minimal channels and one natural full channel were studied and 3% drag reduction difference was observed between the smallest minimal channel of 39.1% and the natural full channel of 36.2%. In the near-wall region, however, the minimal channel flow did not exhibit significant difference in the mean velocity profiles and other lower-order statistics. Finally, from this systematical study, it showed that the accuracy of DNS depends more on the spanwise resolution, and it also confirmed that a minimal channel model is able to catch key structures of turbulence in the near-wall region but is much less expensive.


Author(s):  
Redha Wahidi ◽  
Walid Chakroun ◽  
Sami Al-Fahad

Turbulent boundary layer flows over a flat plate with multiple transverse square grooves spaced 10 element widths apart were investigated. Mean velocity profiles, turbulence intensity profiles, and the distributions of the skin-friction coefficients (Cf) and the integral parameters are presented for two grooved walls. The two transverse square groove sizes investigated are 5mm and 2.5mm. Laser-Doppler Anemometer (LDA) was used for the mean velocity and turbulence intensity measurements. The skin-friction coefficient was determined from the gradient of the mean velocity profiles in the viscous sublayer. Distribution of Cf in the first grooved-wall case (5mm) shows that Cf overshoots downstream of the groove and then oscillates within the uncertainty range and never shows the expected undershoot in Cf. The same overshoot is seen in the second grooved-wall case (2.5mm), however, Cf continues to oscillate above the uncertainty range and never returns to the smooth-wall value. The mean velocity profiles clearly represent the behavior of Cf where a downward shift is seen in the Cf overshoot region and no upward shift is seen in these profiles. The results show that the smaller grooves exhibit larger effects on Cf, however, the boundary layer responses to these effects in a slower rate than to those of the larger grooves.


1976 ◽  
Vol 77 (3) ◽  
pp. 473-497 ◽  
Author(s):  
L. J. S. Bradbury

This paper describes an investigation into the response of both the pulsed-wire anemometer and the hot-wire anemometer in a highly turbulent flow. The first part of the paper is concerned with a theoretical study of some aspects of the response of these instruments in a highly turbulent flow. It is shown that, under normal operating conditions, the pulsed-wire anemometer should give mean velocity and longitudinal turbulent intensity estimates to an accuracy of better than 10% without any restriction on turbulence level. However, to attain this accuracy in measurements of turbulent intensities normal to the mean flow direction, there is a lower limit on the turbulent intensity of about 50%. An analysis is then carried out of the behaviour of the hot-wire anemometer in a highly turbulent flow. It is found that the large errors that are known to develop are very sensitive to the precise structure of the turbulence, so that even qualitative use of hot-wire data in such flows is not feasible. Some brief comments on the possibility of improving the accuracy of the hot-wire anemometer are then given.The second half of the paper describes some comparative measurements in the highly turbulent flow immediately downstream of a normal flat plate. It is shown that, although it is not possible to interpret the hot-wire results on their own, it is possible to calculate the hot-wire response with a surprising degree of accuracy using the results from the pulsed-wire anemometer. This provides a rather indirect but none the less welcome check on the accuracy of the pulsed-wire results, which, in this very highly turbulent flow, have a certain interest in their own right.


Sign in / Sign up

Export Citation Format

Share Document