Temporal and Spatial Flow Variations over a Movable Scour Hole Downstream of a Grade-Control Structure with a PIV System

Weirs or grade-control structures (GCSs) are frequently adopted to protect bridges or control riverbed degradation. Scour holes may develop downstream of these hydraulic structures. Laboratory experiments have been performed in this study, using sophisticated equipment and newly developed procedures. The purpose was to investigate important characteristics of the turbulent flow in the movable scour hole. The results of these experiments demonstrated the significance of instantaneous shear stress in the scouring process. The measured Reynolds stress can be fitted with the theoretical equation reasonably well. Furthermore, the results revealed that the normalized mean vertical velocity profiles in the diffusion region of the scour hole can be fitted with a Gaussian curve. An analysis of the turbulence intensity measurements showed that the turbulent flow is anisotropic in the scour hole. The turbulence intensities also decreased with time as the scour hole gradually approached equilibrium.


Introduction
Many rivers in Taiwan have steep slopes and rapid flows.The flow stage often varies rapidly during typhoon seasons.This can cause severe degradation to the riverbed and levee foundation, and also bridge pile exposures.A grade-control structures (GCS) is frequently adopted to reduce riverbed degradation.For steep rivers, several consecutive low-head GCSs are used.The top of the GCS is flush with the channel bed to reduce possible downstream local scour if required.
Wu and Rajaratnam [15] divided the submerged flow into two broad classes: Impinging jet and surface flow regimes.The latter includes the breaking wave (or surface jump), surface wave, and surface jet.The boundaries of the two regimes were based on the change in direction of the tailwater.On the basis of the experimental data from the literature and also new data (Guan et al. [16]), Guan et al. [17] developed a new transition regime boundary equation as a function of the upstream Froude number and the ratio of the weir height to the tailwater depth.
Lu et al. [1] conducted laboratory experiments with dual GCSs and both steady and unsteady conditions.The first and the second scour holes were under clear-water and live-bed conditions, respectively.It was discovered that three distinct phases can be identified for the evolution of a scour hole downstream of a GCS.These are the initial, developing, and equilibrium phases.A scour evolution simulation scheme was also proposed.The proposed scheme can reasonably well predict the temporal variations of maximum scour depth for unsteady flows with both single and multiple peaks.However, detailed measurements of the turbulent flows were not performed.Guan et al. [18] performed a laboratory study to investigate the flow patterns and turbulence structures in a scour hole downstream of a submerged weir.The study was under conditions of low aspect ratios (less than about 3) and relatively flat slope gradients.The study aims to investigate the temporal and spatial variations of the flow in the scour hole downstream of a GCS for a wide, sloped channel.A sophisticated PIV (particle image velocimetry) system is used to further clarify the scouring phenomenon.The two-dimensional wide channel (high aspect ratio) was selected, as it is closer to the conditions for most of the natural rivers and reduces the flow measuring cost.

Experimental Design and Flume
The experiments were performed in a 17.5 m long, 0.6 m wide, and 0.6 m deep laboratory flume with glass sidewalls (Figure 1).Two horizontal wooden plates with lengths of 11.48 m and 2.85 m were placed upstream and downstream of the sediment recess.A vertical wooden plate was also installed at the upstream end of the sediment recess to represent a GCS.The top of the GCS was flush with the surface of the movable bed.The working length of the flume L which represented the sediment recess was 2.67 m, selected to sufficiently guarantee full scour development without any geometrical interference.In this study, uniform gravel with a median size (D 50 ) of 2.7 mm and a geometric standard deviation (σ g = √ D 84 /D 16 ) of 1.4 was used.In general, ripple is not expected in a laboratory flume with a gravel bed (Gaudio et al. [19] and Ferraro et al. [20]).The clear-water condition (no upstream sediment supply) was considered to slightly simplify the experiments.Table 1 summarizes the experimental conditions.The water level measurements were conducted by using an ultrasonic transducer Keyence, Osaka, Japan, with a resolution ± 1 mm.hole downstream of a GCS.These are the initial, developing, and equilibrium phases.A scour evolution simulation scheme was also proposed.The proposed scheme can reasonably well predict the temporal variations of maximum scour depth for unsteady flows with both single and multiple peaks.However, detailed measurements of the turbulent flows were not performed.Guan et al. [18] performed a laboratory study to investigate the flow patterns and turbulence structures in a scour hole downstream of a submerged weir.The study was under conditions of low aspect ratios (less than about 3) and relatively flat slope gradients.The study aims to investigate the temporal and spatial variations of the flow in the scour hole downstream of a GCS for a wide, sloped channel.A sophisticated PIV (particle image velocimetry) system is used to further clarify the scouring phenomenon.The two-dimensional wide channel (high aspect ratio) was selected, as it is closer to the conditions for most of the natural rivers and reduces the flow measuring cost.

