Next Article in Journal
Numerical Analysis of an Electroless Plating Problem in Gas–Liquid Two-Phase Flow
Next Article in Special Issue
Toward a Better Understanding of Sediment Dynamics as a Basis for Maintenance Dredging in Nagan Raya Port, Indonesia
Previous Article in Journal
Velocity Profile Representation for Fully Developed Turbulent Flows in Pipes: A Modified Power Law
Previous Article in Special Issue
The Hydrodynamic Moment of a Floating Structure in Finite Flowing Water
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Characterization of the Flow Field around Oblong Bridge Piers

1
Water Resources and Hydraulic Structures Division, Hydraulics and Environment Department, National Laboratory of Civil Engineering, 1700-066 Lisboa, Portugal
2
Hydraulics, Water Resources and Environmental Division, Civil Engineering Department, Faculty of Engineering of the University of Porto, 4400-465 Porto, Portugal
3
Interdisciplinary Centre of Marine and Environmental Research of the University of Porto, 4450-208 Matosinhos, Portugal
*
Author to whom correspondence should be addressed.
Current address: Civil Engineering Department, ISISE, School of Engineering, University of Minho, 4800-058 Guimarães, Portugal.
Fluids 2021, 6(11), 370; https://doi.org/10.3390/fluids6110370
Submission received: 7 September 2021 / Revised: 10 October 2021 / Accepted: 12 October 2021 / Published: 20 October 2021
(This article belongs to the Special Issue Environmental Hydraulics, Turbulence and Sediment Transport)

Abstract

:
The prediction of scour evolution at bridge foundations is of utmost importance for engineering design and infrastructures’ safety. The complexity of the scouring inherent flow field is the result of separation and generation of multiple vortices and further magnified due to the dynamic interaction between the flow and the movable bed throughout the development of a scour hole. In experimental environments, the current approaches for scour characterization rely mainly on measurements of the evolution of movable beds rather than on flow field characterization. This paper investigates the turbulent flow field around oblong bridge pier models in a well-controlled laboratory environment, for understanding the mechanisms of flow responsible for current-induced scour. This study was based on an experimental campaign planned for velocity measurements of the flow around oblong bridge pier models, of different widths, carried out in a large-scale tilting flume. Measurements of stream-wise, cross-wise and vertical velocity distributions, as well as of the Reynolds shear stresses, were performed at both the flat and eroded bed stages of scouring development with a high-resolution acoustic velocimeter. The time-averaged values of velocity and shear stress are larger in the presence of a developed scour hole than in the corresponding flat bed configuration.

Graphical Abstract

1. Introduction

The presence of obstacles such as bridge foundations in rivers is often associated with local morphological changes as a result of altering the sediment–flow equilibrium at the vicinity of those hydraulic infrastructures [1,2]. The modification of local geomorphological bed characteristics, known as local scour, should be accounted for when analysing the safety of bridges due to the severe threat that the uncovering of the foundation can pose to its integrity. Local scour around piers is one of the major threats to bridge structures founded in riverbeds and is the result of the complex flow patterns around those foundations [3,4,5]. It is therefore important to be acquainted with the flow structure and the related scour mechanisms around bridge piers.
The formation of highly turbulent vortex systems leads to increased local velocities and bottom shear stress, which results in the formation of scour holes. The sediment entrainment and transport are then controlled by the turbulent flow field therein developed [6,7]. Once formed, the scour hole interferes significantly with the turbulent flow field and the hydrodynamic characteristics [8,9]. For instance, the bed shear stress and turbulence quantities are strengthened, which induces further sediment transport and scour development in the pier vicinity.
Besides the formulation of scour depth predictors, the understanding of complex turbulence patterns in the vicinity of bridge piers is also crucial in planning and designing phases of foundations’ orientation, shape and geometry [6]. Several measuring techniques have been used before to characterize the flow field around piers. They began with the advent of propeller gauges, electromagnetic gauges, hot-film anemometers [10] and five-hole pitot spheres [11], and then developed to more sophisticated measuring equipment, such as Acoustic Doppler Velocimeter (ADV) [8,9,12,13], Laser Doppler Velocimetry (LDV) [14] and Particle Image Velocimetry (PIV) [15,16,17].
Ataie-Ashtiani and Aslani-Kordkandi [18] investigated the ADV measurements of turbulent flow field around single, side-by-side and tandem pile groups with and without a scour hole. Their results corroborated that the in-depth understanding of the flow and turbulence influenced by the scour hole is essential to improving the prediction of scour depth, see [19,20,21] and references therein.
The main flow processes that take place around a circular pier exposed to a steady current on a plane bed have been the subject of several flow field studies, both experimental and numerically, see [6], and references therein. They comprise the horseshoe vortices formed immediately upstream of the pier, the lee-wake vortices (usually in the form of vortex shedding), contracted streamlines at the side edges of the pier, and the downflow process due to the deceleration of the flow at the pier front. Dargahi [7] reported the result of several studies that focused on the visualization of the vertical structures around a bridge pier placed in an incoming fully turbulent flow at different stages of the scouring process. Graf and Yulistiyanto [22] measured the mean velocity components and turbulent Reynolds stresses using an ADV profiler in several planes around a circular pier for equilibrium scour conditions. By contrast, Roulund et al. [23] performed experimental and numerically visualizations of the flow in the scour hole during the evolution of the scour process and measured flow field and bed shear stresses for circular bridge piers on a flat bed. Dey and Raikar [8] have studied experimentally the structure of the turbulent horseshoe vortex flow within the developing scour holes at cylindrical piers by measuring instantaneous 3D velocities using an ADV. Notwithstanding, the flow in the rear part of the pier was not studied; they concluded that the flow and turbulence intensities in the horseshoe vortex remain rather similar during the development of the scour hole. While these flow structures are well-established for circular pier shape, few works have been devoted to the flow field study around bridge foundations of other shapes and geometries.
Recent studies have been focused on the characteristics of turbulence causing scour around circular, tandem, staggered arrangements and oblong piers [24,25], which did not account for the presence of scour holes. In Vijayasree et al. [24], the turbulent kinetic energy, turbulence intensities and Reynolds shear stresses were weaker for oblong piers, suggesting that this shape is a better alternative for the construction of bridge piers. The experimental work conducted by Pasupuleti et al. [25] was performed under clear water as the majority of previous research, which does not totally reproduce the field conditions. Therefore, the relatively limited number of flow field studies induced by the presence of an oblong bridge pier, most common shape of Portuguese bridge foundations of the 19th and 20th centuries, when compared with circular piers, motivated this investigation.
Therefore, the current study aims to identify and analyse the physics of the flow disturbed by the presence of a developed scour hole around oblong bridge piers. For that purpose, flow and turbulence characteristics are analysed, on the basis of an experimental campaign, concerning the change in pier effective width, approach hydraulic conditions and bed morphology stages, at the initial flat and the final eroded bed stages. Six fixed bed experiments, involving models of oblong bridge piers, were carried out in a recirculating tilting flume located in the experimental facilities of the National Laboratory of Civil Engineering (LNEC), in Lisbon, Portugal. Based on velocity measurements, the flow fields at the flume mid-plane surface were analysed, in terms of the time-averaged velocities and Reynolds shear stresses. Furthermore, space-averaged quantities were also obtained to clarify the influence of scour development on the flow and turbulence characteristics.
Following this introduction, the flow field experiments are characterized in Section 2. The description of the experimental facility and bridge pier models, the devices, measured variables and instrumentation, experimental procedure and conditions are given in Section 2. A particular description of the approach used for the characterization of the turbulent flow field is included in Section 3. The results and their discussion are presented in Section 4. Finally, the main conclusions of this study are provided in Section 5.

