Experimental investigation of turbulent flows through a boulder array placed on a permeable bed

Glass beads were used to model permeable beds and boulders (simulated by plastic spherical balls) placed on the centre section of the bed. Flume experiments were conducted to investigate the hydrodynamics through a boulder array over impermeable and permeable beds (i.e. IMPB and PB). For background reference, hydrodynamics investigation was made over smooth beds (SB) with the boulder array. Through measuring the instantaneous velocity field, the major flow characteristics such as mean flow velocity, turbulence intensity, turbulent kinetic energy (TKE) and instantaneous Reynolds stresses (through quadrant analysis) were presented. The results show that the increase in bed permeability through decreasing the exposure height of boulders has little impact on the magnitude of streamwise velocity, but tends to decrease the near-bed velocity gradient, thus affecting the bed shear-stress. For turbulence, similar to the previous studies, the bed permeability is identified to enable a downward shift of the peak of turbulence intensity. The TKE budget analysis shows that bed permeability tends to inhibit the transport and diffusion processes of TKE generation. Finally, the quadrant analysis of turbulence structure clearly shows that the ejections (Q2) and sweeps (Q4) with and without the boulder array are dominated by turbulence structure of different


INTRODUCTION
In mountain rivers, immobile boulders are regarded as large bed obstacles, largely influencing the spatial and temporal variability in flow and topographic dynamics (Euler et

EXPERIMENTS AND METHODS
Experiments were carried out in a rectangular flume equipped with a recirculation system.The flume was 0.5 m wide, 0.6 m deep, and 16 m long with a streamwise bed slope of 0.1%.The flume commenced with a head tank, capably filtering large-scale eddies and providing a smooth flow.The incoming flow was adjusted by an inlet valve to produce various flow configurations and the flow discharge was measured by a sharp-crested weir.The smooth bed (SB) was configured with a layer of glass beads with a constant diameter of d b ¼ 14 mm to mimic the gravel bed configuration, which is commonly observed in mountain river systems.Using glass beads packed in multiple layers to function as a permeable bed can refer to the setup of Manes et al.
(), in which the bead diameter is 12 mm.Because of the one-layer set up, the bed with one-layer of beads is referred to as the impermeable bed (IMPB).To simply construct the permeable configuration, the second layer of beads was well packed on the IMPB, which was able to enhance the vertical motion of flows.Thereafter, it is referred to as the permeable bed (PB) in the current study.points can be recorded so that the turbulent characteristics around a boulder can be well measured.The ADV used was able to measure flow velocities in a range of 0.001 m/s-5 m/s, with an accuracy of ±1%.To improve the data quality, those data with a ratio of signal to noise less than 70% were removed.In addition, high instantaneous velocity spikes were filtered to further refine the quality of the measured velocity data, using the phase-space thresholding technique proposed by Goring & Nikora ().Due to the employment of the one-point-measurement ADV, the measurement efficiency was not as high as with other velocity profilers such as ADV profiler and particle image velocimetry (PIV).For each run, measurements were carried out at every 1 cm across the entire water depth.This In the longitudinal direction, the test section started at a distance of 2 m from the most upstream test section and extended up to 2.205 cm (the ending vertical is located 5.125D downstream of the centre of the boulder).Along the distance of 5.125D, five verticals (0.625D, 1.75D, 2.8D, 4D and 5.125D) were measured.For the high-rate case, after the turbulence intensity reaches its maximum above the boulder, the magnitude decays more rapidly with the vertical distance than that for the PB.This results in a greater deviation between the measured profiles and the original exponential decay function.
The dissipation term (e) is not determined from the velocity measurement since this term involves the secondorder correlation of the fluctuating velocities.An alternative approach to estimating e is by fitting the spectra of the longitudinal velocity in the inertial sub-range according to Kolmogorov's À5/3 theory as below (Baki et al. ): where S( f ) is the power spectra of the streamwise velocity in the frequency domain, f is the frequency, and C is a nondimensional constant with a value of 0.52.Among the above terms, the instantaneous pressure cannot be directly measured by the velocity field, so the pressure transport can be obtained by the residual of the other total terms, yielding:   For the high-rate case, every component to determine the TKE is greater than that for the low-rate case despite the presence of bed permeability.The bed permeability seems to have negative impacts on all terms of the TKE budget.
Specifically, the magnitude of all terms is smaller for the PB than for the IMPB.However, the cause might be uncertain since the bed permeability (increasing bead layers) leading to a decrease in exposure height of the boulders might also decrease turbulent activity near the boulders.
Another interesting observation is that the presence of the bed permeability tends to shift the maximum values of some terms a small distance away from the boulder top.
This is consistent with the distribution of turbulence intensities (Figure 4).

Quadrant analysis of turbulence structure
To analyse the characteristics and contributions of both the burst and sweep events on the turbulence, we used quadrant analysis to examine the conditional statistics of the velocity fluctuations.As suggested by Lu & Willmarth (), a holesize parameter H was defined in the u 0 w 0 plane to determine the dominant instantaneous Reynolds stress, enclosed by four curves ju 0 w 0 j ¼ Hu 0 w 0 .Depending on the selected hole size, this relationship can clearly distinguish between strong and the weak events.The contributions of a particular quadrant to the Reynolds shear-stress is given as: where T is the sampling time and I i,H is the detection function giving: Four events can be described by the above equations, which are the events in the quadrants that were outward interactions, Q1 (i ¼ 1, u 0 > 0, w 0 > 0); ejections, Q2 (i ¼ 2, u 0 < 0, w 0 > 0); inward interactions, Q3 (i ¼ 3, u 0 < 0, w 0 < 0); and sweeps, Q4 (i ¼ 4, u 0 > 0, w 0 < 0).Quadrant analysis was made for H ¼ 0 and 2, respectively.This aims to investigate the contribution of high-frequency small eddies (for H ¼ 0) and low-frequency large vortices (for H ¼ 2) (Balachandar & Bhuiyan ).
Since the Q2 and Q4 represent the dominant events of turbulence and carry the major mass, momentum and energy during the velocity fluctuations, hereby we only present the results of the Q2 and Q4 as shown in For comparison, the results for bed configurations without the boulder array are firstly presented before that for boulder bed.For H ¼ 0, without a boulder array, the near-bed Q4 dominates the entire turbulent activity compared with the Q2 (Figure 7(a) and 7(b)).This indicates that the downward sweeps are the major process in transporting the turbulence momentum and energy from the outer flow to the inner flow.However, the presence of the bed permeability seems to have negligible impact on the magnitude of Q2 and Q4, except for a short downward shifting of the maximum of Q2 and Q4.This identity is somewhat consistent with the distribution of turbulence intensity (see Figure 4).The short shifting might be from the bed porosity providing the extra space so that the Q4 event of turbulence can penetrate the deeper layer.
As the boulder array occurs, the near-bed Q2 behaves comparably to the Q4, suggesting that the turbulence momentum and energy exchange are driven equally by the ejections and sweeps.Even for the SB, the Q2 prevails over the Q4.
Clearly, the presence of the glass beads appears to inhibit the generation of the Q2 and Q4 events.The magnitude of both turbulent events for the SB approximates 5-7 times those for the IMPB and PB.The increase in the glass bead layer also holds a similar but slight inhibiting effect.This suggests that the bed porosity may absorb the turbulent energy of the near-bed energy.However, the decrease in exposure height of the boulders may be the second cause of the weak turbulent activity when the glass bead layers exist.
The presence of the bead layer somehow decreases the roughness of the bed configuration, thus weakening the turbulence.
As some small random fluctuating flow motions are filtered out (by setting H ¼ 2), the Q2 for cases without the boulder array becomes significant, even exceeding the effect of Q4, and the Q4 is highly reduced compared with the values with respect to H ¼ 0. This may be because the ejections (Q2) are dominated by low-frequency large eddies while the sweeps (Q4) are mainly assembled by high-frequency small eddies.With the filtering of small eddies, the large eddies take the major role.However, this change is not applicable to al. ; Afzalimehr et al. ).The investigation of the effects of boulders on flow structures is difficult to perform due to the complexity of the bed surface geometry configuration.Many laboratory studies focused on the mean and turbulent flow characteristics in their vicinity and the complex eddy structures around boulders (e.g.Shamloo et al. ; Sadeque et al. ; Dey et al. ).Regarding an isolated boulder, it is found that three different regions can be identified along the centre line of the boulder, including downstream of the boulder, an intermediate region and upstream of the boulder.For the downstream region of a boulder, the velocity profile appears to be complex over the entire flow depth and a near-wake region forms when flow separates at the edge of the boulder, resulting in downstream flow recirculation.Meanwhile, a far-wake region beyond the near-wake region enables flow to return to a normal state near the bed as upstream.In real situations, boulders in mountain streams commonly occur in arrays.Particularly, multiple arrays of isolated boulders in a regular staggered pattern can be observed in mountain streams (e.g.Lawless & Robert ; Papanicolaou et al. ).In terms of turbulent flow characteristics solely past a boulder array, researchers point out that flows around an individual boulder in an array, different from those around an isolated one, are impacted by flow retardation and backwater effects upstream of the array.Therefore, flow structure around an individual boulder would be largely modified in comparison with the isolated-boulder configuration.Baki et al. (, ) investigated the flow velocity and turbulence characteristics in a staggered-boulder array in a typical fish pass, finding that the boulder array tended to increase the water depth and reduce the flow velocity.Compared with the flow around an isolated boulder, the wake region behind a boulder decayed slowly in the longitudinal direction and had the higher turbulent intensity and turbulent kinetic energy affected by neighbouring boulders.To further understand the effects of a boulder array, Baki et al. () using a three-dimensional numerical model investigated the impacts of variations in boulder size, spacing, and pattern on flow characteristics.Sarkar et al. () used double-averaging methodology (DAM) to successfully resolve the spatial heterogeneity in flows over a boulder array on a gravel bed.Sarkar () found that turbulence characteristics vary insignificantly above the crest level of the array in comparison with those below the crest level.Also, macro-bed-roughness can result in the variability of the curve shape of the velocity profiles (Wang et al. ).Yager et al. () found that boulders in a staggered pattern magnified the reach-averaged total bed shear-stress.Papanicolaou et al. () found that the additional form drag created by the array produced a greater increase in reach-averaged total bed shear-stress compared with the isolated roughness conditions.Fang et al. () investigated flows over submerged boulder arrays and found that boulder arrays promoted the local bed shear-stress distribution.Boulders can promote the formation of discrete clusters, which serve as standalone anchor particles atop the finer particles and on the bed.The fine particles around boulders create a rough, permeable bed.Conventionally, turbulent flows over permeable beds are considered similar to those over impermeable beds since the effect of bed permeability on flows is usually considered negligible.However, when the flow is over a permeable bed, the surface and subsurface flow regions for a permeable bed exist above and within the bed, respectively.The flow over a permeable bed penetrates the bed, driving the exchange of mass, momentum, and energy between the bed surface and subsurface (Zippe & Graf ; Choi & Waller ; Liu et al. ).Bed permeability also controls flow mixing, encouraging a more organized near-wall flow structure and simultaneously governing flow physical processes within porous media (e.g.Cheng et al. ; Manes et al. ).It was observed that the von Karman constant over permeable beds is much lower than its standard value (¼0.41) for terms of permeable smooth and rough beds (Breugem et al. ).Manes et al. (, ) found that bed permeability could enhance the wall-normal fluctuating velocity.This enhancement led to an increase in turbulent shear-stresses and thus further induced higher friction factors over permeable beds.The fact is that numerous studies have reported hydrodynamics of boulder flows, and turbulent flow structure over permeable beds has been extensively investigated.However, the impact of the bed permeability on flows deflected by boulders has rarely been documented.This may lead to a knowledge gap between existing open channel flow theory and realistic flow configurations.Therefore, the present study designs a novel complex porous bed configuration, which describes flow propagating through a boulder array placed on a permeable bed.This study aims to investigate the impacts of the coexistence of boulders and bed permeability on the hydrodynamics.The investigated parameters are mean flow velocity distribution, Reynolds stresses and turbulent kinetic energy (TKE).To deepen the understanding of the interaction mechanism between flow, boulder and bed permeability, we further conduct a quadrant analysis of the turbulence structure, indicating the features of coherent structure over the permeable boulder-array bed.Finally, the budget of TKE is shown to describe the transfer of TKE by different mechanisms.
Figure 1(c).The state of boulders, therefore, can be completely exposed, slightly buried and largely buried.The exposure height, correspondingly, is D 1 ¼ D, DÀd b and DÀ2d b .The flow velocities and turbulence quantities were obtained from the measurement of 3D instantaneous flow velocities using a SonTek ® 10 MHz type Acoustic Doppler Velocimeter (ADV).The sampling frequency was set at measurement resolution is similar to the setup of Nepf & Vinoni (), who used a one-point-measurement ADV to measure flow velocities controlled by near-canopy-top coherent structure.Under this resolution, a good analysis of the TKE budget depending on the spatial derivatives was made.Two flow discharges, i.e.Q ¼ 28 L/s and 48 L/s, were specified and a constant flow depth was maintained by an adjustable tailgate located at the exit of the flume, thus resulting in two Reynolds numbers (Re ¼ 51,638 and 88,523).Therefore, the experiments are referred to as the low-rate and high-rate flows, respectively (Table 1).In experiments, the effective depth h ¼ h o þ z o , the distance from the water surface to the zero-velocity plane, where z o ¼ 0.25d is the displacement position and h o is the distance from the water surface to the top level of the glass beads.As shown in Figure 1(a), the test section has a length of 4 m located at 5 m downstream of the flume inlet to ensure a fully developed flow before entering the boulder array, and this can be confirmed by pre-experiments (examining the water depth).This is also similar to the setup of Papanicolaou et al. () using similar geometry.The velocity measurement was conducted in the region where the flow became developed after moving through the boulder array.

