The effect of riblets on laminar to turbulent transition

Experiments conducted on the effect of riblets on the laminar-to-turbulent transition of a flat plate in a water tunnel are reported. Transition was determined using a Laser Doppler Velocimeter (LDV). A smooth reference surface was compared to five riblet surfaces for a range of Reynolds numbers. Smooth surface transition Reynolds number was about 2.75 × 10⁶. All of the five tested riblet surfaces had lower transition Reynolds numbers. A critical roughness Reynolds number of about 6 was determined for one of the riblet surfaces. This is much lower than the generally accepted value of 25, considered safe for distributed roughness.     

Riblet flow calculation with a low Reynolds number κ - ε model

A low Reynolds number κ - ε model has been used to calculate the turbulent boundary layer over riblets. Calculated mean velocity, Reynolds shear stress and kinetic energy distributions are generally in good agreement with available experimental data. The comparison between these distributions and those in a corner flow points to a significant difference between the two flows and the unlikelihood of counter-rotating vortices within the riblet grooves. One shortcoming of the present κ - ε model is the relatively slow return to a two-dimensional turbulence state as the distance from the riblet surface increases.     

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

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

固体颗粒对沟槽湍流边界层影响的实验研究

本文采用粒子图像测速技术(particles image velocimetry, PIV)研究固体颗粒对放置在平板湍流边界层中的平壁和沟槽壁面减阻效果的影响. 实验对清水和加入粒径为155 μm聚苯乙烯颗粒的流法向二维速度场信息进行采集, 对不同工况下的平均速度剖面、雷诺应力和湍流度等统计量进行对比, 分析流体在边界层中的行为. 运用空间局部平均结构函数提取了不同工况湍流边界层喷射?扫掠行为的空间拓扑结构并进行比较. 结果发现, 在不同的壁面条件下, 粒子加入后的对数律区中无量纲速度均略大于清水组, 雷诺切应力有所降低, 湍流度有所减弱. 对于不同流场速度下的沟槽而言, 颗粒的加入均降低了壁面附近的阻力, 而颗粒单独作用于光滑壁面的减阻效果并不明显. 加入粒子后的相干结构数目有所增加, 法向脉动速度下降. 沟槽壁面附近的相干结构数目有所增加, 法向脉动速度在自由来流速度较大时有所上升, 在速度较小时有所下降. 这表明不同减阻状况下的沟槽均能将大涡破碎成更多的涡, 并且粒子的加入强化了这种破碎作用.      

Riblet drag reduction in mild adverse pressure gradients: A numerical investigation

Riblet films are a passive method of turbulent boundary layer control that can reduce viscous drag. They have been studied with great detail for over 30 years. Although common riblet applications include flows with Adverse Pressure Gradients (APG), nearly all research thus far has been performed in channel flows. Recent research has provided motivation to study riblets in more complicated turbulent flows with claims that riblet drag reduction can double in mild APG common to airfoils at moderate angles of attack. Therefore, in this study, we compare drag reduction by scalloped riblet films between riblets in a zero pressure gradient and those in a mild APG using high-resolution large eddy simulations. In order to gain a fundamental understanding of the relationship between drag reduction and pressure gradient, we simulated several different riblet sizes that encompassed a broad range of s⁺ (riblet width in wall units), similarly to many previously published experimental studies. We found that there was only a slight improvement in drag reduction for riblets in the mild APG. We also observed that peak values of streamwise turbulence intensity, turbulent kinetic energy, and streamwise vorticity scale with riblet width. Primary Reynolds shear stresses and turbulence kinetic energy production however scale with the ability of the riblet to reduce skin-friction.     

脊状表面减阻机理研究

