An Improved Analytical Model for Vertical Velocity Distribution of Vegetated Channel Flows

The existence of vegetation plays an important role to protect the ecosystem and water environment in natural rivers and wetlands, but it alters the velocity field of flow, consequently influencing the transport of pollutant and biomass. As a pre-requisite for the analysis of environmental capacity in a channel, the vertical velocity distribution of flows has attracted much research attention; however, there is yet lack of a good prediction model available. For the channel with submerged vegetation, the vertical velocity distribution in the lower vegetation layer will be different from that in the upper flow layer of non-vegetation. In this paper, after review on the most recent two-layer model proposed by Baptist et al., the author has proposed an improved two-layer analytical model by introducing a different mixing length scale (λ). The proposed model is based on the momentum equation of flow with the turbulent eddy viscosity assumed as a linear relationship with the local velocity. The proposed model is compared with the Baptist model for different datasets published in the literature, which shows that the proposed analytical model can improve the vertical velocity distribution prediction well compared with the Baptist model for a range of data. This study reveals that the λ is well related with the submergence of vegetation (H/h), as suggested by 1 h H h λ β = − . When the constant β is taken as 3/100, the proposed model shows good agreement with a wide range of datasets studied: flow depth (H)/vegetation height (h) in 1.25 to 3.33, different vegetation densities of a in 1.1 to 18.5 m (a defined as the frontal area of the vegetation per unit volume), and bed slopes in (1.38 4.0) × 10.


Introduction
In many natural rivers and wetlands, various types of vegetation exist like bushes, trees, herbs or any other kinds of plant.Sometime they are deliberately arranged to meet engineering or ecological requirement in river or wetland planning and management.The presence of vegetation will affect the hydrodynamic structure of flow by exerting drag force on the flow.In a river or channel with submerged vegetation, its flow is significantly affected by the discontinuous vegetation drag, resulting in a shear layer generated around the top of vegetation.Many researchers have studied the interaction between vegetation and flow, which includes the flow resistance (Kouwen et al., 1969;Temple, 1986;Stone & Shen, 2002;Cheng, 2015), velocity, Reynolds stress and turbulence intensity (Nepf, 1999;Nepf & Vivoni, 2000;Lopez & Garcia, 2001;Ghisalberti & Nepf, 2006;Poggi et al., 2004;Stoesser et al., 2010;Rahimi et al., 2018;Tang et al., 2018).
In the open-channel flow with single-layer submerged vegetation, a strong moment exchange exits between the fast moving flow in the surface layer and the slowly moving flow in the lower vegetation layer, consequently generating a shear region at the interface near the top of vegetation (Poggi et al., 2004;Nezu & Sanjou, 2008).The flow structure in this region is governed through the mixing length of eddy and vertical turbulent transport of momentum from the surface flow, with the negligible contribution from the pressure gradient.The effect of the shear length scale can be related to the vegetation height and flow depth (Klopstra et al., 1997;Defina & Bixio, 2005;Tang, 2018aTang, , 2018bTang, , 2019a;;Singh et al., 2019).
A review on the studies on the vertical velocity distribution of flow with submerged rigid vegetation shows that different models are used to predict the vertical velocity distribution through and above submerged vegetation (Tang & Ali, X. Tang Journal of Geoscience and Environment Protection 2013; Tang, 2019a).These models are so-called two-layer or three-layer methods.In the layer of submerged vegetation, the analytical model of velocity distribution is based on the concept of mixing-layer analogy of eddy (Klopstra et al., 1997;Meijer & Van Velzen, 1999;Ghisalberti & Nepf, 2004;Defina & Bixio, 2005;Baptist et al., 2007;Huai, et al., 2009Huai, et al., , 2014;;Yang & Choi, 2010;Nepf, 2012;Tang, 2018a).
For instance, in the two-layer model by Klopstra et al. (1997) the vertical velocity distribution in the vegetation layer was obtained through the analogy of turbulent shear stress.The turbulent eddy viscosity is approximated as the product of the local velocity (u) and a mixing length (λ), which is related to the submergence ratio (H/h) via an empirical formula, see Figure 1.This empirical formula of λ was obtained by fitting limited experimental data.Afterwards, further studies on this approach were undertaken, which suggested that λ can be approximated by 0.0144 Hh after more data evaluation (Meijer & Van Vel-  zen, 1999; Defina & Bixion, 2005).Similarly, Baptist et al. (2007) established their analytical solution of vertical velocity distribution in the vegetation layer, and they suggested that λ is approximated by 0.05(H − h).Recently, Nepf (2012) proposed an alternative model based on some empirical equations to describe the vertical velocity distribution.
In the evaluation of the above models and other models, Tang & Ali (2013) and Tang (2019a) showed that the predicted velocity distribution by the above mentioned models is acceptable only for certain experimental datasets, but they do not agree with other datasets.In addition, the predictions by these models show significant differences even for the same dataset.This is not surprising because different constants in their models of λ were used.Moreover, none of these models appear to predict well against a wide range of data tested (Tang & Ali, 2013;Tang, 2018b).
In the present paper, firstly review and discussion is carrying out on the most recent two-layer analytical model by Baptist et al. (2007).Then, based on the mixing length concept of turbulent eddy, the author proposed an alternative im- X. Tang , where β is a constant.Finally, a wide range of experimental data from the literature were used to evaluate the proposed analytical model along with the model by Baptist et al. (2007).It was found that when β = 3/100 the proposed model agrees well with a range of datasets, which are given by Dunn et al. (1996), Meijer & Van Velzen (1999), Nguyen (2012) and Hao et al. (2014).However, the Baptist model only suits for certain cases with lower submergence and vegetation density.