Experimental Design and Flume
The experiments were performed in a 17.5 m long, 0.6 m wide, and 0.6 m deep laboratory flume with glass sidewalls (Figure 1).Two horizontal wooden plates with lengths of 11.48 m and 2.85 m were placed upstream and downstream of the sediment recess.A vertical wooden plate was also installed at the upstream end of the sediment recess to represent a GCS.The top of the GCS was flush with the surface of the movable bed.The working length of the flume L which represented the sediment recess was 2.67 m, selected to sufficiently guarantee full scour development without any geometrical interference.In this study, uniform gravel with a median size ( D 50 ) of 2.7 mm and a geometric standard deviation ( g σ = D D 84 16 ) of 1.4 was used.In general, ripple is not expected in a laboratory flume with a gravel bed (Gaudio et al. [19] and Ferraro et al. [20]).The clear-water condition (no upstream sediment supply) was considered to slightly simplify the experiments.Table 1 summarizes the experimental conditions.The water level measurements were conducted by using an ultrasonic transducer Keyence, Osaka, Japan,with a resolution ± 1 mm.The experiments were conducted with two independent variables: Unit flow discharge q and bed slope S b .The flow discharges were selected on the basis of Froude law, with a length ratio (L r ) of 1/50 for the 5-year and 20-year flows of the Tou-Chien River, Hsin-Chu County, Taiwan.The slope range of 1-1.5% was typical for gravel bed rivers in Taiwan.The average approach flow depth h 0 was measured at x 0 = 8.48 m, where the flow was fully developed.The low relative submergences (h 0 /D 50 , ranged from 6 to 9) were mostly due to the gradients of steep slopes.Powell [21] reviewed the studies of Bathurst et al. [22] and Lawrence [23] on the effect of relative submergence.Figure 2 shows the evolution of maximum scour depth (d m ) for a preliminary test (Run S-2).The result indicated a very high scouring rate for the first hour, and the scouring phenomenon almost reached equilibrium after 5 h.The turbidity of the flow was extremely high for the first 15 min.Thus, tests were conducted at Water 2018, 10, 1002 3 of 14 15 min, 1 h, and 5 h.The experimental duration of 5 h in this study represents the peak flow duration for an unsteady flood event.Note: S b = bed slope; h 0 = average approach flow depth; q = unit flow discharge; U 0 = average approach flow velocity; U * = average approach flow shear velocity; F r = Froude number; R e = Reynolds number; B/h 0 = aspect ratio; d m,e = maximum scour depth at the equilibrium stage; x m,e = position of maximum scour depth; l m,e = maximum scour hole length at the equilibrium stage.
Water 2018, 10, x FOR PEER REVIEW 3 of 14 Reynolds number; B h 0 = aspect ratio; d m,e = maximum scour depth at the equilibrium stage; x m,e = position of maximum scour depth; l m,e = maximum scour hole length at the equilibrium stage.
The experiments were conducted with two independent variables: Unit flow discharge q and bed slope b S .The flow discharges were selected on the basis of Froude law, with a length ratio ( r L ) of 1/50 for the 5-year and 20-year flows of the Tou-Chien River, Hsin-Chu County, Taiwan.The slope range of 1-1.5% was typical for gravel bed rivers in Taiwan.The average approach flow depth h 0 was measured at 0 x = 8.48 m, where the flow was fully developed.The low relative submergences ( h D 0 50 , ranged from 6 to 9) were mostly due to the gradients of steep slopes.Powell [21] reviewed the studies of Bathurst et al. [22] and Lawrence [23] on the effect of relative submergence.Figure 2 shows the evolution of maximum scour depth ( d m ) for a preliminary test (Run S-2).The result indicated a very high scouring rate for the first hour, and the scouring phenomenon almost reached equilibrium after 5 h.The turbidity of the flow was extremely high for the first 15 min.Thus, tests were conducted at 15 min, 1 h, and 5 h.The experimental duration of 5 h in this study represents the peak flow duration for an unsteady flood event.) and the aspect ratio ranged from 24-38.The aspect ratios ( B h 0 ) were greater than 10, making it reasonable to assume that the resistance of the side walls was negligible and the flows were two-dimensional (2-D).As a result, a two-dimensional measuring system was adopted in the current study.

