Visualization and laser measurements on the flow field and sand movement on sand dunes with porous fences

The installation of windbreak sand fences around sand dunes is one of the most promising methods to suppress windblown sand movement. In the study reported in this paper, we investigated the influence and validity of a small fence mounted on a model sand dune, in order to understand the fence’s suppression mechanism on the sand movement. The flow field around the dune and the process of sand-dune erosion were measured using LDV, PIV, and laser-sheet visualization techniques. A non-porous fence was found to suppress sand movements in its upstream area, but to enhance erosion downstream of the fence. This intensive erosion was caused by separated shear flow from the leading edge of the fence. In this study, four levels of porosity rate of the fence were tested. The fence-porosity dependences of the turbulent flow field and the erosion were discussed. The shapes of eroded sand dunes were found to depend on the porosity rate. The relationship between the sand-dune erosion and the flow field around the dune was illustrated with schematic diagrams. We concluded that the most desirable fence porosity should be 30% in order to avoid dune erosion if installed at a middle height on the stoss surface of a dune. This porosity provides a mean velocity reduction with avoiding a separated flow, although the flow bleeding through the porous fence is accompanied by grid turbulence and induces serious erosion in a narrow space behind the fence. Furthermore, we confirmed that the empirical correlation of the critical friction velocity can be applied to sand movements influenced by a fence.     

Flow field investigation in rotating rib-roughened channel by means of particle image velocimetry

The turbulent velocity field over the rib-roughened wall of an orthogonally rotating channel is investigated by means of two-dimensional particle image velocimetry (PIV). The flow direction is outward, with a bulk Reynolds number of 1.5 × 10⁴ and a rotation number ranging from 0.3 to 0.38. The measurements are obtained along the wall-normal/streamwise plane at mid-span. The PIV system rotates with the channel, allowing to measure directly the relative flow velocity with high spatial resolution. Coriolis forces affect the stability of the boundary layer and free shear layer. Due to the different levels of shear layer entrainment, the reattachment point is moved downstream (upstream) under stabilizing (destabilizing) rotation, with respect to the stationary case. Further increase in rotation number pushes further the reattachment point in stabilizing rotation, but does not change the recirculation length in destabilizing rotation. Turbulent activity is inhibited along the leading wall, both in the boundary layer and in the separated shear layer; the opposite is true along the trailing wall. Coriolis forces affect indirectly the production of turbulent kinetic energy via the Reynolds shear stresses and the mean shear. Two-point correlation is used to characterize the coherent motion of the separated shear layer. Destabilizing rotation is found to promote large-scale coherent motions and accordingly leads to larger integral length scales; on the other hand, the spanwise vortices created in the separating shear layer downstream of the rib are less organized and tend to be disrupted by the three-dimensional turbulence promoted by the rotation. The latter observation is consistent with the distributions of span-wise vortices detected in instantaneous flow realizations.     

On the Interaction of Turbulent Shear Layers with Harmonic Perturbations

The problem of coherent perturbations in a turbulent shear layer is considered for the purpose of developing a mathematical model based on a triple decomposition that extracts the coherent components of random fluctuations. The governing equations for the mean and the coherent parts of flow are derived, assuming the eddy-viscosity equivalence for the random part of flow, and solved by iterations to provide a coupled solution of the problem as a whole. Calculations agree well with experimental data in the upstream part of the layer where the mean–coherent flow interaction is the most important. In this region, the interaction changes the mean flow velocity distribution in such a manner that the neutral stability curve is shifted upstream relative to its position in the undisturbed layer and the perturbation intensity decreases further downstream. Experiments show that the coherent waves suppress the turbulent Reynolds stress production downstream of this region, but the model fails to predict the layer spreading correctly probably due to an inadequate turbulence closure of the mean flow. For the case of a turbulent mixing layer, we suggest a new closure relation that takes into account this coherent-random interaction.     

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.     

Experimental Investigation of Three-Dimensional Vortex Structures Downstream of Vortex Generators Over a Backward-Facing Step