Figure 2
Figure2shows time-averaged streamwise and vertical velocities normalized by the shear velocity (u/u * ) for an IMPB (one-layer beads) and PB (two-layer beads) without a boulder array for the high-rate flow.The presentation firstly provides a basic understanding of the flow behaviour over an IMPB and PB.It is observed that the difference in the streamwise velocity was relatively small (Figure2).For a PB of the high-rate case, the time-averaged velocity

Figure 3
Figure3shows the vertical distribution of streamwise velocity in regards to different bed configurations.Firstly, the time-averaged velocity at different locations downstream of a boulder is shown in Figure3(a)-3(d).Generally compared with the IMPB, the averaged flow velocity above the PB downstream of the boulder (x/D ¼ 0.625-5.125) is greater by roughly 0.4 ∼ 1u * .The increase of the friction factor, however, was 5% and 7% for Q ¼ 48 L/s and Q ¼ 28 L/s, respectively (Table1).Due to the different layers of glass beads laid on the bed, the magnitude of the streamwise velocities can be observed shifted upward as the layer number increases for both low-rate and high-rate flows (Figure3(e) and 3(f)).However, the maximum velocities for the three bed configurations hold no significant difference.This indicates that the exposure height of boulders might have only a small impact on the mean streamwise velocities.It should be noted that as the bead layer number decreases, the variation rate in velocity near the boulders becomes larger.This is attributed to the boulders deflecting the near-bed velocity more as the exposure height increases.Interestingly, the deflection in velocity is more significant for the low-rate case.The possible cause is that for the high-rate case, vortices generated around the boulders might facilitate the exchange of flow streamwise (Figure 4(e) and 4(f)) and vertical (Figure 4(g) and 4(h)) turbulence intensities at x/D ¼ 0.625 are illustrated.