2. Flow Field Experiments

2.1. Experimental Facility and Bridge Pier Models

The experiments were undertaken in the tilting flume of LNEC (Figure 1). The tilting flume, herein termed CIV, is 40.7 m long, 2.0 m wide and 1.0 m deep, with glass wall sides and a concrete base.
Figure 2 presents the schematic setup of the flume. Two working sections were created within the base of the tilting flume, each being approximately 5.0 m long and 0.4 m deep. These sections were filled with uniform quartz sand characterized by a median diameter, D 50 , of 0.86 mm, a density, ρ s , of 2650 kgm 3 and a standard deviation of sediment particle sizes, σ D , of 1.28. The working sections, herein also referred to as “sand recess boxes”, are preceded by 2.0 m long accelerating ramps and 7.0 m long fixed bed reaches, and are followed by 3.0 m long fixed bed reaches. The accelerating ramps (“Acc. ramp” in Figure 2), placed at flume inlet and at the stilling recess boxes, conduce to uniform flow distributions in the two working sections. The stilling recess area between the two working sections was conceived to restore the flow regime for the downstream working section, and to store the sand bed material transported by the flowing water from the upstream working section during local scouring experiments. In the upstream boundary of the working sections, fine gravel mattresses, properly levelled by the respective concrete base, allow for smooth and progressive transitions, as well as ensure a turbulent flow development.
Two oblong bridge pier models, P i e r 11 and P i e r 14, were placed in the upstream and downstream working sections, respectively. The piers differ in width (W), length (L), and semi-cylindrical surface ratio (R), schematically identified in Figure 3. The dimensions of P i e r 11 are W = 0.11 m, L = 0.433 m and R = 0.065 m. P i e r 14 is slightly larger, with W = 0.14 m, L = 0.463 m, and R = 0.08 m. The bridge pier models were placed vertically on the flume mid-plane at ≈ 2 m from the working sections’ upstream borders (Figure 2). The contraction and wall effects were negligible to the scouring process, since the associated flume “width-to-pier width” ratios ( W / B , B being the flume width) were 5.5% and 7.0% (for P i e r 11 and P i e r 14, respectively), both lower than the limit value of 10% indicated by Chiew and Melville [26].
The present study comprises the hydrodynamic characterization of six experiments divided into fixed flat bed and eroded bed stages of three local scouring experiments, namely E x p . 1U, E x p . 1D and E x p . 2U. The letters U and D denote the pier model, being U for P i e r 11 and D for P i e r 14. These local scouring experiments are characterized in detail in Bento [27,28]. For simplicity, a suffix was added to the local scouring experiments’ name. The suffix F b e d stands for the fixed flat bed, whilst E b e d denotes the fixed eroded bed experiments. Therefore, the six flow field experiments were labelled as follows: E x p . 1U F b e d and E x p . 1U E b e d for E x p . 1U; E x p . 1D F b e d and E x p . 1D E b e d for E x p . 1D; and E x p . 2U F b e d and E x p . 2U E b e d for E x p . 2U.

2.2. Devices, Measured Variables and Instrumentation

Several flow variables were measured during the experiments, namely: (i) the flow discharge at the flume entrance; (ii) the flow depth within CIV and at the immediacy of the working sections; and (iii) the instantaneous three velocity components of the flow at the piers’ vicinity. The following paragraphs describe the measuring process of each variable.
A variable frequency drive was used to control the rotational speed of the pump motor in order to obtain the desired flow discharge. The approach flow discharge was then monitored by an electromagnetic flowmeter with a rate accuracy of ±0.25%. At the flume entrance, stilling grids and a flow straightener were provided for damping the flow containing large scale eddies, leading to vortex free uniform flow in CIV. An adjustable sluice gate at the downstream end of the flume enables obtaining the desirable flow depth within the flume.
The approach flow depth was monitored using the acquired signals from acoustic probes, suitably positioned within the flume so as not to interfere with the development of the experiments. The signal of the three acoustic probes was acquired at a sampling frequency of 10 Hz, whereas the signal of the electromagnetic flowmeter was acquired at 25 Hz. Movable carriages were present on the rails of the flume to allow the access for controlling and measuring purposes of both the approach flow depth and the acquisition of the flow velocity components by the connection of the downlooking vectrino to its respective acquisition firmware (Vectrino Plus).
The instantaneous three-dimensional velocity measurements were acquired by using a high-resolution acoustic velocimeter, namely a downlooking vectrino (hereinafter designated as vectrino), from Nortek®, in two different moments: (i) at the beginning of the scouring process (flat bed stage), and (ii) at the final stage of the scouring experiments (eroded bed stage). The velocity flow field was measured throughout the working section around the pier models at relative elevations for the comparison of turbulence and flow characteristics due to the scouring process (before scour vs. after scour). An independent moving carriage system was developed to clamp the vectrino in order to move it along the different velocity measurement grid points to ease the measurement task, up to an accuracy of ±1 mm in the three directions (x, y and z). This issue will be considered in more detail in Section 3. The origin was considered at the pier centre, where x, y and z denote the longitudinal, transversal and vertical axes, respectively (see Figure 3).

2.3. Experimental Procedure

