Influence of Inlet Shear on Structure of Wake behind a Square Cylinder

The force coefficients and the frequency of vortex shedding in the wake of a square cylinder exposed to a uniform shear flow and the flow structure around it were numerically investigated. The Reynolds number defined on the basis of cylinder width was in the range of 250–1,500. The shear parameter, namely the transverse velocity gradient, which is nondimensionalized using the obstacle width and the average incoming velocity, was varied between 0 and 0.2. Analyses were performed for a number of flow parameters using various combinations of Reynolds number and shear parameters. Results show that mean and root-mean-square values of drag coefficient initially decrease up to certain values of the shear rate and then increase with increase in shear parameter. The root-mean-square values of lift coefficient show a similar behavior. The Strouhal number decreases uniformly with increase in shear parameter. At higher shear rates, the von Kármán vortex street comprising alternating vortices breaks, and the far wake shows mainly clockwise vortices.     

Flow Characteristics around a Circular Cylinder Placed Horizontally above a Plane Boundary

Flow characteristics around a circular cylinder positioned near a plane boundary (on which laminar boundary layer flow develops in the absence of circular cylinder), are investigated for Reynolds numbers R ranging from 7.8×102 to 1.15×104. Particle image velocimetry and fiber laser Doppler velocimetry were used to measure the velocity fields and velocity time histories, respectively. Flow structures are particularly revealed using flow visualization technique at R = 7.8×102 for gap ratios G/D (where G is the net gap between the surface of circular cylinder and the plane boundary), varying from 0 to 4. Based on the experimental results, the variation of Strouhal number of shedding vortex (or eddy) with G/D, the mechanism of vortex shedding suppression, and the streamwise velocity profiles of the upper shear layers and gap flows for small G/D are all discussed. Although the regular, alternate vortex shedding is suppressed for G/D<0.5, the periodicity could be detected due to the vortex (or eddy) shedding from the upper shear layer of the circular cylinder. Gap flow switching randomly is found and first put forward to be the main reason of multipeak or broadband spectral characteristics of the shedding event at a certain small gap ratio. It is also found that the streamwise velocity profiles of the upper shear layer, where periodic shedding eddies originate, exhibit well-behaved similarity. In addition, a unique similarity of mean streamwise velocity profiles of the gap flows is demonstrated for G/D ? 0.3. For R<4×103, the S increases as G/D decreases to its maximum around G/D ? 0.5 and then decreases as G/D decreases. For R ≥ 4×103, although most of the previous studies indicate that the S is insensitive to G/D, the present study shows that S still increases as G/D decreases but the variations of S are in a small range (i.e., 0.18 ? S ? 0.22).     

Flow around Cylinders in Open Channels

This paper presents the results of an experimental study of flow around cylindrical objects in an open channel. Cylindrical objects of equal diameter and four heights were tested under similar flow conditions producing four different levels of submergence, including a surface piercing bridge-pier-like cylinder. Different flow elements and their locations were identified using a set of flow visualization tests. Observations made from the flow visualization tests were then verified by measurements of bed-shear stress and deflected flow velocity around the cylinders. Horse-shoe vortex systems were found to appear closer to the submerged cylinders compared to a surface piercing cylinder. The increase in dimensionless bed-shear stress is found to be inversely related to the level of submergence of the cylinders. Bed-shear stress results presented in this paper will be valuable for a qualitative understanding of the scour potential of flow around submerged cylinders. Mean velocity profiles in the deflected flow region were analyzed in terms of the theories of three-dimensional turbulent boundary layer. Submergence of a cylinder has been found to suppress alternate vortex shedding and produce stronger three-dimensional flows in the downstream wake. Perry and Joubert’s model was found to be sufficiently accurate to predict the deflected velocity magnitudes around submerged cylinders. Overall, the present study will provide valuable knowledge of hydraulics of flow around submerged structures (e.g., simple fish habitat structures).     

Turbulence Structures in Flow over Two-Dimensional Dunes