An experimental investigation of vortex generators has been carried out in turbulent backward-facing step (BFS) flow. The Reynolds number, based on a freestream velocity U₀ = 10 m/s and a step height h = 30 mm, was Re_h = 2.0 × 10⁴. Low-profile wedge-type vortex generators (VGs) were implemented on the horizontal surface upstream of the step. High-resolution planar particle image velocimetry (2D-2C PIV) was used to measure the separated shear layer, recirculation region and reattachment area downstream of the BFS in a single field of view. Besides, time-resolved tomographic particle image velocimetry (TR-Tomo-PIV) was also employed to measure the flow flied of the turbulent shear layer downstream of the BFS within a three-dimensional volume of 50 × 50 × 10 mm³ at a sampling frequency of 1 kHz. The flow control result shows that time-averaged reattachment length downstream of the BFS is reduced by 29.1 % due to the application of the VGs. Meanwhile, the Reynolds shear stress downstream of the VGs is considerably increased. Proper Orthogonal Decomposition (POD) and Dynamic Mode Decomposition (DMD) have been applied to the 3D velocity vector fields to analyze the complex vortex structures in the spatial and temporal approaches, respectively. A coherent bandwidth of Strouhal number 0.3 < St_h < 0.6 is found in the VG-induced vortices, and moreover, Λ-shaped three-dimensional vortex structures at St_h = 0.37 are revealed in the energy and dynamic approaches complementarily.     

“Barchan” dunes in the labDunes de type « Barchanes » au laboratoire

We demonstrate the feasibility of studying dunes in a laboratory experiment. It is shown that an initial sand pile, under a wind flow carrying sand, flattens and gets a shape recalling barchan dunes. An evolution law is proposed for the profile and the summit of the dune. The dune dynamics is shown to be shape invariant. The invariant shape, the “dune function” is isolated. To cite this article: O. Dauchot et al., C. R. Mecanique 330 (2002) 185–191.     

Non-oscillating/Oscillating Shear Layer over a Large Deep Cavity at Low-Subsonic Speeds

The flow over a deep cavity at low subsonic velocity is considered in the present paper. The cavity length-to-depth aspect ratio is L/H = 0.2. Single hot-wire measurements characterized the incident turbulent boundary layer and show the influence of the cavity on the streamwise statistic components just downstream from the cavity. The streamwise mean and fluctuating velocity profiles are affected by the cavity. PIV measurements reveal the presence for ejection-like events responsible of local perturbations of the skewness and the flatness coefficients. Time-resolved PIV technic is also used to characterize phase properties of shear layer oscillating cycle. It is shown that for deep cavity with first Rossiter mode, only one vortical structure is formed at the cavity leading edge. Then, it grows while convecting downstream along the shear layer. A well-defined ejection process begins after the vortex impact near the cavity downstream corner. A cylinder device placed spanwisely near the cavity leading edge eliminates the resonance and highly modifies the behavior of the shear layer flow. In fact, the shear layer could be divided into upper and lower parts with different structure aspects.     

Simultaneous measurements of velocity and deformation in flows through compliant diaphragms

Flow through a circular orifice in a deformable diaphragm mounted in a pipe was studied experimentally as a simple yet suitable case for validating numerical fluid/structure interaction (FSI) codes including structures with significant deformation and strain. The flow was characterized using pressure taps, particle image velocimetry (PIV), and hot-film anemometry while deformation of the compliant diaphragm was determined directly from PIV images. The diaphragm material properties were measured independently by a uniaxial tensile testing machine. The diaphragm material modulus, orifice diameter, and pipe Reynolds number were varied over ranges appropriate for simulations of flows through heart valves. Pipe Reynolds numbers ranged from 600 (laminar upstream condition) to 8800 (turbulent upstream condition). The pressure drop across the diaphragm resulted in a concave deformation for all cases studied. For the range of Reynolds number tested, the Euler number decreased with increasing Reynolds number as a result of orifice expansion. The flow immediately downstream of compliant diaphragms was jet-like with strong inward radial velocity components and vena contracta. Laminar low Reynolds number flow (Re=600) through both rigid and compliant diaphragms yielded early and regular roll up of coherent vortex rings at a fixed frequency in contrast to turbulent higher Reynolds number flow (Re=3900), which yielded a broad range of vortex passage frequencies. Expansion of the compliant orifice for Re=3900 resulted in an initially broader slower jet with delayed shear layer development compared with the equivalent rigid case.     

