Next Article in Journal
Study on the Sensitivity of the Streamwise Location of MVG on SWBLI in MVG-Based Supersonic Flow Control
Previous Article in Journal
Aerodynamic Drag Reduction of Railroad Tank Wagons
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Field Study of Three–Parameter Flow Resistance Model in Rivers with Vegetation Patch

1
Faculty of Civil Engineering, Iran University of Science and Technology, Tehran 16846-13114, Iran
2
Faculty of Civil Engineering, Tehran University, Tehran 14176-13131, Iran
3
Department of Mathematics and Computer Science, West Virginia State University, Institute, WV 25112, USA
4
Department of Civil and Environmental Engineering and Construction, University of Nevada, Las Vegas, NV 89154, USA
*
Author to whom correspondence should be addressed.
Submission received: 15 July 2022 / Revised: 13 August 2022 / Accepted: 18 August 2022 / Published: 22 August 2022

Abstract

:
Bed shear stress in coarse–bed rivers with vegetation patches is one of the challenging parameters in hydraulic engineering, mechanical engineering, fluvial morphology, and environmental studies. Based on this necessity, in this study, the values of bed shear stress in four reaches of rivers in Iran were estimated and compared using the methods of boundary layer characteristics, logarithmic law, and Darcy–Weisbach. Data collection in this study started in February 2021 and ended in April 2021. Estimation of flow resistance is a key factor in many numerical and physical models. In order to obtain a reasonable evaluation of this factor, it is necessary to measure and calculate the key variables of resistance to flow. Accordingly, the experimental design in this study includes surveying operations, velocity measurement, and sampling of bed sediments. The results show that due to bed forms, vegetation patches, and variations of flow depth and grain size in the river, the universal velocity distribution law (the log law) may not be suitable to estimate the shear velocity, which is a key parameter of flow resistance. This calls for more justifiable methods which are not affected by near–the–bed conditions. Accordingly, a three–parameter flow resistance model is presented, which shows an average error of 17%, indicating the accuracy of the model. The investigation of 71 measured velocity profiles shows the occurrence of the Dip phenomenon in the velocity profiles near the vegetation patches. However, by moving away from the vegetation patches, the effect of this phenomenon is decreased, and the profiles illustrate an S–shaped distribution. The results show that the relative differences between the logarithmic law and Darcy–Weisbach methods compared to the boundary layer characteristics method (BLCM) are equal to 87% and 39%, respectively, indicating a more reasonable agreement between the Darcy–Weisbach method and the boundary layer characteristics method. This is due to the application of key parameters of the boundary layer theory to calculate shear velocity by BLCM. However, to simplify data collection in the field, the application of the Darcy–Weisbach method is suggested.

1. Introduction