This paper presents a large eddy simulation (LES) of turbulent open channel flow over two-dimensional periodic dunes. The Reynolds number R based on the bulk velocity U(bulk) and the maximum flow depth h, is approximately 25,000. The instantaneous flow field is investigated with special emphasis on the occurrence of coherent structures. Instantaneous vortices were visualized and it is shown that separated vortices are formed downstream of the dune crest due to Kelvin–Helmholtz instabilities. Near the point of reattachment the so-called kolk-boil vortex evolves in form of a hairpin vortex. Also present are previously separated vortices, which are convected along the stoss side of the downstream dune and elevated toward the water surface. The existence of near wall streaks which reform shortly after reattachment is also shown. The spacing between two low-speed streaks is very similar to that observed previously over smooth and rough walls. For validation, profiles of the time-averaged velocities, streamwise, and wall-normal turbulent intensities and the Reynolds shear stress calculated by the LES are presented and compared with available laser Doppler velocimetry measurements and overall good agreement is found.     

Visualization of vortex shedding at the proximal side-to-end artery-graft anastomosis

The motion of guanine particles was recorded by video to visualize transitional flow phenomena in models of a proximal side-to-end anastomosis. Close examination of successive video fields revealed that above a critical Reynolds number, particles were periodically shed into the graft from a vortex situated near the anastomosis heel, and this disturbed the flow patterns in the graft causing vortex shedding to occur near to the toe of the anastomosis. The images clearly demonstrated that periodic flow structures propagated distally along the graft for over 15 tube diameters from the proximal anastomosis. The frequency of the vortex shedding was found to increase with Reynolds number. Under pulsatile conditions, the primary vortex at the heel of the anastomosis became unstable during the deceleration phase of the flow cycle and particles were shed downstream into the graft. Although it was possible briefly to observe the characteristic banded structure in the bypass graft, the flow patterns were highly three-dimensional and were quickly broken up by the accelerating flow. Dynamic flow visualization using guanine particles was found to be a complementary technique to particle tracer flow visualization and was highly effective in identifying transitional flow phenomena and the mass transport mechanisms associated with them.     

Analysis of the axial flow field in stenosed carotid artery bifurcation models--LDA experiments

Laser Doppler anemometer (LDA) experiments were performed to gain quantitative information on the differences between the large-scale flow phenomena in a non-stenosed and a stenosed model of the carotid artery bifurcation. The influence of the presence of the stenosis was compared to the effect of flow pulse variation to evaluate the feasibility of early detection of stenosis in clinical practice. Three-dimensional Plexiglass models of a non-stenosed and a 25% stenosed carotid artery bifurcation were perfused with a Newtonian fluid. The flow conditions approximated physiological flow. The results of the velocity measurements in the non-stenosed model agreed with the results from previous hydrogen-bubble visualization. A shear layer separated the low-velocity area near the non-divider wall from the high-velocity area near the divider wall. In this shear layer, vortex formation occurred during the deceleration phase of the flow pulse. The instability of this shear layer dictated the flow disturbances. The influences of the mild stenosis, located at the non-divider wall, was mainly limited to the stability of the shear layer. No disturbances were found downstream of the stenosis near the non-divider wall. Using a pulse wave with an increased systolic deceleration time, the velocity distribution showed an extended region with reversed flow, a more pronounced shear layer and increased vortex strength. From these measurements it is obvious that the influence of the presence of a mild stenosis, mainly limited to the stability of the shear layer, can hardly be distinguished from the effects of a variation of the flow pulse. From this it can be concluded that methods for detection of mild stenosis, using solely the large-scale flow phenomena, as can be measured by ultrasound or MRI techniques, will hardly have any clinical relevance.     

Coupling of U-Tube and Shear Layer Oscillations in an Interconnected Flow Compartment

Results of experiments conducted in a 2?m high flume at large Reynolds numbers are reported in this paper. The flume was partitioned into two compartments. Flow entered the bottom of the upstream test compartment as a wall jet, at jet Reynolds number ranging from 11,000 to 170,000. Periodic oscillations of the free surface in the two compartments resembling the oscillatory flow in a liquid-filled U-tube, and large coherent structures formed above the potential core of the wall jet were observed. Coupling of the U-tube oscillations and vortex shedding is attributed to fluid-dynamic and fluid-resonant feedback processes. For test compartment length, Lc = 0.8?m, fluid-resonant feedback was found to be dominant, and the shear layer was observed to oscillate at the natural frequency of the two-compartment, U-tube system. The observed U-tube oscillations are initiated by the oscillations of the shear layer at a frequency equal to the subharmonic component for the U-tube. The flow oscillations were generally weaker for Lc = 1.2 and 2.0?m with oscillation frequencies governed by fluid-dynamic feedback, verified from a comparison with the results from a previously reported study.     

Detached Eddy Simulation Investigation of Turbulence at a Circular Pier with Scour Hole