Figure 2 |
Figure 2 | The time-averaged streamwise velocities normalized by the shear velocity (u/u * ) for the IMPB and PB without a boulder array.

Figure 3 |
Figure 3 | Distribution of streamwise velocity.(a) Low-rate flow and (b) high-rate flow for the PB with a boulder array, and for (c) low-rate flow and (d) high-rate flow for the IMPB with a boulder array.Streamwise velocity for different bed configurations for (e) low-rate flow and (f) high-rate flow.Vertical velocity over different bed configurations for (g) low-rate flow and (h) high-rate flow.

Figure 5
Figure 5 depicts S( f ) in the near-bed region (z/h ¼ 0.2) and compares the boulder wake at this point for the highrate flow downstream of the boulder (x/D ¼ 0.625 and 4).

Figure 5 |
Figure 5 | Power spectral density of the streamwise velocity at z/h ¼ 0.2 (a) at different x/D and (b) for different bed configurations.

Figure 6
Figure 6 shows the TKE budget immediately downstream of the boulder (x/D ¼ 0.625) for the IMPB and PB.Similar to the distribution of the turbulence intensity

Figure 6 |
Figure 6 | TKE budget at x/D ¼ 0.625 over (a) the IMPB for the low-rate flow, (b) the IMPB for the high-rate flow, (c) the PB for the low-rate flow and (d) the PB for the high-rate flow.

H.
Figure 7 | Contribution of different quadrant events to the Reynolds stress for flow over different bed configurations for the high-rate flow case: (a-d) for H ¼ 0 and (e-h) for H ¼ 2.

Table 1 |
Hydraulic conditions for the experiments Notes: H ¼ the effective depth; Re is the bulk Reynolds number; bed shear velocity uÃ ¼ ffiffiffiffiffiffiffiffi ffi gHS p , U d ¼ depth-averaged velocity; friction factor, flow resistance coefficient f ¼ (uÃ=U d ) 2 .decreases by ∼0.8u * and bulk flow resistance increases by ∼4% (Table 1).The performance of streamwise velocity and bulk flow resistance for IMPB and PB is similar to the observation of Zippe & Graf () and Manes et al. ().