Vegetation patches in riverbeds affect the main flow parameters, such as shear stress and velocity, change in the flow structure, and increase in the overall resistance of the flow. The most apparent effect of vegetation patches in water transfer canals is to reduce the shear stress of the bed, increase the flow resistance, and reduce the water transfer capacity. Vegetation patches are widely observed in rivers, and their presence in the flow path can change the velocity profile and shear stress. These patches are different from other obstacles such as bridge piers due to their permeability [1]. In the past, it was supposed that vegetation only increases the flow resistance; for this reason, to increase the flow transmission capacity, the vegetation at front of the stream was removed; while today, we know that it has several benefits. For instance, the patches can create a suitable environment for aquatic life, growth, and spawning [2,3] and reduce erosion by decreasing shear stress [4,5,6,7].
Bed shear stress is one of the main parameters in the river and hydraulic engineering, the estimation of which is always accompanied by an error in river engineering. Up to the present date, a complete model for coarse–bed rivers based on the theoretical foundations of boundary layer and flow velocity distribution is unavailable. Most models are entirely experimental, and the effect of flow structure is not considered in them. Bed shear stress is a significant factor in estimating hydraulic parameters, design of hydraulic structures, and river engineering studies. In rivers, the bed shear stress is essential in determining the local erosion of the bed, the roughness coefficient, the sediment parameters of the bed load, and the suspended load. Accurate estimation of bed shear stress is difficult, especially in mountain rivers. In laboratory studies, Prandtl–Pitot tubes or shear plates are used to measure the bed shear stress, while in field studies, direct measurements are rarely done, and various indirect methods are used. There has been numerous research on the interaction of vegetation and flow in the presence of artificial vegetation and the laboratory flume (e.g., [8,9,10,11]); however, few studies have been conducted on natural vegetation in coarse–bed rivers (e.g., [12]).
Thoman and Niezgoda (2008) investigated the effect of vegetation on bed shear stress, which is a key factor in determining the pattern of sedimentation and erosion. They showed that the presence of vegetation reduces the velocity in the regions near the bed and thus reduces shear stress and erosion [13].
Depending on whether the vegetation is flexible or rigid, as well as submerged or emergent, models have been developed to estimate flow resistance in the presence of vegetation [14,15]. However, so far, a comprehensive model for predicting flow resistance has not been found [16]. Since most of these models have been developed under laboratory and controlled conditions, their generalization to river conditions cannot be easily done. Therefore, it is necessary to develop a comprehensive model for river conditions.
Due to the drought phenomenon, the river flow is not established in all seasons, which results in vegetation growth in the river beds. Nosrati et al. (2022) focused their study on the effect of irregular vegetation on flow resistance in a sandy river. Their results showed that the effect of vegetation cover on the total flow resistance is greater than the effect of the bed form of the river. Therefore, in order to correctly estimate the flow resistance in mountain rivers and its application in hydraulic models, engineers and designers must pay special attention to the contribution of vegetation in the estimation of flow resistance [17]. In another study, Afzalimehr et al. (2019) studied the flow on a sandy bed with walls covered with vegetation. They showed that the presence of vegetation causes the flow to change from uniform to non–uniform despite the constant depth of flow and constant slope. They showed that the reason for this change in the flow structure is the non–uniform distribution of the Reynolds stress and the occurrence of the maximum velocity below the water surface due to the presence of vegetation on the channel wall [18].
The most common method for estimating the average bed shear stress in a reach is defined as follows:
τ 0 ρ = g R S f
where τ 0   is bed shear stress, ρ is the specific gravity of water, g is the gravitational acceleration, R   is the hydraulic radius, and   S f   is the energy slope. This method is not suitable for local estimation of bed shear stress. Assuming that the hydraulic radius is equal to the depth of the water flow, and by substituting the slope of the energy   S f   from the Manning equation, Equation (1) becomes:
τ 0 ρ = g n 2 U 2 h 1 3  
where n is the Manning coefficient, U is the flow velocity, and h is the flow depth. Equation (2) represents the most common approach adopted in most hydrodynamic software packages. The Manning equation is not a reliable method in rivers; because the Manning roughness coefficient considers the water depth and flow velocity as constant, which is different from the actual conditions due to the non–uniformity of the flow and the presence of obstacles in rivers [19]. Afzalimehr et al. (1999) showed that applying the Manning coefficient with the assumption of uniform flow in a river is accompanied by an error of approximately 500% [20].
Bellos et al. (2018) proposed a new three–parameter flow resistance model for open channels in which the bed shear stress is proportional to the flow velocity and inversely to the water depth as the following [21]:
τ 0 ρ = U A B h C
where A, B, and C are the three parameters of the flow resistance model. For A = 2, B = 1 / ( g n ) 2 , and C = 1/3, this model is the same as the classic flow resistance model (Equation (2)), which is based on the Manning equation. Therefore, it can be concluded that the classic flow resistance model is a particular case of the proposed resistance model that includes three parameters instead of one parameter. This model is valid for all flow and roughness properties, and it is proposed to replace the Manning equation in flood simulation models. By combining Equations (1) and (3), the flow velocity is calculated as the following:
U = ( g B R C + 1 S f ) 1 A      
Equation (4) estimates the flow velocity by considering the effect of vegetation patch and roughness size through coefficients A, B, and C, to better evaluate the flow resistance. The three parameters of the mentioned flow resistance model (A, B, C) are calibrated based on random data from a uniform probability distribution of water depth, flow velocity, and roughness height; however, it has not been calibrated based on laboratory and field data. In addition, the effect of vegetation has not been considered for calibration of these parameters, which must be considered. This research aims to calibrate this model based on different bed shear stress methods by considering vegetation patches. There are various methods for calculating bed shear stress. In this study, three widely used methods have been implemented to calibrate the three–parameter flow resistance model. The descriptions of these three methods are given below.
Shear velocity and local bed shear stress can be calculated using the relationship based on the boundary layer characteristics as follows [22]:
u * = ( δ * θ ) u m a x C δ *
τ 0 = ρ u * 2    
δ * = 0 h ( 1 u u m a x ) d y
θ = 0 h u u m a x ( 1 u u m a x ) d y
where δ *   and θ are the thickness of the boundary layer displacement and the thickness of the momentum, respectively; u m a x is the maximum velocity in the velocity profile, and C is a constant value obtained in coarse–bed rivers which is equal to 4.4 [23]. The researchers have used boundary layer characteristics for sandy and gravely beds in the last 20 years [18,24,25]. This method is not applicable in boulder–bed rivers.
The velocity profile in the inner layer where the viscosity is dominant follows a logarithmic law and is as follows:
u u * = 1 k l n ( z k s ) + B r
where z is the distance from the bed, ks is the equivalent roughness height, k is the Von Kármán constant, and Br is the integration parameter. The logarithmic law method for determining the bed shear stress is one of the most valid methods. The following points are significant about this method:
  • Emadzadeh et al. (2010) and Fazlollahi et al. (2007) showed that the Von Kármán constant is less than 0.4 in the conditions of sediment transfer and more than 0.4 in the presence of the bed form; and in the absence of these cases, it is equal to 0.4. In this study, the value of 0.4 has been used [26,27]. Naderi et al. (2021) found a value of 0.4 for reaches with vegetation cover when there is no sediment transport, and the effect of bedforms is insignificant [28].
  • Different researchers have found different values for ks; For example, Song and Chiew (2001) presented the value of 0.25   d 50   and Alonso et al. (2009) presented the value of 2 . 4   d 90 . The amount suggested by Alonso et al. (2009) has been used in this study [29,30]. In this study, to show the effect of larger particles on surface roughness in coarse rivers, ks = 2.4 d90 was considered.
  • In coarse–bed rivers, there may be a flow between the coarse–grained particles and below the measured points; therefore, the flow has a velocity at zero depth and cannot be considered zero. Therefore, to solve this problem, a concept called hypothetical line is used. A hypothetical line is a line on which the velocity is zero. This line is defined as follows:
    z = z m e a s + α d p      
Song and Chiew (2001) suggested a value of 0.25 for α , which is also considered in this study. Different researchers have proposed different values for p . The value of d84 is usually used to show the greater effect of larger particles on surface roughness in coarse–bed rivers [29]. According to the hypothetical line definition, the final relationship of the logarithmic law method that is valid in the inner layer of the flow is as follows:
u = u * k l n ( z + 0.25 d 84 d 84 ) + B r u *
To determine the shear velocity and bed shear stress using a logarithmic law, the validity range of this law must first be calculated. The coefficient of determination index (R2) is used to determine the validity range of the logarithmic law. Depth with the highest coefficient of determination is considered the upper limit of the logarithmic law. After regression between u and l n ( z + 0.25 d 84 d 84 ) and calculation of the slope of the regression line (m) and y–intercept (b), the following equations are used to calculate the shear velocity, bed shear stress, and integration constant:
u * = m k   τ 0 = ρ u * 2    
B r = b u *
Another method for estimating shear stress is the Darcy–Weisbach method, which is defined as follows:
τ 0 = f 8 ρ U 2
where f is the Darcy–Weisbach friction factor. In this study, the method proposed by Bellos et al. [21] was used to calculate the Darcy–Weisbach friction factor (Equation (15)).
The main difference between this research and other research presented in the literature is the use of field data to validate different methods of bed shear stress in the presence of vegetation patches. The originality of this research is the use of field measurements for calibration and validation of the flow resistance model.
To achieve the research objectives, the research assumptions include the following:
  • The selected reaches have a direct path without the presence of obstacles such as boulders.
  • The transfer of sediment from the bed is insignificant.
  • The wind velocity in the area does not affect the velocity profiles.
  • There is no flexibility in vegetation patches and deformation during data collection.
In many seasonal rivers in dry regions, there is no flow in rivers. This leads to growing vegetation patches through rivers. However, no research has reported the influence of these patches on flow resistance estimation. This study applies different methods of shear stress estimation and reveals that the most suitable method to predict flow resistance is the three–parameter flow resistance model.
This paper consists of four main sections. In the first section, the theoretical foundations and various relationships for calculating the shear stress of the bed are discussed. The materials and methods section explains how to measure the data and calculate the parameters used in this study. In the results section, different methods of bed shear stress in the presence of vegetation patches are compared and analyzed. Finally, the main findings of the research are presented.

2. Materials and Methods

2.1. Reaches Location