针对脊状表面流场的特点,通过实验测量和数值模拟的方法对脊状表面微观流场进行了深入研究,获得了脊状表面湍流边界层的时均速度分布曲线、湍流度分布曲线和微观流场结构.为了得到脊状结构对壁面物性的影响,对脊状表面进行了疏水性测试,获得了液滴在脊状表面上的表观接触角,并通过水洞试验验证了脊状表面的减阻效果.研究表明,与光滑表面相比,脊状表面微观流场结构中存在"二次涡",近壁区的黏性底层厚度比平板的要厚得多,湍流度显著降低,且脊状表面表现出明显的疏水性.由此提出了基于壁面隔离效应、增大湍流阻尼效应和改变壁面物性效应的减阻机理.     

脊状表面减阻特性的风洞试验研究

利用热线风速仪,对光滑表面和多个脊状表面在低速风洞中进行了表面流场测试。基于测得的边界层速度分布数据,利用对数律区速度分布公式,编程分别计算出光滑表面和脊状表面的壁面摩擦速度和虚拟原点。研究发现,脊状表面最大减阻量达13．5％;有减阻效果的脊状表面使边界层速度曲线上移、湍流强度下降;与光滑表面相比,脊状表面的位移厚度和动量损失厚度明显减小,也表明脊状表面具有减阻效果;位移厚度和动量损失厚度减少量随槽间距s^＋的增加呈现先变大后变小的趋势,在S^＋=12时达到最大。     

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     

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.     

Large-scale spanwise periodicity in a turbulent boundary layer induced by highly ordered and directional surface roughness

The effect of converging–diverging riblet-type surface roughness (riblets arranged in a ‘herringbone’ pattern) are investigated experimentally in a zero pressure gradient turbulent boundary layer. For this initial parametric investigation three different parameters of the surface roughness are analysed in detail; the converging–diverging riblet yaw angle α, the streamwise fetch or development length over the rough surface F_x and the viscous-scaled riblet height h⁺. It is observed that this highly directional surface roughness pattern induces a large-scale spanwise periodicity onto the boundary layer, resulting in a pronounced spanwise modification of the boundary layer thickness. Hot-wire measurements reveal that above the diverging region, the local mean velocity increases while the turbulent intensity decreases, resulting in a thinner overall boundary layer thickness in these locations. The opposite situation occurs over the converging region, where the local mean velocity is decreased and the turbulent intensity increases, producing a locally thicker boundary layer. Increasing the converging–diverging angle or the viscous-scaled riblet height results in stronger spanwise perturbations. For the strongest convergent–divergent angle, the spanwise variation of the boundary layer thickness between the diverging and converging region is almost a factor of two. Such a large variation is remarkable considering that the riblet height is only 1% of the unperturbed boundary layer thickness. Increasing the fetch seems to cause the perturbations to grow further from the surface, while the overall strength of the induced high and low speed regions remain relatively unaltered. Further analysis of the pre-multiplied energy spectra suggests that the surface roughness has modified or redistributed the largest scale energetic structures.     

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.     

Flow over convergent and divergent wall riblets

Fast swimming sharks have small riblets on their skin, which are assumed to improve the swimming performance of the fish. Fluid dynamic experiments in water as well as in air confirm this assumption. With riblet surfaces as compared to smooth surfaces, drag reductions up to about 10% were measured. The overall riblet pattern on sharks shows parallel riblets directed from head to tail, but besides this overall pattern fast swimming sharks have also small areas with converging riblets and others with diverging riblets. In the present study the velocity field over convergent and divergent riblet patterns is investigated by hot-wire measurements in turbulent pipe flow. Significant changes in the near wall velocity field were found.     

Spalding公式在脊状表面湍壁摩擦力测量中的应用

在低速风洞中来流速度一定的情况下使用IFA300恒温热线风速仪测量了光滑表面和两种不同尺寸的脊状表面湍流边界层平均速度分布剖面,并验证了试验段湍流发展的充分性;通过应用Spalding壁面公式使用最小二乘法精准拟合了实验测量的边界层内层速度分布曲线,得到了湍流边界层壁面摩擦速度并进一步求得湍流壁面摩擦应力,较准确地计算出脊状表面的虚拟原点位置,并通过与对数律公式拟合结果比较分析,证实了该方法更加准确有效. 最后分别计算了3种实验模型的湍流边界层动量损失厚度. 通过对比脊状表面与光滑表面动量损失厚度和壁面摩擦应力,反映了动量损失厚度的大小与壁面摩擦应力的大小具有一致性,充分证实了脊状表面在湍流中具有一定的减阻效果.      