This paper uses results from detached eddy simulation to reveal the dynamics of large-scale coherent eddies in the flow around a circular pier with an equilibrium scour hole. This is important for the sediment transport because the local scour process is controlled to a large extent by the large-scale coherent structures present in the near-bed region. The present paper investigates the dynamics of these coherent structures, their interactions and their role in entraining sediment in the later stages of the scour process when the horseshoe vortex system is stabilized by the presence of a large scour hole. The pier Reynolds number was 2.06×105, outside the range of well-resolved large-eddy simulation (LES). Additionally, scale effects are investigated based on comparison with LES results obtained at a much lower Reynolds number of 16,000 in a previous investigation. The paper provides a detailed study of the dynamics of the main necklace vortices of the horseshoe vortex system, including an investigation of the bimodal oscillations, their effect on the amplification of the turbulence within the scour hole and the interactions of the necklace vortices with the downflow. Several mechanisms for the growth of the downstream part of the scour hole in the later stages of the scour process are discussed. Similar to the low-Reynolds-number simulation, and consistent with experimental observations, the presence of strong upwelling motions near the symmetry plane resulted in the suppression of the large-scale vortex shedding in the wake. The fact that the nondimensional values of the turbulent kinetic energy and pressure RMS fluctuations in the higher Reynolds number simulation were consistently lower inside the regions of high turbulence amplification associated with the main necklace vortex, the separated shear layers and the near-wake shows that changes in the flow and turbulence due to the Reynolds number and scour hole geometry can be quantitatively significant over Reynolds numbers between 104 and 105.     

Characteristics of Shear Layer Structure in Skimming Flow over a Vertical Drop Pool

The characteristics of shear layer structure between the sliding jet and the pool for skimming flows over a vertical drop pool were investigated experimentally, using flow visualization technique and high speed particle image velocimetry. Four series of experiments having different end sill ratios (h/H = 0.12, 0.43, 0.71 and 1.0, where h=end sill height and H=drop height) with various approaching flow discharges were performed to measure the detailed quantitative velocity fields of the shear layer. The mean velocities and turbulence properties were obtained by ensemble averaging the repeated measurements. From the velocity profiles, it is found that the growth of the shear layer in the downward direction as the jet slides down the pool represents the momentum exchange. Analyzing the distribution of measured velocity, the similarity profile of the mean velocity at different cross sections along the shear layer was obtained. The proposed characteristic scales provided unique similarity profiles having promising regression coefficient. The selection of these characteristic scales is also discussed. Further, the spatial variations of mean velocity profiles, turbulence intensities, in-plane turbulent kinetic energy, and Reynolds shear stress were also elucidated in detail. The imperative observation is that the Reynolds shear stress dominates the major part along the shear layer as compared to the viscous shear stress. The study also provides an insight into the flow phenomena through the velocity and turbulent characteristics.     

Experimental Investigation of Flow Past a Square Cylinder at an Angle of Incidence

Flow past a square cylinder placed at an angle to the incoming flow is experimentally investigated using particle image velocimetry, hot wire anemometry, and flow visualization. The Reynolds number based on cylinder size and the average incoming velocity is set equal to 410. Data for four cylinder orientations (θ = 0, 22.5, 30, and 45°) and two aspect ratios [AR = 16 and 28] are reported. Results are presented in terms of drag coefficient, Strouhal number, time averaged velocity, stream traces, turbulence intensity, power spectra, and vorticity field. In addition, flow visualization images in the near wake of the cylinder are discussed. The shape and size of the recirculation bubble downstream of the cylinder are strong functions of orientation. A minimum in drag coefficient and maximum in Strouhal number is observed at a cylinder orientation of 22.5°. The v-velocity profile and time-average stream traces show that the wake and the separation process are asymmetric at orientations of 22.5 and 30°. The corresponding power spectra show additional peaks related to secondary vortical structures that arise from nonlinear interaction between the Karman vortices. The flow visualization images show the streamwise separation distance between the alternating vortices to be a function of cylinder orientation. Further, the flow approaches three dimensionality early, i.e., closer to the cylinder surface for the 22.5° orientation. The drag coefficient decreases with an increase in aspect ratio, while the Strouhal number is seen to increase with aspect ratio. The turbulence intensity is higher for the large aspect ratio cylinder and the maximum turbulence intensity appears at an earlier streamwise location. The overall dependence of drag coefficient and Strouhal number on orientation is preserved for the two aspect ratios studied.     