PIV Measuring System
The images analyzed were captured using a high-resolution digital camera (Phantom Miro ex4) with a resolution of  800 600 pixels and highest frequency array 1265 Hz.In the conventional PIV measurement, lasers are predominantly adopted in PIV setups.The conventional methods are difficult to apply in our movable bed experiments due to the scattering of particles.
To reduce the particle scattering and produce reasonably uniform light, a sheet of soft light paper was attached on the back sidewall and halogen lamps were adopted in this study.Black Acrylonitrile-Butadiene-Styrene (ABS) particles were used as the seeding particles.The particle size ranged from 0.25-0.3mm and the density was 1.04 g/cm 3 .The image velocimetry was then adopted to analyze the moving trajectories of the black particles projected on the soft light paper.

Verification of the PIV Measuring System
A preliminary test was performed in the fully developed zone of the channel to compare the results of the PIV measuring system with those of Nezu and Rodi's [24] empirical formula, which The Froude number F r ranged from 2.29-2.63(supercritical flows, F r = U 0 / gh 0 ) and the aspect ratio ranged from 24-38.The aspect ratios (B/h 0 ) were greater than 10, making it reasonable to assume that the resistance of the side walls was negligible and the flows were two-dimensional (2-D).As a result, a two-dimensional measuring system was adopted in the current study.

PIV Measuring System
The images analyzed were captured using a high-resolution digital camera (Phantom Miro ex4) with a resolution of 800 × 600 pixels and highest frequency array 1265 Hz.In the conventional PIV measurement, lasers are predominantly adopted in PIV setups.The conventional methods are difficult to apply in our movable bed experiments due to the scattering of particles.
To reduce the particle scattering and produce reasonably uniform light, a sheet of soft light paper was attached on the back sidewall and halogen lamps were adopted in this study.Black Acrylonitrile-Butadiene-Styrene (ABS) particles were used as the seeding particles.The particle size ranged from 0.25-0.3mm and the density was 1.04 g/cm 3 .The image velocimetry was then adopted to analyze the moving trajectories of the black particles projected on the soft light paper.

Verification of the PIV Measuring System
A preliminary test was performed in the fully developed zone of the channel to compare the results of the PIV measuring system with those of Nezu and Rodi's [24] empirical formula, which had been verified with 2-D fiber-optic laser doppler velocimetry (Lu et al. [25]).Figure 3 shows three measured dimensionless longitudinal turbulence intensity profiles in the fully developed verification region (center located at x 0 = 8.48 m) with Nezu and Rodi's [24] empirical curves, in which Y + = U * y i ρ/µ and y i is the elevation from the fixed channel bed, ρ is the water density, µ is the dynamic viscosity of water, and u = u 2 (u = longitudinal fluctuating velocity).The grid size for the analysis of PIV data was 7 mm.The results indicate that generally the measured profiles are in good agreement with the curves presented by Nezu and Rodi [24], except for the region very close to the water surface.In this region, light may have been slightly affected by the outdoor environment.As our primary focus was on the flow field near the channel boundary (e.g., boundary shear stress), the error near the water surface was ignored.
Water 2018, 10, x FOR PEER REVIEW 4 of 14 had been verified with 2-D fiber-optic laser doppler velocimetry (Lu et al. [25]).Figure 3 shows three measured dimensionless longitudinal turbulence intensity profiles in the fully developed verification region (center located at 0 x = 8.48 m) with Nezu and Rodi's [24] empirical curves, in which and i y is the elevation from the fixed channel bed, ρ is the water density, μ is the dynamic viscosity of water, and   = u 2 u (  u = longitudinal fluctuating velocity).The grid size for the analysis of PIV data was 7 mm.The results indicate that generally the measured profiles are in good agreement with the curves presented by Nezu and Rodi [24], except for the region very close to the water surface.In this region, light may have been slightly affected by the outdoor environment.As our primary focus was on the flow field near the channel boundary (e.g., boundary shear stress), the error near the water surface was ignored.