Description of the Analytical Model
For the channel flow with vegetation, the bed and wall boundary stress are both assumed to be negligible compared with the drag force on the vegetation.Under the steady uniform flow, the governing momentum equation for fully-developed 1-D vegetated flow can be described by where τ is the shear stress, ρ the density of water, g the gravity, z the vertical coordinate above the bed, S o the bed slope and F v is the drag force per unit weight generated by the vegetation, see Figure 1.The drag force F v is given by: where u is the stream wise time-averaged velocity, h is the height of vegetation, C D is the drag coefficient, a is the density of vegetation, i.e. the frontal area of vegetation (A v ) per unit volume.
To solve the momentum Equation (1) analytically, Baptist et al. (2007) introduced Boussinesq's eddy viscosity concept to approximate the turbulent shear stress (τ) as follows: ( ) ( ) where ν t is the total eddy viscosity of vegetated flow, l is the length scale of vegetation, and c p is the parameter of turbulence intensity.
If the length scale ( l ) and the turbulence intensity parameter (c p ) are assumed constant within the vegetation, Equation (1) becomes For given density (a) of vegetation and drag coefficient (C D ), an analytical solution for u 2 in Equation ( 4) can be obtained: where a v and b v are the integration constants, which can be determined by appropriate boundary conditions, and Afterwards, by ignoring the last term in the right hand of Equation ( 5), Baptist et al. (2007) simplified Equation (5) to Equation ( 8): where a v is described by the following expression: Equation ( 8) describes the velocity profile in the vegetation layer.
For the surface layer above the vegetation (i.e.non-vegetation layer), where the velocity can be described by the well-known logarithmic profile, as used by Klopstra et al. (1997), Nepf & Vivoni (2000), Tang (2018a): ( ) 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 from the top of vegetation to 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 ( ) Furthermore, h s and z o are expressed respectively by: ( ) ( ) Baptist et al. (2007) validated their method using the data of Nepf & Vivoni (2000) and recommended that ( ) Baptist et al. (2007) proposed a simplified analytical model of velocity profile, which mainly aims to introduce an equation for calculating flow resistance in a channel with vegetation.This will help to evaluate the roughness and modelling.

Discussions on the Analytical Model
However, some assumptions have been made in the development of their analytical model, where some aspects of the assumptions are open for further discussions, which are given as follows: In the abovementioned analytical model, the analytical solution ( 8) is an approximated, simplified solution to the full solution (5) by ignoring the term of ( )  , which is usually not zero because b v is not zero.When applying the boundary condition at the top of the vegetation (z = h), where the boundary shear stress is given by ( ) , one can obtain b v = −a v .There-X.Tang fore, the solution of Equation ( 5) may result in some errors, which could not be negligible in certain cases.Thus, it is expected that the solution of Equation ( 8) may over-predict the velocity near the channel bed in the vegetation layer.
Moreover, at the bed (z = 0) of vegetated flow, where the bed shear stress is neglected compared with the drag force of vegetation, the flow velocity can be described by u o , given by Equation ( 6), based on the local equilibrium between gravity force and vegetation drag.However, the solution of Equation ( 8) does not lead to u o when z approaches to zero, because ( ) which is usually not zero but a small limit value.In this respect, the approximate solution of Equation ( 8) may not properly represent the velocity profile very close to the bed, thus consequently affecting the vertical velocity distribution in both layers.
Finally, since only limited data were used for the calibration of the parameters (c p and l) in the model of Baptist et al. (2007), the recommended relationship of ( ) may not suit for a wide range of data.
Based on the discussion above, a full solution, rather than a simplified solution, to the governing Equation ( 1) should be developed in the next section.