Vortex Flow Pattern in a Slab Continuous Casting Mold with Argon Gas Injection

An Eulerian-Eulerian two-fluid model was developed to study the vortex flow inside a slab continuous casting mold with argon gas injection. Interracial momentum transfer that accommodated various interracial forces including drag force, lift force, virtual mass force, and turbulent dispersion force was considered. Predicted results agree well vaith experimental measurements of the water model in two-phase flow pattern and vortex flow structures. Three typical flow patterns with different argon steel ratios （ASRs） have been obtained： ＂double roll＂, ＂three roll＂, and ＂single roll＂. The flow pattern inside the mold alternates among the three types or it may attain some intermedi ate condition. With increasing ASR, the positions of vortices move from the submerged entry nozzle to the narrow face of the mold, and the sizes of vortices are reduced gradually. The rotating directions of vortices are all from high velocity area to low velocity area. Two mechanisms of vortex formation on the top surface have been suggested, i. e. , congruous shear flow and incongruous shear flow.     

Reynolds Number Effects on the Near-Exit Region of Turbulent Jets

In this paper, attention has been focused on the near-exit region of a round turbulent free jet to study the large-scale coherent structures and to document the signatures of the vortices over a range of Reynolds numbers. Particle image velocimeter measurements were conducted at three jet exit Reynolds numbers of 10,000, 30,000, and 55,000. The large-scale structures in the near field (X/D<12) were investigated by performing a proper orthogonal decomposition analysis of the velocity fields. A vortex identification algorithm was complemented by swirling strength maps to identify the vortex centers, rotational sense, size, and circulation of the vortices. The influence of the Reynolds number on the distribution of the number, size, and circulation of the identified vortices was studied. Proper orthogonal decomposition of the velocity fields revealed that Reynolds number has a strong influence on the mean circulation of vortices. The present results show that the axial location where vortices first appear and the number of vortices close to the nozzle exit (X/D<5) are dependent on the Reynolds number.     

Stability Analysis of Velocity Profiles in Water-Hammer Flows

This paper performs linear stability analysis of base flow velocity profiles for laminar and turbulent water-hammer flows. These base flow velocity profiles are determined analytically, where the transient is generated by an instantaneous reduction in flow rate at the downstream end of a simple pipe system. The presence of inflection points in the base flow velocity profile and the large velocity gradient near the pipe wall are the sources of flow instability. The main parameters that govern the stability behavior of transient flows are the Reynolds number and dimensionless timescale. The stability of the base flow velocity profiles with respect to axisymmetric and asymmetric modes is studied and its results are plotted in the Reynolds number∕timescale parameter space. It is found that the asymmetric mode with azimuthal wave number 1 is the least stable. In addition, the results indicate that the decrease of the velocity gradient at the inflection point with time is a stabilizing mechanism whereas the migration of the inflection point from the pipe wall with time is a destabilizing mechanism. Moreover, it is shown that a higher reduction in flow rate, which results in a larger velocity gradient at the inflection point, promotes flow instability. Furthermore, it is found that the stability results of the laminar and the turbulent velocity profiles are consistent with published experimental data and successfully explain controversial conclusions in the literature. The consistency between stability analysis and experiments provide further confirmation that (1) water-hammer flows can become unstable; (2) the instability is asymmetric; (3) instabilities develop in a short (water-hammer) timescale; and (4) the Reynolds number and the wave timescale are important in the characterization of the stability of water-hammer flows. Physically, flow instabilities change the structure and strength of the turbulence in a pipe, result in strong flow asymmetry, and induce significant fluctuations in wall shear stress. These effects of flow instability are not represented in existing water-hammer models.     

Measurement of Vortices in Model Pump-Intake Bay by PIV

This paper describes experiments aimed at improving the understanding of vortices that form in water-pump intake bays. Experiments were conducted in a laboratory model of a simple intake bay comprising a vertical intake pipe with a bell mouth in a rectangular channel under conditions for which surface tension and free-surface elevation effects were small. The experiments comprised flow visualization followed by detailed measurements of the velocity field with particle-image velocimetry (PIV) in selected sections across the free-surface and wall-attached (subsurface) vortices. All vortices were found to meander in space. Depending on the submergence of the intake pipe and its position relative to the channel walls, some vortices had well-defined mean positions, while others were intermittent or exhibited bimodal behavior. The subsurface vortex strength depended on the distance between the wall or floor to which the vortex was attached and the intake pipe, with smaller distances promoting stronger vortices. The experiments clearly demonstrate the usefulness of PIV in obtaining the quantitative information that is required about the number, location, size, and strength of the vortices to develop and validate numerical models of pump intake flows.     