Natural Low-Frequency Riblets

A turbulent flow at modest Reynolds numbers (26,000–32,000, as defined in the text) over flat erodible sandy beds was investigated. The beds were composed of sand grains of uniform size with a diameter smaller than the viscous sublayer of the flow. When a modest sand transport (1.2 × 10-3–8.5 × 10-3 g/s·cm) was present, the near-wall coherent structures of the flow produced a bed form with a streaky longitudinal pattern. This pattern had a lateral spacing of about 100 viscous units and had less drag than a flat wall.The mechanism causing the drag reduction is due to a low-frequency riblet pattern which allows the flow to sustain higher wall shear stresses without additional erosion. That means during the riblet regime the sediment transport remains quasi-constant.The formation of the sand streaks as well as the drag reducing mechanism are discussed.     

Explicit algebraic Reynolds stress model for shock-dominated flows

Shock waves drastically alter the nature of Reynolds stresses in a turbulent flow, and conventional turbulence models cannot reproduce this effect. In the present study, we employ explicit algebraic Reynolds stress model (EARSM) to predict the Reynolds stress anisotropy generated by a shockwave. The model by Wallin and Johansson (2000) is used as the baseline model. It is found to over-predict the post-shock Reynolds stresses in canonical shock turbulence interaction. The budget of the transport equation of Reynolds stresses computed using linear interaction analysis shows that the unsteady shock distortion mechanism and the pressure–velocity correlations are important. We propose improvement to the baseline model using linear interaction analysis results and redistribute the turbulent kinetic energy between the principle Reynolds stresses. The new model matches DNS data for the amplification of Reynolds stresses across the shock and their post-shock evolution, for a range of Mach numbers. It is applied to oblique shock/boundary-layer interaction at Mach 5. Significant improvements are observed in predicting surface pressure and skin friction coefficient, with respect to experimental measurements.     

Parametric study of separation and transition characteristics over an airfoil at low Reynolds numbers

Time-resolved surface pressure measurements are used to experimentally investigate characteristics of separation and transition over a NACA 0018 airfoil for the relatively wide range of chord Reynolds numbers from 50,000 to 250,000 and angles of attack from 0° to 21°. The results provide a comprehensive data set of characteristic parameters for separated shear layer development and reveal important dependencies of these quantities on flow conditions. Mean surface pressure measurements are used to explore the variation in separation bubble position, edge velocity in the separated shear layer, and lift coefficients with angle of attack and Reynolds number. Consistent with previous studies, the separation bubble is found to move upstream and decrease in length as the Reynolds number and angle of attack increase. Above a certain angle of attack, the proximity of the separation bubble to the location of the suction peak results in a reduced lift slope compared to that observed at lower angles. Simultaneous measurements of the time-varying component of surface pressure at various spatial locations on the model are used to estimate the frequency of shear layer instability, maximum root-mean-square (RMS) surface pressure, spatial amplification rates of RMS surface pressure, and convection speeds of the pressure fluctuations in the separation bubble. A power-law correlation between the shear layer instability frequency and Reynolds number is shown to provide an order of magnitude estimate of the central frequency of disturbance amplification for various airfoil geometries at low Reynolds numbers. Maximum RMS surface pressures are found to agree with values measured in separation bubbles over geometries other than airfoils, when normalized by the dynamic pressure based on edge velocity. Spatial amplification rates in the separation bubble increase with both Reynolds number and angle of attack, causing the accompanying decrease in separation bubble length. Values of the convection speed of pressure fluctuations in the separated shear layer are measured to be between 35 and 50% of the edge velocity, consistent with predictions of linear stability theory for separated shear layers.     

