Modelling the time-dependent flow over riblets in the viscous wall region

The flow over drag reducing riblets is examined computationally using a time-dependent model of the viscous wall region. The flow at the upper bound of the computational domain (y ⁺?40) is described using a streamwise eddy model consisting of two scales. A control-volume finite-element method utilizing triangular meshes is used to exactly fit the riblet cross-sectional geometry. Observations of the transient flow conditions suggest that the riblets limit the lateral spread of fluid inrushes towards the wall and retain low momentum fluid in the riblet valleys effectively isolating much of the wall from such inrushes. The generation of intermittent secondary vortices within the riblet valleys also occurs; however, these appear to be quite weak and fairly short-lived.     

Modelling the time dependent flow over riblets in the viscous wall region

The flow over riblets is examined computationally using a time dependent model of the viscous wall region. This 2 1/2 D model, developed by Hatziavramidis and Hanratty (1979) and modified by Nikolaides (1984) and Chapman and Kuhn (1981, 1986) assumes homogeneity in the streamwise direction so that the flow is solved only in the cross-sectional plane. The flow at the upper boundary of the computational domain (y ⁺ 40) is described using a streamwise eddy model consisting of two scales, one of the streak spacing (⁺ 100), which dominates vertical momentum transport, and a larger scale that accounts for the influence of large outer flow eddies.The protrusion height concept (Bechert and Bartenwerfer, 1989) is used to define ay ⁺=0 location for surfaces with riblets. A control volume finite element method utilizing triangular meshes is used to exactly fit the riblet cross-sectional geometry. Results obtained using fairly large riblets compare well with the limited experimental evidence available. Observations of the transient flow suggest that the riblets interact with the near-wall streamwise vortices, weakening them by the generation of intermittent secondary vortices within the riblet valleys. The riblets also appear to limit the lateral spread of inrushes towards the wall and retain low momentum fluid in the riblet valleys effectively isolating much of the wall from such inrushes.     

Influence of wall riblets on diffuser flow

Diffuser flows have been investigated with different riblet combinations both experimentally and numerically. Wall pressures and velocity profiles have been measured throughout a large working domain. In addition numerical calculations have been performed by accounting for the strong coupling between the inviscid core flow and the viscous boundary layer flows by using a simple model based on the negative roughness for the riblet walls.Nomenclature c _f skin friction coefficient - p local static pressure - x streamwise distance - y distance normal to the wall - u friction velocity - A diffuser cross section - H shape factor - U _e fluid velocity in inviscid flow - W ₁ inlet height of the diffuser - * boundary layer displacement thickness - density of working fluid Institut de Machines Hydrauliques et de Mecanique des Fluides     

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.     

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.     

Coherent structures and riblets

This work deals with the effect of the riblets on the coherent structures near the wall. The emphasis is put on the genesis of the quasi-streamwise vortices in the presence of the riblets. The quasi-streamwise vortices regenerate by the tilting of wall normal vorticity induced by prevailing structures. This requires a mechanism which leads to a temporal streamwise dependence near the elongated flow structures and to a subsequent formation of new wall normal vorticity. It is suggested here that the action of existing quasi-streamwise vortices on the sidewalls of wall normal vorticity may create a local, streamwise dependent spanwise velocity and therefore, a secondary wall normal vorticity field. A preliminary analysis of the set-up and the time and space development of this secondary three-dimensional flow associated with the regeneration mechanism, is given. An attempt is made, in order to explain the drag reduction performed by the riblets through an intermittent model, based on the protrusion height. Logical estimates of the amount of drag reduction are obtained. The differences between the mechanism suggested here and those based on forced control experiments are also discussed.     

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

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

Flow of power-law fluids over a moving wedge surface with wall mass injection