Improved Analytical Model
In this paper, like Klopstra et al. (1997), Defina & Bixio (2005) and Tang (2018a), the turbulent shear stress (τ) in Equation ( 1) is approximated by Boussinesq hypothesis through a mixing length concept as: where λ is a mixing length of eddy.
For steady flow, inserting Equations ( 2) & ( 13) into (1) produces: ( ) If vegetation density (a), drag coefficient (C D ) and the mixing length (λ) are given, Equation ( 14) becomes an ordinary differential equation in terms of u 2 .The solution of u 2 can be obtained in certain boundary conditions.The boundary conditions are considered as follows: At the bed (z = 0), where the bed shear stress is negligible in comparison to the drag force of vegetation, the local equilibrium between vegetation drag and gravity force leads to: At the top of vegetation (z = h), the turbulent shear stress is given by: Based on the boundary conditions of Equations ( 15) & ( 16), the solution of Equation ( 15) to u gives the velocity in the vegetation layer, as expressed by where λ/h is related to the submergence H/h.In this study, it is assumed as For the surface layer (i.e.non-vegetation layer), the velocity takes the same form as Equation ( 10), where the unknown parameters h s and z o can be obtained by applying the continuity conditions of the velocity and velocity gradient (or the shear stress) at the interface (i.e. at z = h) between the surface layer and the vegetation layer.Thus, h s and z o can respectively be described as:

Data Used for This Study
To evaluate both the proposed new model and the Baptist model describe above, the author used a range of different experimental data of submerged rigid vegetation published in the literature.A total of 17 datasets used in this study covers various submergence (H/h), which is 1.25 -3.33, the vegetation densities a = 1.1 -18.5 m −1 , and bed slope S o ranging from 1.38 × 10 −3 to 4 × 10 −3 .The summary of the datasets are shown in Table 1, where the experimental data of Nguyen (2012) include various vegetation densities and different diameter sizes of vegetation.The values of C D in the dataset of emergent vegetation were taken from the original papers if available or assumed according to White (1974) if not.

Results
As shown in Section 3, the characteristic length λ reflects the mixing strength of eddy of submerged vegetative flow.The λ is closely affected by the submergence (H/h), as recommended as 1 H h β − in the present study.The larger β value, the bigger λ becomes, i.e. the stronger the eddy near the top of vegetation.This implies that the velocity deceases as increasing β, as demonstrated in Figure 2 as an example.Figure 2 demonstrates that the value of β has much greater impact on the predictive velocity in the surface layer than in the vegetation layer.Especially, the velocity close to the channel bed is almost not affected by the constant β.In other words, λ will have great influence on the velocity distribution in the upper region of depth.After extensive tests on β values in the proposed λ model, the optimal value of β = 3/100 was found for the proposed analytical model of u.This value of β is used in the subsequent figures.
For the simplicity, in the subsequent figures, the model denotes the proposed new model with β = 3/100 in this study, whereas the Baptist model denotes the model proposed by Baptist et al. (2007), as described in Section 2.

X. Tang
Figures 3-7 show the results of comparison between the models (both the new and Baptist models) and the data for the cases of the submergence of H/h ≥ 1.5.The comparison shows that the predicted velocity profiles by the new model generally agree well with all the experimental data tested, whereas the Baptist model does not predict well for most the cases when the flows have relatively high submergence (e.g. 2 H h  ), see Figure 3(a), Figure 4(c), Figure 5( a)   and Figure 6(b).In these high submergence cases, the Baptist model significantly under-predicts the velocity in the surface layer, indicating that the length scale (c p l) is over-estimated because of a large constant of (H-h) used in the model.
Note that in Figure 7 for the case of Nguyen (2012), both models agree the data reasonably well when the density of vegetation is not too high (e.g.ah < 1.8) under moderate submergence (H/h = 1.5).However, if the density of vegetation is very large, e.g.
1.8 ah  in case C30 -15, it appears that the predictive velocities by the two models are over-estimated in both layers, which indicates the stronger mixing of eddy due to the high density of vegetation.In such a case, because the wake effect between the vegetation of finite height becomes strong, it is expected that the drag coefficient (C D ) will be affected by the density and arrangement formation of vegetation.Thus, for the case of C30 -15, as the drag

