Improved methods for thin, surface boundary layer investigations

 New techniques are developed to improve the velocity flow-field measurement capability within a free-surface boundary layer region on which progressive capillary-gravity waves are present. Due to the extremely thin but rather vortical characteristics of the aforementioned boundary layer, conventional particle image velocimetry (PIV) methods fail to estimate velocity (and vorticity) vectors at an acceptable detection rate. This failure is a direct consequence of optimal PIV parameters that are difficult to achieve in practice for such flow situations. A new technique, Sub-pattern PIV, is developed. This method has features similar to both the super-resolution PIV (Keane et al. 1995) and the particle image distortion (PID) technique (Huang et al. 1993), but is predicated upon a very differential philosophy. Another difficulty that arises in experiments to investigate surface boundary layer flows is that the oscillating and deforming air–water interface has a mirror-like behavior that affects the images, and generates very noisy data. An alternative experimental setup that utilizes the Brewster angle phenomenon is adopted and the specular effects of the free-surface are removed successfully. This Brewster angle imaging, along with the Sub-pattern PIV technique, is used for the target application – a free-surface boundary layer investigation. It proved to be very effective. The methodology of both techniques is discussed, and the modified PIV procedure is validated by numerical probabilistic simulations. Application to the capillary-gravity wave boundary layer is presented in a subsequent paper. Received: 31 July 1997/Accepted: 4 February 1998     

Velocity Measurements of a Free-Surface Turbulent Flow Penetrating a Porous Medium Composed of Uniform-Size Spheres

We present a laboratory experimental investigation of the interaction between the turbulent flow in an open channel and the turbulent flow within its very permeable bed. The bed was composed of uniform-size spheres packed in a cubic pattern. Fluid velocities were measured by Particle Image Velocimetry (PIV), which allowed us to investigate the spatial distribution of the time-averaged flow properties in the zone where they are strongly affected by the geometry of the rough bed. We investigate the effect of bed porosity on these flow properties by comparing the results of two experimental configurations: one with an impermeable bed composed of a single layer of spheres and another with a permeable bed composed of five layers. For the latter case, PIV measurements of velocities were also carried out inside two pores adjacent to the bed surface. This data provides an insight into the mechanisms of momentum transfer between the turbulent open channel flow and the turbulent flow within its very permeable bed.     

Measurements in an urban-type boundary layer

Wind tunnel measurements of the boundary layer flow over a very rough surface comprising a staggered array of cubes are presented and discussed. Attention is concentrated on the near-wall region, including the canopy region below the tops of the roughness elements. Particle image velocimetry (PIV) and laser Doppler anemometry were used to identify the dominant features of the mean and turbulent flow and these are compared with the better-known features of the flow above the roughness. Spatial correlation data, extracted from the PIV images, are used to provide information about eddy structures and it is shown that these differ in some crucial respects from those typical of more classical boundary layers. The implications of the results are discussed in terms of their relevance to flows within the urban environment.     

PIV for granular flows

 Particle image velocimetry (PIV) has been adapted for use in measuring particle displacement and velocity fields in granular flows. “Seeding” is achieved by using light and dark particles. The granular flow adjacent to a clear bounding wall is illuminated with a strobe, and the recorded images are analyzed using standard PIV techniques. The application is demonstrated by measuring convection rolls in a granular bed undergoing vertical oscillations. The PIV measured displacement is consistent with displacement of a marked layer of particles. Received: 29 January 1998/Accepted: 8 April 1999     

Heat transfer enhancement in a turbulent natural convection boundary layer along a vertical flat plate

An experimental study on heat transfer enhancement for a turbulent natural convection boundary layer in air along a vertical flat plate has been performed by inserting a long flat plate in the spanwise direction (simple heat transfer promoter) and short flat plates aligned in the spanwise direction (split heat transfer promoter) with clearances into the near-wall region of the boundary layer. For a simple heat transfer promoter, the heat transfer coefficients increase by a peak value of approximately 37% in the downstream region of the promoter compared with those in the usual turbulent natural convection boundary layer. It is found from flow visualization and simultaneous measurements of the flow and thermal fields with hot- and cold-wires that such increase of heat transfer coefficients is mainly caused by the deflection of flows toward the outer region of the boundary layer and the invasion of low-temperature fluids from the outer region to the near-wall region with large-scale vortex motions riding out the promoter. However, heat transfer coefficients for a split heat transfer promoter exhibit an increase in peak value of approximately 60% in the downstream region of the promoter. Flow visualization and PIV measurements show that such remarkable heat transfer enhancement is attributed to longitudinal vortices generated by flows passing through the clearances of the promoter in addition to large-scale vortex motions riding out the promoter. Consequently, it is concluded that heat transfer enhancement of the turbulent natural convection boundary layer can be substantially achieved in a wide area of the turbulent natural convection boundary layer by employing multiple column split heat transfer promoters. It may be expected that the heat transfer enhancement in excess of approximately 40% can be accomplished by inserting such promoters.     

