Model help:TopoFlow-Channels-Dynamic Wave

TopoFlow-Channels-Dynamic Wave

The module is used to compute flow routing in a D8-based, spatial hydrologic model with dynamic wave process.

Model introduction

This process component is part of a spatially-distributed hydrologic model called TopoFlow, but it can now be used as a stand-alone model. The dynamic wave method is the most complete and complex method for modeling flow in open channels. This method retains all of the terms in the full, 1D momentum equation, including the gravity, friction and pressure gradient terms (as used by the diffusive wave method) as well as local and convective acceleration (or momentum flux) terms. This full equation is known as the St. Venant equation. In the current version of TopoFlow it is assumed that the flow directions are static and given by a D8 flow grid. In this case, integral vs. differential forms of the conservation equations for mass and momentum can be used.

Model parameters

Input 1

Input 2

Output 1

Output 2

Parameter	Description	Unit
Component status	Enabled / Disabled	[-]
Input directory	The location of the input files	[-]
Output directory	The location for the output files	[-]
Site prefix		[-]
Case prefix		[-]
Number of steps	Number of simulation steps	[-]
Time step		[sec]
D8 flow code file	grid of D8 flow codes in binary file	[-]
D8 slope file	grid of D8 flow slopes in binary file	[-]
Manning flag	Option to use Manning'n for roughness	[-]
Law of Wall flag	Option to use Law of Wall for roughness	[-]
Manning N type	grid of D8 flow slopes in binary file (Allowed input types: Scalar / Grid /Time_Series /Grid_Sequence )	[-]
Manning N	Manning'n value	[m / s^1/3]
Roughness z0 type	Allowed input types: Scalar / Grid /Time_Series /Grid_Sequence	[m]
Roughness z0	Law of Wall roughness value	[m]

Parameter	Description	Unit
Bed width type	Allowed input types: Scalar / Grid /Time_Series /Grid_Sequence	[-]
Bank width	bed width of trapezoid cross-section	[m]
Bank angle type	Allowed input types: Scalar / Grid /Time_Series /Grid_Sequence	[-]
Bank angle	bank angle of trapezoid cross-section	[degree]
Init. depth type	Allowed input types: Scalar / Grid /Time_Series /Grid_Sequence	[-]
Init. depth	initiate flow depth (If scalar, use 0.0)	[m]
Sinuosity type	Allowed input types: Scalar / Grid /Time_Series /Grid_Sequence	[-]
Sinuosity	absolute channel sinuosity	[m / m]

Parameter	Description	Unit
Save grid timestep	time interval between saved grids	[sec]
Save Q grids toggle	Option to save computed Q grids	[-]
Save Q grids file	file name for Q grid stack	[m^3 / s]
Save u grids toggle	Option to save computed u grids	[-]
Save u grids file	file name for u grid stack	[m / s]
Save d grids toggle	Option to save computed d grids	[-]
Save d grids file	file name for d grid stack	[m]
Save f grids toggle	Option to save computed f grids	[-]
Save f grids file	file name for f grid stack

Parameter	Description	Unit
Save pixels timestep	time interval between time series vales	[sec]
Save Q pixels toggle	Option to save computed Q time series	[-]
Save Q pixels file	file name for Q time series	[m^3 / s]
Save u pixels toggle	Option to save computed u time series	[-]
Save u pixels file	file name for u time series	[m / s]
Save d pixels toggle	Option to save computed d time series	[-]
Save d pixels file	file name for d time series	[m]
Save f pixels toggle	Option to save computed f time series	[-]
Save f pixels file	file name for f time series	[-]

Uses ports

• Meteorology
• Snow (Snowmelt)
• Evap (Evaporation)
• Infil (Infiltration)
• Satzone (Subsurface flow in saturated zone)
• Ice (Icemelt)
• Diversions (sources, sinks, canals)

Provides ports

• Channels (surface water flow in a network of channels)
• Configure (tabbed dialog GUI to change settings)
• Run (only if used as the Driver)

Main equations

• Mass conservation equation

[math]\displaystyle{ \Delta V \left (i,t \right)=\Delta t \ast [ R \left (i,t \right) \Delta x \Delta y -Q \left (i,t \right) +\Sigma_{k} Q \left (k,t \right) ] }[/math]

(1)

• Mean water depth in channel segment (if θ > 0 )

[math]\displaystyle{ d=\{[w^2 + 4 \tan \left (\theta\right) V / L] ^{\frac{1}{2}} -w\} / [ 2 \tan \left (\theta\right)] }[/math]

(2)

• Mean water depth in channel segment (if θ = 0)

[math]\displaystyle{ d= V / [w \ast L ] }[/math]

(3)

• Momentum conservation equation

[math]\displaystyle{ \Delta v \left (i,t \right)=\Delta t \ast \left ( T_{1} + T_{2} + T_{3} + T_{4} + T_{5} \right) / [d \left (i,t \right) \ast A_{w}] }[/math]

(4)

• Efflux term

[math]\displaystyle{ T_{1}=v \left (i,t \right) \ast Q \left (i,t \right) \ast \left ( C - 1 \right) }[/math]

(5)

• Influx term

[math]\displaystyle{ T_{2}=\Sigma_{k} [ v \left (k,t \right) - v \left (i,t \right) \ast C ] \ast Q \left (k,t \right) }[/math]

(6)

• "New mass" momentum term

[math]\displaystyle{ T_{3}=- v \left (i,t \right) \ast C \ast R \left (i,t \right) \ast \Delta x \ast \Delta y }[/math]

(7)

• Gravity term

[math]\displaystyle{ T_{4}=A_{w} \ast [g \ast d \left (i,t \right) \ast S \left (i,t \right)] }[/math]

(8)

• Friction term

[math]\displaystyle{ T_{5}=- A_{w} \ast [ f\left (i,t \right) \ast v \left (i,t \right)^2 ] }[/math]

(9)

• Discharge of water

[math]\displaystyle{ Q=v \ast A_{w} }[/math]

(10)

• Friction factor (for law of the wall)

[math]\displaystyle{ f \left (i,t \right)= [\kappa / LN \left (a \ast d \left (i,t \right) / z_{0}\right)]^2 }[/math]

(11)

• Friction factor for Manning's equation

[math]\displaystyle{ f \left (i,t \right)=g \ast n^2 / R_{h}\left (i,t \right)^{\frac{1}{3}} }[/math]

(12)

• Area ratio appearing

[math]\displaystyle{ C=A_{w} / A_{t} }[/math]

(13)

• Top surface area of a channel segment

[math]\displaystyle{ A_{t}=w_{t} \ast L }[/math]

(14)

• Top width of a wetted trapezoidal cross-section

[math]\displaystyle{ w_{t}= w + [ 2 \ast d \ast \tan \left ( \theta \right ) ] }[/math]

(15)

• hydraulic radius

[math]\displaystyle{ R_{h}= A_{w} /P_{w} }[/math]

(16)

• wetted cross-sectional area of a trapezoid

[math]\displaystyle{ A_{w}= d \ast [ w + \left (d \ast \tan \left (\theta\right)\right)] }[/math]

(17)

• Wetted perimeter of a trapezoid

[math]\displaystyle{ P_{w}= w + [ 2 \ast d / cos\left (\theta\right)] }[/math]

(18)

• wetted volume of a trapezoidal channel

[math]\displaystyle{ V_{w}=d^2 \ast [ L \ast \tan \left (\theta\right)] +d \ast [ L \ast w ] }[/math]

(19)

Nomenclature

Symbol	Description	Unit
ΔV	change in water volume	m3
Δt	time step	sec
R	effective rain rate / excess rainrate, represents the sum of all vertical contributions to a grid cell's mass balance	m / s
Δx	distance change in x direction	m
Δy	distance change in y direction	m
w	width of channel	m
θ	bank angle for trapezoid	deg
V	water volume	m3
L	channel length	m
T1	efflux term	m4 / s2
T2	influx term	m4 / s2
T3	"new mass" momentum term	m4 / s2
T4	gravity term	m4 / s2
T5	friction term	m4 / s2
Aw	wetted cross-sectional area of a trapezoid	m2
n	Manning's n	s / m1/3
Pw	wetted perimeter of a trapezoid	m2
S	bed slope	m / m
g	gravity acceleration	m / s2
z0	roughness length	m
κ	Von Karman's constant, equals to 0.41	-
a	constant	-
C	constant	-
At	top surface area of a channel segment	m2
Wt	top width of a wetted trapezoidal cross section	m
Vw	wetted volume of a trapezoidal channel	m

Output

Symbol	Description	Unit
Q	discharge of water	m3 / s
v	flow velocity	m / s
d	mean channel flow depth	m
f	friction factor	-
Rh	hydraulic radius	m
S_free	free-surface slope	m / m

Notes

Notes on Input Parameters

The input variables for the dynamic wave method should usually be specified as grids, except in special cases.

Single grids and grid sequences are assumed to be stored as RTG and RTS files, respectively. Time series are assumed to be stored as text files, with one value per line. For a time series or grid sequence, the time between values must coincide with the timestep provided.

Flow directions are determined by a grid of D8 flow codes. All grids are assumed to be stored as RTG (RiverTools Grid) files and flow codes are assumed to follow the Jenson (1984) convention (i.e. [NE,E,SE,S,SW,W,NW,N] → [1,2,4,8,16,32,64,128]) that is also used for RiverTools D8 flow grids. Flow grids and slope grids can be created by RiverTools or a similar program and the other grids can be created using tools in the TopoFlow Create menu.

Bed slope, S, can be computed from a DEM by using the Create → Profile-smoothed DEM dialog or by using hydrologic GIS software. The current version assumes that all channels have trapezoidal cross-sections (see Notes below) but allows bottom-width and bank angle to vary spatially as grids. TopoFlow has pre-processing tools in the Create menu for creating grids of bed width, bank angle and bed roughness. The Create → Channel Geometry Grids → With Area Grid tool allows you to parameterize these variables as power-law functions of contributing area. The Create → Channel Geometry Grids → With HS Order Grid tool allows you to assign values based on Horton-Strahler order.

Each pixel is classified as either a hillslope pixel (overland flow) or a channel pixel (channelized flow) and appropriate parameters must be used for each. For overland flow, w >> d, Rh → d, and bank angle drops out. Overland flow can then be modeled with a large value of Manning's n, such as 0.3. For channelized flow, the variation of n with bed grain size can be modeled with Strickler's equation as explained in the Notes below.

If a sinuosity greater than 1 is specified, then bed slopes are reduced by dividing them by this value. As with the other variables, it is most appropriate to specify a grid in this case.

It is physically unrealistic to specify a spatially uniform initial flow depth by entering a scalar value greater than zero for init_depth. This will result in a very large peak in the hydrograph and may cause TopoFlow to crash. The Create → RTG File for Initial Depth tool can be used to create a grid of initial flow depths that varies spatially and is in steady-state equilibrium with a specified baseflow recharge rate.

Notes on the Equations

All variables and their units can be seen by expanding the Nomenclature section above.

The dynamic wave method is the most complete and complex method for modeling flow in open channels. This method retains all of the terms in the full, 1D momentum equation, including the gravity, friction and pressure gradient terms (as used by the diffusive wave method) as well as local and convective acceleration (or momentum flux) terms. This full equation is known as the St. Venant equation. In the current version of TopoFlow it is assumed that the flow directions are static and given by a D8 flow grid. In this case, integral vs. differential forms of the conservation equations for mass and momentum can be used.

Conservation of mass, in integral form, is represented by the equations above. The quantity, R, that appears in the first equation is known as the effective rainrate or excess rainrate and represents the sum of all vertical contributions to a grid cell's mass balance. R is computed as R = (P + M + G) - (I + E), where P = precipitation, M = snowmelt, G = seepage from subsurface, I = infiltration and E = evapotranspiration. (Note that R is technically not the same as the runoff, since runoff includes horizontal fluxes.) The summation sign in the first equation adds up all horizontal inflows to a grid cell from its neighbor grid cells. Mean channel flow depth, d, is then computed from channel geometry and the water volume that is computed for the corresponding grid cell. Note that channel length depends on distance between grid cell centers and sinuosity, while cross-sections are trapezoidal. When the bank angle, θ is greater than zero, the flow depth required to accomodate the water volume is computed by solving the last equation (a quadratic) for d to get the second equation.

Conservation of momentum, in integral form, is represented by the remaining equations above. The term T3 results from using the equation for mass conservation to simplify the one for momentum conservation. It represents the change in a grid cell's momentum balance due to the addition or subtraction of mass within the cell from vertical fluxes. Recall that momentum represents a product of mass and velocity, so changes in momentum can result from changes in mass or changes in velocity (chain rule).

Notice that boundary-layer theory (which leads to the law of the wall) allows us to compute the loss of momentum due to friction directly as a separate process and isn't really restricted to the special case of steady, uniform flow. By contrast, Manning's equation is empirical and was found by studying steady, uniform flows. It does not allow us to cleanly separate the effects of friction and gravity; it assumes that they are in balance. However, if we use the fact that Manning's equation is really just a power-law approximation to the law of the wall, we can back-calculate a Manning's-equation version of the friction term (and friction factor) which has been freed from the assumption of steady, uniform flow.

In the current version of TopoFlow (1.5 beta), water-surface slopes are set to zero if they ever become negative (implying upstream flow).

Note on the current version

It has also been assumed here that the water is incompressible; this allows us to take the density of water, ρ, to be a constant.

Examples

An example run with input parameters, BLD files, as well as a figure / movie of the output

Follow the next steps to include images / movies of simulations:

• Upload file: https://csdms.colorado.edu/wiki/Special:Upload
• Create link to the file on your page: [[Image:<file name>]].

See also: Help:Images or Help:Movies

Developer(s)

Scott Peckham

References

Peckham, S.D. (2008) Geomorphometry and spatial hydrologic modeling, In: Geomorphometry, chapter 25, Elsevier, New York, in press.

Peckham, S.D. and J.D. Smith (2008) Manning's equation and the best power-law approximation to the logarithmic law of the wall, in preparation.

Schlicting, H. (1960) Boundary Layer Theory, 4th ed., McGraw-Hill, New York, 647 pp.

Links

• Model help:TopoFlow-Channels-Kinematic Wave
• Model help:TopoFlow-Channels-Diffusive Wave
• Model help:TopoFlow-Diversions

• Model:TopoFlow-Channels-Dynamic_Wave (Model metadata)
• Model:TopoFlow