Several factors are considered in choosing the appropriate reach in the river, the most important of which is the safety of the researcher and the objectives of the research. Accordingly, the flow velocity and the average depth of the selected reach should be such that the data collection is done with minimal errors and fluctuations and the researcher’s safety is not risked. According to these points, finally, two reaches of Shapur River (Fars province), one reach of Fahlyan River (Fars province), and one reach of Dalaki River (Bushehr province) were selected. All three selected rivers are located in Iran. Figure 1 shows the geographical location of the studied reaches.
The selected reaches had a straight path without physical obstacles such as boulders. Data collection time started in February 2021 and ended in April 2021. A summary of the general information of the selected reaches is presented in Table 1. In the table, Lr is the reach length, D is the average hydraulic depth, W is the average width, S0 is the longitudinal slope, and Q is the discharge. According to the table, Shapur1 and Shapur2 have the highest cross–section width and depth of the flow, respectively, and the largest discharge belongs to Shapur1. Brierley and Fryirs (2005) stated that the longitudinal slope for sandy and gravely rivers is in the range of 0.05% to 0.5%, which in this study is in the same range for the longitudinal slope as well [31].

2.2. Vegetation Patches in Selected Reaches

The four studied reaches were selected with different amounts of vegetation patches. Shapur1 and Dalaki reaches on their left bank (looking downstream) have emerged vegetation patches. There are scattered submerged vegetation patches in the bed of these two reaches with different dimensions. Selected reaches in the Fahlyan and Shapu2 rivers have submerged vegetation patches in the bed and emerged vegetation patches on their right bank. The average height of submerged vegetation patches in Shapur1, Shapur2, Fahlyan, and Dalaki reaches were 18, 15, 12, and 9 cm, respectively. Figure 2 shows some vegetation patches in the selected reaches. Table 2 presents the characteristics of the vegetation patches of the reaches, and average relative submergence (ratio of vegetation patch height to flow height), cross–sectional blockage factor (area of the cross–section occupied by the vegetation patch), surface blockage factor (area of the study reach occupied by the vegetation patch), and volumetric blockage factor (ratio of volume covered by vegetation patch to canal volume) are h p h , BX, BA, and BV, respectively.

2.3. Measured Data

Measurement data in this study include velocity measurement, surveying operations, and sampling of bed sediments, respectively. The order of data measurement in river engineering is important; for example, if sampling of bed sediments is done before surveying operations, changes in riverbed level may cause the surveying operation data to be inaccurate.

2.3.1. Velocity Measurement

Velocity measurement at different points and depths of the stream is one of the essential tasks in river engineering studies. Various velocimeters are used to measure flow velocity. The most common velocimeter for measuring point velocity in field studies is the butterfly current meter (BCM), which was also used in this study. To reduce the error in velocity measurement, in this study, the measurement time was considered 20 seconds, and the number of repetitions in velocity measurement was considered 3 to 5 times. In velocity profiles, the number of measured points near the bed (inner layer) was considered more so that the distance between the measuring points was 0.5 cm. In each reach, 4 to 5 cross–sections, and in each cross–section, 3 to 4 velocity profiles, one of which was located in the middle of the cross–section, were measured. Moreover, 16 to 22 points (an average of 19 points) were considered in each profile. In addition, to calculate the discharge in each cross–section more accurately, the flow velocity was measured at 0.2 and 0.8 depths from the bed at a distance of 1 m in each cross–section. For example, for the first cross–section of Shapur1, which had a width of 28.4 m, four velocity profiles with an average of 19 points per profile and 28 two–point velocity profiles were measured. Figure 3 shows the measurement of velocity at selected reaches. Table 3 shows the number of velocity measurement points at the four selected reaches. In the Table, N p is the number of velocity measurement points in the 19–point profile, N T P is the number of velocity measurement points in the two–point profile, and Fr is the frequency of velocity measurements by considering 3 to 5 repetitions at each point. According to the table, the total velocity of 2174 points has been measured in these four reaches, which by considering 3 to 5 repetitions for each point, the total number of velocity measurements is 8184 times.

2.3.2. Surveying Operations

In order to calculate the area, wet perimeter of cross–sections, longitudinal slope, and the dimensions of vegetation patches in the studied reaches, in this study, Leica TS02 Total Station Camera was used. For surveying, the selected reaches were networked by ropes with dimensions of 1 m in the longitudinal direction and 0.5 m in the cross–section direction of the river, and then the surveying was performed. Figure 4 shows the surveying in the studied reaches. Table 4 presents the number of surveyed points in the studied reaches. In the table, NB is the number of surveyed points from the riverbed, and NV is the number of surveyed points from the vegetation patches. DB and DV are the density of surveyed points for riverbed and vegetation, respectively. According to Table 4, a total of 7272 points (including bed and vegetation patches) have been surveyed.
Figure 5 shows images of the plan of selected reaches along with the position of the measured profiles (including two–point and 19–point profiles) and surveyed points. In this figure, the black dots indicate the places where the surveying was done, and in the red dots, the velocity profile was also measured in addition to the surveying. Cross–sections were selected to measure velocity and mapping as perpendicular to the flow direction.

2.3.3. Sampling of Bed Sediments

Using Wolman’s method, 100 particles were randomly sampled in each cross–section in the central, left bank, and right bank axes to present grain size distribution (Figure 6 and Figure 7) [32].

2.4. Calculation of Darcy–Weisbach Friction Factor

Bellos (2018) presented a relationship for estimating the Darcy–Weisbach friction factor that is valid for all flow regimes (from laminar to turbulent) and depends on the flow characteristics (water depth and flow velocity) and the roughness height as follows [21]:
f = ( 24 R h ) α × [ 0.86 e W ( 1.35 R h ) R h ] 2 ( 1 α ) b × { 1.34 [ l n ( 12.21 h k s ) ] 2 } ( 1 α ) ( 1 b )
where Rh is the Reynolds number of the flow and ks is the roughness height. In this study, the value of ks = 2.4 d90 has been used [30]. The parameters α , b, and W(1.35 Rh) are defined as follows:
α = 1 1 + ( R h 678 ) 8.4      
b = 1 1 + ( R h 150 h k s ) 1.8
W ( 1.35 R h ) = l n ( 1.35 R h ) l n [ l n ( 1.35 R h ) ] + l n [ l n ( 1.35 R h ) ] l n ( 1.35 R h ) + l n [ l n ( 1.35 R h ) ] 2 2 l n [ l n ( 1.35 R h ) ] 2 [ l n ( 1.35 R h ) ] 2

2.5. The General Framework of This Research

In this research, first, the bed shear stress is calculated using logarithmic law methods, boundary layer characteristics, and the Darcy–Weisbach method. Then, these methods are compared with each other. Finally, the three–parameter flow resistance model is calibrated and validated based on these three methods. Figure 8 shows a summary of the overall research framework.

3. Results and Discussion

3.1. Sediment and Hydraulic Parameters