Investigation of the unsteadiness of a shock-reflection interaction with time-resolved particle image velocimetry

The spatio-temporal dynamics of an impinging shock/boundary layer interaction at Mach 2 and under incipient separation conditions, has been investigated experimentally by means of high-speed particle image velocimetry (PIV). The available PIV acquisition rate of up to 20 kHz permits a time-resolved characterization of the interaction. The dynamics of different flow regions—notably the separation region and the reflected shock—were quantified by means of temporal auto-correlation fields and pseudo-spectral analysis. The PIV data further enable to investigate the relationship between spatially extended flow features, such as shock position and bubble size, as well as the influence of the upstream boundary layer. The results confirm earlier studies that there is an important upstream effect on the present incipient interaction.     

Application of an acoustic analogy to PIV data from rectangular cavity flows

The present paper describes a method to derive information about the acoustic emission of a flow using particle image velocimetry (PIV) data. The advantage of the method is that it allows studying sound sources, the related flow phenomena and their acoustic radiation into the far field, simultaneously. In a first step the time history of two-dimensional instantaneous pressure fields is derived from planar PIV data. In a successive step the Curle’s acoustic analogy is applied to the pressure data to obtain the acoustic radiation of the flow. The test cases studied here are two rectangular cavity flows at very low Mach number with different aspect ratios L/H. The main sound source is located at the cavity trailing edge and it is due to the impingement of vortices shed in the shear layer. It is shown that the flow emits sound with a main directivity in the upstream direction for the smaller aspect ratio and the directivity is more uniform for the larger aspect ratio. In the latter case the acoustic pressure spectra has a broader character due to the impact of the downstream recirculation zone onto the shear layer instabilities, destroying their regular pattern and alternating the main sound source.     

Time-resolved and volumetric PIV measurements of a transitional separation bubble on an SD7003 airfoil

To comprehensively understand the effects of Kelvin–Helmholtz instabilities on a transitional separation bubble on the suction side of an airfoil regarding as to flapping of the bubble and its impact on the airfoil performance, the temporal and spatial structure of the vortices occurring at the downstream end of the separation bubble is investigated. Since the bubble variation leads to a change of the pressure distribution, the investigation of the instantaneous velocity field is essential to understand the details of the overall airfoil performance. This vortex formation in the reattachment region on the upper surface of an SD7003 airfoil is analyzed in detail at different angles of attack. At a Reynolds number Re _c < 100,000 the laminar boundary layer separates at angles of attack >4°. Due to transition processes, turbulent reattachment of the separated shear layer occurs enclosing a locally confined recirculation region. To identify the location of the separation bubble and to describe the dynamics of the reattachment, a time-resolved PIV measurement in a single light-sheet is performed. To elucidate the spatial structure of the flow patterns in the reattachment region in time and space, a stereo scanning PIV set-up is applied. The flow field is recorded in at least ten successive light-sheet planes with two high-speed cameras enclosing a viewing angle of 65° to detect all three velocity components within a light-sheet leading to a time-resolved volumetric measurement due to a high scanning speed. The measurements evidence the development of quasi-periodic vortex structures. The temporal dynamics of the vortex roll-up, initialized by the Kelvin–Helmholtz (KH) instability, is shown as well as the spatial development of the vortex roll-up process. Based on these measurements a model for the evolving vortex structure consisting of the formation of c-shape vortices and their transformation into screwdriver vortices is introduced.     

The Flow Structure Produced by Pulsed-jet Vortex Generators in a Turbulent Boundary Layer in an Adverse Pressure Gradient

The effect of pulsed jet vortex generators on the structure of an adverse pressure gradient turbulent boundary layer flow was investigated. Two geometrically optimised vortex generator configurations were used, co-rotating and counter-rotating. The duty cycle and pulse frequency were both varied and measurements of the skin friction (using hot films) and flow structure (using stereo PIV) were performed downstream of the actuators. The augmentation of the mean wall shear stress was found to be dependent on the net mass flow injected by the actuators. A quasi steady flow structure was found to develop far downstream of the injection location for the highest pulse frequency tested. The actuator near field flow structure was observed to respond very quickly to variations in the jet exit velocity.     

