Evaluating Two-Layer Models for Velocity Profiles in Open-Channels with Submerged Vegetation

For submerged vegetated flow, the velocity profile has two distinctive distributions in the vegetation layer in the lower region and the surface layer in the upper non-vegetated region. Based on a mixing-layer analogy, different analytical models have been proposed for the velocity profile in the two layers. This paper evaluates the four analytical models of Klopstra et al., Defina & Bixio, Yang et al. and Nepf against a wide range of independent experimental data available in the literature. To test the applicability and robust of the models, the author used the 19 datasets with various relative depths of submergence, different vegetation densities and bed slopes (1.8 × 10 4.0 × 10). This study shows that none of the models can predict the velocity profiles well for all datasets. The three models except Yang’s model performed reasonably well in certain cases, but Yang’s model failed in most the cases studied. It was also found that the Defina model is almost the same as the Klopstra model, if the same mixing length scale of eddies (λ) is used. Finally, close examination of the mixing length scale of eddies (λ) in the Defina model showed that when λ/h = 1/40(H/h), this model can predict velocity profiles well for all the datasets used.


Introduction
Vegetation exists in many natural rivers and woodlands.No matter what type the vegetation is, it will alter the velocity field of flow, consequently increasing X.Tang DOI: 10.4236/gep.2019.7100669 Journal of Geoscience and Environment Protection figflow resistance.The influence of vegetation on the velocity of flow depends on its type (rigid or flexible) and condition (submerged or emergent).As a prerequisite for analysis of flow resistance, pollutant mixing process, etc., the velocity profile of vegetated flows has drawn the interest of many researchers.For example, Temple (1986) studied velocity distribution in channels with grass; turbulent structure of vegetated layer was studied (Shimzu & Tsujimoto, 1994;Nepf & Koch, 1999;Nepf & Vivoni, 2000;Lopez & Garcia, 2001;Dimitris & Panayotis, 2011).There are extensive studies on the velocity distribution via experimental and analytical methods (Tsujimoto & Kitamur 1990, Klopstra et al., 1997;Meijer & Van Velzen, 1999;Nepf & Koch, 1999;Ghisalberti & Nepf, 2004;Defina & Bixio, 2005;Baptist et al., 2007;Kubrak, et al. 2008;Huai et al. 2009;Yang & Choi, 2010;Nepf, 2012;Tang, 2018).However, due to different types of vegetation and flow conditions, various methods are proposed for predicting velocity profiles.For the flows with submerged vegetation, the most common method is a two-layer approach, in which different analytical models are applied in the lower vegetation layer and the upper surface layer, using a mixing-layer analogy (Klopstra et al., 1997;Meijer & Van Velzen, 1999;Defina & Bixio, 2005;Baptist, et al., 2007;Huai, et al. 2009;Yang & Choi, 2010;Nepf, 2012;Singh et al. 2019;Tang, 2018).Klopstra et al. (1997) proposed a two-layer model of velocity: one in the vegetation layer and one above it called the surface layer.In the vegetation layer, the turbulent eddy viscosity is expressed by the product of a characteristic length (λ) and velocity, following Boussinesq hypothesis.λ was found to empirically relate to H/h (H is the flow depth and h denotes the height of vegetation).After further investigation on this model through more data, Meijer & Van Velzen (1999) found that λ could be approximated by 0.0144 Hh .Similarly, Defina & Bixio (2005) estab- lished a different form of analytical solution for velocity in the vegetation layer.Yang & Choi (2010) and Nepf (2012) proposed different empirical models for velocity in the vegetated layer.
Considering only certain data used in the validation of the above-mentioned models, the present study aims at comparing the four models of Klopstra et al. (1997), Defina & Bixio (2005), Yang & Choi (2010) and Nepf (2012) for the prediction of velocity profiles in submerged rigid vegetation against a wide range of independent datasets.Thus it is to demonstrate the capability and robust of each model, consequently making some recommendation for the future application.In this study, the 19 datasets used cover different submergence ratios (H/h) ranging from 1.25 to 3.4, different vegetation densities (defined as a, the frontal area of the vegetation per unit volume) (a = 1.1 -10 m −1 ) and bed slopes S o (1.8 × 10 −6 -4.0 × 10 −3 ).

