Next Article in Journal
The Use of Phytoplankton in the Assessment of Water Quality in the Lower Section of Poland’s Largest River
Next Article in Special Issue
Local Scour around Tandem Double Piers under an Ice Cover
Previous Article in Journal
Combined Effects of Hillslope-Concentrated Flows and Riverine Stream Waves on Soil Erosion in the Reservoir Riparian Zone
Previous Article in Special Issue
Investigation of the Effect of Vegetation on Flow Structures and Turbulence Anisotropy around Semi-Elliptical Abutment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characteristics of Turbulence in the Downstream Region of a Vegetation Patch

1
Faculty of Civil Engineering, Iran University of Science and Technology, Tehran 16846-13114, Iran
2
School of Engineering, University of Northern British Columbia, Prince George, BC V2N 4Z9, Canada
*
Author to whom correspondence should be addressed.
Water 2021, 13(23), 3468; https://doi.org/10.3390/w13233468
Submission received: 18 November 2021 / Revised: 29 November 2021 / Accepted: 1 December 2021 / Published: 6 December 2021
(This article belongs to the Special Issue Fluvial Hydraulics Affected by River Ice and Hydraulic Structures)

Abstract

:
In presence of vegetation patches in a channel bed, different flow–morphology interactions in the river will result. The investigation of the nature and intensity of these structures is a crucial part of the research works of river engineering. In this experimental study, the characteristics of turbulence in the non-developed region downstream of a vegetation patch suffering from a gradual fade have been investigated. The changes in turbulent structure were tracked in sequential patterns by reducing the patch size. The model vegetation was selected carefully to simulate the aquatic vegetation patches in natural rivers. Velocity profile, TKE (Turbulent Kinetic Energy), turbulent power spectra and quadrant analysis have been used to investigate the behavior and intensity of the turbulent structures. The results of the velocity profile and TKE indicate that there are three different flow layers in the region downstream of the vegetation patch, including the wake layer, mixing layer and shear layer. When the vegetation patch is wide enough (Dv/Dc > 0.5, termed as the patch width ratio, where Dv is the width of a vegetation patch and Dc is the width of the channel), highly intermittent anisotropic turbulent events appear in the mixing layer at the depth of z/Hv = 0.7~1.1 and distance of x/Hv = 8~12 (where x is streamwise distance from the patch edge, z is vertical distance from channel bed and Hv is the height of a vegetation patch). The results of quadrant analysis show that these structures are associated with the dominance of the outward interactions (Q1). Moreover, these structures accompany large coherent Reynolds shear stresses, anomalies in streamwise velocity, increases in the standard deviation of TKE and increases in intermittent Turbulent Kinetic Energy (TKEi). The intensity and extents of these structures fade with the decrease in the size of a vegetation patch. On the other hand, as the size of the vegetation patch decreases, von Karman vortexes appear in the wake layer and form the dominant flow structures in the downstream region of a vegetation patch.

1. Introduction

Submerged vegetation patches are fundamental participants of aquatic ecosystems, providing important habitats for fauna, including fishes [1,2,3]. As the reduction in the aquatic vegetation is a global concern, researchers have conducted studies to investigate the concept of a large-scale environmental issue [4,5,6]. However, while most of the current hydraulic-orientated research is focused on the flow characteristics in the canopy region of fully vegetated channels, finite dimensional vegetation patches with no fully developed region downstream of vegetation patches are very common in natural rivers. These sparse vegetation patches have a substantial effect on flow structures depending on the dimensions of the vegetation patch [7,8,9]. The effect of vegetation patch dimensions on flow structures has significant consequences for the ecosystem through various mechanisms. For example, an intermediate coverage (10–40%) of vegetation may promote high species richness due to the conditions provided by the optimal behavior of flow [10,11,12]. In addition, flow structures around a vegetation patch have a significant impact on the sediment transport, and consequently on the topography of the river bed [13,14,15]. Overall, the interaction between flow and vegetation in natural channels has attracted a lot of attention from hydraulic researchers. To mention only a few, the most recent research works have been conducted by researchers [16,17,18], in which various aspects of this interaction have been investigated.
Given the importance of the research topic, this present study aims to assess the effects of different sizes of vegetation patch on the flow characteristics in the downstream region of vegetation patches. The sizes of vegetation patches range from full coverage across the entire channel to a narrow vegetation patch located in the middle of the channel.

2. Flow Structures behind Vegetation Patches