Sediment and hydraulic parameters of the flow in four selected reaches are presented in Table 5. In the table, D is the hydraulic depth of the flow, Re is the Reynolds number, F r = U / ( g D ) 1 / 2 is the Froude number, and σ g = d 84 d 16 is the geometric standard deviation where d16 and d84 are the grain sizes that are larger than 16% and 84% of the bed material’s mass, respectively. The difference in the amount of discharge in different cross–sections of each reach is due to water supply for agricultural lands and urban irrigation downstream of the selected reaches. The Froude number presented in Table 5 shows that in all cross–sections, this number is less than 1; therefore, the flow in all four studied reaches is sub–critical. The study of more than 100 reaches in Iranian rivers showed that the range of Froude number in mountainous rivers is 0.12 to 0.56, confirming this statement [28]. Moreover, the Reynolds number values show that the flow is turbulent in all four selected reaches. According to Table 5, the particle size distribution of the bed in 74% of the cross–sections is non–uniform where the geometric standard deviation is greater than 1.4 [33].

3.2. Velocity Profiles

As mentioned, 71 velocity profiles (19 points) were measured in 19 cross–sections in this study. Figure 9 shows the velocity profiles measured in cross–sections number 2 and 4 of the reaches of Dalaki and Shapur2, respectively. In this figure, the unit of velocity is meters per second (m/s). To better investigate the effect of vegetation patches on velocity distribution, more data were collected close to the vegetation patch (See Figure 9A,B). According to Figure 9A, in the profile located on the left bank near the patch of vegetation, the dip phenomenon occurred due to the presence of vegetation and the generation of secondary flows. In this figure, by moving away from the vegetation patch to the right bank, the effect of this phenomenon is decreased, and the velocity distribution illustrates an S–shaped form partly due to variations in the bed sediment. Other researchers have reported S–shaped profiles in other sand and gravel bed streams [34,35].

3.3. Validity of Logarithmic Law

Figure 10 shows the method of determining the validity of the logarithmic law in the profiles of the central axis in the four selected reaches, where the black triangles indicate the fitness of the logarithmic law to the inner layer data and green triangles indicate deviation from the logarithmic law in the outer layer. In the figure, the profiles are named as (reach name-cross-section number-distance from the right bank). For example, Sh1_1_14.2 shows the existing profile in Shapur1 reach and the cross-section of No. 1 at a distance of 14.2 m from the right bank. According to Figure 10, the coefficient of determination values (R2) in all profiles are more than 90%, indicating the reasonable agreement of the logarithmic law with the measured data in the inner layer. Investigation of the logarithmic law in velocity profiles showed that the profiles located on the crest of the bedform (e.g., Sh1-3-16 and F-4-10) are more irregular than the profiles of the concave region of the bedform (e.g., Sh1-1-14.2 and D-2-5.25). Moreover, the coefficient of determination values (R2) related to the logarithmic law fit in these profiles was less than the profiles of the concave region. This could be due to a sudden increase in velocity on the crest of the bedform and the occurrence of the flow separation zone downstream of the crest [36]. In profiles where the outer layer data are below the regression line, the favorable pressure gradient ( p / x < 0 ) is dominant, and if the deviation is above the regression line, the adverse pressure gradient ( p / x > 0 ) is dominant. According to this figure, in most profiles, the deviation of the outer layer data is below the regression line, indicating the dominance of the adverse pressure gradient ( p / x > 0 ) due to increasing turbulence [37].

3.4. Calculation of Shear Velocity

River engineers use indirect methods to estimate the bed shear stress. Three methods of boundary layer (BL) characteristics, logarithmic law, and Darcy–Weisbach (DW) were used to estimate the bed shear in this study. Table 6 shows the bed shear stress values calculated using these three methods. The values of the Darcy–Weisbach friction factor in the table were obtained using the method of Bellos et al. (Equations (15)–(18)). Due to the large Reynolds number in the selected reaches, the values of α and b in Equations (16) and (17) are almost equal to zero, which makes the first and second terms of Equation (15) approximately equal to one. Therefore, the Darcy–Weisbach friction factor depends only on the flow depth and roughness height. For this reason, in Table 6, profiles with the same cross–section and depth (for example, profiles Sh1_2_10 and Sh1_2_15) have an equal friction factor.
Comparing the shear velocity calculated with three different methods in Table 6 shows that:
  • The logarithmic law method reveals the highest range of changes in estimating the shear velocity due to the sensitivity of this method to adjusting the reference level and using the velocity variation near the bed. However, the Darcy–Weisbach method shows the lowest range of changes in shear velocity estimation because it uses all data measured in a velocity profile.
  • The boundary layer characteristics method is reliable in coarse–grained rivers for calculating shear velocity [23] because it considers all velocity data in each profile. It was found that the relative differences between the logarithmic law and Darcy–Weisbach methods compared to the boundary layer characteristics method are equal to 87% and 39%, respectively, which indicates more reasonable agreement between Darcy–Weisbach method and the boundary layer characteristics method.
In this research, the relative difference in the percentage of the logarithmic law and Darcy–Weisbach methods was calculated and compared to the boundary layer characteristics using the following equation:
E = i = 1 n 1 n | u * L o g   o r   D w u * B L u * B L | × 100  
where n is the total number of data.

3.5. Calibration of Three–Parameter Flow Resistance Model

The three–parameter flow resistance model presented in Equation (3) has two independent variables (depth h and flow velocity v ) and three constant parameters ( A ,   B ,   C ). After calculating the shear velocity using the three methods, the constant parameters (A, B, C) of the three–parameter flow resistance model were determined through the best fit between the results obtained from each method and the results of the three–parameter flow resistance model. This was done using the MATLAB software curve fitting toolbox.
In this study, 90% of the profiles (64 profiles) were used for calibration, and 10% (7 profiles) were used to validate the three–parameter flow resistance model. Table 7 shows the calibration results of the mentioned model using three methods of boundary layer characteristics, logarithmic law, and Darcy–Weisbach. According to the table, the Darcy–Weisbach method has the highest coefficient of determination (R2) and the lowest root mean square error (RMSE), suggesting the ability of the Darcy–Weisbach method to calibrate the three–parameter flow resistance model compared to other methods. In addition, the logarithmic law method has the lowest coefficient of determination (R2) and the highest root mean square error (RMSE), revealing the low accuracy of the logarithmic law method for calibrating the three–parameter flow resistance model. The value of coefficient A (velocity power) obtained using Darcy–Weisbach method and the boundary layer characteristics are close to each other due to the use of all points of the velocity profile. Figure 11 shows the fitted surface on the data obtained using the three methods.

3.6. Validation of Three–Parameter Flow Resistance Model