Simultaneous Particle Image Velocimetry and Laser Doppler Velocimetry Measurements of Periodical Oscillatory Horseshoe Vortex System near Square Cylinder-Base Plate Juncture

This paper presents simultaneous measurements using particle image velocimetry (PIV) and laser Doppler velocimetry (LDV) techniques on the study of a horseshoe vortex system. The horseshoe vortex system is generated near the juncture of a vertical square cylinder and a horizontal base plate. The combination of PIV and LDV not only gives the spatial distribution and time history of velocity near the juncture for spatial and time domain analyses, it also allows phase averaging the PIV velocity data to reduce noise and, in a turbulent flow, result in turbulence statistics. A flow visualization technique displaying particle streaklines has also been used to help the classification of the vortex system and visualize the flow motion and vortex evolution. The classification of the horseshoe vortex was briefly categorized as steady, periodical oscillatory, and turbulence-like chaotic vortex systems through the use of the flow visualization technique and time-domain spectral analysis. Phase-averaged flow characteristics of the periodical oscillatory vortex system with a Reynolds number of 2,250 are presented in detail through the use of PIV and LDV as well as the flow visualization technique.     

Characteristics of Steady Horseshoe Vortex System near Junction of Square Cylinder and Base Plate

This paper presents an experimental investigation on the characteristics of a horseshoe vortex system near the juncture of a square cylinder and a horizontal base plate, using particle image velocimetry and flow visualization technique. Experiments were conducted for Reynolds numbers (based on the free stream velocity and the width of square cylinder) ranging from 2.0×102 to 6.0×103. The flow patterns are first classified into four major regimes: Steady horseshoe vortex system, periodic oscillation vortex system with small displacement, periodic breakaway vortex system, and irregular vortex system. The classifications can be demonstrated as a figure of Reynolds number versus the ratio of the height of square cylinder to undisturbed boundary layer thickness. The study then mainly focused on the characteristics of steady horseshoe vortex system (corresponding to Reynolds numbers ranging from 2.0×102 to 2.5×103). The nondimensional characteristics, including the horizontal and vertical distances from the primary vortex core to frontal face of the vertical square cylinder and bottom boundary of the base plate, respectively, the height of stagnation point at frontal face of the square cylinder, and the down-flow discharge as well as circulation of the primary vortex, all increase with increase of the ratio of the height of square cylinder to undisturbed boundary layer thickness. However, they all decrease with the increase of the aspect ratio (i.e., the height-to-width ratio) of the square cylinder. The study provides essential properties of a steady horseshoe vortex system and gives an insight for related engineering applications. It can be served as a basis for more complicated horseshoe vortex systems occurring at high Reynolds numbers.     

Effects of Bed Roughness on Flow around Bed-Mounted Cylinders in Open Channels

This paper presents the results of an experimental study of flow around cylindrical objects on a rough bed in an open channel. This is an extension of a previous study of flow around cylinders on a smooth bed. The purpose of this study is to explore the effects of bed roughness on the characteristics of the deflected flow around cylindrical objects and the resulting bed-shear stress distributions. Similar to the previous study cylindrical objects of equal diameter and four heights were tested under similar flow conditions producing four different levels of submergence. Bed shear stress and deflected flow velocities were measured by a thin yaw-type Preston probe after a set of flow visualization tests. Flow visualization tests showed that the horse-shoe vortex systems on the rough bed occupy a relatively greater width compared to the smooth bed. Unlike smooth bed observations, the flow separation point upstream of the cylinder was not dependent on the level of submergence as the separation points were found to appear within a short range of x = ?1D to ?1.2D. Bed shear stress has been found to increase significantly near the shoulder of the cylinders, and its ratio with respect to the approach bed-shear stress was twice as large compared to the smooth bed case. Mean velocity profiles were analyzed in terms of three-dimensional turbulent boundary layer theories. Bed roughness was found to oppose the effect of the lateral pressure gradient that causes skewing in the boundary layer. Perry and Joubert’s model has been found to be equally accurate on smooth and rough beds for predicting the deflected velocity magnitudes around cylinders. The present study will enhance the knowledge of hydraulics of flow around bed-mounted objects (e.g. fish-rocks) in natural streams.     