The boundary layer problem of a power-law fluid flow with fluid injection on a wedge whose surface is moving with a constant velocity in the opposite direction to that of the uniform mainstream is analyzed. The free stream velocity, the injection velocity at the surface, moving velocity of the wedge surface, the wedge angle and the power law index of non-Newtonian fluid are assumed variables. The fourth order Runge–Kutta method modified by Gill is used to solve the non-dimensional boundary layer equations for non-Newtonian flow field. Without fluid injection, for every angle of wedge β, a limiting value for velocity ratio λ _cr (velocity of the wedge surface/velocity of the uniform flow) is found for each power-law index n. The value of λ _cr increases with the increasing wedge angle β. The value of wedge angle also restricts the physical characteristics of the fluid to be used. The effects of the different parameters on velocity profile and on skin friction are studied and the drag reduction is discussed. In case of C = 2.5 and velocity ratio λ = 0.2 for wedge angle β = 0.5 with the fluid with power law-index n = 0.5, 48.8% drag reduction is obtained.     

Experimental study of developing turbulent flow and heat transfer in ribbed convergent/divergent square ducts

The local heat transfer and pressure drop characteristics of developing turbulent flows of air in three stationary ribbed square ducts have been investigated experimentally. These are: ribbed square duct with constant cross-section (straight duct), ribbed divergent square duct and ribbed convergent square duct. The convergent/divergent duct has an inclination angle of 1°. The measurement was conducted within the range of Reynolds numbers from 10 000 to 77 000. The heat transfer performance of the divergent/convergent ducts is compared with the ribbed straight duct under three constraints: identical mass flow rate, identical pumping power and identical pressure drop. Because of the streamwise flow acceleration or deceleration, the local heat transfer characteristics of the divergent and convergent ducts are quite different from those of the straight duct. In the straight duct, the fluid flow and heat transfer become fully developed after 2–3 ribs, while in the divergent and convergent ducts there is no such trend. The comparison shows that among the three ducts, the divergent duct has the highest heat transfer performance, the convergent duct has the lowest, while the straight duct locates somewhere in between.     

On the relationship between drag and vertical velocity fluctuations in flow over riblets and liquid infused surfaces

Direct numerical simulations (DNS) of flow over triangular and rectangular riblets in a wide range of size and Reynolds number have been carried out. The flow within the grooves is directly resolved by exploiting the immersed-boundary method. It is found that the drag reduction property is primarily associated with the capability of inhibiting vertical velocity fluctuations at the plane of the crests, as in liquid-infused surfaces (LIS) devices. This is mimicked in DNS through artificial suppression of the vertical velocity component, which yields large drag decrease, proportionate to the riblets size. A parametrization of the drag reduction effect in terms of the vertical velocity variance is found to be quite successful in accounting for variation of the controlling parameters. A Moody-like friction diagram is thus introduced which incorporates the effect of slip velocity and a single, geometry-dependent parameter. Reduced drag-reduction efficiency of LIS-like riblets is found as compared to cases with artificially imposed slip velocity. Last, we find that simple wall models of riblets and LIS-like devices are unlikely to provide accurate prediction of the flow phenomenon, and direct resolution of flow within the grooves in necessary.     

Experimental and numerical study of developing turbulent flow and heat transfer in convergent/divergent square ducts

 The work reported in this paper is a systematic experimental and numerical study of friction and heat transfer characteristics of divergent/convergent square ducts with an inclination angle of 1^∘ in the two direction at cross section. The ratio of duct length to average hydraulic diameter is 10. For the comparison purpose, measurement and simulation are also conducted for a square duct with constant cross section area, which equals to the average cross section area of the convergent/divergent duct. In the numerical simulation the flow is modeled as being three-dimensional and fully elliptic by using the body-fitted finite volume method and the kɛ turbulence model. The uniform heat flux boundary condition is specified to simulate the electrical heating used in the experiments. The heat transfer performance of the divergent/convergent ducts is compared with the duct with uniform cross section under three constraints (identical mass flow rate, pumping power and pressure drop). The agreement of the experimental and numerical results is quite good except at the duct inlet. Results show that for the three ducts studied there is a weak secondary flow at the cross section, and the circumference distribution of the local heat transfer coefficient is not uniform, with an appreciable reduction in the four corner regions. In addition, the acceleration/deceleration caused by the cross section variation has a profound effect on the turbulent heat transfer: compared with the duct of constant cross section area, the divergent duct generally shows enhanced heat transfer behavior, while the convergent duct has an appreciable reduction in heat transfer performance. Received on 18 September 2000 / Published online: 29 November 2001     

Hybrid nanomaterial flow and heat transport in a stretchable convergent/divergent channel:a Darcy-Forchheimer model

