1. Introduction
Local scour around bridge piers is resulted from the interaction of local flow field with riverbed sediment, which is one of concerned issues in hydraulic engineering [
1]. The severe scour will undermine the structure and even give rise to the failure of bridge [
2]. It is therefore essential to accurately estimate the scour depth around piers to guarantee the safety of bridges.
In the past several decades, numerous studies on pier scour have been conducted focusing on the equilibrium scour depth, that is the maximum scour depth when it reaches to the equilibrium state [
3,
4]. However, in reality, natural flood duration is usually far shorter than what is needed to reach the equilibrium scour depth [
5]. Hence, in order to predict the scour depth reasonably and accurately from the viewpoints of both the safe and economical design of bridge pier, it is necessary to study and identify the variation of scour depth with time during the pier scouring process.
In engineering practice, simple models (e.g., empirical equations) are normally used to predict the scour process. Compared with large number of equations to predict the equilibrium scour depth, fewer predictive equations are available on the temporal variation of scour depth (see the review [
6,
7]). The widely-used predictor for the temporal scour depth is normally given by an exponential expression [
7,
8]. However, there are two empirical coefficients
c1 and
c2 included in the equations, of which the value is usually ad-hoc and needs to be further clearly determined based on a solid physical explanation. For example, Sumer and Fredsøe (2002) [
7] simply used
c1 =
c2 = 1 for convenience. Simarro and Martin (2004) [
9] noted that the value of
c1 controls the scour rate, that is the curve of scour development is shifted to the left by increasing
c1 [
9]. The value of
c2 is usually considered as a constant [
9], while it can vary under different flow conditions. Cheng et al. (2016) [
8] pointed out that
c1 and
c2 are dependent on sediment coarseness. Ben and Mossa (2020, 2006) [
10,
11] applied this exponential function into predicting the temporal scour evolution downstream of a sloped grade-control structure. They found that compared to the scour evolution at the initial and equilibrium phases, the scour evolution during the intermediate phase under different conditions shows different trend and specifically different value of the two empirical coefficients
c1 and
c2. This indicates that the scour evolution at the intermediate stage is more influenced by the structural and hydraulic conditions. In the current study, the two adjustable coefficients were thoroughly discussed from both theoretical and experimental points of view.
The pier scour is a time-dependent process, and the flow fields around the pier is correspondingly changed with time during the bed scouring process [
12,
13,
14]. The presence of scour hole radically alters the turbulent flow field, and the hydrodynamic characteristics (e.g., bed shear stress and turbulence quantities) are strengthened, which induces further sediment transport and scour development around pier [
15]. For instance, the horseshoe vortex as well as the turbulent energy are enhanced with the development of scour hole [
16,
17], indicating that the flow mechanism effected by the presence of scour hole is different from the flatbed cases without scour hole, which needs to be further explored to deepen the understanding of the pier scouring process. Recently, Ben et al. (2020) [
18] experimentally investigate the flow field within the scour hole downstream of bed sills, and further obtain a new scaling of the equilibrium scour depth based on the phenomenological theory of turbulence. Their results also indicates that the in depth understanding of the flow and turbulence influenced by the scour hole is essential to improve the prediction of scour depth.
Particle image velocimetry (PIV) is a type of non-intrusive velocity measurement technique which is capable of measuring the instantaneous velocity field in the defined flow region. Some scholars have utilized PIV to study the flow field for the pier on the fixed flatbed, e.g., Lin et al. (2003) [
19], Huang et al. (2014, 2018) [
20,
21], Schanderl et al. (2014) [
22], Ozturk et al. (2007) [
23]. Recently, Li et al. (2018) [
24] used the high-resolution PIV system to study the flow characteristics induced by the horseshoe vortex in front of the pier. Similar experimental studies can be found by Chen et al. (2019) [
25] who systematically analyzed the transient evolution process of horseshoe vortex. The previous relevant PIV studies are mainly confined on the flow of a pier mounted on a flatbed without scour hole, which can be regarded as the initial scour state. The limited and valuable study on the flow affected by the pier scour hole can be found from Unger and Hager (2007) [
26] who investigated the temporal evolution of the vertical deflected flow in front of the pier with increasing scour hole by PIV. Overall, the study on the flow field influenced by the development of scour hole at a pier by the high-resolution PIV measurement is rather limited so far.
To our knowledge, a comprehensive understanding of the flow field around piers affected by the scour development is still deficient. The purpose of the present study is twofold: (1) Obtain a simple practical model to predict the temporal scour depth; and (2) study the effect of scour development on turbulent flow field based on a high-resolution PIV system. The paper is organized as follows: The experimental setup and PIV measurement are described in
Section 2. In
Section 3 we discuss the scour development around the pier, involving the scour hole evolution and the temporal variation of scour depth, and a simple practical model to predict the temporal scour depth is developed and validated. In
Section 4 the turbulent flow fields in front of the pier were investigated in terms of the time-averaged velocities, turbulence intensities, Reynolds shear stress and turbulent kinetic energy. Finally, summary and conclusions are drawn in
Section 5. The present experimental results provide not only new insight into the complex flow fields around piers but also valuable experimental data to validate and verify numerical models.
2. Experimental Setup and Methodology
The present experiment was conducted in a water-recirculating and glass-walled flume with 6 m length, 0.25 m width and 0.25 m depth, at Bejing Jiaotong University. The Flow straightener was set at the entrance of the flume to calm the turbulence of flow, and a tailgate was set at the outlet of the flume to adjust the water depth. As depicted in
Figure 1, a circular cylinder with diameter
D = 4 cm was used to model the bridge pier, and mounted vertically on the bed at 4 m downstream of the flume entrance to ensure the fully developed flow occurred [
27]. The sediment bed was set up on the bottom with median diameter
d50 = 0.6 mm and specific mass
ρs = 2650 kg/m
3. The blockage ratio (i.e., flume width to pier diameter) in the present test is 6.25, which is in the proposed acceptable range [
12,
28,
29]. The origin of the coordinate system is defined at the bottom center of the cylinder, with the longitudinal
x-axis parallel to the flume bed, the
y-axis normal to the side wall of the flume, and the
z-axis vertical to the flume bottom.
The experiment was performed under a steady-flow condition, i.e., the incoming discharge from the inlet of the flume is constant. The flow depth
h = 5 cm and the approaching depth-averaged flow velocity
U0 = 0.233 m/s, the Reynolds number based on pier diameter (
ReD) was 9400, and the Froude number based on water depth (
Fr) was 0.28. The present Shields number was calculated as Θ
0 =
u∗2/[(
ρs − ρ)/
ρ]
gd50 = 0.029, where
u∗ (=0.016 m/s) is the bed friction velocity obtained from the best fitting of a log-law profile of the approaching velocity, and
ρ is the fluid density. In this test, the critical Shields number for incipient sediment motion on flat planar bed Θ
cr can be determined by the well-known diagram of Shields curve. Here, for convenience, Θ
cr = 0.033 is given by the method outlined in Soulsby and Whitehouse [
30], and thus the test was performed in a clear-water scour condition (i.e., Θ
0 < Θ
cr).
During the scouring process, the temporal variation of scour depth around the pier has been measured. A scale label was marked on the surface of the pier to facilitate the direct observation of the scour depths. The scour hole morphology was measured by a bed profiler, with space resolution of 2 mm. Since the geometry of the scour hole was observed symmetric about the symmetry plane of the pier, only half of the scour hole was measured.
The two-dimensional instantaneous flow fields in the symmetry plane in front of the pier were measured using a time-resolved PIV system, which consists of a high-speed camera, a pulse laser, tracer particles and post-processing software. It is worth noting that the scour holes that typically develop are, in fact, reasonably symmetric about the symmetry plane, and thus the time-averaged flow field on the symmetry plane can be regarded as two-dimensional in principle. The laser was placed directly above the flume and produced a 1 mm-thick light sheet illuminating the field of view. In order to perform PIV measurements, the flow was seeded with micro-spheres with mean diameters of 10 μm and specific gravity of 1.06. A high-speed digital CMOS camera with a Nikkor AF 50 mm f/1.4D lens was placed on one side of the flume to capture the high-resolution particle images. The camera has a maximum frame rate of 500 fps at full spatial resolution of 1280 × 1024 pixels, and the exposure time between two consecutive frames is 2 ms. The size of the camera viewing area is 52 × 42 mm such that the spatial resolution is 24.5 pixels/mm. Due to the limited memory of camera (2 GB), a sequence of 2000 frames was captured for each case, and we used 2000 frames to calculate the time-averaged flow field. It should be noted that since the cameras were positioned parallel to the test section and due to the presence of the sand bed, the field of view of the cameras was blocked. Hence, the measurements inside the scour hole greater than 1 cm could not be acquired.
The solution of instantaneous velocity vectors between two consecutive frames was achieved by using the iterative multigrid image deformation algorithm [
31]. For the present application, the final interrogation window size was 32 × 32 pixels which have an overlap of 50%, that is the adjacent velocity vector interval is 16 pixels. Based on the flow velocity in our experiment, the interval time of 2 ms assured that the maximum particle displacement between two consecutive frames is 5 pixels, which satisfies the one-quarter rule for PIV correlation analysis [
32,
33]. The accuracy of the velocity vectors is approximately 0.1 pixels, which corresponds to approximately 0.7% of the full-scale velocity in the present measurement. For more specific details on the PIV setup, the readers can refer the previously published paper by the authors [
24,
34]. The validation of the PIV measurement can be seen as the
Appendix A.
4. Turbulent Flow Fields
To identify the physics of the flow disturbed by the generation of the scour hole around the pier, the flow and turbulence characteristics at different scour stages (t* = 2, 11, 133, 266, 931) were analyzed and presented in this section. The flow fields at the symmetry plane in front of the pier were obtained as well as the typical flow and turbulence characteristics are described in detail, in terms of the time-averaged velocities, turbulence intensities, Reynolds shear stress, and turbulent kinetic energy. Moreover, the space-averaged quantities (can be regarded as the representative of the focused area) were obtained to clarify the influence of scour development on the flow and turbulence characteristics.
The instantaneous velocity components were denoted as
u for the streamwise and
v for the vertical direction along the
x- and
z-axes, respectively. Such instantaneous flow velocities could be separated into time-averaged (
U,
V) and fluctuating (
u’,
v’) components according to Reynolds decomposition (i.e.,
u =
U +
u’). The time-averaged velocity is obtained by averaging the instantaneous velocity. The streamwise and vertical turbulence intensities are calculated by
urms (
) and
vrms (
), respectively. The Reynolds shear stress (RSS) is calculated as
, and the turbulent kinetic energy (TKE) is calculated as
[
41].
4.1. Time-Averaged Velocities
Figure 9 shows the velocity vectors superimposed on the velocity contours at different time during the scouring process. The velocity components, including the streamwise and vertical velocities are made non-dimensional with the friction velocity (
U+ = U/u*,
V+ = V/u*). It is worth noting that the results within the scour hole (i.e., below the initial bed) was not available because the sand bed at the lateral side blocks the optical access of the field view inside the scour hole. The vectors clearly illustrate the downward movement of flow (downflow) in front of the pier, and the downflow is more visible with the scour hole development. The streamwise velocity decreases when flow approaching to the pier due to its blocking effect. The contours of
V+ show that the area with negative value of
V+ in front of the pier expands, indicating that the downflow extends with developing scour hole. Moreover, the region with large value of
V+ moves downward along the vertical direction as the scour depth increases.
To further compare the variation of velocity under different scour depth, the space-averaged velocities were obtained and plotted in
Figure 10. The left subplot is the velocity variation along the center line for which the value was obtained by averaging the velocities over vertical direction. To clearly show the influence of scoured bed on the velocity, the corresponding space-averaged results (i.e., by averaging the data in the measured area) versus scour depth were presented in the right subplot. That can be regarded as the double-averaged results obtained by both time and space averaging the instantaneous flow field data. The distribution of
U+ (
Figure 10a) along
x direction exhibits that it is reduced significantly close to the upstream face of the pier. The magnitude of
U+ generally increases with increasing scour depth. 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 is correspondingly increased.
Figure 10b indicates that the value of negative
V+ increases near the pier due to the formation of downflow, and the magnitude of
V+ increases with increasing scour depth, especially for the early scouring stage with high scour rate. Combining the velocity contour analyzed above, it further indicates that with the increase of scour depth the downflow in front of the pier becomes more prominent, and both the extent and strength of the downflow increase with the development of scour hole. This finding is consistent with the results from the acoustic Doppler velocimeter (ADV) measurement by Dey and Raikar (2007) [
12].
4.2. Turbulence Intensities
The contours of turbulence intensities at different time during the pier scouring process are depicted as
Figure 11, in which the streamwise and vertical turbulence intensities (
urms and
vrms) are normalized by the friction velocity.
Figure 11a shows that at the initial scour stage (e.g.,
t* = 2) the region with high value of
urms is near the bed outside of the scour hole. With the development of scour, the region with large
urms moves towards the inside of the scour hole meanwhile the extent with large
urms becomes more significant. This is related to the horseshoe vortex within the scour hole, which was also observed by Dey and Raikar (2007) [
12]. As the scour hole size increases, the strength and size of the horseshoe vortex increase, resulting in the associated turbulence intensity increases correspondingly.
Figure 11b shows that the maximum value of
vrms is located close to the upstream face of the pier, which can be further illustrated by the following
Figure 12.
Figure 12 presents the distributions of turbulence intensities along the center line at different time during the pier scouring process, and the corresponding space-averaged value versus dimensionless pier scour depth is also plotted in the right subplot. It can be seen that the streamwise turbulence intensity is decreased close to the pier, while the vertical turbulence intensity is increased and reached the maximum just near the upstream face of the pier. This variation tendency is reasonably similar to the distribution of flow velocity as shown in the left subplot of
Figure 10. The space-averaged value for both
urms and
vrms increase with increasing scour depth, and
vrms is comparatively smaller than
urms. Overall, the flow in front of the pier becomes more turbulent with the development of scour due to the effect of horseshoe vortex. This implies that the influence of turbulence on sediment transport becomes more prominent, which should be considered in the numerical computation of pier scour [
24].
4.3. Reynolds Shear Stress and Turbulent Kinetic Energy
Figure 13 presents the contours of Reynolds shear stress (RSS) and turbulent kinetic energy (TKE) at different scour stage, in which the value is normalized by the friction velocity. It indicates that the region with large RSS expands with the scour hole development. For instance, the core of the high RSS is located at just the front of the pier at the initial scour stage (
t* = 2). When
t* = 11, the maximum RSS is occurred at the upstream edge of the scour hole (turning point) as a result of flow separation. As the scour evolves upstream, the flow with high RSS expands to the outside of the scour hole (e.g., see the red-color region at
t* = 931), which is probably due to the effect of turbulent horseshoe vortex in the scour hole. The distribution of TKE is similar to that of streamwise turbulence intensity, which is consistent with the results of Dey and Raikar (2007) [
12].
Figure 14 further illustrates the RSS increase with developing scour hole, and the large value of RSS enlarges and moves upstream as shown in the right subplot of
Figure 14a. The variation of TKE with scour depth is relatively complicated rather than monotonic increasing or decreasing. It is inferred that the TKE is closely generated by the horseshoe vortex. With the development of scour depth the horseshoe vortex moves towards the bottom of the scour hole. Hence most of the TKE is taken into the scour hole, and the TKE in the region up the scour hole is correspondingly reduced. However, the size and strength of the horseshoe vortex increase at the same time, and the influence may expand to the flow outside of the scour hole, which leads to the increase of TKE as clarified in the right subplot of
Figure 14b. However, more in-depth investigations are needed to confirm this conjecture and it is beyond the scope of this study.
4.4. Production of Turbulent Kinetic Energy
The TKE budget based on the PIV data obtained in present measurement is analyzed. The production terms included in the TKE budget equation are hereby mainly concerned, including the productions from the streamwise stress
Puu and shear stress
Puv, which are denoted as [
42,
43]
where
u’ and
v’ represent streamwise and vertical fluctuating velocity, respectively;
U is the time-averaged streamwise velocity.
The contours of TKE production at different time during the pier scouring process are depicted as the following
Figure 15. It shows that at the initial scour stage (i.e.,
t* = 2) the region with high absolute value of
Puu is near the bed outside of the scour hole, and the patch of large
Puu seems to move towards the front face of the pier. This is similar to the distribution of
urms, which illustrates that the turbulence intensity term dominates the TKE production in front of the pier due to streamwise fluctuation. With the development of scour, the region with large
Puu moves towards the inside of the scour hole meanwhile the absolute value of
Puu is generally decreased because most of turbulent kinetic energy is transported to the scour hole.
Figure 15b clearly shows that the core of the high
Puv is located at just the front of the pier at the initial scour stage (
t* = 2). In this region, the streamlines experiences large curvature and thus there exists a large shear rate in the streamwise velocity component. The location of its maximum corresponds approximately to the peak position of RSS and TKE. When
t* = 11, the maximum
Puv is formed at the upstream edge of the scour hole (turning point) as a result of flow separation, which is consistent to the distribution of RSS. This patch of shear production is probably a main contributor to the turbulent kinetic energy. With the development of scour, the region with large
Puv moves towards the inside of the scour hole, giving rise to the value of
Puv decreased at this region.
Figure 16 shows the variation of double-averaged TKE budget under different scour depth, in which the left subplot is the velocity variation along the center line for which the value was obtained by averaging the velocities over vertical direction. The corresponding space-averaged results (i.e., by averaging the data in the measured area) versus scour depth were presented in the right subplot.
Figure 16a demonstrates that the absolute value of
Puu is decreased closing to the pier and almost zero near the front face of the pier. This variation tendency is reasonably similar to the distribution of
urms. The absolute space-averaged value for
Puu is decreased with increasing scour depth.
Figure 16b shows that the
Puv is small along the
x direction except for the region close to the pier with a significant change to a large absolute value. With the development of scour depth, the space-averaged value of
Puv is reduced because most of turbulent kinetic energy is transported into the scour hole.
5. Summary and Conclusions
This paper experimentally investigated the influence of scour development on turbulent flow field in front of a bridge pier. The scour development, including the scour hole evolution and the temporal variation of scour depth around the pier were thoroughly discussed. The flow fields in front of the pier at different instants during the scour process were obtained by the high-resolution PIV measurement. Furthermore, the hydrodynamic characteristics, in terms of the time-averaged velocities, turbulence intensities, Reynolds shear stress, turbulent kinetic energy, and production of turbulent kinetic energy were analyzed in detail.
During the scour process, the bed slope upstream of the scour hole at different moment remains almost constant, which is approximately equal to the angle of repose. The scour depth at the pier front exceeds that of the pier side at the later scouring stage. The temporal development of scour depth can be well predicted by a simple practical engineering model based on an exponential function with a change in the two adjustable coefficients (c1 and c2), of which c1 represents the scour rate and c2 is found analytically to be related to the scour hole geometry with a value smaller than 1. Comparison with other previously experimental data demonstrates that the present simple model does a good job of predicting the scour depth observed at both pier front and side.
The flow field indicates that with the scour hole development, the downward flow in front of the pier becomes more prominent, and the flow becomes more turbulent due to the effect of turbulent horseshoe vortex within the scour hole. The variation tendency for both velocities and turbulence intensities along streamwise direction in front the pier shows similarity. The Reynolds shear stress increases with developing scour hole, and the region with large value enlarges and moves upstream of the scour hole. The variation of turbulent kinetic energy with scour depth is relatively complicated rather than monotonic increasing or decreasing. The analysis of TKE budget illustrates that the turbulence intensity term dominates the TKE production in front of the pier due to streamwise fluctuation. With the development of scour, most of TKE production is transported into the scour hole.