Since it is difficult to measure the bed shear stress in field research directly, it is impractical to calculate the model accuracy by comparing the bed shear stress calculated with the measured one. Therefore, to validate and calculate the model accuracy in this study, the amount of flow velocity calculated by the model (Equation (4)) was compared with the flow velocity measured by the current meter. Since the Darcy–Weisbach method has the highest accuracy in model calibration, only the calibration results using this method were used to validate the model. As mentioned, 10% of the velocity profiles (7 profiles) were used to validate and calculate the model accuracy. In Figure 12, the velocity value calculated using the three–parameter flow resistance model is compared with the velocity values measured by the current meter. According to this figure, the three–parameter flow resistance model in 6 of 7 profiles (86%) estimates values higher than the measured values. The average error of the model was calculated to be 17%, which indicates the high accuracy of the model in estimating the velocity and consequently the bed shear stress. It should be noted that up to 200% error is acceptable in estimating certain fluvial parameters [38]. Reasons for accepting this amount of error in river engineering are:
  • Difficulty working in the river with rough bed and non–uniform flow;
  • Difficulty working with high–precision measuring instruments in field research;
  • Not recognizing all the factors affecting error production;
  • Assumptions and simplifications performed that generate errors.

4. Conclusions

Resistance to flow is a key parameter in river engineering studies and its estimation plays a significant role in engineering models. Accordingly, in this study, resistance to flow was estimated in four reaches of coarse–bed rivers with vegetation patches using three methods of boundary layer characteristics, logarithmic law, and Darcy–Weisbach. In addition, the coefficients of the three–parameter flow resistance model (A, B, C) were calibrated based on these three methods using 90% of the data. In summary, the following results can be presented:
  • The three–parameter flow resistance model shows that the average error percentage of the model is 17%, indicating the accuracy of the model.
  • The logarithmic law method has the highest range of changes in estimating the bed shear stress compared to the methods of Darcy–Weisbach and the boundary layer characteristics method. This is due to the sensitivity of the logarithmic law method to adjusting the reference surface and using the near bed data of the velocity profile. The Darcy–Weisbach method reveals the lowest range of changes.
  • The relative difference percentages between the logarithmic law and the Darcy–Weisbach methods compared to the boundary layer characteristics method were equal to 87% and 39%, respectively. This indicated a more reasonable agreement between the Darcy–Weisbach method and the boundary layer characteristics method.
  • By investigating 71 measured velocity profiles, it was found that the logarithmic law was well applicable in reaches with vegetation patches. A high coefficient of determination ( R 2 > 0.9 ) in fitting the velocity profile data to the logarithmic law indicated their reasonable agreement with this law.
  • The investigation of measured velocity profiles shows the occurrence of the Dip phenomenon in the velocity profiles near the vegetation patches. However, by moving away from the vegetation patches, the effect of this phenomenon is decreased, and the profiles illustrate an S–shaped distribution. The Dip location plays a significant role in the boundary layer thickness and estimation of key hydraulic parameters.
The results of this study can be used in many environmental models, including fluvial geomorphology, and mountain rivers with vegetation patches in order to improve the flow resistance evaluation and its application in hydraulic models.

Author Contributions