Low Reynolds number axisymmetric turbulent boundary layer on a cylinder

In the present study, an axisymmetric turbulent boundary layer growing on a cylinder is investigated experimentally using hot wire anemometry. The combined effects of transverse curvature as well as low Reynolds number on the mean and turbulent flow quantities are studied. The measurements include the mean velocity, turbulence intensity, skewness and flatness factors in addition to wall shear stress. The results are presented separately for the near wall region and the outer region using dimensionless parameters suitable for each case. They are also compared with the results available in the open literature.The present investigation revealed that the mean velocity in near wall region is similar to other simple turbulent flows (flat plate boundary layer, pipe and channel flows); but it differs in the logarithmic and outer regions. Further, for dimensionless moments of higher orders, such as skewness and flatness factors, the main effects of the low Reynolds number and the transverse curvature are present in the near wall region as well as the outer region.     

SPH modelling of depth‐limited turbulent open channel flows over rough boundaries

A numerical model based on the smoothed particle hydrodynamics method is developed to simulate depth‐limited turbulent open channel flows over hydraulically rough beds. The 2D Lagrangian form of the Navier–Stokes equations is solved, in which a drag‐based formulation is used based on an effective roughness zone near the bed to account for the roughness effect of bed spheres and an improved sub‐particle‐scale model is applied to account for the effect of turbulence. The sub‐particle‐scale model is constructed based on the mixing‐length assumption rather than the standard Smagorinsky approach to compute the eddy‐viscosity. A robust in/out‐flow boundary technique is also proposed to achieve stable uniform flow conditions at the inlet and outlet boundaries where the flow characteristics are unknown. The model is applied to simulate uniform open channel flows over a rough bed composed of regular spheres and validated by experimental velocity data. To investigate the influence of the bed roughness on different flow conditions, data from 12 experimental tests with different bed slopes and uniform water depths are simulated, and a good agreement has been observed between the model and experimental results of the streamwise velocity and turbulent shear stress. This shows that both the roughness effect and flow turbulence should be addressed in order to simulate the correct mechanisms of turbulent flow over a rough bed boundary and that the presented smoothed particle hydrodynamics model accomplishes this successfully. © 2016 The Authors International Journal for Numerical Methods in Fluids Published by John Wiley & Sons Ltd     

Mobile bed thickness in skewed asymmetric oscillatory sheet flows

A new instantaneous mobile bed thickness model is presented for sediment transport in skewed asymmetric oscillatory sheet flows. The proposed model includes a basic bed load part and a suspended load part related to the Shields parameter, and takes into account the effects of mass conservation, phase-lag, and asymmetric boundary layer development, which are important in skewed asymmetric flows but usually absent in classical models. The proposed model is validated by erosion depth and sheet flow layer thickness data in both steady and unsteady flows, and applied to a new instantaneous sediment transport rate formula. With higher accuracy than classical empirical models in steady flows, the new formula can also be used for instantaneous sediment transport rate prediction in skewed asymmetric oscillatory sheet flows.     

Free-surface flow of concentrated suspensions

Free-surface flows of concentrated suspensions exhibit many interesting phenomena such as particle segregation and surface corrugation. In this work the flow structures associated with free-surface has been studied experimentally. The free-surface velocity for neutrally buoyant suspension of uniform spheres in a gravity driven inclined channel flow was determined by particle imaging velocimetry (PIV) technique. Experiments were carried out for concentrated suspensions with particle fractions ? ranging from 0.40 to 0.50. The measured velocities show blunted profile in the channel. The blunting of the velocity profile increases with the particle concentration. The rms velocity fluctuations measured at the free-surface progressively increase with particle fraction ? and are linear in shear rate γ. The surface roughness were characterized by analyzing the power spectral density of the refracted light from the free-surface. The characteristics observed are in support of earlier findings.     

Turbulent flow over a rough backward-facing step

This work characterizes the impacts of the realistic roughness due to deposition of foreign materials on the turbulent flows at surface transition from elevated rough-wall to smooth-wall. High resolution PIV measurements were performed in the streamwise-wall-normal (x–y) planes at two different spanwise positions in both smooth and rough backward-facing step flows. The experiment conditions were set at a Reynolds number of 3450 based on the free stream velocity U_∞ and the mean step height h, expansion ratio of 1.01, and the ratio of incoming boundary layer thickness to the step height, δ/h, of 8. The mean flow structures are observed to be modified by the roughness and they illustrate three-dimensional features in rough backward-facing step flows. The mean reattachment length X_r is significantly reduced by the roughness at one PIV measurement position while is slightly increased by the different roughness topography at the other measurement position. The mean velocity profiles at the reattachment point indicate that the studied roughness weakens the perturbation of the step to the incoming turbulent flow. Comparisons of Reynolds normal and shear stresses, productions of normal stresses, quadrant analysis of the instantaneous shear-stress contributing events, and mean spanwise vorticity reveal that the turbulence in the separated shear layer is reduced by the studied roughness. The results also indicate an earlier separation of the turbulent boundary layer over the current rough step, probably due to the adverse pressure gradient produced by the roughness topography even before the step.     