To facilitate the performance of the experiments on fixed beds, the surface layers of the bed material (Step 1—Figure 4a) were frozen by spraying uniformly a synthetic glue over the bed morphology stages (flat and eroded). This practice is frequently used in flow field studies under the assumption that the results are simpler and more explanatory than those from complex scenarios with a variety of acting processes, such as sediment burial and reappearance processes [29], and references therein.
Once the sand beds were properly stabilized, after approximately 72 h, the same hydraulic conditions of the corresponding local scouring experiments [27,28] were imposed, and the flume was filled up to the desired flow discharge and flow depth (conditions referred later). Since the vectrino uses the Doppler effect to measure current velocity by transmitting short pairs of acoustic waves, they must bounce off from particulate matter in the flow. As the laboratory water was quite clean, an artificial introduction of particles (commonly known as seeding particles) was deemed necessary, and thus a silica powder was added to the flow. An indicative aspect of the quality of the seeded flow is the parameter signal-to-noise ratio (SNR), provided by the Vectrino Plus interface; the higher the SNR, the better the seeding, and consequently, the more reliable the velocity measurements (Step 2—Figure 4b).
Five acoustic probe tips (one transmitter and four receivers, illustrated in Figure 4c) send and receive acoustic information in a fluid volume below the vectrino, referred to as the sampling volume. Since the instrument is inserted into the flow, the sampling volume must be sufficiently far away from the probe tip so as not to disturb the flow and the velocity measurements. Using acoustic waves, the vectrino takes real-time measurements (Step 3—Figure 4c). To evaluate quality of velocity measurements, one of the real-time outputs that vectrino provides is a correlation statistic. This is a value ranging from 0% to 100%; the closer the correlation is to 100%, the less noisy and more reliable a velocity measurement is. A minimum correlation of 70% is strongly recommended by several researchers [12], among others, and thus considered in the present study. In addition, a minimum signal-to-noise ratio ( S N R ) of 15 dB was designated for ensuring the reliability of the instantaneous velocity data [30,31].
Before conducting each experiment, several trials were undertaken to estimate the uncertainty in the measurement sampling duration. To tackle this, an assessment of the signal acquisition time’s influence on the statistics of turbulent quantities was performed in Bento et al. [32]. Thus, signal acquisition times within the range of 120 to 300 s were considered in order to achieve time-independent velocity components, regarding the magnitude of turbulence at each measuring location. The vectrino was operated with an acoustic frequency of 10 MHz and a sampling rate of 200 Hz.
Besides the flow velocity characterization of the six fixed bed experiments by a comprehensive measuring system (resumed later), the flow depth was also continuously monitored for ensuring the predefined hydraulic conditions. After the completion of each fixed bed experiment, the surface layers of the sand beds were removed, the sand box re-levelled and once again fixed by the synthetic glue (to perform the corresponding flat bed experiment). A careful re-establishment of the hydraulic conditions was performed with a gradual filling of the flume.

2.4. Hydraulic Conditions and Time Duration

Table 1 presents the main hydraulic variables of the fixed bed experiments, namely the approach flow discharge (Q), flow depth (h), flow velocity (V), Froude number ( F r = V / g h ), flow Reynolds number ( R e = V h / ν , being ν the kinematic viscosity of the fluid), and the parameter concerning the relation with the pier model: the pier Reynolds number ( R e D = V W / ν ). In Table 1, these parameters are given with a 95% of confidence interval; the estimated error corresponded to the sample error, expressed as 1.96 σ / n s (being σ the standard deviation and n s the sample size).
The flow depth within CIV was considered as the arithmetic average of the acoustic probes data for E x p . 1U F b e d , E x p . 1U E b e d , E x p . 1D F b e d , E x p . 1D E b e d , E x p . 2U F b e d and E x p . 2U E b e d experiments. These analyses constitute a relevant and efficient source of information to account for the uncertainty in the hydraulic conditions of each experiment. The characterization of the experimental error and propagation of uncertainty demonstrated that the variability of the hydraulic variables, particularly flow discharge, flow depth, and flow velocity, had no impact on the final results (Table 1).
The oblong bridge pier models behaved as narrow or transitional [1], depending on the resultant flow shallowness ratios ( h / W ). A flow shallowness of 1.13 was obtained for E x p . 1D F b e d and E x p . 1D E b e d (transitional piers), while values varying between 1.40 and 1.48 were attained for the remaining experiments (narrow pier category).
The flume width-to-flow depth ratios ( B / h ) ensured negligible effects of the secondary currents due to flume side walls on the flow. Values of 12.4, 12.3, 12.6, 12.7, 10.2 and 13.5 were obtained, respectively, for E x p . 1U F b e d , E x p . 1U E b e d , E x p . 1D F b e d , E x p . 1D E b e d , E x p . 2U F b e d and E x p . 2U E b e d (Table 1). In accordance with Yang et al. [33], that effect is solely noticeable for ratios of less than 6. The flow was, thus, free from the secondary currents at the flume mid-plane [34], where the pier models were installed. Fully developed incoming flows were obtained for all experiments, as indicated by the magnitudes of the flow and pier Reynolds number ( R e and R e D ) in Table 1. The values of Froude number ( F r ) lower than unity confirms the subcritical flow for experiments.

3. Flow Field Characterization

In the present research, instantaneous velocity maps were acquired for the fixed bed experiments, before the beginning of the scouring process, at a flat bed, and at the final stage of experiments, at an eroded bed. Figure 5 provides a schematic diagram of the velocity point measurements. It comprises five half-planes, namely A, J, E, G and I. Each half-plane is composed by 4 or 5 vertical profiles, as indicated in Figure 5, where the flow velocity was measured for several points along the water depth. The position of the measuring spots varied for each application; it was dependent on flow depth and pier dimensions. In line with what was performed in the works of Kirkil et al. [35], Chang et al. [36], and Beheshti and Ataie-Ashtiani [13], among others, the measurements in each experiment were taken in the Cartesian coordinate system (x, y, z) also indicated in Figure 5. Further details can be found in Bento [28].
A minimum correlation of 70% and a minimum signal-to-noise ratio ( S N R ) of 15 dB were thus considered for ensuring the reliability of the instantaneous velocity data [30,31]. During the velocity measurements, the downlooking vectrino sampling frequency was chosen to be 200 Hz with a duration within the range of 120 to 300 s, depending on the magnitude of turbulence at each measuring location. These instantaneous three-dimensional velocity measurements allowed obtaining: (i) mean velocities (u, v and w); (ii) turbulent fluctuation velocities ( u , v and w ); and (iii) Reynolds shear stresses ( R S S e s ) including u v ¯ , u w ¯ and v w ¯ . The derivation of such parameters provides insights of the flow magnitude and structure for the particularities of each experiment and bed morphology stage (before v s after scour).
Prior to data analysis, all records were visually inspected to identify potential erroneous data such as spikes, trends or abrupt discontinuities in the velocity time series. The velocity data was then processed, which streamlined the analysis of the turbulent and the time-averaged 3D velocity data, as well as the statistical and the spectral analysis of the downlooking vectrino data by computing the respective skewness and kurtosis measures for each measuring point. Results including mean values, variance, skewness and kurtosis of each of the aforementioned variables, for all experiments, can be found in Database S1 in the Supplementary Materials. For the sake of simplicity, a limited number of the most relevant plots is analysed in Section 4, namely for the different vertical profiles of the half-planes A and J. A total of 680 velocity point measurements were acquired for these two half-planes. For each vertical profile, a relative study on the characteristic flow and turbulence quantities is conducted, regarding the bridge’s effective width, approach flow condition and bed configuration (flat v s eroded), addressed in Section 4.
A characterization of the approach flow was also performed. For fully turbulent flow, over hydraulically rough sand bed surface, the dimensionless stream-wise mean velocity ( u / u * , being u * the shear velocity), was found to follow the log-law as expected, see [37] and references therein. The roughness coefficient value adopted, k s , was 0.00215 m ( k s = 2.5 D 50 ). The adjustment of the log-law to E x p . 1U E b e d velocity profile upstream the P i e r 11’s front, for instance, allowed obtaining the shear velocity, u * , of 0.0174 and constant B s = 9.74, as shown in Figure 6. These values, along with the bed Reynolds number, u * k s / ν ≈ 37, are in agreement with the literature for the hydraulically transitional boundary layer [38].