Investigation of the flow in a circular cavity using stereo and tomographic particle image velocimetry

The turbulent flow over a circular cavity with an aspect ratio of D/H = 2 is investigated by multi-planar stereoscopic particle image velocimetry and with tomographic particle image velocimetry (PIV). The main aim of the study is the flow topology and the turbulent structure of the asymmetrical flow pattern that forms inside the cavity at these specific conditions. The flow field is measured in the vertical symmetry plane to describe the overall recirculation pattern in the cavity and the turbulent shear layer developing from the separation point. In this specific regime the shear layer fluctuations are recognized as those caused by instabilities together with the effect of the incoming boundary layer turbulence. Additional observations performed at several wall-parallel planes at different height inside the cavity allow to further evaluate the secondary flow circulation generated by this asymmetric regime. The observed flow pattern consists of a steady vortex, occupying the entire cavity volume and placed diagonally inside the cavity such to entrain the external flow from one side, capture it into a circulatory motion and eject it from the opposite side of the cavity. The spatial distribution of the turbulent fluctuations also reveals the same structure. The tomographic PIV measurement returns a visual inspection to the instantaneous three-dimensional structure of the turbulent fluctuations, which at the investigated height exhibit a low level of coherence with slightly elongated vortices in the recirculating flow inside the cavity.     

In vitro post-stenotic flow quantification and validation using echo particle image velocimetry (Echo PIV)

Echo particle image velocimetry (Echo PIV) presents itself as an attractive in vivo flow quantification technique to traditional approaches. Promising results have been acquired; however, limited quantification and validation is available for post-stenotic flows. We focus here on the comprehensive evaluation of in vitro downstream stenotic flow quantified by Echo PIV and validated in relation to digital particle image velocimetry (DPIV). A Newtonian blood analog was circulated through a closed flow loop and quantified immediately downstream of a 50 % axisymmetric blockage at two Reynolds numbers (Re) using time-averaged Echo PIV and DPIV. Centerline velocities were in good agreement at all Re; however, Echo PIV measurements presented with elevated standard deviation (SD) at all measurements points. SD was improved using increased line density (LD); however, frame rate or field of view (FOV) is compromised. Radial velocity profiles showed close agreement with DPIV with the largest disparity in the shear layer and near-wall recirculation. Downstream recirculation zones were resolved by Echo PIV at both Re; however, magnitude and spatial coverage was reduced compared to DPIV that coincided with reduced contrast agent penetration beyond the shear layer. Our findings support the use of increased LD at a cost to FOV and highlight reduced microbubble penetration beyond the shear layer. High local SD at near-wall measurements suggests that further refinement is required before proceeding to in vivo quantification studies of wall shear stress in complex flow environments.     

Three-dimensional particle-tracking velocimetry measurement of turbulence statistics and energy budget in a backward-facing step flow

Using a three-dimensional (3-D) particle-tracking velocimeter, detailed turbulent flow measurements were made in a plane channel with a one-sided 50% abrupt expansion, which acted as a backward-facing step. The turbulent channel flow reached a fully developed state well upstream of the step. The Reynolds number based on the upstream centerline velocity and the step height H was 5540. With the mean reattachment point located at 6.51H downstream of the step, the measurement region ranged from −2H upstream to 12H downstream of the step. Various turbulent statistics and the energy budget were calculated from numerous instantaneous vector distributions. As in previous experimental investigations, the Reynolds normal and shear stresses had maximum values upstream of the reattachment. The stress anisotropy tensor revealed a peculiar phenomenon near the reattachment wall, wherein the spanwise normal stress was the largest among the three normal stresses. The triple velocity correlations indicated large values in the separating shear layer, and hence the turbulent diffusion was a major term in the energy budget. Comparison was made between the present results and those of the direct numerical simulation (DNS) of Le et al. (1993), and it was found that the mean and fluctuating velocities, the Reynolds shear stress, and the turbulent energy budget were in excellent agreement, although there was a considerable difference in the inflow conditions.     