The flow behavior in non-parallel walls is an important factor of any physical model including cavity flow and canals, which is applicable for diverging/converging channel. The present communication explains that the flow of the hybrid nanomaterial subjected to the convergent/divergent channel has non-parallel walls. It is assumed that the hybrid nanomaterial movement is in the porous region. A Darcy-Forchheimer medium of porosity is considered to interpret the porosity features. A useful simila...     

Flow and heat transfer over a generalized stretching/shrinking wall problem—Exact solutions of the Navier-Stokes equations

In this paper, we investigate the steady momentum and heat transfer of a viscous fluid flow over a stretching/shrinking sheet. Exact solutions are presented for the Navier-Stokes equations. The new solutions provide a more general formulation including the linearly stretching and shrinking wall problems as well as the asymptotic suction velocity profiles over a moving plate. Interesting non-linear phenomena are observed in the current results including both exponentially decaying solution and algebraically decaying solution, multiple solutions with infinite number of solutions for the flow field, and velocity overshoot. The energy equation ignoring viscous dissipation is solved exactly and the effects of the mass transfer parameter, the Prandtl number, and the wall stretching/shrinking strength on the temperature profiles and wall heat flux are also presented and discussed. The exact solution of this general flow configuration is a rare case for the Navier-Stokes equation.     

Effect of heat transfer on the flow of a second‐grade fluid in divergent/convergent channel

This paper studies the effects of a second‐grade fluid on the flow and heat transfer characteristics in a divergent/convergent channel. The momentum and energy equations are first given in a nondimensional form and then solved analytically using the method of homotopy analysis method. Convergence of derived series solutions is shown. Graphical results for the velocity and the temperatures are presented and discussed for various emerging parameters. Copyright © 2009 John Wiley & Sons, Ltd.     

Turbulent wall pressure fluctuations over wavy surfaces

Measurements of the spectral characteristics of the wall pressure fluctuations produced by a turbulent boundary layer flow over solid sinusoidal surfaces of moderate wave amplitude to wave-length ratios have been obtained. The wave amplitudes were sufficiently small so that the flow remained attatched. The results show that the root mean square pressure level reaches a maximum on the adverse pressure gradient side of the wave at a position somewhat before the trough. Spectral analysis of the pressure fluctuations in narrow frequency bands reveals considerable differences in low and high frequency behavior. At low frequencies, the peak fluctuation amplitude was found at the trough whereas at high frequencies, the peak occurs just after the crest and a minimum is found at the trough. Pressure fluctuations having streamwise correlation lengths on the order of or larger than the wavelength of the surface do not return to their equilibrium (crest) amplitudes as they travel the length of a wave. Pressure fluctuations having streamwise correlation lengths about one order of magnitude less than a wavelength return exactly to their equilibrium amplitudes. Two-point correlation measurements show a decrease in longitudinal coherence on the adverse pressure gradient side of the wave at low frequencies and a considerable increase over a broad frequency range on the positive pressure gradient side. No change is found in the lateral coherence.List of symbols C _f skin friction coefficient - C _p pressure coefficient - C _n Fourier amplitudes of the pressure coefficient - C _dp pressure drag coefficient - d pinhole diameter - f frequency - h half the crest to trough distance - h ₊ nondimensional wave amplitude = - k _n wavenumber = - k fundamental wavenumber = - l _p pressure correlation length - p _s mean surface pressure - P ambient pressure - p fluctuating pressure - p ² mean square pressure - q dynamic head = 1/2 U ² - R space-time correlation - P Reynolds number based on wavelength = - R Reynolds number based on momentum thickness = - t time - R free stream velocity - U mean streamwise velocity - U _e streamwise velocity at the edge of the boundary layer - u ^* friction velocity = - x streamwise coordinate - y wall-normal coordinate - z spanwise coordinate - ⁺ non-dimensional wavelength = ^*) - phase of the cross-spectral density - ^* boundary layer displacement thickness - _long longitudinal coherency - _lat lateral coherency - wavelength of wavy surface - v kinematic viscosity - radian frequency = 2 f - spectral or cross-spectral density - _n phase of the Fourier series - density - time delay - _w wall shear stress - boundary layer momentum thickness