Two-Dimensional Velocity Distribution
Analysis of the mean velocity profiles for the conditions investigated (Table 1) revealed that the flow in the scour hole can be qualitatively divided into three regions shown in Figure 4a: Jet and diffusion, transition, and acceleration regions (Run S-2, b S = 0.015, q = 0.0283 m 2 /s, r F = 2.59).In Figure 4, U is the vectorial mean velocity.For the purpose of demonstration, Run S-2 was selected as a representative case as it has the highest flow intensity (high q and b S values) and the largest scour hole.
As the flow passes the submerged GCS, the supercritical high-velocity fluid behaves like a submerged jet, oscillating occasionally.When the jet deflection angle is large (the early stage of the scouring process), the impinging jet (Wu and Rajaratnam [15]) may occur.This may entrain a significant amount of bed sediment, causing a "bursting phenomenon" (Lu et al. [1]).A surface jump is usually formed slightly downstream.Severe mixing of turbulence and energy dissipation occur in the jump.Point b is the highest point of the surface jump, while point a is the lowest point of the water surface profile before the surface jump.As depicted in Figure 5, the normalized flow velocity profiles can be approximated as a normal distribution near zone ab, indicating the diffusion characteristic of the submerged jet.As a result, the first region from the entrance (GCS) to point b is defined as the jet and diffusion region.A recirculation zone is also formed when the scour hole is deep enough.The reattachment point downstream of the recirculation zone usually continues to move slowly downstream with time.

Two-Dimensional Velocity Distribution
Analysis of the mean velocity profiles for the conditions investigated (Table 1) revealed that the flow in the scour hole can be qualitatively divided into three regions shown in Figure 4a: Jet and diffusion, transition, and acceleration regions (Run S-2, S b = 0.015, q = 0.0283 m 2 /s, F r = 2.59).In Figure 4, U is the vectorial mean velocity.For the purpose of demonstration, Run S-2 was selected as a representative case as it has the highest flow intensity (high q and S b values) and the largest scour hole.
As the flow passes the submerged GCS, the supercritical high-velocity fluid behaves like a submerged jet, oscillating occasionally.When the jet deflection angle is large (the early stage of the scouring process), the impinging jet (Wu and Rajaratnam [15]) may occur.This may entrain a significant amount of bed sediment, causing a "bursting phenomenon" (Lu et al. [1]).A surface jump is usually formed slightly downstream.Severe mixing of turbulence and energy dissipation occur in the jump.Point b is the highest point of the surface jump, while point a is the lowest point of the water surface profile before the surface jump.As depicted in Figure 5, the normalized flow velocity profiles can be approximated as a normal distribution near zone ab, indicating the diffusion characteristic of the submerged jet.As a result, the first region from the entrance (GCS) to point b is defined as the jet and diffusion region.A recirculation zone is also formed when the scour hole is deep enough.The reattachment point downstream of the recirculation zone usually continues to move slowly downstream with time.The normalized velocity profiles can be approximated as a normal distribution in zone ab above the recirculation zone:   The normalized velocity profiles can be approximated as a normal distribution in zone ab above the recirculation zone:  The normalized velocity profiles can be approximated as a normal distribution in zone ab above the recirculation zone: where u, h, u max and h max are the mean longitudinal velocity, vertical distance measured above the recirculation zone, maximum u value and h value at u max , respectively.The upper envelope curve of the recirculation zone in Figure 4 was obtained by interpolation.Based on the regression analysis, the standard deviation σu was 0.370.The result indicates that although the maximum velocity of the submerged jet decreases as the flow expands downstream, the similarity criterion can be applied to zone ab.The last region within the scour hole is called acceleration region, occurring when the channel bed elevation increases rapidly near the downstream end of the scour hole.The flow accelerates due to the decrease of the flow depth, and gradually approaches the incoming flow condition.The region between the jet and diffusion region and the acceleration region is called transition region.This region elongates with time, and the elevation of the channel bottom flattens out (or only slightly increases) downstream, as shown in Figure 4c.
The maximum scour depth d m,e , position of maximum scour depth x m,e and maximum scour hole length l m,e at the equilibrium stage (5 h) are also summarized in Table 1. Figure 4c indicates that after the vertical evolution of the maximum scour depth (refer to Figure 2), the scour hole tends to gradually evolve in the downstream direction, revealing that the length of the sediment recess (L in Figure 1) is long enough.