Description of Models
In open-channel flow with vegetation, the bed and wall boundary stress are both negligible compared with the drag force on the vegetation (Nepf & Vivoni, 2000;Stone & Shen, 2002).Thus, the momentum equation of steady fully-developed X. Tang 1-D vegetated flow may be described as (Klopstra et al., 1997;Ghisalberti & Nepf, 2004;Defina & Bixio, 2005;Kubrak et al., 2008;Tang, 2018): where τ is the shear stress, g the gravity, z the vertical coordinate above the bed, S o the bed slope and F v is the drag force per unit mass generated by the vegetation, see Figure 1.The drag force F v is given by: where h is the height of vegetation, u the local stream wise velocity, C D the drag coefficient, a the frontal area of vegetation (A v ) per unit volume, representing the density of vegetation, and m is the number of vegetation per unit area.
Using Boussinesq's hypothesis on eddy viscosity combined with a mixing-length concept (Defina & Bixio, 2005;Baptist et al. 2007), the Reynolds stress (τ) is described as: where ν T is the total eddy viscosity of vegetated flow, and λ is a mixing length.
Under steady flow condition, inserting Equations ( 2) & (3) into (1) gives: ( ) For given a and C D , one can obtain an analytical solution for u 2 in Equation (4) under appropriate boundary conditions (e.g.Klopstra et al., 1997;Defina & Bixio, 2005).The boundary conditions are as follows: At the bed (z = 0), where the bed shear stress is neglected compared with drag force of vegetation, the local equilibrium between vegetation drag and gravity force leads to: At the top of the vegetation (z = h), the Reynolds stress is described as: Figure 1.Velocity profile in a channel with submerged vegetation.
X. Tang

Model by Klopstra et al. (1997)
For the vegetation layer, Equation (4) has the following analytical solution for the velocity under boundary conditions ( 5) & ( 6): 0.0793 ln 0.0009 when 0.001 For the surface layer, the velocity is described by the well-known logarithmic profile: ( ) where κ is von Karman's constant (0.40), z m (=h − h s ) is the zero-plane displacement of the logarithmic profile, h s is the distance between the top of vegetation and the virtual bed of the surface layer (see Figure 1), z o is the equivalent height of bed roughness, and u * is the shear velocity, given by: ( ) where

Model by Defina and Bixio (2005)
Based on Equation (4) with similar assumptions, and using the mixing length concept, Defina & Bixio (2005) proposed another form of velocity profile by introducing a parameter β, which depends on λ/h.The velocity profile is described as follows: For the vegetation layer: where λ is recommended to be 0.0144 Hh .
For the surface layer, the velocity is described by the same Equation ( 11), where h s and z o are respectively given by:

Model by Yang and Choi (2010)
Yang and Choi (2010) proposed a two-layer model that predicts the velocity profile for flows with both flexible and rigid submerged vegetation.The flow through the vegetation is assumed to be driven by gravity and the shear force generated at the top of the vegetation.This is a result of the co-flowing between the low speed vegetation layer and relatively higher speed surface layer, together with the drag force of vegetation.Thus, the flow velocity in the vegetation layer is given by: The flow in the upper non-vegetated layer is assumed to be driven by gravity force and Reynolds stresses, leading to a logarithmic profile as: where c u is recommended as 2 when a > 5; otherwise it is 1.

Model by Nepf (2012)
Nepf ( 2012) proposed a model using an empirical equation to describe the velocity profile in the vegetation layer.Although the model is supposed to be applicable to both flexible and rigid vegetation, most of its parameters are only proposed for rigid vegetation.Therefore, only the parameters for the rigid vegetation are introduced here.For the surface layer, the velocity profile is assumed as Equation ( 11), where z m and z 0 are given by: where z 0 is recommended for ah > 0.1; otherwise, z 0 is estimated as 0.1h (Nepf & Ghisalberti, 2008).

X. Tang
The vegetation layer is divided into two zones, which depends on the penetration of the turbulence stresses from the surface layer into the vegetation layer.An empirical Equation ( 24) is proposed for estimating the depth of turbulent shear penetration (δ) for (ah) values of (0.2 -3): 0.23 0.6 Therefore in the upper zone of the vegetation layer (h − δ < z < h), the flow is driven by gravity and turbulent stresses and balanced by vegetation drag.The velocity profile is assumed to be exponential: where u h = u at the top of the vegetation, and u 2 = velocity in the lower region of vegetation (z < h − δ), which is described by Equation (5).