PIV measurements of turbulent flow in planar mixing layer

A turbulent mixing layer consists of two different flow types, i.e. shear layer (shear-flow turbulence) and free stream regions (nearly homogeneous turbulence). The inherent non-uniform seeding tracer distributions observed around the interfaces between the shear layer and two free stream regions usually lead to a difficulty in particle image velocimetry (PIV) measurements. A parametric study on the application of PIV to the measurement of velocity field in a planar mixing layer is made by means of six factors, including interrogation window size, aspect ratio of interrogation window, interrogation window offset, threshold of data validation, sharpening spatial filters (Prewitt and Sobel masks), and smoothing spatial filter (median mask). The objective of this study is to obtain accurate turbulent measurements in both mean and fluctuating velocities using PIV under an appropriate parametric setting. The optimal levels, which are trade-off in between the accuracy and fine spatial resolution of velocity field measurements, are determined with the aid of the Taguchi method. It is shown that the PIV measurements made with this optimal set of parameters are in good agreement with the measurements made by a two-component hot-wire anemometer. Case independency of the proposed optimal set of parameters on the flow condition of the mixing layer is validated through the applications to two additional tests under the different experimental conditions in changing solely either velocity ratio of high-speed to low-speed free stream velocities or Reynolds number.     

Scanning PIV investigation of the laminar separation bubble on a SD7003 airfoil

A laminar separation bubble occurs on the suction side of the SD7003 airfoil at an angle of attack α =  4–8° and a low Reynolds number less than 100,000, which brings about a significant adverse aerodynamic effect. The spatial and temporal structure of the laminar separation bubble was studied using the scanning PIV method at α =  4° and Re = 60,000 and 20,000. Of particular interest are the dynamic vortex behavior in transition process and the subsequent vortex evolution in the turbulent boundary layer. The flow was continuously sampled in a stack of parallel illuminated planes from two orthogonal views with a frequency of hundreds Hz, and PIV cross-correlation was performed to obtain the 2D velocity field in each plane. Results of both the single-sliced and the volumetric presentations of the laminar separation bubble reveal vortex shedding in transition near the reattachment region at Re = 60,000. In a relatively long distance vortices characterized by paired wall-normal vorticity packets retain their identities in the reattached turbulent boundary layer, though vortices interact through tearing, stretching and tilting. Compared with the restricted LSB at Re = 60,000, the flow at Re = 20,000 presents an earlier separation and a significantly increased reversed flow region followed by “huge” vortical structures.     

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.     

Vortex shedding induced by a solitary wave propagating over a submerged vertical plate

Experimental study was conducted on the vortex shedding process induced by the interaction between a solitary wave and a submerged vertical plate. Particle image velocimetry (PIV) was used for quantitative velocity measurement while a particle tracing technique was used for qualitative flow visualization. Vortices are generated at the tip of each side of the plate. The largest vortices at each side of the plate eventually grow to the size of the water depth. Although the fluid motion under the solitary wave is only translatory, vortices are shed in both the upstream and downstream directions due to the interaction of the generated vortices as well as the vortices with the plate and the bottom. The process can be divided into four phases: the formation of a separated shear layer, the generation and shedding of vortices, the formation of a vertical jet, and the impingement of the jet onto the free surface. Similarity velocity profiles were found both in the separated shear layer and in the vertical jet.     

Correlation study in shock wave�Cturbulent boundary layer interaction

Shock wave–turbulent boundary layer interaction is a critical problem in aircraft design. Therefore, a thorough understanding of the processes occurring in such flows is necessary. The most important task is to study the unsteady phenomena, in particular, the low-frequency ones, for this interaction. An experimental study of separated flow has been performed in the zone of interaction of the incident oblique shock wave with a turbulent boundary layer at Mach 2. Two-point correlation data in the separation zone and the upstream flow were obtained and showed that low-frequency oscillations of the reflected shock waves are related to pulsations in the inflow turbulent boundary layer.