Compared to a rough bed without a vegetation patch, the presence of a vegetation patch in a channel results in changes in the flow structures and momentum transfer. Therefore, flow resistance and turbulence characteristics are expected to differ in the downstream regions of vegetation patches [19]. To assess the characteristics of turbulent flow, the following methodologies are commonly used: spectral analysis of velocity time series, determination of Reynolds stresses, and Turbulent Kinetic Energy (TKE), as well as quadrant analysis [20,21,22]. The analysis of the coherent Reynolds shear stress is another method that is used to determine coherent turbulent structures. However, this approach has hardly been used in hydraulic-oriented works due to the difficulties in calculating the phase velocity.
In order to make clear the terminologies used in the above-mentioned approaches, an introduction is presented in this section alongside the literature review. In the classic terminology of turbulences, for a sufficiently long duration of velocimetry (T), and for the velocity component of ui(t), the average velocity ( u ¯ i ) is defined as:
u ¯ i = lim T 1 T 0 T u i ( t ) d t
Considering i = 1, 2, 3 as the x, y and z directions in the Cartesian coordinate system, u ¯ ,   v ¯ and w ¯ are defined as the averages of velocities in the x, y, and z directions, respectively. Consequently, the velocity fluctuations are defined by:
u = u ( t ) u ¯ v = v ( t ) v ¯ w = w ( t ) w ¯
Then, the Turbulent Kinetic Energy (TKE) is also calculated as:
TKE = 0.5 ( u 2 ¯ + v 2 ¯ + w 2 ¯ )
In general, as was reported by researchers, a sharp increase in TKE   is linked to the presence of vortexes [23,24,25,26]. In partially covered channels with emergent vegetation patches, where the flow is almost 2D, the flow structure in the downstream region of vegetation patches is dominated by the horizontal von Karman vortexes developed by instabilities in the trailing edges of the vegetation patch [27,28]. The distance between the trailing edge of the barrier and the formation of von Karman vortexes is known as Lkv. According to reported research based on TKE analysis, for a relatively dense vegetation patch with a diameter of Dv, in which the leakage velocity from the dense vegetation patch is negligible, Lkv ≈ 2.5Dv [27,28]. On the other hand, the 3D flow structures in the downstream region of a submerged vegetation patch are more complicated. In addition to the constant presence of a vertical recirculation zone beyond the vegetation patch, the presence and magnitude of the von Karman vortexes in the horizontal plane depends on the geometry of the vegetation patch [28,29,30]. For instance, the presence of von Karman vortexes at a depth of Hv/H > 0.55~0.7, where Hv is the height of the vegetation patch and H is the flow depth [28,30]. On the other hand, results indicated that the dominant structure of a wide submerged path occurs in the vertical plane, as vortexes forming a vertical recirculation zone [28,29]. The distance between the trailing edge of the barrier and the forming of the vertical vortex is known as Lkv. For a solid submerged barrier, there is no gap between the trailing edge of the barrier and the forming vortex, and the length of wake behind the barrier is zero (Lvv = 0) [31]. However, for a porous vegetation patch, the velocity of flow passing through the vegetation patch may delay the formation of the recirculating region (Lvv ≠ 0). The reported value of Lvv ranges between 1Hv and 5Hv depending on the porosity and diameter of the vegetation patch [28,29]. As reported by Liu et al. (2018), the value tends to be about 1Hv in a relatively wide vegetation patch, and increases up to 5Hv in vegetation patches with a low blockage ratio [28]. These results indicate that Lvv and Lkv are affected by the vegetation patch dimensions of h and D, respectively. In addition, the bleed velocity that penetrates into the wake region can also increase these values considerably [23,31].
In addition to the above-mentioned classic variables, coherent turbulent variables have become fundamental terms of flow analysis in recent years [32,33,34]. According to the triple decomposition approach, the instantaneous velocity can be written as [35]:
u ( x , t ) = u ¯ ( x ) + u ˜ c ( x , t ) + u r ( x , t )
where u ( x , t ) is the instantaneous velocity, ur is the deflection of velocity (known as u’ in the classic approach), and u ¯ ( x ) is the classic average of velocity during time period T. In addition to the classic averaging method, the periodic phase average can be defined as:
u ( x , t ) = lim N 1 T P 1 N u i ( x , t + i T P )
where Tp is the period of occurrence of a coherent structure and is equal to 1/fd, where fd is the dominant frequency of the occurrence. Generally, when a coherent event occurs in a particular region of flow, fd will reach an obvious peak in the power density spectrum of the related time series in the same region. For a time series with a length of T, N is the number of cycles with a period of Tp, which can be calculated as N = T/Tp. Following these definitions, the coherent velocity deflection is:
u ˜ c ( x , t ) = u ( x , t ) u ¯ ( x )
Considering the same approach for other components of velocity, the coherent Reynolds shear stress can be written as:
τ c = ρ u ˜ c w ˜ c
The total coherent and non-coherent Reynolds shear stress is defined as follows:
τ r = ρ u ¯ r w ¯ r ¯
Although it is very rare to implement τ c in studies of river hydraulics with the presence of vegetation patches, a few research works have been conducted to study the variation in τ r in river hydraulics in the presence of vegetation patches. For instance, in the downstream region of a submerged vegetation patch, there is a high gradient pattern of Reynolds shear stress in the vertical plane, which increases toward the surface of the flow [36]. This phenomenon is induced by the high-speed flow that passes above the submerged vegetation patch and causes a considerable velocity gradient when compared with the low velocity flow leaking through the vegetation patch. The region with a high velocity gradient is known as the shear layer. This region of flow has particular turbulence characteristics that will be studied in this paper.
As another method of analysis, the decomposition of bursting events (known as quadrant analysis in 2D space and octant analysis in 3D space) is also widely applied to determine the dominant turbulent events in the presence of both emerged and submerged vegetation patches [22,37,38,39]. The occurrence probability of event k is calculated as the normalized occurrence frequency, fk, for a particular class of events related to different classes of events [40]:
f k = n k N with N = 1 4 n k   k = 1 , 2 , 3 ,   4
where   n k is the number of events belonging to class k, and N is the total number of events.
For the quadrant analysis, a “hole” region has been proposed in the majority of previous studies. In a “hole” region, the event must be filtered and not be considered. In octant analysis, this threshold can be written as:
| u ( t ) w ( t ) | > C H | u w ¯ |
For the threshold parameter CH, the low-intensity events below a certain limit were omitted, which were scaled by the average of the velocity fluctuations. The high value of CH implies the selection of the strongest events, but the total number of instantaneous [u’(t) w’(t)] decreases so much that the contribution region in each quadrant becomes meaningless. For quadrant analysis, the threshold level CH = 1 is suggested for use to reach a good compromise between the clear identification of the events and the preservation of a number of instantaneous events of each class [40,41]. In this paper, CH is set as 1. The classes of events are also defined as:
Q1—occurrence of the outward interaction when u’ > 0 and w’ > 0;
Q2—occurrence of the ejection when u’ < 0 and w’ > 0;
Q3—occurrence of the inward interaction when u’ < 0 and w’ < 0;
Q4—occurrence of the sweep when u’ > 0 and w’ < 0.
To date, most of the reported research has applied quadrant analysis to study coherent flow structures within or above a vegetation patch. However, there are a few published studies in which quadrant analysis has been conducted to assess the coherent flow structures in the downstream region of a vegetation patch. Okamoto and Nezu (2013) point out that in the region immediately behind the submerged vegetation patch, there is a small margin of ejection dominancy around the top edge [7]. Similarly, based on data collected from experiments performed on dryland vegetation in a wind tunnel, Mayaud et al. (2016) revealed that there were elevated frequencies of Q2 (ejection) and Q4 (sweep) events in the immediate toe of the vegetation patch [42]. In contrast, it is reported the dominance of outward and inward interactions in the shear layer induced by the flow passing above the vegetation patch, which is a notable characteristic of the flow in the downstream region of a submerged vegetation patch in a more distant region [38]. It must be mentioned that the relation between quadrant occurrences and velocity structures is an interesting research topic that has received attention from researchers recently. Wang et al. (2019) established interconnections between the classic definition of vortex groups and quadrant occurrences, which are used in the next section of this paper [43].
In conclusion, the introduced techniques have been successfully used in the hydrodynamics analysis in channels with the presence of vegetation patches. However, the reported studies on the hydrodynamics of channels with submerged vegetation patches are limited to those with the flow either above the patches or immediately downstream of the patches. For example, the flow structures was traced up to a distance of 8Hv from the toe of the patch towards the downstream [36]. Liu et al. (2018) studied the flow structures to a distance of just 5~6Dv from the patch toe in the downstream direction [28]. The quadrant analysis has been conducted in the field; here, the distance from the patch toe was even more limited, to ~3Dv [38]. Although this distance from the patch toe is enough to address the characteristics of patch-induced turbulence in the near field, it is not enough to cover the more distant events of turbulence. For example, the investigation of the characteristics of the shear layer generated above the patch requires an extended range of velocimetry in the downstream direction. In addition, there is a need to combine the above-mentioned methods to evaluate the turbulent characteristics of flow behind the submerged vegetation patches. This requires a particular emphasis on the implementation of τ c in the turbulence analysis, which has hardly been studied in recent works. In order to fill the above-mentioned gaps in the previous research, in the present study, a combined approach to turbulent analysis is used to investigate the turbulent characteristics of flow in the region downstream of a vegetation patch with an extended range up to 17Hv.

3. Experimental Setup

The experiments were conducted using a glass flume at the Laboratory of Hydraulics of the Iran University of Science and Technology. The flume is 14 m long, 0.9 m wide and 0.6 m deep. The discharge was controlled by an electromagnetic flow meter installed at the entrance of the flume and was set for 31 L/s in the present study. The water level in the flume was adjusted by a tailgate located at the end of the flume and was set for a depth of 18.5 ± 0.3 cm in the present study. The distance between the flume entrance and vegetated flume zone was 6 m to ensure a fully developed flow in the region upstream of the vegetation patch.
The velocity measurements were conducted after the flow reached the steady state condition. The velocity profiles were measured using an acoustic Doppler velocimeter (ADV), placed at the centerline of each row of vegetation patch. There were 19–26 measuring points along each vertical line for velocity measurements, and the vertical distance between two adjacent measuring points was 4–10 mm. The sampling frequency and the measuring time of the ADV were 200 Hz and 120 s, respectively, resulting in 24,000 instantaneous velocities for point measurement.
To date, most previous studies have employed artificial simple rods of regular shapes to simulate natural vegetation patches. As real vegetation is flexible and irregular, this method might not represent the nature of vegetation’s behavior [36]. On the other hand, the experiments showed that natural vegetation would lose its stiffness and develop a long-lasting curvature in the flow direction after couple of days. Therefore, a well-shaped synthetic plant was used to represent the natural vegetation, which was selected on the basis of a real world sample of patches in a gravel bed river.
Each model plant had three branches. Each branch had 12 leaves, and the diameter of the branch trunk was approximately 3 mm, as shown in Figure 1. The average height of the vegetation patch was 105 ± 5 mm, and the lateral and longitudinal spread widths of the leaves were approximately 9–19 and 11–22 mm, respectively. The model plants had a certain degree of flexibility and could swing in a flowing current similar to vegetation in a natural river. The vegetation patch was attached to a perforated board in a staggered arrangement. Four different layouts of vegetation patches were used to simulate both fully covered and non-fully covered channel beds. The dimensions of these vegetation patches, namely, length (Lv) × width (Dv), in this study were 120 × 90 cm, 90 × 60 cm, 60 × 45 cm and 40 × 30 cm, respectively. Table 1 summarizes the experimental runs. The material on the flume bed was a mixture of natural gravel similar to that in a natural gravel bed river (the Marbor River, Zagros Mountains region, Iran). The equivalent particle diameter, which 90% of the total particles were smaller than (d90), was 18.8 mm, and d50 = 14 mm. The layout of the experimental device and employed materials are shown in Figure 1.