Drag measurements on planar riblet surfaces at high subsonic speeds

This paper presents the results of an investigation of riblet performance at high subsonic Mach numbers, and Reynolds numbers of about 20 000 based on the momentum thickness, in both zero and adverse pressure gradient boundary layers. The maximum length Reynolds number of the ribbed section was 3.4×10⁷ so the results were directly relevant to flight applications on the engine nacelles of civil airliners. Seven different sizes of riblets with heights h (equal to spacing s) ranging from 0.0007 (0.0178mm) to 0.006 (0.1524 mm) have been studied, covering a range of h+, s+ from 10 to 106. The maximum percentage skin friction reduction, as deduced from velocity profiles measured at the downstream end of the riblet surfaces, under nominally zero pressure gradient conditions was 5.5±1; rather less than that recorded in low speed studies, but consistent with a recent theoretical analysis of the effect of Reynolds number. The values of h+ required for maximum and zero skin friction reduction agreed closely with other data. In addition subsequent floating element drag balance measurements revealed little effect of yaw angles up to 15°, again in line with other findings, and also suggested that the extent of the initial development length on, and recovery length behind, the riblets was approximately 5. The adverse pressure gradient studies indicated that riblet performance was essentially unaffected by mild gradients (=0.25), but diminished to zero in a more severe gradient (=0.5).     

Numerical and experimental investigation of the laminar boundary layer over riblets

One of the main aims of this work is to show to what extent drag reduction in a turbulent boundary layer can be ascribed to a purely viscous effect. A numerical and experimental study is performed in a laminar boundary layer over triangular riblets. The 2-D parabolic equations of motion are integrated using an x marching method and the discretised system is solved with the MSI algorithm. The influence of the riblet geometrical parameters and of the number of grid points is studied. Measurements are carried out in a water tunnel with forward scatter and backscatter laser-Doppler velocimetry extending within the riblets. The longitudinal velocity component measurements and computations are practically identical. Numerical results presented herein show that a slight drag reduction is obtained for s/h=1.2. It appears that, as far as friction is concerned, the wetted area is not the surface to be considered. Thus, the boundary layer over riblets would behave like a boundary layer on an equivalent smooth plate located beneath the crest plane. The numerical study in terms of the riblet height h shows best results are for h tending to zero, with the ratio s/h being equal to 1.2.     

Effects of riblet bends on pipe flows

Four riblet bends were tested to investigate the effects of riblets on pipe flows including the secondary flow on the Reynolds numbers; Re _D =6×10³–4×10⁴. The pressure gradients on the smooth pipe downstream from the riblet bends were measured, and also the pressure losses of the bends only were measured. All riblet bends reduced the pressure gradient on the smooth pipe downstream from them, which means a drag reduction. Two of the riblet bends showed the maximum drag reduction of about 4 percent at Re _D = 6500; this reduction rate was significant considering the uncertainty of the present experiments. Since the pressure losses of these two riblet bends were almost identical to that of the smooth bend at Re _D = 6500, they could cause a net drag reduction of about 4 percent on the piping system including these bends at that Reynolds number. Furthermore, the velocity profiles measured by LDV indicated that the secondary flow becomes weaker downstream from the riblet bends when a drag reduction is recognized there.Nomenclature D pipe diameter - D ₀ the distance from the valley to the valley passing through the pipe center - H height of groove - P nondimensional static pressure (p/it/(U ₀ ² ):p is gauge pressure) - dP/dX nondimensional pressure gradient - Rc curvature of bend - Re _D Reynolds number based on bulk velocity and pipe diameter - s spacing of groove - U mean streamwise velocity along the horizontal diameter - U ₀ bulk velocity - V mean vertical velocity along the horizontal diameter - x streamwise direction along the pipe axis (see Fig. 1) - X nondimensionalx (=x/D) - y radial direction in the horizontal plane which is perpendicular to the plane including the bend (see Fig. 1) - yUV swirl intensity (nondimensional swirl intensity:yUV/(DU ₀ ² ))