Data for Study
To compare the four models in Section 2, we used a wide range of different experimental data from the literature for submerged rigid vegetation.summarizes a total of 19 datasets used in this study.These datasets include different submergence ratio, i.e.H/h = 1.25 -3.4,various vegetation densities a = 1.1 -10 m −1 and bed slopes S o between 1.8 × 10 −6 and 4.0 × 10 −3 .Note that the data of Dunn et al. (1996) and runs 34 and 36 of Meijer & Van Velzen (1999) are extracted from Dimitris & Panayotis (2011).The values of C D in the dataset of emergent vegetation were estimated for all runs, taken from the original papers.

Comparison of Results between the Models
The comparisons between the four models are shown in Figures 2-4, which demonstrate that all the models, except the Yang model, can predict profiles reasonably well for the lower zone of vegetation.However, none of them can predict the velocity profile well in the surface layer for all the data tested.For the    Meijer & Van Velzen (1999) as the same value of λ was used in the Defina model.It appears that both models performs reasonably well when H/h is not very large (e.g.H/h < 2.0).
Despite the Klopstra and Defina models having similar assumptions, the Defina model performs better than the Klopstra model in most cases.The reason appears from the different λ models used.Close examination shows that if the same λ model is used, the predicted velocity profiles are almost the same although z 0 and z m proposed by the two models are different.Note that the Klopstra model did not consider the effect of the 1st term on the right side of Equation ( 18).
The Nepf model does not perform well in some cases.This is mainly due to the insufficiency of the formulas of velocity profile in the vegetation layer.The proposed formula for calculating the velocity in the upper zone of vegetation depends on two velocity predictions: one (u 2 ) at the top of vegetation and the other (u o ) in the lower part of the vegetation.Therefore, if these two velocities deviate from the measured ones significantly then the formula for the upper zone of vegetation will underperform.The proposed formula for the lower zone of vegetation neglects the dispersive, turbulent and viscous stresses and, as mentioned before, in the above experiments the vegetation is deeply penetrated by the turbulent stress and in such cases the dispersive stress might become significant.Bearing in mind that this model is based on the empirical equations for estimation of z m and z o , which is dependent on C D and vegetation density (a), further study is needed about the empirical Equations ( 22) & ( 23).
The predicted velocity by the Yang model significantly deviates from the measured ones.It can be argued that assuming a uniform distribution of velocity profile over the whole vegetation depth is not a realistic assumption.The measured velocity profile in the vegetation in most of the cases tested herein is exponential for a large part of the vegetation depth.Moreover, in some cases illustrated here, this model has predicted the velocity very close to the measured one but the trend is not well captured because the assumption of a uniform velocity does not reflect reality.The inaccuracy of the logarithmic part might be explained from the following two aspects.First, the equation for the surface layer depends on the predicted velocity in the vegetation layer.Therefore, the accuracy of the model in the surface layer is partly dependent on the precision of prediction in the vegetation layer.Second, examining the measured data indicates that the penetration of turbulent stresses into the vegetation layer is significant, so the effective roughness should be much smaller than the vegetation height, which is taken as the effective roughness in this model.Therefore, in all cases the velocity in the upper region is under-estimated.

Further Discussion on Klopstra & Defina Models
Further examination of the Defina and Klopstra models shows that their differences are very small if the same λ value is used.Because both z 0 and h s are im-X.Tang portant to describe the velocity in the surface layer, they relate to the λ value in the models.Since both flow depth and vegetation height will have impact on λ (characteristic length of eddy), λ may be represented by k Hh λ = , where k is a constant.The predicted velocity by both models decreases as increasing k value, because the larger the k value, the bigger λ becomes, i.e. the stronger the eddy, indicating a smaller velocity.Close studies on all the data tested suggest that when k = 1/40 both models have good agreement with the data, as shown in Figures 5-8.

Conclusion
In a river with submerged vegetation, the velocity profile in the lower vegetation layer is significantly different from the upper non-vegetation surface layer.
Comparison has been made on four commonly used two-layer models of velocity profile: Klopstra et al. (1997), Defina & Bixio (2005), Yang & Choi (2010) and Nepf (2012).This study has demonstrated that all the models except the Yang model can predict the velocity reasonably well in the vegetation layer near the bed.However, in the surface layer none of the models using their default parameter    Close examination of λ parameters in the Defina and Klopstra models shows that the two models are very close when the same parameter λ value is used.This study shows that when λ is described by k Hh with k having an optimum value of 1/40, the Defina and Klopstra models can both predict the velocity profiles well for a wide range of vegetated flows.More data are needed for examining the new recommended value of k in future research.

Figure 8 .
Figure 8.Comparison for the experimental data of Hao et al. (2014).

Table 1 Table 1 .
The dataset used for evaluating the models of submerged rigid vegetation.