4. Results and Discussion

The results are presented as follows: First, both the velocity profile and the TKE have been investigated to determine the flow layers that formed behind the vegetation patch. Then, spectral analysis and coherent Reynolds shear stress analysis have been carried out to determine the nature of turbulence and coherent Reynolds shear stresses in different flow layers. Consequently, the results of the quadrant analysis of bursting events are provided to compare and evaluate the coherent occurrences in the framework of the time domain. Then, by considering the particular characteristics of turbulence in the mixing layer, the temporal characteristics of turbulence in the mixing layer are investigated. Finally, by considering the results of previous sections, the transformation of coherent structures in the downstream region of a vegetation patch is provided.

4.1. Velocity Profile and TKE

For all experimental runs, three layers of flow can be observed based on the inflection points of the velocity profiles (Figure 2). The associated regions of these layers are known as the wake zone, the mixing layer and the Log–Law shear zone. These three regions resemble the results of trough patch velocimetry reported by other researchers [44,45,46,47]. For instance, in case 1 and for a distance from the vegetation edge x/Hv = 8, the wake zone formed at the water depth of z/Hv = 0~0.7 (z/H = 0~0.56), the mixing layer formed at the water depth of z/Hv = 0.7~1.1 (z/H = 0.56~0.62) and the Log–Law shear zone formed at the water depth of z/Hv > 1.1 (z/H > 0.62). The extent and thickness of the mixing layer reduced as the vegetation patch decayed. In this layer, a severe deflection is detectable in the velocity profile, which resembles an adverse pressure gradient. For case 1, this anomalous phenomenon covered a distance (from the vegetation edge) between 8 < x/Hv and x/Hv ≥ 17. In case 2, this phenomenon covered a distance of 3 < x/Hv < 8, and it was only visible around x/Hv = 12 in case 3. However, this phenomenon was absent in case 4. The wake zone was formed in the lower region beyond the patch. For case 1, the wake zone occurred at a distance from the vegetation edge of x/Hv = 1, with the upper limit of this layer at z/H = 0.45 or z/Hv = 0.8, which increased up to z/Hv = 1 for the smallest patch (case 4). In this layer, the velocity profiles tended to adopt a Log–Law shape beyond the distance of x/Hv = 12 for all cases. For a channel bed partially covered by a vegetation patch, the thickness of the wake layer varied along the flow direction, and increased slightly in such a way that it reached the lower boundary of the Log–Law shear region. Consequently, the mixing layer became narrower and even disappeared in the presence of smaller patches. However, for the case of a channel bed fully covered with a vegetation patch (case 1), the thickness of the mixing layer increased considerably at a distance (from the vegetation edge) of x/Hv = 8, and seemed to be effectively present at a distance of x/Hv >> 17.
While the vertical vortexes that normally generate behind a barrier have attracted a lot of attention from researchers, the development of horizontal vortexes (the von Karman Vortex Street) needs to be investigated too. Generally, the formation of von Karman vortexes is associated with the occurrence of peaks in the TKE values [23]. Thus, for each case, the TKE values have been calculated for the relative flow depths of z/Hv = 0.5 and z/Hv = 1, and also for the near-bed region behind the vegetation patch along the centerline of the flume (Figure 3). The variation in the streamwise velocity along the flow direction is also shown in this figure. In case 1, there was only one peak in the graph of TKE/(U0)2 around x/Dv = 0.9 and x/Hv = 7.5~8, which is associated with severe deflection in the velocity profile. However, no flow escaped from the sides of the vegetation patch (since the channel bed was fully covered with vegetation). This peak in TKE/(U0)2 clearly indicates the presence of a vertical vortex behind the vegetation patch and near the water’s surface. While a sharp increase in TKE occured after a distance of x/Hv = 4.5~5, the length of the wake behind the patch can be considered as Lvv = 5Hv. Importantly, the anomalies in the mixing layer appeared in the same point as was reported [28] for a length of wake of Lvv = (3.5~5)Hv with a wide submerged vegetation patch. However, Folkard (2005) reported a slightly lower range for a wide but highly submerged vegetation patch [29]. Overall, these results are comparable to those of Liu et al. (2018) for a wide and partially vegetation-covered bed with slight submergence [28].
For experiment cases 2, 3 and 4, the location of the minimum velocity was shifted towards the flow direction. However, since the patch width ratio of Dv/Dc was relatively high for case 2, Lvv was only shifted to the near-bed region. For cases 3 and 4, the shift of wake was also detectable at depths of z/Hv = 0.5 and z/Hv = 1. In these cases, the shift extent of the upper layers exceeded that of the near-bed zone. For all four cases, the minimum velocity was negative in the near-bed region, implying the occurrence of a vertical vortex. Considering the velocity values in the upper zone, it was observed that the center of the vortex (or rotation) moved upward if Dv/Dc decreased. In case 1, the center of vortex was located in the near-bed region and z/Hv = 0.5, while it was located between z/Hv = 0.5 and z/Hv = 1 in case 4 (smallest vegetation patch). This elevated center of rotation in case 4 pushed the larger part of flow down into the lower zones. Subsequently, the velocity pattern in the middle region became similar to that of the near-bed region (one can observe this by comparing the minimum velocity located between the near-bed region and z/Hv = 0.5). In addition, the larger radius of the rotation in the channel with a smaller patch width ratio of Dv/Dc caused the minimum streamwise velocity to be in the near-bed region rather than in the mixing layer. The reason for the upward movement of the rotation center in the case of a smaller patch in the channel may be the higher through-patch velocity. Regarding the larger values of Lvv due to the through-patch flow, this statement is consistent with the findings of other researchers [27,28,31]. In addition, for experiment cases 2, 3 and 4, the flow from the sides of the vegetation patches created the second peak in the TKE curve at the depth of z/Hv = 1. However, for cases 2 and 3, while the Dv/Dc values were still relatively high, the second peak in the TKE curve did not occur in the middle and near-bed region. Additionally, the location of the single peak in the TKE curve in the lower zones was shifted forward.
Overall, the results of both the TKE and velocity analyses show that the three layers (wake, mixing and shear) were present in all cases of vegetation patches. However, the wake and shear layers were considerably affected by the flow passing from the sides of the vegetation patches, and as a consequence, different flow structures and associated length scales emerged, dependent on the upward movement of the center of the vertical rotation beyond the patch.

4.2. Spectral Analysis and Coherent Reynolds Shear Stress