PIV study on a shock-induced separation in a transonic flow

A transonic interaction between a steady shock wave and a turbulent boundary layer in a Mach 1.4 channel flow is experimentally investigated by means of particle image velocimetry (PIV). In the test section, the lower wall is equipped with a contour profile shaped as a bump allowing flow separation. The transonic interaction, characterized by the existence in the outer flow of a lambda shock pattern, causes the separation of the boundary layer, and a low-speed recirculating bubble is observed downstream of the shock foot. Two-component PIV velocity measurements have been performed using an iterative gradient-based cross-correlation algorithm, providing high-speed and flexible calculations, instead of the classic multi-pass processing with FFT-based cross-correlation. The experiments are performed discussing all the hypotheses linked to the experimental set-up and the technique of investigation such as the two-dimensionality assumption of the flow, the particle response assessment, the seeding system, and the PIV correlation uncertainty. Mean velocity fields are presented for the whole interaction with particular attention for the recirculating bubble downstream of the detachment, especially in the mixing layer zone where the effects of the shear stress are most relevant. Turbulence is discussed in details, the results are compared to previous study, and new results are given for the turbulent production term and the return to isotropy mechanism. Finally, using different camera lens, a zoom in the vicinity of the wall presents mean and turbulent velocity fields for the incoming boundary layer.     

Turbulence Characteristics in Flows Over Solid and Porous Square Ribs Mounted on Porous Walls

PIV measurements have been performed for turbulent flows in a rib-mounted channel whose bottom wall is made of a porous layer. The ratio of the rib and channel heights is fixed at 0.5. The effects of the wall and rib permeability are investigated focusing on the separating and reattaching flows at the bulk Reynolds number of 10³???10⁴. Three kinds of foamed ceramics are employed as the porous media. They have the same porosity of 0.8 but each permeability is different from the others. Its normalized values by the rib height are 0.89 × 10^???4, 1.47 × 10^???4 and 3.87 × 10^???4. Two kinds of square cylinder ribs: an impermeable smooth solid rib or a permeable porous rib which is made of the same porous medium as that for the bottom wall are used. The obtained turbulent velocity fields of the solid rib flows indicate that the turbulent intensity behind the rib becomes weak and the recirculation bubble in the clear channel tends to vanish as the the wall permeability increases. In the porous rib flow, the recirculation and the reattachment point shift downstream and turbulence becomes weaker due to the bleeding flow through the rib. In the higher permeability cases, the recirculation bubble hardly exists due to the flows through not only the bottom wall but also the porous rib. From the measurements, it is suggested that in the solid rib flows, a reverse flow region exists inside the porous wall whereas in porous rib flows, such reverse flow does not exist at higher permeability.     

Experimental study on swirling flow of dilute surfactant solution with deformed free-surface

An experimental investigation was performed on a swirling flow of dilute surfactant solution with deformed free-surface in a cylindrical container driven by the constantly rotating bottom wall. The purpose of the experiment was to estimate weak viscoelasticity in the tested surfactant solutions as well as to investigate the flow characteristics. The tested fluid was an aqueous solution of CTAC (CTAC: cetyltrimethyl ammonium chloride), which is a cationic surfactant. Water, 40 ppm, 60 ppm and 200 ppm CTAC solution flows were tested at Froude numbers ranging from 2.59 to 16.3. Particle image velocimetry (PIV) was used to measure the secondary velocity field in the meridional plane. The deformed free-surface level was extracted from the PIV images. At a similar Froude number, the depth of the dip formed at the center region of the free surface was decreased for CTAC solution flow compared with water flow. The inertia-driven vortex at the up-right corner in the meridional plane becomes more and more weakened with increase of the solution concentration or viscoelasticity. Through analyzing the overall force balance compared with water flow, the first normal stress difference characterizing the viscoelasticity was estimated for the dilute CTAC solution flows. The result supports the viscoelasticity-based turbulent drag-reduction mechanism of surfactant solution flow.     

Measurements of turbulence characteristics in an open-channel flow over a transitionally-rough bed using particle image velocimetry

The particle image velocimetry technique was used to measure characteristics of a turbulent flow over a transitionally-rough fixed bed in an open-channel flow. These conditions are typical of flows encountered in sediment transport problems. Measurements obtained with this technique were used to investigate the distributions of velocities, turbulence intensities, Reynolds stress, and third- and fourth-order moments in a region above y ⁺ = 10. The present results are in good agreement to those previously obtained on smooth walls and provide further evidence that PIV can be applied successfully to investigate turbulence in open-channel flows over a rough bed.     

Reconstruction of velocity fields from wall pressure measurements in a shock wave/turbulent boundary layer interaction

We present here experimental results in a shock wave/turbulent boundary layer interaction at Mach number of 2.3 impinged by an oblique shock wave, with a deflection angle of 9.5°, as installed in the supersonic wind tunnel of the IUSTI laboratory, France. For such a shock intensity, strong unsteadiness are developing inside the separated zone involving very low frequencies associated with reflected shock motions.The present work consists in simultaneous PIV velocity fields and unsteady wall pressure measurements. The wall pressure and PIV measurements were used to characterize the pressure distribution at the wall in an axial direction, and the flow field associated. These results give access for the first time to the spatial-time correlation between wall pressure and velocity in a shock wave turbulent boundary layer interaction and show the feasibility of such coupling techniques in compressible flows. Linear Stochastic Estimation (LSE) coupled with Proper Orthogonal Decomposition (POD) has been applied to these measurements, and first results are presented here, showing the ability of these techniques to reproduce both the unsteady breathing of the recirculating bubble at low frequency and the Kelvin–Helmholtz instabilities developing at moderate frequency.     

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

Feasibility of particle image velocimetry in vegetative fire spread experiments

This study is part of an ongoing effort to improve the understanding of mechanisms that control the spread of fires with a focus on the turbulent flow modified by the flame front. A large-scale PIV system was used to measure the flow field inside and in the vicinity of a flame front spreading across a bed of fuel in an open environment. The vegetative fuel consisted of a 10-m-long and 5-m-wide bed of excelsior (1?kg/m2 fuel load) leading to a nearly 1.5-m-high flame front. The velocity field was investigated in a measurement region about 1.5?m high and 2?m long. In such a configuration, a 450-mJ laser source was used to generate the light sheet, and the flow was seeded using zirconium oxide particles (ZrO₂). The PIV measurements in the presence of flame were improved by the use of a liquid crystal shutter in front of the PIV camera, allowing very short exposure times and eliminating the flame trace in the tomographic pictures. Despite the variability of the external conditions, leading to a difficult seeding over the whole PIV area, the present study shows the feasibility of the optical method of fluid visualization in the field. The measurements of the velocity fields show some features of the dynamics of fire plumes. This preliminary study demonstrates the feasibility of the method in the open, but some strong efforts to improve the seeding of the flow must be made.     

Reynolds Stress Budgets in Couette and Boundary Layer Flows

Reynolds stress budgets for both Couette and boundary layer flows are evaluated and presented. Data are taken from direct numerical simulations of rotating and non-rotating plane turbulent Couette flow and turbulent boundary layer with and without adverse pressure gradient. Comparison of the total shear stress for the two types of flows suggests that the Couette case may be regarded as the high Reynolds number limit for the boundary layer flow close to the wall. The limit values of turbulence statistics close to the wall for the boundary layer for increasing Reynolds number approach the corresponding Couette flow values. The direction of rotation is chosen so that it has a stabilizing effect, whereas the adverse pressure gradient is destabilizing. The pressure-strain rate tensor in the Couette flow case is presented for a split into slow, rapid and Stokes terms. Most of the influence from rotation is located to the region close to the wall, and both the slow and rapid parts are affected. The anisotropy for the boundary layer decreases for higher Reynolds number, reflecting the larger separation of scales, and becomes close to that for Couette flow. The adverse pressure gradient has a strong weakening effect on the anisotropy. All of the data presented here are available on the web [36]. This revised version was published online in July 2006 with corrections to the Cover Date.     

POD investigation of the unsteady turbulent boundary layer developing over porous moving flexible fishing net structure

Particle image velocimetry (PIV) measurements are made to investigate the boundary layer developing over a modeled bottom trawl. The random motion of the fishing net structure as well as the flexibility and the porosity of this structure means that it is not enable to access the main characteristics of such a flow, using classical post-processing mathematical tools. An innovative post-treatment tool based on proper orthogonal decomposition (POD) is then developed to extract the mean velocity flow field from each available PIV instantaneous unsteady velocity field. In order to do so, the whole available velocity database is used to compute POD eigenfunctions and the first POD modes are identified as representing the mean flow field. It is then possible to deduce the mean boundary layer flow field for each position of the fishing net structure during PIV measurements. It is then observed that the mean flow field strongly depends on multiple parameters such as surface curvature, structure porosity, random motion of the structure. Streamwise evolution of classical thicknesses of boundary layer flow are also analyzed. The present work also provides benchmark PIV data of the unsteady flow developing on fishing net porous structures, which helps the progress in unsteady numerical codes for this investigation.