4. Results and Discussion

4.1. Mean Velocities

Figure 7, Figure 8 and Figure 9 depict the profiles of the velocity components at the eight mid-plane vertical profiles—A and J half-planes (Figure 5)—for the three flat bed experiments: E x p . 1U F b e d , E x p . 1D F b e d and E x p . 2U F b e d . The velocity points were scaled by a reference velocity magnitude of 0.30 ms 1 .
In the upstream area, the stream-wise component of the velocity (u) increases with the distance from the bed to an elevation near the flow height. Along the mid-plane surface, it decreases with proximity to the oblong bridge pier model, due to the effect of velocity reduction close to the pier’s blocking effect (Figure 7a, Figure 8a and Figure 9a). The flat bed experiments, E x p . 1U F b e d , E x p . 1D F b e d and E x p . 2U F b e d , exhibit nearly the same pattern. The u vertical profile is diminished for E x p . 1D F b e d when compared with the E x p . 1U F b e d experiment. This can be attributed to large flow shallowness of E x p . 1U F b e d ( h / W = 1.45) when compared with the correspondent value for E x p . 1D F b e d ( h / W = 1.14). The higher mean approach flow velocity of E x p . 2U F b e d (in the order of 20%) justifies an increment of approximately 34% of the mean values of each of the four vertical profiles upstream the bridge pier model ( P i e r 11).
In E x p . 1U F b e d , the maximum stream-wise component of velocity, u, occurred for an elevation (z) equal to 0.097 m in all the vertical profiles at the upstream side of the oblong pier ( P i e r 11); the magnitude reached a maximum value of 0.377 ms 1 (skewness and kurtosis coefficients of −0.70 and 1.02, respectively) for the vertical profile “1” of the A half-plane (situated 0.5 m upstream from the pier front edge, Figure 4), whereas a magnitude of 0.313 ms 1 (skewness and kurtosis coefficients of 2.55 and 39.43, respectively) was attained for the vertical profile “4” of the A half-plane (at roughly 0.09 m upstream from the pier front edge. For E x p . 2U F b e d , those values registered a maximum of 0.509 ms 1 and 0.422 ms 1 for the vertical profiles “1” and “4” of A, respectively. The corresponding computed skewness and kurtosis coefficients were −0.94 and 2.31 for vertical profile “1” and −1.99 and 17.42 for vertical profile “4”).
Statistical moments, such as skewness and kurtosis coefficients, differing from zero means that the temporal distributions of u (for the previously mentioned example) were slightly skewed and sharpened, especially for the points measured at the vertical profile “4”. For the corresponding velocity-point measurements, the v statistical moments were closer to zero (see Database S1: Flow field database). Even though observed deviations from the Gaussian distribution [39,40] are relatively small, they may reflect significant features of the large-scale turbulence dynamics. According to Nikora et al. [39] and references therein, large-scale vortices or coherent structures enhance the kurtosis coefficient, which is also corroborated by the aforementioned statistical results. Further comparisons can be undertaken for the remaining experiments.
The stream-wise component of velocity in the downstream region of the pier models presents a similar pattern between the three experiments: a general increase of u as the vertical profiles move away from the oblong bridge pier. However, it should be noted that for profile “1” of the J half-plane (see Figure 7a, Figure 8a and Figure 9a), in all three flat bed experiments, alternate positive and negative values were observed along the vertical profile, which can be attributed to the effect of the shedding vortexes therein developed. In line with Figure 5, the vertical profiles “1” of J half-planes are situated at 0.087 m and 0.093 m from the downstream end of P i e r 11 and P i e r 14, respectively.
The time average v values at the symmetry plane (y = 0) are nearly zero, as expected (Figure 7b) and reported in previous research, such as Beheshti and Ataie-Ashtiani [13] and Pandey et al. [41]. For that reason, the v vertical distributions from the remaining flat bed experiments ( E x p . 1D F b e d and E x p . 2U F b e d ) are not herein depicted.
The vertical profiles of the w—component velocity, represented in Figure 7c, Figure 8b and Figure 9b, illustrate the downward movement of flow (downflow) in front of the respective oblong bridge pier models. The w component is highest for E x p . 2U F b e d , as compared to the other two experiments ( E x p . 1U F b e d and E x p . 1D F b e d , respectively). In those, the smaller values are attributed to the areas with less scouring near the pier as the stream bed surface did not permit the flow acceleration in the downward direction of the pier model. These results are in line with previous studies such as Dey and Raikar [8].
Figure 10, Figure 11 and Figure 12 illustrate the distribution of the velocity components at the same eight mid-plane vertical profiles, from the A and J half-planes (as previously performed), for the three eroded bed experiments: E x p . 1U E b e d , E x p . 1D E b e d and E x p . 2U E b e d . Similarly, the velocity points were scaled by a reference velocity magnitude of 0.30 ms 1 , except for E x p . 2U E b e d , where a value of 0.40 ms 1 was considered.
Unlike in other flow field experimental studies, see [21] and references therein, the velocity measurements were also allowed within the scour hole of the eroded bed experiments, which explains the negative values for elevations z below the initial flat bed level (see Figure 10, Figure 11 and Figure 12). As can be observed, vertical profile “1” of A is situated outside the upstream border of the scour hole in all experiments. The remaining vertical profiles of A and J half-planes are placed within the eroded scour holes.
The stream-wise velocity (u) profiles are marked by negative velocity within the developed scour holes, which indicates the occurrence of flow separation and recirculation, as can be seen in Figure 10a, Figure 11a and Figure 12a. Besides the roughly same approach hydraulic conditions of E x p . 1U E b e d and E x p . 1D E b e d , the corresponding flow shallowness values fit into distinct pier categories; the pier models behave as narrow pier for E x p . 1U E b e d and transitional pier category for E x p . 1D E b e d . As initially stated by Raudkivi [42] and corroborated by several other researchers ([43], and references therein), the parameter h / W affects the vorticity, which is associated with the generation of scour. For decreasing flow shallowness values, the influence of the surface roller on the bed increases, damping the vortices in front of the pier and consequently explaining the slightly lower u values.
A region of reversal flow was noticeable at vertical profile “4” of the A half-plane (immediately upstream from the pier), then it moved further upstream from the pier towards the scour hole edge. The u magnitudes generally increase with larger scour hole dimensions; thus, higher values were obtained for E x p . 2U E b e d . This can be explained that as the scour evolves upstream, the bed slope increases due to the formation of the scour hole, and the flow velocity increases correspondingly [21]. In Figure 10a, Figure 11a and Figure 12a, the downstream vertical profiles of u (from J half-plane in Figure 5) are characterized by increasing magnitudes as they move away from the oblong bridge pier models.
As performed for E x p . 1U F b e d , the v profile of E x p . 1U E b e d (shown in Figure 10b) was analysed and a similar behaviour (showing negligible changes) was observed. Hence, the remaining v profiles were not represented.
The vertical profiles of the w-component velocity illustrate the downward movement of flow (downflow) in front of the pier (attaining its maximum value for vertical profile "4" of A half-plane in all eroded bed experiments). The downward vertical velocity increased with the scour hole dimension, registering the maximum value for E x p . 2U E b e d . This finding is consistent with the results from the measurement of Dey and Raikar [8] and Li et al. [21]. Downstream of the oblong bridge pier models, w decreases with x / h for all experiments ( E x p . 1U E b e d , E x p . 1D E b e d and E x p . 2U E b e d ), thus reducing the lifting action. This pattern is more pronounced for the eroded bed experiments (with developed scour holes) than for the corresponding flat bed configurations.
To sum up, the mean flow velocities for the six fixed experiments indicate that with the developed scour hole, the downward flow in front of the oblong bridge pier model becomes more prominent, and the flow becomes more turbulent due to the effect of a turbulent horseshoe vortex within the scour hole, as documented by previous researchers [21,43], among others.