The dominant frequency of the streamwise Reynolds shear stress ( f d u w ) was the peak of the power spectrum diagram calculated via u′w′ fluctuations (Suw). Therefore, Suw is calculated at different points for all experimental cases (Figure 4), and the results are used to calculate the coherent Reynolds shear stresses ( u ˜ c w ˜ c ). The results of the spectral analysis can be placed in the following three categories:
(1)
The first category is the points where the turbulence was well matched with the von Karman isotropic–homogeneous turbulence, but no dominant frequency on Suw was observable. At these points, the power spectra reached the slope of −5/3 in the inertial range. As there was no dominant frequency, f d u w was too low and no coherent events occurred in the streamwise Reynolds shear stresses (the blank area in Figure 5 represents these points);
(2)
The second category is the points where the turbulence was well matched with the von Karman isotropic–homogeneous turbulence. The dominant frequency on Suw was observable in the form of an obvious peak around 0.07~0.3 Hz. At these points, the power spectra reached the slope of −5/3 in the inertial range. Where there was an obvious dominant frequency, f d u w was used to calculate the streamwise coherent Reynolds shear stresses. Excluding points at the depth of z/Hv ≈ 1, all other points for the coherent Reynolds shear stress belong to this group;
(3)
The third category is the points where the turbulence was not compatible with the von Karman isotropic–homogeneous turbulence. The dominant frequency on Suw was observable in the form of an obvious peak around 0.07~2 Hz. At these points, the power spectra reached the slope of −1 in the inertial range. According to the analysis provided by Tchen (1953), this slope is associated with the anisotropic turbulence characterized by a large vorticity and strong resonance [47]. These points were observed at a depth of z/Hv ≈ 1 (mixing layer) for cases 1, 2 and 3, while they were absent in case 4. This finding confirms the upward movement of the rotation center in the presence of smaller patches in the channel (case 4), because the downward velocity beyond the patch was too weak to penetrate into the lower zones and the center of rotation developed in the upper zone. Such a weak vortex cannot produce the strong vertical vortexes required to alter the isotropic turbulence. A Matlab@ code was used to calculate the phase velocity and its deflection for the dominant frequency at each point. The coherent vertical momentum transfer triggered by coherent vertical vortexes continued to appear in both the near-canopy and near-bed regions in cases 1 to 3. However, near-bed coherent shear stress did not occur in case 4. It can be inferred that the flow regime in the near bed region of case 4 was dominated by the flow through the vegetation patch and the wall effect of the rough bed. The presence of the vegetation patch had no coherent effect on this region. In addition, for cases 1, 2 and 3, the majority of coherent occurrences were accompanied by a high spatial gradient of the Reynolds shear stress, which is associated with the boundary layer separation zone and large scale vortexes—a characteristic that was described by Lian (1990) [48]. However, in case 1, the occurrence of coherent Reynolds shear stresses was associated with very low Reynolds shear stresses with a low spatial gradient. Consequently, it can be referred that the coherent Reynolds shear stresses can be classified into two categories according to the origins of the coherence. The first and most prevalent category is the strong coherent Reynolds stress, which was prevalent in the mixing layers of cases 1, 2 and 3. The second category is the weak coherent Reynolds shear stress that occurred in the wake layer. Based on the results of the spectra analysis, the strong coherent shear stresses were associated with peaks of the anisotropic turbulent spectra. In contrast, the weak form of the coherent shear stress occurred in the presence of the peak isotropic turbulence.
On the other hand, for the patch width ratio of Dv/Dc = 0.5, a new type of coherent structure formed in the wake zone, and extended as the patch size reduced. With the patch width ratio of Dv/Dc = 0.33 (case 4), the greatest part of wake zone was covered by this structure. While, in this region, the turbulent spectra of instantaneous velocity resembled isotropic turbulence, the effect of the interaction between the mixing layer and the shear layer seems to be insignificant. On the other hand, the presence of strong flow from the patch sides could form strong von Karman vortexes in the horizontal plane. As reported by Siddique et al. (2008) and Liu et al. (2018), von Karman vortexes appeared in partially covered patches with Hv/H > 0.55~0.7 [28,30]. Consequently, for these regions, the frequency of vortex shedding was calculated as 0.07~0.11 Hz, which is equivalent to an average Strouhal number of about 0.25. These results are comparable with those reported by other researchers [27,28].

4.3. Quadrant Analysis of Bursting Events

