Model help:AgDegNormalGravMixHyd: Difference between revisions
mNo edit summary |
mNo edit summary |
||
Line 153: | Line 153: | ||
| Q<sub>w</sub> | | Q<sub>w</sub> | ||
| constant water discharge for the Wth flood | | constant water discharge for the Wth flood | ||
| | | L <sup>3</sup> / T | ||
|- | |- | ||
| Q<sub>p</sub> | | Q<sub>p</sub> | ||
| discharge for each time step p (p =1...P) | | discharge for each time step p (p =1...P) | ||
| | | L <sup>3</sup> / T | ||
|- | |- | ||
| p | | p | ||
Line 165: | Line 165: | ||
| Δt<sub>f</sub> | | Δt<sub>f</sub> | ||
| a short flood time step | | a short flood time step | ||
| | | T | ||
|- | |- | ||
| n<sub>step,w</sub> | | n<sub>step,w</sub> | ||
| an integer number of the flood time step | | an integer number of the flood time step | ||
| | | T | ||
|- | |- | ||
| x | | x | ||
| streamwise coordinate | | streamwise coordinate | ||
| | | T | ||
|- | |- | ||
| η | | η | ||
| river bed elevation | | river bed elevation | ||
| | | T | ||
|- | |- | ||
| t | | t | ||
| time step | | time step | ||
| | | T | ||
|- | |- | ||
| B | | B | ||
| river width | | river width | ||
| | | L | ||
|- | |- | ||
| D | | D | ||
| grain size of the bed sediment | | grain size of the bed sediment | ||
| | | L | ||
|- | |- | ||
| D<sub>bi</sub> | | D<sub>bi</sub> | ||
| bound diameter | | bound diameter | ||
| | | L | ||
|- | |- | ||
| λ<sub>p</sub> | | λ<sub>p</sub> | ||
Line 221: | Line 221: | ||
| D<sub>bnew</sub> | | D<sub>bnew</sub> | ||
| bound of the distribution such that the fraction of sediment finer than D<sub>bnew</sub> is equal to 0 | | bound of the distribution such that the fraction of sediment finer than D<sub>bnew</sub> is equal to 0 | ||
| | | L | ||
|- | |- | ||
| F<sub>subfli</sub> | | F<sub>subfli</sub> | ||
Line 241: | Line 241: | ||
| D<sub>s90</sub> | | D<sub>s90</sub> | ||
| the diameter of the bed surface such that the 90% of the sediment is finer | | the diameter of the bed surface such that the 90% of the sediment is finer | ||
| | | L | ||
|- | |- | ||
| n<sub>a</sub> | | n<sub>a</sub> | ||
Line 253: | Line 253: | ||
| η<sub>d</sub> | | η<sub>d</sub> | ||
| fixed bed elevation at the downstream end of the modeled reach | | fixed bed elevation at the downstream end of the modeled reach | ||
| | | L | ||
|- | |- | ||
| R | | R | ||
Line 261: | Line 261: | ||
| ξ<sub>d</sub> | | ξ<sub>d</sub> | ||
| downstream water surface elevation | | downstream water surface elevation | ||
| | | L | ||
|- | |- | ||
| q<sub>w</sub> | | q<sub>w</sub> | ||
| water discharge per unit width | | water discharge per unit width | ||
| | | L<sup>2</sup> / T | ||
|- | |- | ||
| k<sub>c</sub> | | k<sub>c</sub> | ||
| composite roughness height | | composite roughness height | ||
| | | L | ||
|- | |- | ||
| G | | G | ||
| imposed annual sediment transfer rate from upstream | | imposed annual sediment transfer rate from upstream | ||
| | | M / T | ||
|- | |- | ||
| G<sub>tf</sub> | | G<sub>tf</sub> | ||
Line 281: | Line 281: | ||
| ξ<sub>d</sub> | | ξ<sub>d</sub> | ||
| downstream water surface elevation | | downstream water surface elevation | ||
| | | L | ||
|- | |- | ||
| L | | L | ||
| length of reach under consideration | | length of reach under consideration | ||
| | | L | ||
|- | |- | ||
| q<sub>w</sub> | | q<sub>w</sub> | ||
| water discharge per unit width | | water discharge per unit width | ||
| | | L<sup>2</sup> / T | ||
|- | |- | ||
| i | | i | ||
Line 317: | Line 317: | ||
| U | | U | ||
| flow velocity | | flow velocity | ||
| | | L / T | ||
|- | |- | ||
| C<sub>f</sub> | | C<sub>f</sub> | ||
Line 325: | Line 325: | ||
| g | | g | ||
| acceleration of gravity | | acceleration of gravity | ||
| | | L / T<sup>2</sup> | ||
|- | |- | ||
| α<sub>r</sub> | | α<sub>r</sub> | ||
Line 333: | Line 333: | ||
| k<sub>s</sub> | | k<sub>s</sub> | ||
| grain roughness | | grain roughness | ||
| | | L | ||
|- | |- | ||
| n<sub>k</sub> | | n<sub>k</sub> | ||
Line 345: | Line 345: | ||
| ρ | | ρ | ||
| fluid density | | fluid density | ||
| | | M / L<sup>3</sup> | ||
|- | |- | ||
| ρ<sub>s</sub> | | ρ<sub>s</sub> | ||
| sediment density | | sediment density | ||
| | | M / L<sup>3</sup> | ||
|- | |- | ||
| τ<sub>c</sub> | | τ<sub>c</sub> | ||
Line 365: | Line 365: | ||
| q<sub>t</sub> | | q<sub>t</sub> | ||
| volume sediment transport rate per unit width | | volume sediment transport rate per unit width | ||
| | | L<sup>2</sup> / T | ||
|- | |- | ||
| I<sub>f</sub> | | I<sub>f</sub> | ||
Line 373: | Line 373: | ||
| t<sub>f</sub> | | t<sub>f</sub> | ||
| cumulative time the river has been in flood | | cumulative time the river has been in flood | ||
| | | T | ||
|- | |- | ||
| G<sub>t</sub> | | G<sub>t</sub> | ||
| the annual sediment yield | | the annual sediment yield | ||
| | | M / T | ||
|- | |- | ||
| t<sub>a</sub> | | t<sub>a</sub> | ||
Line 385: | Line 385: | ||
| Q<sub>f</sub> | | Q<sub>f</sub> | ||
| sediment transport rate during flood discharge | | sediment transport rate during flood discharge | ||
| L<sup>2</sup> / T | |||
|- | |- | ||
| α<sub>t</sub> | | α<sub>t</sub> | ||
Line 403: | Line 404: | ||
| C<sub>Z</sub> | | C<sub>Z</sub> | ||
| dimensionless Chezy resistance coefficient. | | dimensionless Chezy resistance coefficient. | ||
| - | |||
|- | |- | ||
| S<sub>l</sub> | | S<sub>l</sub> | ||
Line 410: | Line 412: | ||
| η<sub>i</sub> | | η<sub>i</sub> | ||
| initial bed elevation | | initial bed elevation | ||
| | | L | ||
|- | |- | ||
| D<sub>sub50</sub> | | D<sub>sub50</sub> | ||
| median size of the substrate layer | | median size of the substrate layer | ||
| | | L | ||
|- | |- | ||
| D<sub>subg</sub> | | D<sub>subg</sub> | ||
| geometric mean size of the substrate layer | | geometric mean size of the substrate layer | ||
| | | L | ||
|- | |- | ||
| L<sub>a</sub> | | L<sub>a</sub> | ||
| thickness of the active layer | | thickness of the active layer | ||
| | | L | ||
|- | |- | ||
| x | | x | ||
| downstream coordinate | | downstream coordinate | ||
| | | L | ||
|- | |- | ||
| τ | | τ | ||
| shear stress on bed surface | | shear stress on bed surface | ||
| | | - | ||
|- | |- | ||
| q<sub>b</sub> | | q<sub>b</sub> | ||
| bed material load | | bed material load | ||
| | | M / T | ||
|- | |- | ||
| Δx | | Δx | ||
| spatial step length, equals to L / M | | spatial step length, equals to L / M | ||
| | | L | ||
|- | |- | ||
| Q<sub>w</sub> | | Q<sub>w</sub> | ||
| flood discharge | | flood discharge | ||
| | | L<sup>3</sup> / T | ||
|- | |- | ||
| Δt | | Δt | ||
| time step | | time step | ||
| | | T | ||
|- | |- | ||
| Ntoprint | | Ntoprint | ||
Line 472: | Line 474: | ||
| η | | η | ||
| bed surface elevatioon | | bed surface elevatioon | ||
| | | L | ||
|- | |- | ||
| H | | H | ||
| water depth | | water depth | ||
| | | L | ||
|- | |- | ||
| ξ | | ξ | ||
| water surface elevation | | water surface elevation | ||
| | | L | ||
|- | |- | ||
| τ<sub>b</sub> | | τ<sub>b</sub> | ||
| bed shear stress | | bed shear stress | ||
| | | M / (T<sup>2</sup> L) | ||
|- | |- | ||
| S | | S | ||
Line 492: | Line 494: | ||
| q<sub>t</sub> | | q<sub>t</sub> | ||
| total bed material load | | total bed material load | ||
| | | L<sup>2</sup> / T | ||
|- | |- | ||
|} | |} |
Revision as of 14:01, 26 May 2011
AgDegNormalGravMixHyd
This is a module that calculates the evolution of a gravel bed river under an imposed cycled hydrograph.
Model introduction
This program is a close relative of AgDegNormGravMixP. It computes aggradation and degradation in gravel-bed river subject to a repeated hydrograph. The sediment is modeled as mixture of different grain sizes and the bedload formulation is that of Parker (1990) that was derived to compute the transport of gravel only.
Model parameters
Uses ports
This will be something that the CSDMS facility will add
Provides ports
This will be something that the CSDMS facility will add
Main equations
- time duration of a short flood
[math]\displaystyle{ \Delta t_{w} = n_{step,w} \Delta t_{f} }[/math] (1)
- Total time for flood
[math]\displaystyle{ P = \sum\limits_{w=1}^W n_{step,w} }[/math] (2)
Symbol | Description | Unit |
---|---|---|
Qw | constant water discharge for the Wth flood | L 3 / T |
Qp | discharge for each time step p (p =1...P) | L 3 / T |
p | time step in the hydrograph | |
Δtf | a short flood time step | T |
nstep,w | an integer number of the flood time step | T |
x | streamwise coordinate | T |
η | river bed elevation | T |
t | time step | T |
B | river width | L |
D | grain size of the bed sediment | L |
Dbi | bound diameter | L |
λp | bed porosity | - |
α | the parameter that governs the grain size distribution of the sediment at the active layer-substrate interface during bed aggredation | - |
Ffi | grain size distribution of the active layer for initial condition | - |
Fi | fraction of sediment in the ith grain size range in the active layer | - |
fsubfi | fraction of sediment in the ith grain size range in the substrate layer for initial condition | - |
npp | number of grain sizes for which a percent finer is specified, it must have the value between 2 and 9 | - |
Dbnew | bound of the distribution such that the fraction of sediment finer than Dbnew is equal to 0 | L |
Fsubfli | percent finer than ith grain size range for the substrate layer for initial condition | - |
fli | fraction of sediment in the ith grain size range in the active-layer substrate interface | - |
pi | fraction of sediment in the ith grain size range in the bedload | - |
Fsub,i | fraction of substrate material in the ith size range | - |
Ds90 | the diameter of the bed surface such that the 90% of the sediment is finer | L |
na | user specified order-one non dimensional constant | - |
pffi | the percent that is finer than the ith size range for upstream boundary conditon | - |
ηd | fixed bed elevation at the downstream end of the modeled reach | L |
R | submerged specific gravity | - |
ξd | downstream water surface elevation | L |
qw | water discharge per unit width | L2 / T |
kc | composite roughness height | L |
G | imposed annual sediment transfer rate from upstream | M / T |
Gtf | upstream sediment feed rate | - |
ξd | downstream water surface elevation | L |
L | length of reach under consideration | L |
qw | water discharge per unit width | L2 / T |
i | number of time steps per printout | - |
p | number of printouts desired | - |
M | number of spatial intervals | - |
R | submerged specific gravity of sediment | - |
Sf | friction slope | - |
Fr | Froude number | - |
U | flow velocity | L / T |
Cf | bed friction coefficient | - |
g | acceleration of gravity | L / T2 |
αr | coefficient in Manning-Stricker, dimensionless coefficient between 8 and 9 | - |
ks | grain roughness | L |
nk | dimensionless coefficient typically between 2 and 5 | - |
τ* | Shield number | - |
ρ | fluid density | M / L3 |
ρs | sediment density | M / L3 |
τc | critical Shields number for the onset of sediment motion | - |
ψs | the fraction of bed shear stress | - |
qt * | Einstein number | - |
qt | volume sediment transport rate per unit width | L2 / T |
If | flood intermittency | - |
tf | cumulative time the river has been in flood | T |
Gt | the annual sediment yield | M / T |
ta | the number of seconds in a year | - |
Qf | sediment transport rate during flood discharge | L2 / T |
αt | dimensionless coefficient in the sediment transport equation, equals to 8 | - |
nt | exponent in sediment transport relation, equals to 1.5 | - |
τc * | reference Shields number in sediment transport relation, equals to 0.047 | |
Cf | bed friction coefficient, equals to τb / (ρ U2 ) | - |
CZ | dimensionless Chezy resistance coefficient. | - |
Sl | initial bed slope of the river | - |
ηi | initial bed elevation | L |
Dsub50 | median size of the substrate layer | L |
Dsubg | geometric mean size of the substrate layer | L |
La | thickness of the active layer | L |
x | downstream coordinate | L |
τ | shear stress on bed surface | - |
qb | bed material load | M / T |
Δx | spatial step length, equals to L / M | L |
Qw | flood discharge | L3 / T |
Δt | time step | T |
Ntoprint | number of time steps to printout | - |
Nprint | number of printouts | - |
aU | upwinding coefficient (1=full upwind, 0.5=central difference) | - |
αs | coefficient in sediment transport relation | - |
Output
Symbol | Description | Unit |
---|---|---|
η | bed surface elevatioon | L |
H | water depth | L |
ξ | water surface elevation | L |
τb | bed shear stress | M / (T2 L) |
S | bed slope | - |
qt | total bed material load | L2 / T |
Notes
This program computes the time evolution of the long profile of a river of constant width carrying a mixture of gravel sizes, the downstream end of which has a prescribed elevation. In particular, the program computes the time evolution of the spatial profiles of bed elevation, total gravel bedload transport rate and grain size distribution of the surface (active) layer of the bed.
The flow is specified in terms of a hydrograph repeated a specified number of times annually, rather than a constant flood discharge and an intermittency.
The river has constant width. The upstream point, at which sediment is fed, is fixed in the horizontal to be at x = 0. The vertical elevation of the upstream point may change freely as the bed aggrades or degrades.In the program it is assumed that all gravel reaching the topset-foreset break is captured in the delta.
The reach has constant length L, so that the downstream point is fixed in the horizontal at x = L. This downstream point has a user-specified initial elevation ηdI.
Gravel bedload transport of mixtures is computed using the Parker (1990) surface-based formulation. Sand and finer material must be excluded from the grain size distributions before implementing this relation.
The grain size distributions of the sediment feed, initial surface material and substrate material must be specified. It is assumed that the grain size distribution of the sediment feed rate does not change in time, the initial grain size distribution of the surface material is the same at every node, the grain size distribution of the substrate is the same at every node and does not vary in the vertical.
The program does not store the vertical and streamwise structure of the new substrate created as the bed aggrades. As a result, is cannot capture the case of aggradation followed by degradation. Again, the constraint is easy to relax, but at the price of increased memory requirements for storing the newly-created substrate.
To perform a numerical calculation with a flow hydrograph, the actual hydrograph must be specified in terms of W constant water discharges Qw, where w = 1 … W, each extending for time duration Δtw. The river is assumed to be morphologically inactive when it is not in flood. The morphodynamic evolution is computed solving the equation of sediment continuity (i.e. Exner equation).
- Note on model running
The program may take a few seconds to calculate, depending on the user’s inputs—this is a calculation intensive program—please give the program some time to make the calculations.
There are no output values at time=0 for the geometric mean diameter of the load and the ration between the total bedload transport rate and the feed rate, because these are calculated in the time loop; they are not initial values.
Due to the many variables that are dependent upon user inputted values in this function the “ReadIn” function is called in the main, instead of in the Initialize portion of the code.
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: http://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)
References
- Parker, G., 1990, Surface-based bedload transport relation for gravel rivers, Journal of Hydraulic Research, 28(4): 417-436.
- Wilcock, P. R., and Crowe, J. C., 2003, Surface-based transport model for mixed-size sediment, Journal of Hydraulic Engineering, 129(2), 120-128.