MUSIC includes a model for grass swales, which models the storage-discharge relationship using Manning’s equation. A review of relevant literature shows, however, that Manning’s roughness coefficient, n, varies with vegetation type and height (relative to flow depth), as well as slope. An algorithm is adopted to incorporate this variation in roughness into the modelled discharge.
The following algorithm for modelling flow through a channel with flexible vegetation is based on work described by Kouwen (1988), Kouwen and Li (1980) and Kouwen and Unny (1973). The algorithm is based on determining the roughness of flexible vegetation for a defined flow depth.
The vegetation model is based on an empirical model, which uses an empirically defined vegetation parameter MEI to model the density and stiffness characteristics of the vegetation. The aim of the model is to reproduce the bending of the vegetation at higher discharges, and the consequent reduction of the bed roughness. The algorithm presented is a trial and error procedure to calculate the depth in uniform flow in the grass-lined channel.
The following equations are used in the algorithm:
n Manning roughness,
R Hydraulic radius of channel,
g Acceleration due to gravity,
k Deflected roughness height of vegetation, defined by Eq. 2
a, b Relative roughness parameters, defined in Table 1.
MEI Empirical vegetation stiffness parameter,
h Length of vegetation.
Bed Shear stress, defined by Eq. 3.
So Channel bed slope.
Darby (1999) provides equations for the specification of MEI for growing and dormant grasses after a study by Temple (1987).
For green growing grasses:(4)
Table 1: Values of a and b.
|v*/v*crit ≤ 1.0|
|1.0 < v*/v*crit ≤ 1.5|
|1.5 < v*/v*crit ≤ 2.5|
|2.5 < v*/v*crit|
v* Shear velocity, defined by Eq. 5,
v*crit Critical shear velocity, defined by the minimum of Eq. 6 and Eq. 7.
Application in MUSIC
MUSIC undertakes the following computations:
- Step 1 Calculate MEI for the grass based on Eq. 4.
- Step 2 Calculate the critical shear velocity, v*crit, from the minimum of Eq. 6 & Eq. 7.
- Step 3 Estimates an initial depth of flow in the channel
- Step 4 Calculate the cross-sectional properties of wetted perimeter, area and hydraulic radius (based on the guessed depth, and the swale design parameters provided by the user).
- Step 5 Calculate the bed shear stress, , from Eq. 3.
- Step 6 Calculate the deflected roughness height of the vegetation, k, from Eq. 2. Set the deflected roughness height k equal to the depth when k, calculated from Eq. 2, is greater than the depth.
- Step 7 Calculate the shear velocity, v*, from Eq. 5.
- Step 8 Determine the parameters a and b from Table 1.
- Step 9 Calculate the Manning’s n from Eq. 1.
- Step 10 Calculate the velocity and discharge in the channel from Manning Equation.
- Step 11 Compare the calculated discharge, Qcalc, with the required discharge, Qreq. Reduce the depth for Qcalc > Qreq and increase depth for Qcalc < Qreq. Return to Step 4 and repeat process until the calculated discharge, Qcalc, equals the required discharge, Qreq.
Darby, S. E., (1999), "Effect of Riparian Vegetation on Flow Resistance and Flood Potential", Journal of Hydraulic Engineering, ASCE, Vol. 125, No. 5, May 1999, pp. 443 - 454.
Kouwen, N. and Unny, T. E., (1973), "Flexible Roughness in Open Channels", Journal of the Hydraulic Division, Proceedings of ACSE, Vol 99, No. HY5, May 1973, pp. 713 - 728.
Kouwen N. and Li, R. M., (1980), "Biomechanics of Vegetative Channel Linings", Journal of the Hydraulic Division, Proceedings of ACSE, Vol 106, No. HY6, June 1980, pp. 1085 - 1103.
Kouwen, N., (1988), "Field Estimation of the Biomechanical Properties of Grass", Journal of Hydraulics Research, Vol. 26, 1988, No. 5., pp. 559 - 568.
Temple, D. M., (1987), "Closure of ‘Velocity Distribution Coefficients for Grass-Lined Channels’", Journal of Hydraulic Engineering, ASCE, Vol 113(9), 1987 pp. 1224 - 1226.