4.2. Reynolds Shear Stresses

The time-averaged Reynolds shear stresses, u v ¯ , u w ¯ and v w ¯ , are depicted in Figure 13, Figure 14 and Figure 15 for experiments E x p . 1U F b e d , E x p . 1D F b e d and E x p . 2U F b e d , respectively.
The comparative study of u v ¯ , u w ¯ and v w ¯ along the pier’s mid-plane for the eroded bed experiments was also performed, as illustrated in Figure 16, Figure 17 and Figure 18 for E x p . 1U E b e d , E x p . 1D E b e d and E x p . 2U E b e d , respectively.
Bearing in mind that the Reynolds shear stresses are linked to the capacity of the stream to transport sediments and indicate the exchange of turbulent momentum between layers of a stream [44,45], it is understandable that such capacity increases with the developing scour holes, and the region with large value enlarges and moves upstream of the correspondent scour hole.
According to Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18, no significant changes were observed in the vertical profiles of Reynolds shear stresses for different bed configurations (flat v s eroded). In the downstream region, the magnitudes of u v ¯ are greater than the corresponding points at the upstream side of the bridge pier model at both bed configurations (flat v s eroded). This behaviour can be attributed to the turbulent mix of fluid as a result of vortical flow in the rear side of the pier. For vertical profiles of the J half-plane (Figure 5), the complex flow is distinctly identified by the saw-toothed trend of the Reynolds shear stress u w ¯ , magnified for the eroded fixed bed experiments. The values of v w ¯ are generally less pronounced than the other components of the Reynolds shear stress, u v ¯ and u w ¯ . These results, from Figure 13, Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18, demonstrate that the usual assumption of isotropic turbulence is not exactly true, although it is most used as the basis of much of the theoretical analysis of turbulent flow.
Therefore, the location of the regions of the maximum Reynolds shear stresses ( u w ¯ vertical profiles) occurs at the interfaces of the negative and positive velocity components, meaning that the velocity in this region can be very unstable. This pattern was also observed in Guan et al. [46], where the turbulent flow field, in a developing clear water (with a flow intensity of 0.53), was quantified at a circular pier.
In spite of several studies discussing the individual trends of each component of the Reynolds shear stresses, there are few authors that argue that on the mid-plane (i.e., y = 0), the Reynolds shear stresses involving the lateral fluctuations v (i.e., u v ¯ and v w ¯ ) must be zero for symmetry reasons. In the present study, those components differed from zero mainly for the vertical profiles of the J half-plane. Nevertheless, this issue has not been yet thoroughly understood in the literature. It seems that theoretically the non-zero values are solely related with a signal acquiring time, of the instantaneous velocity components, not long enough to obtain reliable statistics. Further and detailed investigations regarding that issue are still needed.

5. Conclusions

This paper investigates the turbulent flow field around oblong bridge pier models in a well-controlled laboratory environment, for understanding the mechanisms of flow responsible for current-induced scour. Two different bed morphology configurations, namely (i) at the beginning of the scouring process (flat bed) and (ii) at the final stage of the scouring experiments (eroded bed) of the three local scouring experiments were herein considered. A total of six fixed bed experiments at the vicinity of the oblong bridge pier models of distinctive effective width and corresponding flow shallowness were carried out.
Measurements of the instantaneous stream-wise, cross-wise and vertical components of the flow velocity were performed by using a high-resolution acoustic velocimeter (downlooking vectrino) at eight mid-plane vertical profiles. The results showed that the presence of the developed scour hole, by creating extra flow resistance and augmenting the velocity difference within the four vertical profiles at the upstream side of the pier, changes the velocity distribution pattern and the Reynolds shear stresses. At the rear part of the pier, significant differences in flow field patterns, due to the presence of the developed scour holes, were also observed. Their magnitudes were in accordance with the approach flow conditions, such as the corresponding depth and velocity.
An important contribution was made to the understanding of the physics of the flow and turbulence characteristics with respect to the change in pier effective width, approach hydraulic conditions and bed morphology stages. The derivation of mean velocities and Reynolds shear stresses, around the oblong bridge pier models, contributes to a comprehensive description of the flow magnitude and structure for the particularities of each experiment and bed morphology stage (flat v s eroded). The collected data may also constitute a basis for the calibration and validation of numerical models, frequently used to characterize the scouring phenomenon, such as Computational Fluid Dynamics (CFD) tools.
In future work, a spectra analysis of the velocity measurements, and corresponding Reynolds shear stresses, at several points, both at the upstream and downstream sides of the oblong bridge pier models, can be performed. Such analysis aids the evaluation of the shedding frequency of the large scale coherent structures and further understands the turbulent dynamics, and the inertial sub-ranges of turbulence due to the scouring phenomena.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/fluids6110370/s1, Database S1: Flow field database.