The occurrence probabilities of the quadrant classes are shown in Figure 6. One can see from Figure 6 that in the presence of smaller patches, the ejection-dominated zone was shifted to a higher elevation at a shorter distance from the patch. This phenomenon was accompanied with lower Reynolds shear stress in the near-bed region. It was also associated with a marginally thicker wake layer in the presence of smaller patches, which was described in the previous section.
In contrast, the sweep-dominated zone, observed above the top of the patch in case 1, faded with the decrease in the patch dimensions. As regards the Reynolds shear stresses of this region, the occurrences of the sweep were associated with the high Reynolds shear stress of the Log–Law layer in the downstream region of the patches. The same pattern was reported in the downstream region of dryland vegetation patches facing wind flow [42].
As strong coherent shear stresses were prevalent in the boundary between the mixing layer and the shear Log–Law layer, the quadrant analysis of the bursting events of this zone are outlined in Figure 7. The propinquity between the outward-dominated regions and regions with coherent Reynolds shear stresses is another notable result of the quadrant analysis. For all vegetation patches, the coherent shear stress in the mixing layer with a frequency of f d u w appeared in exactly the same region in which the outward class of quadrant occurrences were dominant, where fk=1 > 0.25. This was a one-way relationship, and there were some points with fk=1 > 0.25 at which no coherent shear stress was observable. To assess this important outcome, for point x and quadrant class   i , two conditional functions can be defined, as below:
A ( x , i ) = { 1   i f   f k = i = max { f k = 1 .. 4 } 0   i f   f k = i max { f k = 1 .. 4 }
B ( x ) = { 1   i f   f d ( u w ) 0 0   i f   f k = i = 0
Subsequently, ρ A , B   (the Pearson correlation coefficient between A ( x , i ) and ( x ) ) can be used to determine which class of quadrant occurrences is well-matched with the coherent Reynolds shear stress. For a sample including 278 points for all four cases of patch setup, ρ A , B is 0.68 for Q1 (outward interaction). In comparison, this value is −0.45, 0.26 and −0.38 for Q2, Q3 and Q4, respectively. This result confirms that the outward interactions are the representative quadrant class of the coherent Reynolds shear stresses in the downstream region of vegetation patches, particularly for strong coherent occurrences. Figure 8 illustrates the values of f k and f d u w for the 278 points beyond the vegetation patches that were investigated in this research. According to the definition of the outward quadrant, the upward momentum flux within the mixing layer toward the shear layer must be considered (see Figure 7). This upward flux is also detectable in the velocity profiles (Figure 2), where anomalies in the mixing layer are exacerbated; the streamwise velocity was reduced in the mixing layer and increased in the bottom part of the shear layer. This phenomenon illustrates a particular behavior in the velocity time series, which will be discussed in the next section.

4.4. Temporal Characteristics of Turbulence in the Mixing Layer

As was described in the previous sections, the mixing layer beyond the vegetation patch is prone to anisotropic turbulence, characterized by outward interaction events and the existence of coherent Reynolds stresses. The deviation from an ideal isotropic turbulence originated from the turbulent intermittency. The turbulent intermittency was characterized by a preference for turbulence with large velocity gradients, which is reflected in the strongly non-Gaussian tails of the probability density functions of velocity differences. These tails are determined by extreme events [49]. Intermittency is also defined as an abnormality in turbulent flow initiated by the interaction between turbulent regions or the interaction between a turbulent region and a vicinal laminar behavior region. The second mechanism is the common source of intermittency in shear flows. Moreover, the occurrence of a highly intermittent region in the near-boundary regions was associated with the prevalence of coherent vortices [50]. The mechanism generates quasi-laminar intervals in the velocity time series of a turbulent flow, which postpones the non-viscous intermediate sub-range in the energy cascade. Thus, there was a slight slope in the energy spectra. Consequently, these intervals persist up to arbitrarily large Reynolds numbers through the energy cascade sequences, which produce extreme events in velocity time series [51]. The determination of a particular intermittency index is a challenging issue and requires the spatial analysis of the velocity time series measured at different locations at the same time. Clearly, there is a need for appropriate sensors, such as PIV and hot wires.
However, temporal analyses of velocity time series can also provide a proper description of the quantity and intensity of intermittent occurrences. As the intermittency is linked with the non-Gaussian tails of the probability density functions of velocity differences, the velocity time series for the mixing layer were subjected to a normality test. The traditional empirical rule of normality (or a simple normality test) states that, in a normally distributed data set, 68%, 95%, and 99.7% of the values lie within one, two, and three standard deviations (σ) of the mean (μ), respectively. Most of the points in the wake area could barely pass these restrictions; however, the normality test of Kolmogorov–Smirnov identified almost all points as non-normally distributed data sets. Regarding the methodology and application of the Kolmogorov–Smirnov normality test, please refer to Massey (1951) [52]. The results of the simple and Kolmogorov–Smirnov normality tests for case 1 are shown in Table 2. As the Kolmogorov–Smirnov normality test is very sensitive to non-Gaussian deviations, it even identifies semi-normal velocity time series for the wake layer as non-normal distributions. However, according to the concept of the Kolmogorov–Smirnov normality test, the intensity of non-normality can be observed on the basis of the deviation from the normal cumulative distribution. In this way, an error index can be defined as:
E r r o r C D F = 1 C D F E C D F N
where CDFE refers to the cumulative distribution of the empirical velocity time series, and CDFN is the cumulative distribution of a normal data set with the same μ and σ as the empirical data set. Figure 9 shows the variation in ErrorCDF for different depths at a distance of x/Hv = 8, where the mixing layer is well developed. One can see that the values of ErrorCDF are considerably higher at a depth of 0.8 < z/Hv < 1. However, the most notable issue is that the max values reduce as the vegetation patch disappears. Comparing the results of a simple normality test to ErrorCDF, it can be concluded that a threshold around ErrorCDF = 0.2 can be assumed for the lower limit of intermittency. This threshold is also confirmed by the results of the intermittent turbulence and coherent Reynolds shear stresses at a depth of 0.8 < z/Hv < 1. Correspondingly, in case 4, where ErrorCDF < 0.2 at all depths, no anisotropy was observed in the turbulent spectra. In contrast, the anisotropic turbulence was limited to the mixing layer in cases 1, 2 and 3. In the smaller z/Hv, the no intermittency is observed for all cases and the effect of bed roughness is the determining factor in the near bed region. Generally, in the near bed region of gravel bed channels, the effect of bed roughness exceeds the patch-produced vortexes [53,54,55].
According to the results of this research, the intermittency and anisotropy of the turbulence are associated with the non-Gaussian distribution of instantaneous velocity. This relationship was also identified by researchers, and is known as a fundamental characteristic of the anisotropic intermittent turbulence [49,50,51]. However, in the present research, this characteristic was detected in the mixing layer beyond vegetation patches with a relatively high blockage ratio (Dv/Dc > 0.5).
In addition, the formation of intermittent turbulence in the downstream region of the patch has been investigated in this research. One can see from Table 1, by comparing to a normally distributed data set, that the non-Gaussian data sets are dependent on a greater concentration of data between μ − σ and μ + σ. In addition, the values of the normalized standard deviation are higher in these velocity time series. Consequently, the variation in the probability density functions (PDF) of the instantaneous velocity was considered along the flow direction. The results have been compared with the turbulent spectra of the corresponding points. For the mixing layer of case 1, these results are shown in Figure 10. As can be seen in Figure 10, the intermittent anisotropic turbulence formed at a depth of 8 < x/Hv < 12, which is the same range as the occurrence of the outward quadrant dominancy, leading to strong coherent Reynolds shear stresses and anomalies in streamwise velocity. This event can also be inferred from the flattened shape of PDF, which is associated with high standard deviations.
Furthermore, the intensity of the intermittency can also be evaluated in the energy domain. Researchers used an analysis method based on the filtering of TKE time series via a certain threshold [56,57]. These thresholds include MTKE − aσTKE and M + aσTKE, where MTKE is the mean, “a” is a constant multiplayer (commonly 1 and 3) and σTKE is the standard deviation of the TKE time series in 2D space. The values that exceed this range are considered as strong occurrences characterizing intermittent events. Although the ratio of filtered samples describes the quantity of intermittent events, it cannot provide enough information about the intensity of the intermittency. It is reported that the TKE of extreme events is a suitable indicator of the intermittency [58]. Since extreme events are associated with coherent occurrences in the mixing layer, the Turbulent Kinetic Energy of these extreme events (TKEi) can be used for this propose. This approach is used in this paper for calculating TKE and TKEi in a 3D space. Overall, Figure 11 summarizes the results regarding the formation of intermittent turbulence in the mixing layer behind the vegetation patch by means of demonstrating the filtered time series of the instantaneous velocity and the variation in the normalized standard deviation of TKE and normalized TKEi.

4.5. Transformation of Coherent Structures beyond a Vegetation Patch

Considering the results of the previous sections, the transformation of coherent structures beyond a vegetation patch can be summarized as below:
(1)
In the downstream region of a fully channel-spanning vegetation patch, the coherent structures are observable just behind the patch. These structures originate from the stem-scale vortexes that are formed by the leaking flow passing through the vegetation patch [24,28]. As the patch width ratio of Dv/Dc reduces to 0.66, the leaking flow increases. Consequently, the stem-scale coherent structures are observed in a larger area (Figure 5). With the patch width ratio of Dv/Dc = 0.5, these structures spread to the near-bed region and create a large area of coherent shear stress, in which a Q1-dominant core is surrounded with a sweep ejection-dominant region. However, with a patch width ratio of Dv/Dc = 0.33, the area of stem-scale coherent structures reduces suddenly, and is limited to x/Dv < 3 in the top and middle zones of the wake layer behind the patch;
(2)
In the mixing layer of a fully channel-spanning vegetation patch, coherent shear stresses are associated with largely intermittent fluctuations in the instantaneous velocity forming anisotropic turbulence at 8Dv < x < 12Dv. In addition, this structure is accompanied by an outward interaction. The intensity and extent of this structure are reduced in the presence of smaller patches, and disappear with a patch width ratio of Dv/Dc = 0.33. However, Q1 dominancy is also detectable in the mixing regions of small patches.
(3)
With a patch width ratio of Dv/Dc = 0.5, a new type of coherent shear stresses emerges in the wake zone. These structures grow as the dimensions of the vegetation patch reduce, and they come to cover most area of the wake zone with a patch width ratio of Dv/Dc = 0.33.
In summary, the flow in the downstream region of a decaying vegetation patch alters as the size of vegetation patch decreases. This change includes not only the intensity of occurrences, but also the nature of structure. This means that the dominant coherent structures are transferred from the intermittent anisotropic turbulent fluctuations in the mixing layer of a fully covered patch into the isotropic turbulence associated with the von Karman vortexes on the horizontal plane of the wake layer in the presence of the smallest patch.

5. Conclusions

In the present study, the changes in turbulent structures in the regions downstream of vegetation patches have been investigated. The model vegetation was selected carefully to simulate aquatic vegetation patches in natural rivers. Velocity profile, TKE, turbulent power spectra and quadrant analysis have been used via various approaches to investigate the features and intensity of the turbulent structures. Three different flow layers were detected in the downstream regions of vegetation patches, including the wake layer, the mixing layer and the shear layer. Overall, the results of the TKE and velocity analysis show that those three flow layers (wake, mixing and shear) occurred in all cases of patch setup. However, the wake and shear layers were considerably affected by the flow passing around the sides of the patches. Consequently, different flow structures and associated length scales can be formed, which are dependent on the upward movement of the center of the vertical rotation beyond the patch. The results of spectral analysis indicate that the strong coherent shear stresses are associated with the peaks of the anisotropic turbulent spectra. In contrast, the weak form of coherent shear stresses occurs in the presence of peak isotropic turbulence. The presence of von Karman vortexes near partially covered patches at a depth of Hv/H > 0.55~0.7 shows that the frequency of vortex shedding was around 0.07–0.11 Hz, which is equivalent to an average Strouhal number of about 0.25. The main characteristics of turbulent structures beyond a wide patch are associated with the highly intermittent anisotropic turbulent events in the mixing layer, and appear at a depth of z/Hv = 0.7~1.1 and distance of x/Hv = 8~12. The outward interactions are the representative quadrant class of coherent Reynolds shear stresses in the regions downstream of vegetation patches—particularly the strong coherent occurrences. The intensity and extent of these structures decrease as the size of the patch reduces. In addition, the intermittency and anisotropy of the turbulence are associated with the non-Gaussian distribution of the instantaneous velocity detected in the mixing layer beyond patches with a relatively high patch width ratio of Dv/Dc > 0.5. Finally, when the size of the patch reduces, von Karman vortexes appear in the wake layer, and form the dominant flow structures in the downstream region of the vegetation patch. It should be noted that the results of this research are limited by our use of a certain type of vegetation, and so it is strongly suggested that other types of the vegetation patches often found in natural rivers be considered in future research.

Author Contributions

M.K.; laboratory works, methodology, software and writing—original draft preparation, H.A.; laboratory supervisory, methodology and validation, J.S.; methodology and writing reviewing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study did not involve humans or animals.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available upon request.

Acknowledgments

Special thanks to G. Goodrati Amiri for his contribution in providing the laboratory devices and Esmaeel Dodangeh for his contribution in laboratory setup.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carpenter, S.R.; Lodge, D.M. Effects of submersed macrophytes on ecosystem processes. Aquat. Bot. 1986, 26, 341–370. [Google Scholar] [CrossRef]
  2. Martin, C.W.; Valentine, J.F. Eurasian milfoil invasion in estuaries: Physical disturbance can reduce the proliferation of an aquatic nuisance species. Mar. Ecol. Prog. Ser. 2012, 449, 109–119. [Google Scholar] [CrossRef]
  3. Short, F.T.; Kosten, S.; Morgan, P.A.; Malone, S.; Moore, G.E. Impacts of climate change on submerged and emergent wetland plants. Aquat. Bot. 2016, 135, 3–17. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Jeppesen, E.; Liu, X.; Qin, B.; Shi, K.; Zhou, Y.; Thomaz, S.M.; Deng, J. Global loss of aquatic vegetation in lakes. Earth-Sci. Rev. 2017, 173, 259–265. [Google Scholar] [CrossRef]
  5. Gillard, M.; Thiébaut, G.; Deleu, C.; Leroy, B. Present and future distribution of three aquatic plants taxa across the world: Decrease in native and increase in invasive ranges. Biol. Invasions 2017, 19, 2159–2170. [Google Scholar] [CrossRef]
  6. O’Hare, M.T.; Aguiar, F.C.; Asaeda, T.; Bakker, E.; Chambers, P.A.; Clayton, J.S.; Elger, A.; Ferreira, T.M.; Gross, E.M.; Gunn, I.D.M.; et al. Plants in aquatic ecosystems: Current trends and future directions. Hydrobiologia 2018, 812, 1–11. [Google Scholar] [CrossRef]
  7. Okamoto, T.-A.; Nezu, I. Spatial evolution of coherent motions in finite-length vegetation patch flow. Environ. Fluid Mech. 2013, 13, 417–434. [Google Scholar] [CrossRef]
  8. Afzalimehr, H.; Riazi, P.; Jahadi, M.; Singh, V.P. Effect of vegetation patches on flow structures and the estimation of friction factor. ISH J. Hydraul. Eng. 2019, 1–11. [Google Scholar] [CrossRef]
  9. Kazem, M.; Afzalimehr, H.; Sui, J. Formation of Coherent Flow Structures beyond Vegetation Patches in Channel. Water 2021, 13, 2812. [Google Scholar] [CrossRef]
  10. Killgore, K.J.; Dibble, E.D.; Hoover, J.J. Relationships between Fish and Aquatic Plants: A Plan of Study; Miscellaneous Paper; U.S. Army Corps of Engineers: Vicksburg, MI, USA, 1993; pp. A-91–A-93.
  11. Lauridsen, T.; Buenk, I. Diel changes in the horizontal distribution of zooplankton in the littoral zone of two shallow eutrophic lakes. Archiv. Hydrobiol. 1996, 137, 161–176. [Google Scholar] [CrossRef]
  12. Looby, A.; Reynolds, L.K.; Adams, C.R.; Martin, C.W. Submerged Aquatic Vegetation Patch Size Affects Fish Communities in a Turbid-Algal Lake. Front. Conserv. Sci. 2021, 2, 12. [Google Scholar] [CrossRef]
  13. Neary, V.S.; Constantinescu, S.G.; Bennett, S.J.; Diplas, P. Effects of Vegetation on Turbulence, Sediment Transport, and Stream Morphology. J. Hydraul. Eng. 2012, 138, 765–776. [Google Scholar] [CrossRef]
  14. Follett, E.; Nepf, H.M. Sediment patterns near a model patch of reedy emergent vegetation. Geomorphology 2012, 179, 141–151. [Google Scholar] [CrossRef] [Green Version]
  15. Huai, W.-X.; Li, S.; Katul, G.G.; Liu, M.-Y.; Yang, Z.-H. Flow dynamics and sediment transport in vegetated rivers: A review. J. Hydrodyn. 2021, 33, 400–420. [Google Scholar] [CrossRef]
  16. Zhang, J.; Wang, W.-J.; Shi, H.; Wang, W.; Li, Z.; Tang, X.; Xia, Z. An analytical two-layer model for velocity distribution in open-channel flows with submerged flexible canopies considering multiply fluids mechanics. J. Hydrol. 2021, 603, 127102. [Google Scholar] [CrossRef]
  17. Marjoribanks, T.I.; Paul, M. Modelling flow-induced reconfiguration of variable rigidity aquatic vegetation. J. Hydraul. Res. 2021, 1–16. [Google Scholar] [CrossRef]
  18. D’Ippolito, A.; Calomino, F.; Alfonsi, G.; Lauria, A. Flow Resistance in Open Channel Due to Vegetation at Reach Scale: A Review. Water 2021, 13, 116. [Google Scholar] [CrossRef]
  19. Li, S.; Shi, H.; Xiong, Z.; Huai, W.; Cheng, N. New formulation for the effective relative roughness height of open channel flows with submerged vegetation. Adv. Water Resour. 2015, 86, 46–57. [Google Scholar] [CrossRef]
  20. Zhao, F.; Huai, W. Hydrodynamics of discontinuous rigid submerged vegetation patches in open-channel flow. HydroResearch 2016, 12, 148–160. [Google Scholar] [CrossRef]
  21. Hopkinson, L.C. Flow through vegetation patches with multiple layers: A quadrant analysis. In Crossing Boundaries, Proceedings of the World Environmental and Water Resources Congress, Albuquerque, NM, USA, 20–24 May 2012; American Society of Civil Engineers: Reston, VI, USA, 2012; pp. 1406–1415. [Google Scholar]
  22. Devi, T.B.; Kumar, B. Channel Hydrodynamics of Submerged, Flexible Vegetation with Seepage. J. Hydraul. Eng. 2016, 142, 04016053. [Google Scholar] [CrossRef]
  23. Chen, Z.; Ortiz, A.; Zong, L.; Nepf, H. The wake structure behind a porous obstruction and its implications for deposition near a finite patch of emergent vegetation. Water Resour. Res. 2012, 48. [Google Scholar] [CrossRef]
  24. Liu, C.; Nepf, H. Sediment deposition within and around a finite patch of model vegetation over a range of channel velocity. Water Resour. Res. 2016, 52, 600–612. [Google Scholar] [CrossRef] [Green Version]
  25. Ortiz, A.C.; Ashton, A.; Nepf, H. Mean and turbulent velocity fields near rigid and flexible plants and the implications for deposition. J. Geophys. Res. Earth Surf. 2013, 118, 2585–2599. [Google Scholar] [CrossRef]
  26. Zong, L.; Nepf, H. Flow and deposition in and around a finite patch of vegetation. Geomorphology 2010, 116, 363–372. [Google Scholar] [CrossRef]
  27. Zong, L.; Nepf, H. Vortex development behind a finite porous obstruction in a channel. J. Fluid Mech. 2012, 691, 368–391. [Google Scholar] [CrossRef]
  28. Liu, C.; Hu, Z.; Lei, J.; Nepf, H. Vortex structure and sediment deposition in the wake behind a finite patch of model sub-merged vegetation. J. Hydraul. Eng. 2018, 144, 04017065. [Google Scholar] [CrossRef]
  29. Folkard, A.M. Hydrodynamics of model Posidonia oceanica patches in shallow water. Limnol. Oceanogr. 2005, 50, 1592–1600. [Google Scholar] [CrossRef]
  30. Sadeque, M.A.; Rajaratnam, N.; Loewen, M.R. Effects of Bed Roughness on Flow around Bed-Mounted Cylinders in Open Channels. J. Eng. Mech. 2009, 135, 100–110. [Google Scholar] [CrossRef]
  31. Perera, M.D.A.E.S. Shelter behind two-dimensional solid and porous fences. J. Wind Eng. Ind. Aerodyn. 1981, 8, 93–104. [Google Scholar] [CrossRef]
  32. Sui, J.; Faruque, M.A.; Balachandar, R. Local Scour Caused by Submerged Square Jets under Model Ice Cover. J. Hydraul. Eng. 2009, 135, 316–319. [Google Scholar] [CrossRef]
  33. Jafari, R.; Sui, J. Velocity Field and Turbulence Structure around Spur Dikes with Different Angles of Orientation under Ice Covered Flow Conditions. Water 2021, 13, 1844. [Google Scholar] [CrossRef]
  34. Kabiri, F.; Afzalimehr, H.; Sui, J. Flow structure over a wavy bed with vegetation cover. Int. J. Sediment Res. 2017, 32, 186–194. [Google Scholar] [CrossRef]
  35. Hussain, A.F. Role of coherent structures in turbulent shear flows. Proc. Indian Acad. Sci. Sect. C Eng. Sci. 1981, 4, 129–175. [Google Scholar]
  36. Shahmohammadi, R.; Afzalimehr, H.; Sui, J. Impacts of turbulent flow over a channel bed with a vegetation patch on the incipient motion of sediment. Can. J. Civ. Eng. 2018, 45, 803–816. [Google Scholar] [CrossRef]
  37. Afzalimehr, H.; Moghbel, R.; Gallichand, J.; Sui, J. Investigation of turbulence characteristics in channel with dense vegetation. Int. J. Sediment Res. 2011, 26, 255–282. [Google Scholar] [CrossRef]
  38. Przyborowski, Ł.; Łoboda, A.M.; Bialik, R.J. Effect of two distinct patches of Myriophyllum species on downstream turbulence in a natural river. Acta Geophys. 2019, 67, 987–997. [Google Scholar] [CrossRef] [Green Version]
  39. Huai, W.-X.; Zhang, J.; Katul, G.G.; Cheng, Y.-G.; Tang, X.; Wang, W.-J. The structure of turbulent flow through submerged flexible vegetation. J. Hydrodyn. 2019, 31, 274–292. [Google Scholar] [CrossRef]
  40. Termini, D. Experimental Analysis of Horizontal Turbulence of Flow over Flat and Deformed Beds. Arch. Hydro-Eng. Environ. Mech. 2015, 62, 77–99. [Google Scholar] [CrossRef] [Green Version]
  41. Termini, D.; Sammartano, V. Experimental observation of horizontal coherent turbulent structures in a straight flume. In Proceedings of the River, Coastal and Estuarine Morphodynamics—RCEM, Santa Fe City, Argentina, 21–29 September 2009. [Google Scholar]
  42. Mayaud, J.R.; Wiggs, G.F.; Bailey, R.M. Dynamics of skimming flow in the wake of a vegetation patch. Aeolian Res. 2016, 22, 141–151. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, H.; Peng, G.; Chen, M.; Fan, J. Analysis of the Interconnections between Classic Vortex Models of Coherent Structures Based on DNS Data. Water 2019, 11, 2005. [Google Scholar] [CrossRef] [Green Version]
  44. El-Hakim, O.; Salama, M.M. Velocity Distribution Inside and Above Branched Flexible Roughness. J. Irrig. Drain. Eng. 1992, 118, 914–927. [Google Scholar] [CrossRef]
  45. Nezu, I.; Sanjou, M. Turburence structure and coherent motion in vegetated canopy open-channel flows. HydroResearch 2008, 2, 62–90. [Google Scholar] [CrossRef]
  46. Shivpure, V.; Devi, T.B.; Kumar, B. Turbulent characteristics of densely flexible submerged vegetated channel. ISH J. Hydraul. Eng. 2016, 22, 220–226. [Google Scholar] [CrossRef]
  47. Tchen, C. On the spectrum of energy in turbulent shear flow. J. Res. Natl. Inst. Stand. Technol. 1953, 50, 51. [Google Scholar] [CrossRef]
  48. Lian, Q.X. A visual study of the coherent structure of the turbulent boundary layer in flow with adverse pressure gradient. J. Fluid Mech. 1990, 215, 101–124. [Google Scholar] [CrossRef]
  49. Staicu, A.D. Intermittency in Turbulence. Ph.D. Thesis, Technische Universiteit Eindhoven, Eindhoven, The Netherlands, 2002. [Google Scholar] [CrossRef]
  50. Coscarella, F.; Penna, N.; Servidio, S.; Gaudio, R. Turbulence anisotropy and intermittency in open-channel flows on rough beds. Phys. Fluids 2020, 32, 115127. [Google Scholar] [CrossRef]
  51. Avila, M.; Hof, B. Nature of laminar-turbulence intermittency in shear flows. Phys. Rev. E 2013, 87, 063012. [Google Scholar] [CrossRef] [Green Version]
  52. Massey, F.J., Jr. The Kolmogorov-Smirnov test for goodness of fit. J. Am. Stat. Assoc. 1951, 46, 68–78. [Google Scholar] [CrossRef]
  53. Afzalimehr, H.; Sui, J.; Moghbel, R. Hydraulic parameters in channels with wall vegetation and gravel bed–field observations and experimental studies. Int. J. Sediment Res. 2010, 25, 81–90. [Google Scholar] [CrossRef]
  54. Afzalimehr, H.; Barahimi, M.; Sui, J. Non-uniform flow over cobble bed with submerged vegetation strip. Proc. ICE–Water Manag. 2019, 172, 86–101. [Google Scholar] [CrossRef] [Green Version]
  55. Fazlollahi, A.; Afzalimehr, H.; Sui, J. Impacts of pool and vegetated banks on turbulent flow characteristics. Can. J. Civil Eng. 2015, 42, 979–986. [Google Scholar] [CrossRef] [Green Version]
  56. Cox, D.T.; Kobayashi, N. Identification of intense, intermittent coherent motions under shoaling and breaking waves. J. Geophys. Res. Space Phys. 2000, 105, 14223–14236. [Google Scholar] [CrossRef]
  57. Cox, D.T.; Anderson, S.L. Statistics of Intermittent Surf Zone Turbulence and Observations of Large Eddies using PIV. Coast. Eng. J. 2001, 43, 121–131. [Google Scholar] [CrossRef]
  58. McCaffrey, K.; Fox-Kemper, B.; Hamlington, P.E.; Thomson, J. Characterization of turbulence anisotropy, coherence, and intermittency at a prospective tidal energy site: Observational data analysis. Renew. Energy 2015, 76, 441–453. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (A) The layout of the experimental device, (B) a sample of vegetation patches (120 × 90 cm, Dv/Dc = 1), (C) sample of vegetation patch in a natural gravel bed river (Marbor River, Iran), (D) a sample of bed materials, (E) a single synthetic plant used to simulate patch; (F) a sketch of the experimental setup (not to scale).
Figure 1. (A) The layout of the experimental device, (B) a sample of vegetation patches (120 × 90 cm, Dv/Dc = 1), (C) sample of vegetation patch in a natural gravel bed river (Marbor River, Iran), (D) a sample of bed materials, (E) a single synthetic plant used to simulate patch; (F) a sketch of the experimental setup (not to scale).
Water 13 03468 g001
Figure 2. Streamwise velocity profiles along the centerline of the flume; x/Hv is the normalized distance from the downstream edge of the vegetation patch, and H is the depth of flow. The boundaries between zones affected by the wake and over the canopy flow are shown as red lines. The velocity of the Log–Law zone is shown as dashed and is not to scale.
Figure 2. Streamwise velocity profiles along the centerline of the flume; x/Hv is the normalized distance from the downstream edge of the vegetation patch, and H is the depth of flow. The boundaries between zones affected by the wake and over the canopy flow are shown as red lines. The velocity of the Log–Law zone is shown as dashed and is not to scale.
Water 13 03468 g002
Figure 3. Normalized TKE and streamwise velocity behind 4 different patch layouts near channel bed (z/Hv = 0), at the middle height of patch (z/Hv = 0.5) and top of patch (z/Hv = 1).
Figure 3. Normalized TKE and streamwise velocity behind 4 different patch layouts near channel bed (z/Hv = 0), at the middle height of patch (z/Hv = 0.5) and top of patch (z/Hv = 1).
Water 13 03468 g003
Figure 4. Samples of the power spectra behind the patch showing different classes of turbulence: (A) isotropic–homogenious turbulent flow without dominant frequency on Suw, (B) isotropic–homogenious turbulent flow with dominant frequency on Suw, observed below the height of the patch, (C) anisotropic turbulent flow with dominant frequency on Suw, observed alongside the height of patch, (D) Suw of the same anisotropic sample.
Figure 4. Samples of the power spectra behind the patch showing different classes of turbulence: (A) isotropic–homogenious turbulent flow without dominant frequency on Suw, (B) isotropic–homogenious turbulent flow with dominant frequency on Suw, observed below the height of the patch, (C) anisotropic turbulent flow with dominant frequency on Suw, observed alongside the height of patch, (D) Suw of the same anisotropic sample.
Water 13 03468 g004
Figure 5. Values of u w ¯ along the centerline of the channel beyond the vegetation patch (colored figures) beside values of τ c / τ r (graysclae). Note that the blank area in τ c / τ r is where no dominant frequency was detectable in Suw. The end of the patch is located at x = 800.
Figure 5. Values of u w ¯ along the centerline of the channel beyond the vegetation patch (colored figures) beside values of τ c / τ r (graysclae). Note that the blank area in τ c / τ r is where no dominant frequency was detectable in Suw. The end of the patch is located at x = 800.
Water 13 03468 g005
Figure 6. Probability of occurrence of quadrant classes behind the patch; the end of the patch is located at x = 800. The height of the patches is 105 ± 5 mm, and the values are shown up to level of the maximum domain of ADV capability (~140 mm).
Figure 6. Probability of occurrence of quadrant classes behind the patch; the end of the patch is located at x = 800. The height of the patches is 105 ± 5 mm, and the values are shown up to level of the maximum domain of ADV capability (~140 mm).
Water 13 03468 g006
Figure 7. Conceptual illustration of the quadratic bursting events in the boundary between the mixing layer and the shear layer.
Figure 7. Conceptual illustration of the quadratic bursting events in the boundary between the mixing layer and the shear layer.
Water 13 03468 g007
Figure 8. Values of f k ( Q 1 ) and f d u w for 278 points in the downstream region of the vegetation patches.
Figure 8. Values of f k ( Q 1 ) and f d u w for 278 points in the downstream region of the vegetation patches.
Water 13 03468 g008
Figure 9. Values of ErrorCDF for different depths behind the decaying patch at a distance of x/Hv = 8.
Figure 9. Values of ErrorCDF for different depths behind the decaying patch at a distance of x/Hv = 8.
Water 13 03468 g009
Figure 10. Probability density function and the turbulent spectra for the mixing layer of Case 1.
Figure 10. Probability density function and the turbulent spectra for the mixing layer of Case 1.
Water 13 03468 g010
Figure 11. Formation of the intermittent fluctuations in the mixing layer beyond the patch (case 1). The filtered velocity time series is shown as intermittent fluctuations above the original velocity time series. Variations in the normalized standard deviation of TKE and normalized intermittent Turbulent Kinetic Energy are shown at the bottom.
Figure 11. Formation of the intermittent fluctuations in the mixing layer beyond the patch (case 1). The filtered velocity time series is shown as intermittent fluctuations above the original velocity time series. Variations in the normalized standard deviation of TKE and normalized intermittent Turbulent Kinetic Energy are shown at the bottom.
Water 13 03468 g011
Table 1. Experimental parameters for the cases discussed in this paper.
Table 1. Experimental parameters for the cases discussed in this paper.
CaseQ (Discharge, L/s)n/m2 (Number of Veg. per Square Meter)Lv (Length of Patch, cm)Dv (Width of Patch, cm)Hv (Height of Patch, cm)
131 L/s611.11209010
231 L/s611.1906010
331 L/s611.1604510
431 L/s611.1403010
No. Veg.31 L/s----
Table 2. Results of the normality test of velocity time series of case 1 at a distance of x/Hv=8. Note that the unmatched values are highlighted.
Table 2. Results of the normality test of velocity time series of case 1 at a distance of x/Hv=8. Note that the unmatched values are highlighted.
z/Hv% of Sample between μ − σ and μ + σ% of Sample between μ − 2σ and μ + 2σ% of Sample between μ − 3σ and μ + 3σSimple Normality TestKolmogorov–Smirnov Test
Reference Value0.68950997
0.20.6790.9580.997Gaussiannon-Gaussian
0.50.6960.9540.995~Gaussiannon-Gaussian
0.80.7540.9480.981non-Gaussiannon-Gaussian
10.7100.9380.998non-Gaussiannon-Gaussian
1.20.7530.9280.993non-Gaussiannon-Gaussian
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kazem, M.; Afzalimehr, H.; Sui, J. Characteristics of Turbulence in the Downstream Region of a Vegetation Patch. Water 2021, 13, 3468. https://doi.org/10.3390/w13233468

AMA Style

Kazem M, Afzalimehr H, Sui J. Characteristics of Turbulence in the Downstream Region of a Vegetation Patch. Water. 2021; 13(23):3468. https://doi.org/10.3390/w13233468

Chicago/Turabian Style

Kazem, Masoud, Hossein Afzalimehr, and Jueyi Sui. 2021. "Characteristics of Turbulence in the Downstream Region of a Vegetation Patch" Water 13, no. 23: 3468. https://doi.org/10.3390/w13233468

APA Style

Kazem, M., Afzalimehr, H., & Sui, J. (2021). Characteristics of Turbulence in the Downstream Region of a Vegetation Patch. Water, 13(23), 3468. https://doi.org/10.3390/w13233468

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop