Flow Rate and Velocity

Manning’s equation for uniform flow in a channel is used to calculate the rate and velocity of flow in a reach segment for a given time step:

qch=AchRch2/3slpch1/2nq_{ch}=\frac{A_{ch}*R_{ch}^{2/3}*slp_{ch}^{1/2}}{n} 7:1.2.1

vc=Rch2/3slpch1/2nv_c=\frac{R_{ch}^{2/3}*slp_{ch}^{1/2}}{n} 7:1.2.2

where qchq_{ch} is the rate of flow in the channel (m3^3/s), AchA_{ch} is the cross-sectional area of flow in the channel (m2^2), RchR_{ch} is the hydraulic radius for a given depth of flow (m), slpchslp_{ch} is the slope along the channel length (m/m), nn is Manning’s “n” coefficient for the channel, and vcv_c is the flow velocity (m/s).

SWAT+ routes water as a volume. The daily value for cross-sectional area of flow, AchA_{ch}, is calculated by rearranging equation 7:1.1.7 to solve for the area:

Ach=Vch1000LchA_{ch}=\frac{V_{ch}}{1000*L_{ch}} 7:1.2.3

where AchA_{ch} is the cross-sectional area of flow in the channel for a given depth of water (m2^2), VchV_{ch} is the volume of water stored in the channel (m3^3), and LchL_{ch} is the channel length (km). Equation 7:1.1.4 is rearranged to calculate the depth of flow for a given time step:

depth=Achzch+(Wbtm2zch)2Wbtm2zchdepth=\sqrt{\frac{A_{ch}}{z_{ch}}+(\frac{W_{btm}}{2*z_{ch}})^2}-\frac{W_{btm}}{2*z_{ch}} 7:1.2.4

where depthdepth is the depth of flow (m), AchA_{ch} is the cross-sectional area of flow in the channel for a given depth of water (m2^2), WbtmW_{btm} is the bottom width of the channel (m), and zchz_{ch} is the inverse of the channel side slope. Equation 7:1.2.4 is valid only when all water is contained in the channel. If the volume of water in the reach segment has filled the channel and is in the flood plain, the depth is calculated:

depth=depthbnkfull+(AchAch,bnkfull)zfld+(Wbtm,fld2zfld)2Wbtm,fld2zflddepth=depth_{bnkfull}+\sqrt{\frac{(A_{ch}-A_{ch,bnkfull})}{z_{fld}}+(\frac{W_{btm,fld}}{2*z_{fld}})^2}-\frac{W_{btm,fld}}{2*z_{fld}} 7:1.2.5

where depthdepth is the depth of flow (m), depthbnkfulldepth_{bnkfull} is the depth of water in the channel when filled to the top of the bank (m), AchA_{ch} is the cross-sectional area of flow in the channel for a given depth of water (m2^2), Ach,bnkfullA_{ch,bnkfull} is the cross-sectional area of flow in the channel when filled to the top of the bank (m2^2), Wbtm,fldW_{btm,fld} is the bottom width of the flood plain (m), and zfldz_{fld} is the inverse of the flood plain side slope.

Once the depth is known, the wetting perimeter and hydraulic radius are calculated using equations 7:1.1.5 (or 7:1.1.10) and 7:1.1.6. At this point, all values required to calculate the flow rate and velocity are known and equations 7:1.2.1 and 7:1.2.2 can be solved.

Table 7:1-2: SWAT+ input variables that pertain to channel flow calculations.

Variable Name
Definition
File Name

CH_S(2)

slpchslp_{ch}: Average channel slope along channel length (m m1^{-1})

.rte

CH_N(2)

nn: Manning’s “n” value for the main channel

.rte

CH_L(2)

LchL_{ch}: Length of main channel (km)

.rte

Last updated