Author Contributions

Conceptualization, A.M.B., T.V., J.P.P. and L.C.; methodology and design, A.M.B., T.V., J.P.P. and L.C.; formal analysis, A.M.B.; investigation, A.M.B.; writing—original draft preparation, A.M.B.; writing—review and editing, T.V., J.P.P. and L.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Portuguese Foundation for Science and Technology (FCT) grants number PD/BD/127798/2016, UIDB/04423/2020 and UIDP/04423/2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on reasonable request from the authors.

Acknowledgments

The first author thanks the financial support of the Portuguese Foundation for Science and Technology (FCT), through the Ph.D. scholarship PD/BD/127798/2016. This research was co-supported by strategical funding from FCT UIDB/04423/2020 and UIDP/04423/2020. We would also like to thank the team of the Laboratory of Construction and Modelling of the National Laboratory of Civil Engineering (LNEC) for the technical assistance in the experimental work.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BFlume width
B s Constant of log-law profile
DPier Model ( P i e r 14)
D 50 Median grain size diameter
F r Froude number
gGravity acceleration
hFlow depth
K s Roughness coefficient
LOblong bridge pier length
n s Sample size
QFlow discharge
RPier semi-cylindrical surface ratio
R e Flow Reynolds number
R e D Pier Reynolds number
UPier Model ( P i e r 11)
uStream-wise velocity component
u Stream-wise velocity fluctuation
u * Shear velocity
VApproach flow velocity
vCross-wise velocity component
v Cross-wise velocity fluctuation
WOblong bridge pier width
wVertical velocity component
w Vertical velocity fluctuation
xLongitudinal axis
yTransversal axis
zVertical axis
ν Fluid kinematic viscosity
ρ s Sediment density
σ Standard deviation
σ D Standard deviation of sediment particle sizes
A D V Acoustic Doppler Velocimeter
C F D Computational Fluid Dynamics
L N E C National Laboratory of Civil Engineering
L D V Laser Doppler Velocimetry
P I V Particle Image Velocimetry
R S S e s Reynolds Shear Stresses ( u v ¯ , u w ¯ and v w ¯ )
S N R Signal-to-Noise Ratio