Conceptualization, H.A.; methodology, H.A. and A.D.; software, M.N.; validation, H.A., A.D. and N.D.; investigation, M.N.; resources, H.A. and A.D.; data curation, M.N.; writing—original draft preparation, M.N. and H.A.; writing—review and editing, M.N.-S. and M.K.; visualization, M.N.; supervision, H.A., A.D., N.D., M.N.-S. and M.K.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, J.; Shi, F.-Y.; Chen, P.-P.; Wu, P.; Sui, J. Impact of bridge pier on the stability of ice jam. J. Hydrodyn. 2015, 27, 865–871. [Google Scholar] [CrossRef]
  2. Folkard, A.M. Vegetated flows in their environmental context: A review. Proc. Inst. Civ. Eng. Eng. Comput. Mech. 2011, 164, 3–24. [Google Scholar] [CrossRef]
  3. Marion, A.; Nikora, V.; Puijalon, S.; Bouma, T.; Koll, K.; Ballio, F.; Tait, S.; Zaramella, M.; Sukhodolov, A.; O’Hare, M. Aquatic interfaces: A hydrodynamic and ecological perspective. J. Hydraul. Res. 2014, 52, 744–758. [Google Scholar] [CrossRef]
  4. Kothyari, U.C.; Hayashi, K.; Hashimoto, H. Drag coefficient of unsubmerged rigid vegetation stems in open channel flows. J. Hydraul. Res. 2009, 47, 691–699. [Google Scholar] [CrossRef]
  5. Palmer, M.R.; Nepf, H.M.; Pettersson, T.J.; Ackerman, J.D. Observations of particle capture on a cylindrical collector: Implications for particle accumulation and removal in aquatic systems. Limnol. Oceanogr. 2004, 49, 76–85. [Google Scholar] [CrossRef] [Green Version]
  6. López, F.; García, M. Open-channel flow through simulated vegetation: Suspended sediment transport modeling. Water Resour. Res. 1998, 34, 2341–2352. [Google Scholar] [CrossRef]
  7. Cotton, J.; Wharton, G.; Bass, J.; Heppell, C.; Wotton, R. The effects of seasonal changes to in-stream vegetation cover on patterns of flow and accumulation of sediment. Geomorphology 2006, 77, 320–334. [Google Scholar] [CrossRef]
  8. Setayesh, P.; Afzalimehr, H. Effect of Reedy Emergent Side-Vegetation in Gravel-Bed Streams on Bed Shear Stress: Patch Scale Analysis. Iran. J. Sci. Technol. Trans. Civ. Eng. 2021, 46, 1375–1392. [Google Scholar] [CrossRef]
  9. 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]
  10. Liu, X.-D.; Tang, L.-C.; Han, Y.; Chen, J.; Yang, S.-Q. Experimental study on flow resistance over rigid vegetated channel. IEEE Access 2019, 7, 93974–93985. [Google Scholar] [CrossRef]
  11. Kazem, M.; Afzalimehr, H.; Sui, J. Formation of Coherent Flow Structures beyond Vegetation Patches in Channel. Water 2021, 13, 2812. [Google Scholar] [CrossRef]
  12. Green, J.C. Comparison of blockage factors in modelling the resistance of channels containing submerged macrophytes. River Res. Appl. 2005, 21, 671–686. [Google Scholar] [CrossRef]
  13. Thoman, R.W.; Niezgoda, S.L. Determining erodibility, critical shear stress, and allowable discharge estimates for cohesive channels: Case study in the Powder River basin of Wyoming. J. Hydraul. Eng. 2008, 134, 1677–1687. [Google Scholar] [CrossRef]
  14. Box, W.; Järvelä, J.; Västilä, K. Flow resistance of floodplain vegetation mixtures for modelling river flows. J. Hydrol. 2021, 601, 126593. [Google Scholar] [CrossRef]
  15. Van Dijk, W.; Teske, R.; Van de Lageweg, W.; Kleinhans, M. Effects of vegetation distribution on experimental river channel dynamics. Water Resour. Res. 2013, 49, 7558–7574. [Google Scholar] [CrossRef] [Green Version]
  16. Shields, F.D.; Coulton, K.G.; Nepf, H. Representation of Vegetation in Two-Dimensional Hydrodynamic Models. J. Hydraul. Eng. 2017, 143, 8. [Google Scholar] [CrossRef] [Green Version]
  17. Nosrati, K.; Afzalimehr, H.; Sui, J. Drag Coefficient of Submerged Flexible Vegetation Patches in Gravel Bed Rivers. Water 2022, 14, 743. [Google Scholar] [CrossRef]
  18. Afzalimehr, H.; Barahimi, M.; Sui, J. Non-uniform flow over cobble bed with submerged vegetation strip. Proc. Inst. Civ. Eng. Water Manag. 2019, 172, 86–101. [Google Scholar] [CrossRef] [Green Version]
  19. Ferguson, R. Time to abandon the Manning equation? Earth Surf. Processes Landf. 2010, 35, 1873–1876. [Google Scholar] [CrossRef]
  20. Afzalimehr, H.; Levesque, B.; Molls, T.; Zhao, G.; Molls, F. Friction Slope in Depth-Averaged Flow. J. Hydraul. Eng. 1999, 125, 549. [Google Scholar] [CrossRef]
  21. Bellos, V.; Nalbantis, I.; Tsakiris, G. Friction modeling of flood flow simulations. J. Hydraul. Eng. 2018, 144, 04018073. [Google Scholar] [CrossRef] [Green Version]
  22. Afzalimehr, H.; Anctil, F. Accelerating shear velocity in gravel-bed channels. Hydrol. Sci. J. 2000, 45, 113–124. [Google Scholar] [CrossRef] [Green Version]
  23. Afzalimehr, H.; Rennie, C.D. Determination of bed shear stress in gravel-bed rivers using boundary-layer parameters. Hydrol. Sci. J. 2009, 54, 147–159. [Google Scholar] [CrossRef] [Green Version]
  24. Afzalimehr, H.; Maddahi, M.R.; Naziri, D.; Sui, J. Effects of non-submerged boulder on flow characteristics—A field investigation. Int. J. Sediment Res. 2019, 34, 136–143. [Google Scholar] [CrossRef]
  25. Afzalimehr, H.; Maddahi, M.R.; Sui, J.; Rahimpour, M. Impacts of vegetation over bedforms on flow characteristics in gravel-bed rivers. J. Hydrodyn. 2019, 31, 986–998. [Google Scholar] [CrossRef]
  26. Emadzadeh, A.; Chiew, Y.M.; Afzalimehr, H. Effect of accelerating and decelerating flows on incipient motion in sand bed streams. Adv. Water Resour. 2010, 33, 1094–1104. [Google Scholar] [CrossRef]
  27. Fazlollahi, A.; Afzalimehr, H.; Haghighi, A.; Barati, R.; Neyshabouri, S.A.A.S.; Ahmadi, G.; Omidvar, P.; Kheirkhahan, M.; Hosseini, K.; Kabiri-Samani, A. 101. Validation of Spatially Averaging Method for Using the Law of the Wall over Concave Bed Form (pool). SMR 2007, 2, 27–42. [Google Scholar]
  28. Naderi, M.; Afzalimehr, H.; Sohrabi, S. Investigation of three-parameter flow resistance model in coarse-bed rivers (Case study: Deryuk River). Iran. J. Soil Water Res. 2021, 52, 1423–1435. [Google Scholar]
  29. Song, T.; Chiew, Y. Turbulence measurement in nonuniform open-channel flow using acoustic Doppler velocimeter (ADV). J. Eng. Mech. 2001, 127, 219–232. [Google Scholar] [CrossRef]
  30. Alonso, R.L.; Fernández, J.B.; Cugat, M.C. Flow resistance equations for mountain rivers. For. Syst. 2009, 18, 81–91. [Google Scholar]
  31. Brierley, G.J.; Fryirs, K.A. Geomorphology and River Management: Applications of the River Styles Framework; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  32. Wolman, M.G. A method of sampling coarse river-bed material. EOS Trans. Am. Geophys. Union 1954, 35, 951–956. [Google Scholar] [CrossRef]
  33. Okhravi, S.; Gohari, S. Form friction factor of armored riverbeds. Can. J. Civ. Eng. 2020, 47, 1238–1248. [Google Scholar] [CrossRef]
  34. Afzalimehr, H.; Ramazany, R. Validation of Double-Averaged Velocity Method for Decelerating Flow in Coarse-Bed Rivers. J. Hydraul. 2019, 2, 47–67. [Google Scholar]
  35. Afzalimehr, H.; Subhasish, D. Influence of bank vegetation and gravel bed on velocity and Reynolds stress distributions. Int. J. Sediment Res. 2009, 24, 236–246. [Google Scholar] [CrossRef]
  36. Maddahi, M.R.; Afzalimehr, H.; Rowinski, P.M. Flow characteristics over a gravel bedform: Kaj River case study. Acta Geophys. 2016, 64, 1779–1796. [Google Scholar] [CrossRef] [Green Version]
  37. Houra, T.; Tsuji, T.; Nagano, Y. Effects of adverse pressure gradient on quasi-coherent structures in turbulent boundary layer. Int. J. Heat Fluid Flow 2000, 21, 304–311. [Google Scholar] [CrossRef]
  38. Wu, B.; Molinas, A.; Julien, P.Y. Bed-material load computations for nonuniform sediments. J. Hydraul. Eng. 2004, 130, 1002–1012. [Google Scholar] [CrossRef]