Discussion
In the region close to the bed, the predicted velocity by the Baptist model is often larger than the predictive velocity by the proposed model, see Figures 3-9.This is not surprising as discussed in Section 2.2.This is because the simplified solution ( 8) is used after a removal of the term ( ) , where a v is always a positive value.
To further check the robustness of the models with the experimental data, error analyses were carried out.MAPE (Mean Absolute Percentage Error) is a measure of the error percentage for predicted values of model against the measured values.These individual differences are called as residuals for the data sample that is used for estimation, and the residuals are known as estimation errors for the sample (Tang, 2017(Tang, , 2019b)).
The percentage of error in predicted discharge of each flow depth is calculated by.
, , , , X. Tang where E u,i is the error percentage of predicted velocity, and u a,i and u e,i are the predicted and observed velocity at i-th flow depth, respectively.Therefore, the average error percentage of model for an experiment can be calculated by X ( ) The values (E u ) of MAPE for all datasets are shown in Table 2, which shows that the proposed model has less the averaged percentage error (9%) than the Baptist model (about 11.1%).
As demonstrated in this study, the proposed model with λ/h taken as 1 H h β − shows good agreement with a range of the data tested, whereas the Baptist model does not agree well for most of the data when H/h is very high (e.g.

H h 
).The reason may be due to the fact that the simplified solution of Equation ( 8) is used in the Baptist model.When the 3 rd term in the full solution of Equation ( 5) is ignored, the velocity will be over-predicted, consequently affecting the flow velocity in the upper layer (i.e.over-predicted).Another reason may be due to only the flow depth (H-h) of free flow used for the evaluation of mixing length parameter (c p l).For the cases with higher H/h ( 2  ), the value of (H-h) can be relatively large, which could result in a large value of c p l. Thus, a v in Equation ( 9) is under-valued.Therefore, the predicted velocity from Equation (8) becomes smaller.As a result, the velocity by the Baptist model is under-predicted for higher H/h cases, see Figure 3 In the improved analytical model, the analytical solution of Equation ( 17) is a full solution to Equation ( 15).H/h is also directly considered in both the solution Equation ( 17) and the evaluation of mixing length λ.In other words, the proposed model has considered the impact of H/h on the prediction of velocity.
Nevertheless, further study may need to establish the recommended value of constant β = 3/100 for even a wide range of data in the future.

Conclusion
In the open-channel flow of submerged vegetation, the additional resistance of vegetation drag will significantly affect the vertical velocity distribution within the vegetation, which differentiates from that in the upper non-vegetation region.Based on Boussinesq approximation for the mixing length of eddy in the momentum equation, an analytical solution of velocity can be obtained as Equation (17), which describes the velocity profile in the vegetation region.Equation ( 17) meets the u o near the bed while Equation ( 8) of the Baptist model does not meet.Through tests against a range of experimental datasets, it was found that the proposed model can predict the velocity reasonably well against a wide range of data for the cases of moderate submergence (H/h ≥ 1.5) when the normalized mixing length (λ/h) is approximated by 1 H h β − with β having the optimal value of 3/100.In comparison, the Baptist model does not predict well for the cases when the submergence (H/h) is high (e.g.


).The MAPE error analysis shows that the proposed model has overall less error compared with the Baptist model.
the constant β is taken as 3/100, the proposed model shows good agreement with a wide range of datasets studied: flow depth (H)/vegetation height (h) in 1.25 to 3.33, different vegetation densities of a in 1.1 to 18.5 m −1 (a defined as the frontal area of the vegetation per unit volume), and bed slopes in (1.38 -4.0) × 10 −3 .
Figure 1.Vertical velocity distribution of flow through submerged vegetation.
where β is a constant.The velocity of u in Equation (17) will become u o when z approaches to zero (i.e.close to the bed) because 17) does not have such an issue in the Baptist model.

Figure 2 .
Figure 2. The effect of β on the vertical velocity distribution for the data of Dunn et al.(1996).

X
. Tang Journal of Geoscience and Environment Protection

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

Table 2 .
The MAPE errors (in percentage) of models for the dataset.