References

  1. Melville, B.W.; Coleman, S.E. Bridge Scour; Water Resources Publications, LLC: Highlands Ranch, CO, USA, 2000. [Google Scholar]
  2. Arneson, L.; Zevenbergen, L.; Lagasse, P.; Clopper, P. Evaluating Scour at Bridges; Technical Report; Federal Highway Administration: Fort Collins, CO, USA, 2012. [Google Scholar]
  3. Deng, L.; Wang, W.; Yu, Y. State-of-the-Art review on the causes and mechanisms of bridge collapse. J. Perform. Constr. Facil. 2016, 30, 04015005. [Google Scholar] [CrossRef]
  4. Flint, M.M.; Fringer, O.; Billington, S.L.; Freyberg, D.; Diffenbaugh, N.S. Historical analysis of hydraulic bridge collapses in the Continental United States. J. Infrastruct. Syst. 2017, 23, 04017005. [Google Scholar] [CrossRef] [Green Version]
  5. Proske, D. Bridge Collapse Frequencies Versus Failure Probabilities; Springer: Berlin, Germany, 2018. [Google Scholar]
  6. Ettema, R.; Constantinescu, G.; Melville, B.W. Flow-field complexity and design estimation of pier-scour depth: Sixty years since Laursen and Toch. J. Hydraul. Eng. 2017, 143, 03117006. [Google Scholar] [CrossRef]
  7. Dargahi, B. The turbulent flow field around a circular cylinder. Exp. Fluids 1989, 8, 1–12. [Google Scholar] [CrossRef]
  8. Dey, S.; Raikar, R.V. Characteristics of horseshoe vortex in developing scour holes at piers. J. Hydraul. Eng. 2007, 133, 399–413. [Google Scholar] [CrossRef]
  9. Kumar, A.; Kothyari, U.C.; Raju, K.G.R. Flow structure and scour around circular compound bridge piers—A review. J. Hydro-Environ. Res. 2012, 6, 251–265. [Google Scholar] [CrossRef]
  10. Melville, B.W.; Raudkivi, A.J. Flow characteristics in local scour at bridge piers. J. Hydraul. Res. 1977, 15, 373–380. [Google Scholar] [CrossRef]
  11. Dey, S. Three-dimensional vortex flow field around a circular cylinder in a quasi-equilibrium scour hole. Sadhana 1995, 20, 871–885. [Google Scholar] [CrossRef] [Green Version]
  12. Beheshti, A.A.; Ataie-Ashtiani, B. Experimental study of three-dimensional flow field around a complex bridge pier. J. Eng. Mech. 2010, 136, 143–154. [Google Scholar] [CrossRef]
  13. Beheshti, A.A.; Ataie-Ashtiani, B. Scour hole influence on turbulent flow field around complex bridge piers. Flow Turbul. Combust. 2016, 97, 451–474. [Google Scholar] [CrossRef]
  14. Devenport, W.J.; Simpson, R.L. Time-depeiident and time-averaged turbulence structure near the nose of a wing-body junction. J. Fluid Mech. 1990, 210, 23–55. [Google Scholar] [CrossRef]
  15. Said, N.M.; Mhiri, H.; Bournot, H.; Le Palec, G. Experimental and numerical modelling of the three-dimensional incompressible flow behaviour in the near wake of circular cylinders. J. Wind Eng. Ind. Aerodyn. 2008, 96, 471–502. [Google Scholar] [CrossRef]
  16. Apsilidis, N.; Diplas, P.; Dancey, C.L.; Bouratsis, P. Time-resolved flow dynamics and Reynolds number effects at a wall–cylinder junction. J. Fluid Mech. 2015, 776, 475–511. [Google Scholar] [CrossRef]
  17. Jenssen, U.; Manhart, M. Flow around a scoured bridge pier: A stereoscopic PIV analysis. Exp. Fluids 2020, 61, 1–18. [Google Scholar] [CrossRef]
  18. Ataie-Ashtiani, B.; Aslani-Kordkandi, A. Flow field around single and tandem piers. Flow Turbul. Combust. 2013, 90, 471–490. [Google Scholar] [CrossRef]
  19. Link, O.; González, C.; Maldonado, M.; Escauriaza, C. Coherent structure dynamics and sediment particle motion around a cylindrical pier in developing scour holes. Acta Geophys. 2012, 60, 1689–1719. [Google Scholar] [CrossRef]
  20. Kirkil, G.; Constantinescu, G. Effects of cylinder Reynolds number on the turbulent horseshoe vortex system and near wake of a surface-mounted circular cylinder. Phys. Fluids 2015, 27, 075102. [Google Scholar] [CrossRef] [Green Version]
  21. Li, J.; Yang, Y.; Yang, Z. Influence of scour development on turbulent flow field in front of a bridge pier. Water 2020, 12, 2370. [Google Scholar] [CrossRef]
  22. Graf, W.; Yulistiyanto, B. Experiments on flow around a cylinder; the velocity and vorticity fields. J. Hydraul. Res. 1998, 36, 637–654. [Google Scholar] [CrossRef]
  23. Roulund, A.; Sumer, B.M.; Fredsøe, J.; Michelsen, J. Numerical and experimental investigation of flow and scour around a circular pile. J. Fluid Mech. 2005, 534, 351–401. [Google Scholar] [CrossRef]
  24. Vijayasree, B.; Eldho, T.; Mazumder, B. Turbulence statistics of flow causing scour around circular and oblong piers. J. Hydraul. Res. 2020, 58, 673–686. [Google Scholar] [CrossRef]
  25. Pasupuleti, L.N.; Timbadiya, P.V.; Patel, P.L. Flow Field Measurements Around Isolated, Staggered, and Tandem Piers on a Rigid Bed Channel. Int. J. Civ. Eng. 2021, 1–18. [Google Scholar]
  26. Chiew, Y.M.; Melville, B.W. Local scour around bridge piers. J. Hydraul. Res. 1987, 25, 15–26. [Google Scholar] [CrossRef]
  27. Bento, A.M.; Couto, L.; Viseu, T.; Pêgo, J.P. Image-Based Techniques for the Advanced Characterization of Scour around Bridge Piers in Laboratory. Available online: https://www.e3s-conferences.org/articles/e3sconf/abs/2018/15/e3sconf_riverflow2018_05066/e3sconf_riverflow2018_05066.html (accessed on 4 September 2021).
  28. Bento, A.M. Risk-Based Analysis of Bridge Scour Prediction. Ph.D. Thesis, Faculdade de Engenharia da Universidade do Porto, Porto, Portugal, 2021. [Google Scholar]
  29. Lee, S.O.; Hong, S.H. Turbulence characteristics before and after scour upstream of a scaled-down bridge pier model. Water 2019, 11, 1900. [Google Scholar] [CrossRef] [Green Version]
  30. García, C.M.; Cantero, M.I.; Niño, Y.; García, M.H. Turbulence measurements with Acoustic Doppler Velocimeters. J. Hydraul. Eng. 2005, 131, 1062–1073. [Google Scholar] [CrossRef]
  31. Montero, V.G.G.; Romagnoli, M.; García, C.M.; Cantero, M.I.; Scacchi, G. Optimization of ADV sampling strategies using DNS of turbulent flow. J. Hydraul. Res. 2014, 52, 862–869. [Google Scholar] [CrossRef]
  32. Bento, A.M.; Pêgo, J.P.; Couto, L.; Viseu, T. Assessing the flow field around an oblong bridge pier. Vectrino acquisition time sensitivity analysis. Publ. Inst. Geophys. Pol. Acad. Sci. Geophys. Data Bases Process. Instrum. 2021, 434, 131–132. [Google Scholar] [CrossRef]
  33. Yang, S.Q.; Tan, S.K.; Lim, S.Y. Velocity distribution and dip-phenomenon in smooth uniform open channel flows. J. Hydraul. Eng. 2004, 130, 1179–1186. [Google Scholar] [CrossRef]
  34. Cardoso, A.H.; Graf, W.H.; Gust, G. Uniform flow in a smooth open channel. J. Hydraul. Res. 1989, 27, 603–616. [Google Scholar] [CrossRef]
  35. Kirkil, G.; Constantinescu, S.G.; Ettema, R. Coherent structures in the flow field around a circular cylinder with scour hole. J. Hydraul. Eng. 2008, 134, 572–587. [Google Scholar] [CrossRef]
  36. Chang, W.Y.; Constantinescu, G.; Lien, H.C.; Tsai, W.F.; Lai, J.S.; Loh, C.H. Flow structure around bridge piers of varying geometrical complexity. J. Hydraul. Eng. 2013, 139, 812–826. [Google Scholar] [CrossRef]
  37. Schlichting, H.; Gersten, K. Boundary-Layer Theory; Springer: Berlin/Heidelberg, 2017. [Google Scholar] [CrossRef]
  38. Cardoso, A.H.; Cunha, L.V.d. Hidráulica Fluvial; Fundação Calouste Gulbenkian: Lisbon, Portugal, 1998. [Google Scholar]
  39. Nikora, V.; Nokes, R.; Veale, W.; Davidson, M.; Jirka, G. Large-scale turbulent structure of uniform shallow free-surface flows. Environ. Fluid Mech. 2007, 7, 159–172. [Google Scholar] [CrossRef]
  40. Ouro, P.; Juez, C.; Franca, M. Drivers for mass and momentum exchange between the main channel and river bank lateral cavities. Adv. Water Resour. 2020, 137, 103511. [Google Scholar] [CrossRef]
  41. Pandey, M.; Sharma, P.K.; Ahmad, Z.; Singh, U.K.; Karna, N. Three-dimensional velocity measurements around bridge piers in gravel bed. Mar. Georesources Geotechnol. 2017, 36, 663–676. [Google Scholar] [CrossRef]
  42. Raudkivi, A.J. Functional trends of scour at bridge piers. J. Hydraul. Eng. 1986, 112, 1–13. [Google Scholar] [CrossRef]
  43. Pizarro, A.; Manfreda, S.; Tubaldi, E. The science behind scour at bridge foundations: A review. Water 2020, 12, 374. [Google Scholar] [CrossRef] [Green Version]
  44. Bernard, P.S.; Handler, R.A. Reynolds stress and the physics of turbulent momentum transport. J. Fluid Mech. 1990, 220, 99–124. [Google Scholar] [CrossRef] [Green Version]
  45. Carnacina, I.; Leonardi, N.; Pagliara, S. Characteristics of flow structure around cylindrical bridge piers in pressure-flow conditions. Water 2019, 11, 2240. [Google Scholar] [CrossRef] [Green Version]
  46. Guan, D.; Chiew, Y.M.; Wei, M.; Hsieh, S.C. Characterization of horseshoe vortex in a developing scour hole at a cylindrical bridge pier. Int. J. Sediment Res. 2019, 34, 118–124. [Google Scholar] [CrossRef]