Figure 1. The geographical location of the studied reaches. Blue frames show study areas.
Figure 1. The geographical location of the studied reaches. Blue frames show study areas.
Fluids 07 00284 g001
Figure 2. Some vegetation patches in selected reaches.
Figure 2. Some vegetation patches in selected reaches.
Fluids 07 00284 g002
Figure 3. Velocity measurement at selected reaches.
Figure 3. Velocity measurement at selected reaches.
Fluids 07 00284 g003
Figure 4. Surveying in the studied reaches.
Figure 4. Surveying in the studied reaches.
Fluids 07 00284 g004
Figure 5. Plan of the studied reaches along with the position of the measured profiles.
Figure 5. Plan of the studied reaches along with the position of the measured profiles.
Fluids 07 00284 g005aFluids 07 00284 g005b
Figure 6. (A) Sampling of riverbed sediment, (B) measurement of riverbed sediments.
Figure 6. (A) Sampling of riverbed sediment, (B) measurement of riverbed sediments.
Fluids 07 00284 g006
Figure 7. Grain size distribution in the studied reaches.
Figure 7. Grain size distribution in the studied reaches.
Fluids 07 00284 g007
Figure 8. The general framework of this research.
Figure 8. The general framework of this research.
Fluids 07 00284 g008
Figure 9. Cross−section and velocity profiles measured in (A) Section No. 1 from Shapur1; (B) Section No. 2 from Fahlyan reach.
Figure 9. Cross−section and velocity profiles measured in (A) Section No. 1 from Shapur1; (B) Section No. 2 from Fahlyan reach.
Fluids 07 00284 g009
Figure 10. Determining the validity of the logarithmic law in the profiles in the central axis from four selected reaches.
Figure 10. Determining the validity of the logarithmic law in the profiles in the central axis from four selected reaches.
Fluids 07 00284 g010aFluids 07 00284 g010b
Figure 11. Surface fitted on the data obtained using the methods (A) Darcy–Weisbach, (B) boundary layer characteristics, (C) logarithmic law.
Figure 11. Surface fitted on the data obtained using the methods (A) Darcy–Weisbach, (B) boundary layer characteristics, (C) logarithmic law.
Fluids 07 00284 g011
Figure 12. Comparison between the velocity calculated from the model with the measured velocity.
Figure 12. Comparison between the velocity calculated from the model with the measured velocity.
Fluids 07 00284 g012
Table 1. Summary of general information of selected reaches.
Table 1. Summary of general information of selected reaches.
Q
(m3/s)
S0
(%)
D
(m)
W
(m)
Lr
(m)
LatitudeLongitudeReach Name
3.250.250.3430.7443 51 ° 26   N 29 ° 35   E Shapour1
3.000.210.4024.5830 51 ° 26   N 29 ° 35   E Shapour2
0.860.160.3017.5635 51 ° 27   N 30 ° 14   E Fahlyan
1.390.250.2410.1528 51 ° 26   N 29 ° 16   E Dalaki
Table 2. Vegetation patches’ characteristics of the reaches.
Table 2. Vegetation patches’ characteristics of the reaches.
B V B A B X h p h Reach Name
0.0620.1180.0890.59Shapur1
0.0350.0880.1080.37Shapur2
0.0540.0390.0820.48Fahlyan
0.0830.1450.1350.39Dalaki
Table 3. The number of velocity measuring points in the four studied reaches.
Table 3. The number of velocity measuring points in the four studied reaches.
FrNP + NPTNTPNPProfiles Measured NumberCross–Section NumberReach Name
2630698303395205Shapur1
1902516194322164Shapur2
1992560175385205Fahlyan
166040098302155Dalaki
8184217477014047119Sum
Table 4. The number of surveyed points in studied reaches.
Table 4. The number of surveyed points in studied reaches.
DV
(Point/m2)
DB
(Point/m2)
NVNBReach Name
1223582745Shapur1
1222861554Shapur2
1221951302Fahlyan
122228604Dalaki
--10676205Sum
Table 5. Sediment and hydraulic parameters of reaches.
Table 5. Sediment and hydraulic parameters of reaches.
σ g d84
(mm)
d16
(mm)
FrReA
(m2)
D
(m)
Q
(m3/s)
Section NumberReach
1.5949.0019.330.15484,52411.530.4063.441Shpur1
1.4447.0022.670.14425,09611.100.3833.082
1.6752.0018.670.17405,30810.520.3293.243
1.7346.6715.670.24456,3429.360.2843.764
1.7642.3313.670.20370,6118.670.2772.905
1.6653.3319.330.16431,8469.370.3602.811Shpur2
1.5850.6720.330.17512,7779.680.3833.242
1.5943.0017.000.15486,6679.840.4103.923
1.4547.6722.670.15542,12310.060.4373.124
1.4146.6723.330.13225,3394.870.2950.931Fahlyan
1.4238.3319.000.12204,3185.280.3100.872
1.3937.3319.330.14176,8235.540.3100.793
1.5444.0018.670.12186,7345.560.3090.844
1.5241.0017.670.13192,1845.220.2850.885
1.4735.3316.330.26525,5463.130.2981.381Dalaki
1.7243.3314.670.36568,7022.620.2501.492
1.5942.0016.670.48576,1041.930.2091.333
1.5742.0017.000.29444,0002.490.2491.374
1.6846.0016.330.47529,3972.090.1991.395
Table 6. Shear velocity and shear stress of the substrate calculated by different methods.
Table 6. Shear velocity and shear stress of the substrate calculated by different methods.
u * D W
( m s )
τ D W
( N m 2 )
f U
( m s )
τ B L
( N m 2 )
u * B L
( m s )
θ
( c m )
δ *
( c m )
τ L o g
τ L o g
u * L o g
( m s )
α Profile Name
0.0290.840.1220.2351.140.0344.618.472.610.0510.1278Sh1_1_4.5
0.0391.520.0970.3552.010.0455.8610.262.350.0480.1212Sh1_1_14.2
0.0310.990.0810.3131.190.0356.4710.842.190.0470.1169Sh1_1_18.2
0.0391.510.0820.3841.590.0403.996.744.760.0690.1724Sh1_1_24.9
0.0250.630.0880.2401.370.0374.699.121.180.0340.0859Sh1_2_10
0.0351.220.0880.3341.740.0426.5811.511.510.0390.0971Sh1_2_15
0.0381.460.0930.3552.740.0525.7911.509.960.0320.0789Sh1_2_19.5
0.0502.510.1110.4262.180.0473.445.962.020.0450.1124Sh1_2_26
0.0351.220.01280.2761.570.0403.667.003.650.0610.1511Sh1_3_11
0.0462.070.1070.3932.330.0485.228.943.030.0550.1375Sh1_3_16
0.0411.710.1180.3402.270.0484.458.274.170.0650.1615Sh1_3_18.4
0.0492.400.0950.4513.970.0635.129.945.560.0740.1848Sh1_3_26
0.0371.350.1450.2733.440.0593.067.009.90.0990.2488Sh1_4_10.5
0.0351.230.1650.2442.760.0533.537.900.820.0290.0714Sh1_4_14.5
0.0492.350.1230.3922.530.0502.995.554.500.0670.1681Sh1_4_16.5
0.0573.280.1040.5035.710.0765.4510.593.100.0560.1392Sh1_4_25
0.0361.300.1410.2722.200.0473.507.062.660.0520.1290Sh1_5_9.3
0.0512.600.1180.4203.210.0573.566.545.200.0720.1804Sh1_5_15.65
0.0623.810.1000.5515.700.0764.378.238.730.0930.2336Sh1_5_22.3
0.0210.450.1180.1750.710.0265.559.821.470.0380.0959Sh1_5_27.3
0.0240.570.0910.2242.840.0537.9520.420.2250.0150.0375Sh2_1_8
0.0441.950.0780.4472.100.0467.6312.082.360.0490.1217Sh2_1_13
0.0401.600.0730.4192.330.0489.9216.472.040.0450.1130Sh2_1_15
0.0240.560.0910.2232.810.0538.0720.570.280.0170.0419Sh2_1_18
0.0250.630.0830.2461.090.0336.0911.441.890.0430.1087Sh2_2_8
0.0341.130.0760.3461.640.0419.8516.183.150.0560.1404Sh2_2_12.65
0.0381.460.0730.4002.340.04810.2617.231.320.0360.0909Sh2_2_15.4
0.0180.3150.1000.1591.010.0325.4613.110.330.0180.0459Sh2_2_18.4
0.0240.590.0800.2431.080.0336.2511.570.420.0200.0510Sh2_3_7
0.0401.610.0680.4352.930.0547.9314.143.520.0590.1484Sh2_3_12
0.0381.440.0680.4122.190.0478.2513.841.820.0430.1067Sh2_3_14
0.0250.600.0740.2561.340.0378.6016.000.460.0210.0535Sh2_3_17
0.0280.780.0790.2810.810.0286.7510.610.680.0260.0654Sh2_4_6
0.0391.490.0740.4031.840.0436.6110.882.320.0480.1206Sh2_4_10.5
0.0351.240.0710.3742.150.04610.1717.362.250.0470.1185Sh2_4_14
0.0270.720.0840.2621.600.0407.5614.541.060.0330.0814Sh2_4_17
0.0120.150.1230.1000.280.0173.777.751.260.0360.0889F_1_4.5
0.0230.540.1290.1831.260.0363.758.281.220.0350.0874F_1_8.25
0.0351.260.0930.3292.060.0455.8110.571.370.0370.0925F_1_11.5
0.0280.790.0870.2721.470.0388.5714.892.390.0490.1222F_1_13.5
0.0140.210.1290.1130.490.0222.956.690.810.0280.0710F_2_6
0.0240.580.1290.1900.980.0313.146.152.720.0520.1304F_2_8.5
0.0290.810.0870.2732.170.0476.5512.963.240.0570.1423F_2_11
0.0260.650.0810.2540.980.0315.088.931.050.0320.0811F_2_14
0.0100.110.1030.0920.250.0163.467.560.120.0110.0271F_3_5
0.0140.200.1130.1200.740.0273.158.540.340.0180.0460F_3_9
0.0250.620.0770.2551.500.0398.9414.731.400.0370.0934F_3_13
0.0190.0350.0790.1880.580.0247.4812.632.650.0520.1288F_3_15
0.0120.150.1030.1090.310.0184.068.500.670.0260.0646F_4_8
0.0160.260.0990.1440.440.0214.648.832.270.0480.1190F_4_10
0.0260.650.0800.2471.200.0355.8711.021.530.0390.0979F_4_13
0.0220.490.0900.2090.720.0275.8410.191.830.0430.1069F_4_15
0.0150.210.1210.1190.520.0223.137.360.760.0270.0687F_5_5
0.0170.280.1020.1480.700.0264.049.092.220.0470.1177F_5_9.15
0.0280.780.0910.2611.090.0334.007.271.470.0380.0959F_5_12
0.0220.480.0890.2080.800.0284.598.652.250.0480.1188F_5_15
0.0714.990.0770.7205.910.0774.907.998.610.0930.2320D_1_3
0.0462.100.0820.4542.010.0453.605.741.760.0420.1051D_1_5.25
0.0310.960.0870.2962.490.0507.1513.283.280.0570.1432D_1_6
0.0867.370.1150.71710.720.1042.725.5418.000.1340.3362D_2_3
0.0684.600.0970.6166.540.0813.887.177.000.0840.2093D_2_5.25
0.0623.840.1000.5544.040.0643.285.643.560.0600.1492D_2_7.5
0.0694.730.1180.5667.970.0893.577.0212.340.1110.2778D_3_1.55
0.0704.940.1280.5566.560.0813.045.7045.270.2130.5319D_3_4.6
0.11011.990.1280.86615.590.1252.815.349.100.0950.2385D_3_6.05
0.0613.760.0910.5747.220.0857.1112.598.470.0920.2301D_4_1.5
0.0908.010.1320.69711.240.1062.625.1924.940.1580.3948D_4_5
0.0512.570.1240.40710.490.1023.599.5826.180.1620.4045D_4_8
0.0999.770.0910.92718.260.1354.489.0416.010.1270.3163D_5_2
0.13317.610.2110.81820.150.1421.513.3671.180.2670.6670D_5_5.25
0.0765.790.1360.58513.010.1143.247.3410.440.1020.2554D_5_9.5
Table 7. Calibration results of three–parameter resistance model using three different methods.
Table 7. Calibration results of three–parameter resistance model using three different methods.
R M S E R 2 ModelCBAMethod
0.16080.9929 τ 0 ρ = v 1.937 141.80 h 0.5131 0.5131141.801.937Darcy–Weisbach
1.31100.8461 τ 0 ρ = v 1.955 94.33 h 0.5977 0.597794.331.955boundary layer characteristics
4.5870.5938 τ 0 ρ = v 1.625 251.50 h 1.638 1.638251.501.625logarithmic law
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Naderi, M.; Afzalimehr, H.; Dehghan, A.; Darban, N.; Nazari-Sharabian, M.; Karakouzian, M. Field Study of Three–Parameter Flow Resistance Model in Rivers with Vegetation Patch. Fluids 2022, 7, 284. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids7080284

AMA Style

Naderi M, Afzalimehr H, Dehghan A, Darban N, Nazari-Sharabian M, Karakouzian M. Field Study of Three–Parameter Flow Resistance Model in Rivers with Vegetation Patch. Fluids. 2022; 7(8):284. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids7080284

Chicago/Turabian Style

Naderi, Masoud, Hossein Afzalimehr, Ayoub Dehghan, Nader Darban, Mohammad Nazari-Sharabian, and Moses Karakouzian. 2022. "Field Study of Three–Parameter Flow Resistance Model in Rivers with Vegetation Patch" Fluids 7, no. 8: 284. https://0-doi-org.brum.beds.ac.uk/10.3390/fluids7080284

Article Metrics

Back to TopTop