Investigation of Fluid Structures in a Smooth Open-Channel Flow Using Proper Orthogonal Decomposition

This paper reports particle image velocimetry (PIV) measurements of the instantaneous velocity fields in a smooth open-channel flow. The Reynolds number of the flow based on the water depth was 21,000. The instantaneous velocity fields were analyzed using proper orthogonal decomposition (POD) to expose the vortical structures. The velocity fields were reconstructed using different combination of modes; the first 12 modes to expose the energetic structures, and from Modes 13 to 100 to expose the less energetic structures. The first set recovered about 50% of the turbulent kinetic energy while the second group of modes recovered about 33% of the energy. The POD results were further combined with the results from the momentum analysis as well as with the conditional quadrant analysis performed at three different threshold levels. The POD results revealed the existence of hairpin vortices of different sizes and energy levels. Most of the large eddies are elongated and inclined toward the boundaries in the streamwise direction. The results also revealed patterns of strong ejection and sweep events which are common features in wall-bounded flows. Closer to the free surface (y/d>0.6), it was observed that the existence of hairpin vortices with legs possibly extended upward toward the free surface. As well, the distribution of the uniform momentum zones was consistent with the location of the vortices and their induced flow. While POD exposes the large- and small-scale structures based on the amount of turbulent kinetic energy, the quadrant analysis performed on the PIV maps shows the spatial distribution of the events related to the momentum transport.     

Wall-Wake Flows Downstream of a Sphere Placed on a Plane Rough Wall

The time-averaged characteristics of turbulent wall-wake flows downstream of a sphere placed on a rough wall are studied. The profiles of the defect of streamwise velocity, Reynolds shear stress, and turbulence intensities exhibit some degree of similarities when they are scaled by their respective peak defect values. For the velocity defect profiles, the vertical distances are scaled by the height of the location of the half-peak velocity defect. However, for the defect profiles of the Reynolds shear stress and the turbulence intensities, the vertical distances are scaled by the height of the location of the half-peak Reynolds shear stress defect. The magnitudes of the peak defect of all the quantities diminish with the distance downstream of the sphere characterizing the recovery of their undisturbed profiles. Additionally, the theoretical similarity solution for the velocity defect profiles is obtained. The third-order correlations imply that in the inner layer of wall wakes, a streamwise acceleration is prevalent and associated with a downward flux, suggesting sweeps. In contrast, in the outer layer, a streamwise deceleration exists and is associated with an upward flux, suggesting ejections. The profiles of the energy budget show that the turbulent and pressure energy diffusions oppose each other. The turbulent production has a positive peak, and the pressure energy diffusion has a negative peak, indicating a large gain in turbulence production in the wall-wake flows. The quadrant analysis confirms that in wall-wake flows, sweeps are the governing mechanism resulting from an inrush of fluid streaks. The bursting events have shorter duration, but they are more frequent than those in upstream.     

Experimental Study of Three-Dimensional Flow Field around a Complex Bridge Pier

In this paper, three-dimensional turbulent flow field around a complex bridge pier placed on a rough fixed bed is experimentally investigated. The complex pier foundation consists of a column, a pile cap, and a 2×4 pile group. All of the elements are exposed to the approaching flow. An acoustic-Doppler velocimeter was used to measure instantaneously the three components of the velocities at different horizontal and vertical planes. Profiles and contours of time-averaged velocity components, turbulent intensity components, turbulent kinetic energy, and Reynolds stresses, as well as velocity vectors are presented and discussed at different vertical and horizontal planes. The approaching boundary layer at the upstream of the pile cap separated in two vertical directions and induced an upward flow toward the column and a contracted downward flow below the pile cap and toward the piles. The contracted upward flow on the pile cap interacts with downflow in the front of the column and deflects toward the side of the pier, which in return produces a strong downflow along the side of the pile cap. The flow at the rear of the pile cap is very complex. The strong downward flow at the downstream and near the top of the pile cap in interaction with the reverse flow behind the column and upward flow near the bed produce two vortices close to the upper and lower corners of the pile cap with opposite direction of rotation. On the other hand, the back-flow from the wake of the pile cap is forced into the top region resulting in a secondary recirculation at the wake of the column. The contracted flow below the pile cap and toward the piles, a strong downflow along the sides of the pile cap at the upstream region, and a vortex flow behind the pile cap and an amplification of turbulence intensity along the sides of the pile cap at the downstream region are the main features of the flow responsible for the entrainment of bed sediments.