Figure 1. Experimental facility and bridge pier models.
Figure 1. Experimental facility and bridge pier models.
Fluids 06 00370 g001
Figure 2. Schematic layout of the experimental setup (enlarged vertical scale).
Figure 2. Schematic layout of the experimental setup (enlarged vertical scale).
Fluids 06 00370 g002
Figure 3. Geometry of the oblong bridge pier models.
Figure 3. Geometry of the oblong bridge pier models.
Fluids 06 00370 g003
Figure 4. Setup for the flow velocity measurements.
Figure 4. Setup for the flow velocity measurements.
Fluids 06 00370 g004
Figure 5. Schematic diagram of the velocity point measurements. Profiles on the longitudinal half-planes (A and J) are numbered from left to right, while those on the transversal half-planes (E, G, and I) are numbered from the pier outwards.
Figure 5. Schematic diagram of the velocity point measurements. Profiles on the longitudinal half-planes (A and J) are numbered from left to right, while those on the transversal half-planes (E, G, and I) are numbered from the pier outwards.
Fluids 06 00370 g005
Figure 6. Longitudinal velocity profile measured in the flume mid-plane and at the distance 0.50 m upstream the P i e r 11’s front. Experimental data (∘) and theoretical approximation line (dashed line).
Figure 6. Longitudinal velocity profile measured in the flume mid-plane and at the distance 0.50 m upstream the P i e r 11’s front. Experimental data (∘) and theoretical approximation line (dashed line).
Fluids 06 00370 g006
Figure 7. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 1U F b e d . (a) Stream-wise component of velocity (u). (b) Cross-wise component of velocity (v). (c) Vertical component of velocity (w).
Figure 7. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 1U F b e d . (a) Stream-wise component of velocity (u). (b) Cross-wise component of velocity (v). (c) Vertical component of velocity (w).
Fluids 06 00370 g007
Figure 8. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 1D F b e d . (a) Stream-wise component of velocity (u). (b) Vertical component of velocity (w).
Figure 8. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 1D F b e d . (a) Stream-wise component of velocity (u). (b) Vertical component of velocity (w).
Fluids 06 00370 g008
Figure 9. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 2U F b e d . (a) Stream-wise component of velocity (u). (b) Vertical component of velocity (w).
Figure 9. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 2U F b e d . (a) Stream-wise component of velocity (u). (b) Vertical component of velocity (w).
Fluids 06 00370 g009
Figure 10. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 1U E b e d . (a) Stream-wise component of velocity (u). (b) Cross-wise component of velocity (v). (c) Vertical component of velocity (w).
Figure 10. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 1U E b e d . (a) Stream-wise component of velocity (u). (b) Cross-wise component of velocity (v). (c) Vertical component of velocity (w).
Fluids 06 00370 g010
Figure 11. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 1D E b e d . (a) Stream-wise component of velocity (u). (b) Vertical component of velocity (w).
Figure 11. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 1D E b e d . (a) Stream-wise component of velocity (u). (b) Vertical component of velocity (w).
Fluids 06 00370 g011
Figure 12. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 2U E b e d . (a) Stream-wise component of velocity (u). (b) Vertical component of velocity (w).
Figure 12. Vertical profiles of the velocity components along the pier’s mid-plane for E x p . 2U E b e d . (a) Stream-wise component of velocity (u). (b) Vertical component of velocity (w).
Fluids 06 00370 g012
Figure 13. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 1U F b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Figure 13. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 1U F b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Fluids 06 00370 g013
Figure 14. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 1D F b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Figure 14. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 1D F b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Fluids 06 00370 g014
Figure 15. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 2U F b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Figure 15. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 2U F b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Fluids 06 00370 g015
Figure 16. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 1U E b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Figure 16. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 1U E b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Fluids 06 00370 g016
Figure 17. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 1D E b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Figure 17. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 1D E b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Fluids 06 00370 g017
Figure 18. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 2U E b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Figure 18. Vertical profiles of the Reynolds shear stresses along the pier’s mid-plane for E x p . 2U E b e d . (a) Reynolds shear stresses u v ¯ . (b) Reynolds shear stresses u w ¯ . (c) Reynolds shear stresses v w ¯ .
Fluids 06 00370 g018
Table 1. Hydraulic variables of the fixed bed experiments.
Table 1. Hydraulic variables of the fixed bed experiments.
Experiment Pier (-)Q ( m 3 s 1 ) h ( m ) V ( ms 1 ) Fr ( ) Re ( ) Re D ( )
E x p . 1U F b e d P i e r 110.09270.1610 ± 3 × 10 6 0.2882 ± 6 × 10 6 0.1825 ± 8 × 10 6 46,200 ± 0.0631,600 ± 0.7
E x p . 1U E b e d P i e r 110.09250.1628 ± 3 × 10 6 0.2849 ± 4 × 10 6 0.1784 ± 4 × 10 6 46,200 ± 0.0231,200 ± 0.4
E x p . 1D F b e d P i e r 140.09270.1588 ± 5 × 10 6 0.2921 ± 9 × 10 6 0.1875 ± 1 × 10 6 46,200 ± 0.132,000 ± 1.0
E x p . 1D E b e d P i e r 140.09250.1579 ± 3 × 10 6 0.2937 ± 5 × 10 6 0.1896 ± 6 × 10 6 46,200 ± 0.0532,200 ± 0.5
E x p . 2U F b e d P i e r 110.1244 ± 5 × 10 7 0.1952 ± 1 × 10 5 0.3187 ± 2 × 10 5 0.1665 ± 2 × 10 5 61,900 ± 0.234,900 ± 2.0
E x p . 2U E b e d P i e r 110.1237 ± 4 × 10 7 0.1486 ± 5 × 10 7 0.4164 ± 2x10 6 0.2856 ± 2 × 10 6 61,600 ± 0.1845,600 ± 0.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bento, A.M.; Viseu, T.; Pêgo, J.P.; Couto, L. Experimental Characterization of the Flow Field around Oblong Bridge Piers. Fluids 2021, 6, 370. https://doi.org/10.3390/fluids6110370

AMA Style

Bento AM, Viseu T, Pêgo JP, Couto L. Experimental Characterization of the Flow Field around Oblong Bridge Piers. Fluids. 2021; 6(11):370. https://doi.org/10.3390/fluids6110370

Chicago/Turabian Style

Bento, Ana Margarida, Teresa Viseu, João Pedro Pêgo, and Lúcia Couto. 2021. "Experimental Characterization of the Flow Field around Oblong Bridge Piers" Fluids 6, no. 11: 370. https://doi.org/10.3390/fluids6110370

APA Style

Bento, A. M., Viseu, T., Pêgo, J. P., & Couto, L. (2021). Experimental Characterization of the Flow Field around Oblong Bridge Piers. Fluids, 6(11), 370. https://doi.org/10.3390/fluids6110370

Article Metrics

Back to TopTop