Turbulence Intensities
Figures 6 and 7 show the dimensionless longitudinal and vertical turbulence intensity contours along the channel centerline, in which u = u 2 (u = longitudinal fluctuating velocity), and v = v 2 (v = vertical fluctuating velocity), and U * = approach flow shear velocity.The measuring frequency of the PIV system was 1/500 s.As shown in Table 1, all the Froude numbers of the approach flows are greater than 2 (supercritical flows).Figures 6 and 7 reveal that high turbulence intensity values occur near the entrance region and zone ab (defined previously, see Figure 4a), where the surface jump occurs.
In general, the dimensionless longitudinal turbulence intensity is higher than the dimensionless vertical turbulence intensity, indicating the turbulent flow is anisotropic.The distributions of u /U * and v /U * are very similar, as there were two local maxima occurring at the entrance region and the surface jump.Regarding the jet and diffusion region (x ≈ 0.15-0.3m), only the values near the scour hole bottom decrease for T in the range of 15 min to 1 h.For T in the range of 1 to 5 h, u /U * decreases by approximately 25%, while v /U * decreases by approximately 20%; and u /U * ∼ = 1.5v /U * near the scour hole bottom.By contrast, u /U * decreases by approximately 50% near the water surface, while v /U * only decreases slightly, especially for the area near the surface jump.A knowledge of turbulent flow characteristics within the scour hole may help to improve the accuracy of the scour prediction models.For instance, Dodaro et al. [26] and Dodaro et al. [27] employed instantaneous flow velocity data to accurately predict scour evolution downstream of a rigid bed.    Figure 8a shows that the dimensionless mean absolute Reynolds shear stress t * τ U 2 was very low ( < 0.2) along the scour hole boundary, except for the values near the entrance region (  x 0.15 m).This is due to the effect of the control point at GCS ( x = 0).The GCS prevents the erosion of the channel bed and the reduction of the flow intensity at the entrance point.However, the corresponding

Reynolds Shear Stress
Figure 8 shows the variations in the dimensionless mean absolute Reynolds shear stress |τ t |/U * 2 and the exceeding probability P exc along the boundary of the scour hole with time.The exceeding probability is defined as where N |τ t |>τ c is the number of |τ t | values greater than the critical shear stress τ c , and N total is the total number of the instantaneous Reynolds shear stress τ t values collected during the sampling period (8 s).The Reynolds shear stress (τ t = −ρu v ) was obtained based on the Shields diagram, and τ c = 2.1 N/cm 2 for D 50 = 2.7 mm.With consideration of particle size roughness in the movable boundary (D 50 = 2.7 mm), a grid size of 7 mm was adopted in the analysis of PIV data.The lowest point selected for estimating the bed shear stress was therefore set at 3.5 mm above the bed.
Water 2018, 10, 1002 9 of 14 exc P values were greater than 0.17 at T = 15 min, indicating the significance of the instantaneous shear stress, especially at the early stage of the scouring process.Figure 8b,c shows that the exc P value decreases with time as the maximum scour depth gradually approaches equilibrium.In fact, it also quantified Shen and Lu's [28] concept that the turbulent shear stress moves the sediment particles.Additionally, Figure 8a reveals that the exc P value increases downstream in the acceleration region (0.4 < x < 0.55 m), possibly causing scour holes to evolve downstream.Table 2 summarizes the temporal variations in the first four moments of the longitudinal and vertical velocities at the lowest point of the scour hole.The mean velocities and their standard deviations are in m/s.The third moment (skewness) and fourth moment (kurtosis) of the longitudinal velocity were calculated using the following: Figure 8a shows that the dimensionless mean absolute Reynolds shear stress |τ t |/U * 2 was very low (<0.2) along the scour hole boundary, except for the values near the entrance region (x ≤ 0.15 m).This is due to the effect of the control point at GCS (x = 0).The GCS prevents the erosion of the channel bed and the reduction of the flow intensity at the entrance point.However, the corresponding P exc values were greater than 0.17 at T = 15 min, indicating the significance of the instantaneous shear stress, especially at the early stage of the scouring process.Figure 8b,c shows that the P exc value decreases with time as the maximum scour depth gradually approaches equilibrium.In fact, it also quantified Shen and Lu's [28] concept that the turbulent shear stress moves the sediment particles.Additionally, Figure 8a reveals that the P exc value increases downstream in the acceleration region (0.4 < x < 0.55 m), possibly causing scour holes to evolve downstream.
Table 2 summarizes the temporal variations in the first four moments of the longitudinal and vertical velocities at the lowest point of the scour hole.The mean velocities and their standard deviations are in m/s.The third moment (skewness) and fourth moment (kurtosis) of the longitudinal velocity were calculated using the following: Water 2018, 10, 1002 10 of 14 where n represents the number of samples.Similarly, the skewness and kurtosis of the vertical velocity, SK v and KT v were calculated.Table 2 reveals that most of the |SK u | and |SK v | values were less than 1, and most of the |KT u | and |KT v | values were less than 2. As a rule for sizes greater than 300, a continuous distribution is close to a normal distribution if the absolute values of the skewness and excess kurtosis are less than 2 and 4, respectively (Kim [29]).
It is reasonable to assume that both longitudinal and vertical velocities are normally distributed since the sample size is 4410 in the PIV data analysis of this study.Most of the correlation coefficients r in Table 2 are less than 0.2, except that of Run S-1 at T = 15 min (r = −0.40).It is also reasonable to assume that both longitudinal and vertical velocity components are independent.Based on the video recorded by the charge coupled device (CCD) camera, the value of r was relatively high for Run S-1 at T = 15 min due to the rugged channel bottom.If random variables U and V are independently and normally distributed, with U : N 0, σ u 2 , V : N 0, σ v 2 , and let W = UV, Lu et al. [30] theoretically derived the corresponding probability density function of W to be: where σu and σv are the standard deviations for variables U and V.
In this study, U and V are the random variables corresponding to the longitudinal and vertical fluctuating velocities u and v .By proper integration with MATLAB, the probability density function of the normalized Reynolds shear stress −u v /U * 2 can be obtained.
Figure 9 shows the comparisons of the measured longitudinal and vertical flow velocities, and the fitted normal distributions for Run S-2 at T = 15 min.Figure 10 shows the comparisons of the measured probability distributions of −u v /U * 2 and fitted theoretical distributions for Run S-2 at T = 15 min, 1 h, and 5 h.As can be seen in the figure, the theoretical distributions fitted the measured histograms reasonably well.Table 2 reveals that in general σu decreases with time for all the experimental conditions, and σv decreases with time for Run S-1 and S-2.As a result, the range of the dimensionless Reynolds shear stress −u v /U * 2 also decreases with time in Figure 10, causing a reduction in the extreme values.Table 3 summarizes the temporal variations in the mean (τ t ) and standard deviation ( στ t ) of Reynolds stress and the exceeding probability (P exc , defined in Equation ( 2)) at the lowest point of the scour hole.Theoretically, the critical shear stress (τ c ) for a sloping bed needs to be properly adjusted to reflect the gravitational effect (Chiew and Parker [31]).However, in this study we are mainly interested in the τ c value at the lowest point of the scour hole.It is therefore reasonable to assume that the τ c value for a plane bed can be adopted.As expected, the στ t value decreases with time for all the experimental conditions.The P exc value also decreases with time as the maximum scour depth gradually approaches the dynamic equilibrium condition.
In general, Table 3 reveals that the exceeding probability P exc increases with the flow intensity (bed slope S b and unit flow discharge q).On the basis of the regression analysis, the following equation is developed to provide a qualitative trend of the physical process: Figure 11 shows a comparison of the measured and predicted P exc values [(P exc ) P vs. (P exc ) M ], with a determination coefficient r 2 = 0.94.Two dashed lines with discrepancy ratios (DR) of 0.5 and 2.0 are plotted in the figure.The discrepancy ratio is defined as the ratio of the predicted P exc value to the corresponding measured value.In general, Equation ( 6) fitted the data reasonably well.For Run S-1 and Run M-2 at T = 1 h, the deviations are relatively large due to the relatively rugged bed surfaces.
Water 2018, 10, x FOR PEER REVIEW 12 of 14 Theoretically, the critical shear stress ( c τ ) for a sloping bed needs to be properly adjusted to reflect the gravitational effect (Chiew and Parker [31]).However, in this study we are mainly interested in the c τ value at the lowest point of the scour hole.It is therefore reasonable to assume that the c τ value for a plane bed can be adopted.As expected, the ˆt τ σ value decreases with time for all the experimental conditions.The exc P value also decreases with time as the maximum scour depth gradually approaches the dynamic equilibrium condition.In general, Table 3 reveals that the exceeding probability exc P increases with the flow intensity (bed slope b S and unit flow discharge q ).On the basis of the regression analysis, the following equation is developed to provide a qualitative trend of the physical process:

Conclusions
Based on the laboratory investigation of the nearly two-dimensional turbulent flows over a movable scour hole downstream of a GCS with a PIV system, the following conclusions can be drawn: 1.A PIV system with a high-resolution digital camera, halogen lamps, and ABS seeding particles was adopted and employed to successfully measure the temporal and spatial variations in the flow over the movable scour hole.2. According to the measured vectorial mean velocity profiles (Figure 4), the flow field in the scour hole can be classified into three regions-jet and diffusion, transition, and acceleration regions.Additionally, a recirculation zone was noted near the channel bottom.The reattachment point of the recirculation zone usually moves slowly downstream with time.3. The vertical mean velocity profiles in the diffusion region (between points a and b in Figure 4) of the scour hole can be normalized and fitted with a normal curve (Figure 5, Equation ( 1)). 4. In general, the dimensionless longitudinal turbulence intensity is higher than the dimensionless vertical turbulence intensity, indicating the turbulent flow is anisotropic in the jet and diffusion

Conclusions
on the laboratory investigation of the nearly two-dimensional turbulent flows over a movable scour hole downstream of a GCS with a PIV system, the following conclusions can be drawn: 1.
A PIV system with a high-resolution digital camera, halogen lamps, and ABS seeding particles was adopted and employed to successfully measure the temporal and spatial variations in the flow over the movable scour hole.

2.
According to the measured vectorial mean velocity profiles (Figure 4), the flow field in the scour hole can be classified into three regions-jet and diffusion, transition, and acceleration regions.Additionally, a recirculation zone was noted near the channel bottom.The reattachment point of the recirculation zone usually moves slowly downstream with time.

3.
The vertical mean velocity profiles in the diffusion region (between points a and b in Figure 4) of the scour hole can be normalized and fitted with a normal curve (Figure 5, Equation ( 1)).

4.
In general, the dimensionless longitudinal turbulence intensity is higher than the dimensionless vertical turbulence intensity, indicating the turbulent flow is anisotropic in the jet and diffusion Water 2018, 10, 1002 13 of 14 region.In this study, the distributions of u /U * and v /U * were similar because there were two local maxima occurring at the entrance and the surface jump.In general, u /U * and v /U * values decreased with time as the scour hole gradually approached equilibrium.5.
By assuming that the longitudinal mean velocity ( u ) and vertical mean velocity ( v ) are independently and normally distributed, a theoretical Reynolds stress (τ t = −ρu v ) distribution can be derived (Equation ( 5)).The measured Reynolds stress can be fitted with the theoretical equation reasonably well (Figure 10).6.
This study demonstrated the significance of the instantaneous shear stress, especially at the early stage of the scouring process, which also quantified Shen and Lu's [28] findings.Furthermore, the experimental results show that the exceeding probability increased with the unit flow discharge and bed slope, and decreased gradually with time.

Figure 1 .
Figure 1.Plane view of the experimental channel (not to scale).

SL = 2 Figure 1 .
Figure 1.Plane view of the experimental channel (not to scale).

Figure 3 .
Figure 3.Comparison of the measured dimensionless longitudinal turbulence intensity profiles in the fully developed zone with the empirical profiles presented by Nezu and Rodi [24] at: (a) 0 x = 8.41

Figure 3 .
Figure 3.Comparison of the measured dimensionless longitudinal turbulence intensity profiles in the fully developed zone with the empirical profiles presented by Nezu and Rodi [24] at: (a) x 0 = 8.41 m, (b) x 0 = 8.48 m, and (c) x 0 = 8.55 m.

Figure 5 .
Figure 5.Comparison of the normalized mean velocity profiles in zone ab (Figure 4).

Figure 5 .
Figure 5.Comparison of the normalized mean velocity profiles in zone ab (Figure 4).

Figure 5 .
Figure 5.Comparison of the normalized mean velocity profiles in zone ab (Figure 4).

Figure 6 .
Figure 6.Measured dimensionless longitudinal turbulence intensity (  U * u ) contours in the scour

Figure 7 .
Figure 7. Measured dimensionless vertical turbulence intensity (  U * v ) contours in the scour hole

Figure 8
Figure8shows the variations in the dimensionless mean absolute Reynolds shear stress t

Figure 8 .
Figure 8. Variations in the dimensionless Reynolds shear stress ( t

Figure 8 .
Figure 8. Variations in the dimensionless Reynolds shear stress (|τ t |/U * 2 ) and the exceeding probability (P exc ) along the boundary of the scour hole at (a) T = 15 min, (b) T = 1 h, and (c) T = 5 h.

Figure 9 .
Figure 9. Comparisons of the measured flow velocities and fitted normal distributions for Run S-2 ( T = 15 min): (a) u and (b) v .

Figure 9 . 14 Figure 9 .
Figure 9. Comparisons of the measured flow velocities and fitted normal distributions for Run S-2 (T = 15 min): (a) u and (b) v.

Figure 11
Figure 11 shows a comparison of the measured and predicted exc P values [   p exc P

Figure 11 .
Figure 11.Comparison of the measured and predicted exceeding probabilities (DR = Predicted exc P

Figure 11 .
Figure 11.Comparison of the measured and predicted exceeding probabilities (DR = Predicted P exc /Measured P exc ).

Table 1 .
Summary of the experimental conditions.

Table 1 .
Summary of the experimental conditions.

Table 2 .
Summary of the temporal variations in the first four moments and the correlation coefficients of the longitudinal and vertical velocities at the lowest point of the scour hole.

Table 3 .
Summary of the temporal variations in the Reynolds stress and the exceeding probability ( exc P ) at the lowest point of the scour hole.

Table 3 .
Summary of the temporal variations in the Reynolds stress and the exceeding probability ( exc P ) at the lowest point of the scour hole.

Table 3 .
Summary of the temporal variations in the Reynolds stress and the exceeding probability (P exc ) at the lowest point of the scour hole.