Model help:AgDegNormalFault: Difference between revisions
No edit summary |
m (Text replacement - "http://csdms.colorado.edu/wiki/" to "https://csdms.colorado.edu/wiki/") |
||
(7 intermediate revisions by 2 users not shown) | |||
Line 3: | Line 3: | ||
1) Log in to the wiki | 1) Log in to the wiki | ||
2) Create a new page for each model, by using the following URL: | 2) Create a new page for each model, by using the following URL: | ||
* | * https://csdms.colorado.edu/wiki/Model help:<modelname> | ||
* Replace <modelname> with the name of a model | * Replace <modelname> with the name of a model | ||
3) Than follow the link "edit this page" | 3) Than follow the link "edit this page" | ||
Line 172: | Line 172: | ||
| Q | | Q | ||
| flood discharge | | flood discharge | ||
| | | L <sup>3</sup> / 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 | ||
|- | |- | ||
| λ<sub>p</sub> | | λ<sub>p</sub> | ||
| bed porosity | | bed porosity | ||
| - | | - | ||
|- | |- | ||
| 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 224: | Line 208: | ||
| ξ<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 | ||
|- | |- | ||
| i | | i | ||
Line 260: | Line 240: | ||
| U | | U | ||
| flow velocity | | flow velocity | ||
| | | L / T | ||
|- | |- | ||
| C<sub>f</sub> | | C<sub>f</sub> | ||
Line 268: | Line 248: | ||
| g | | g | ||
| acceleration of gravity | | acceleration of gravity | ||
| | | L / T<sup>2</sup> | ||
|- | |- | ||
| α<sub>r</sub> | | α<sub>r</sub> | ||
Line 276: | Line 256: | ||
| k<sub>s</sub> | | k<sub>s</sub> | ||
| grain roughness | | grain roughness | ||
| | | L | ||
|- | |- | ||
| n<sub>k</sub> | | n<sub>k</sub> | ||
Line 288: | Line 268: | ||
| ρ | | ρ | ||
| 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 308: | Line 288: | ||
| 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 315: | Line 295: | ||
|- | |- | ||
| t<sub>f</sub> | | t<sub>f</sub> | ||
| the | | 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 328: | Line 308: | ||
| 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 339: | Line 320: | ||
| τ<sub>c</sub> <sup>*</sup> | | τ<sub>c</sub> <sup>*</sup> | ||
| reference Shields number in sediment transport relation, equals to 0.047 | | reference Shields number in sediment transport relation, equals to 0.047 | ||
|- | |- | ||
| C<sub>Z</sub> | | C<sub>Z</sub> | ||
Line 353: | Line 330: | ||
| η<sub>i</sub> | | η<sub>i</sub> | ||
| initial bed elevation | | initial bed elevation | ||
| | | 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 393: | Line 366: | ||
| α<sub>s</sub> | | α<sub>s</sub> | ||
| coefficient in sediment transport relation | | coefficient in sediment transport relation | ||
| - | |||
|- | |||
| u<sub>*</sub> | |||
| shear velocity | |||
| L / T | |||
|- | |||
| α<sub>r</sub> | |||
| coefficient in Manning-Strickler resistance relation | |||
| - | |||
|- | |||
| τ<sub>b</sub> <sup>*</sup> | |||
| non-dimensional total shear stress | |||
| - | | - | ||
|- | |- | ||
Line 401: | Line 386: | ||
| Δη | | Δη | ||
| the height of faulting | | the height of faulting | ||
| | | L | ||
|- | |- | ||
|} | |} | ||
Line 411: | Line 396: | ||
| η | | η | ||
| 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 | ||
| bed slope | | bed slope | ||
| - | | - | ||
|- | |- | ||
|} | |} | ||
Line 445: | Line 426: | ||
<span class="remove_this_tag">Follow the next steps to include images / movies of simulations:</span> | <span class="remove_this_tag">Follow the next steps to include images / movies of simulations:</span> | ||
* <span class="remove_this_tag">Upload file: | * <span class="remove_this_tag">Upload file: https://csdms.colorado.edu/wiki/Special:Upload</span> | ||
* <span class="remove_this_tag">Create link to the file on your page: <nowiki>[[Image:<file name>]]</nowiki>.</span> | * <span class="remove_this_tag">Create link to the file on your page: <nowiki>[[Image:<file name>]]</nowiki>.</span> | ||
Line 451: | Line 432: | ||
==Developer(s)== | ==Developer(s)== | ||
[[User:Gparker|Gary Parker]] | |||
==References== | ==References== | ||
Line 459: | Line 440: | ||
==Links== | ==Links== | ||
* [[ | * [[Model:AgDegNormalFault]] | ||
* [[ | * [[Model_help:AgDegNormal]] | ||
[[Category:Utility components]] | [[Category:Utility components]] |
Latest revision as of 17:17, 19 February 2018
AgDegNormalFault
This is used to calculate aggradation and degradation of a river reach using the normal flow approximation; with an extension for calculation of the response to a sudden fault along the reach.
Model introduction
This program computes 1D bed variation in rivers due to differential sediment transport in which it is possible to allow the bed to undergo a sudden vertical fault of a specified amount, at a specified place and time. Faulting is realized by moving all notes downstream of the specified point downward by the amount of the faulting. It uses the same principles of AgDegNormal model but with extension for calculation of the response to a sudden fault along the reach.
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
Use the same equations as AgDegNormal model
Symbol | Description | Unit |
---|---|---|
Q | flood discharge | L 3 / T |
t | time step | T |
B | river width | L |
D | grain size of the bed sediment | L |
λp | bed porosity | - |
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 |
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 | |
CZ | dimensionless Chezy resistance coefficient. | |
Sl | initial bed slope of the river | - |
ηi | initial bed elevation | 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 | - |
u* | shear velocity | L / T |
αr | coefficient in Manning-Strickler resistance relation | - |
τb * | non-dimensional total shear stress | - |
rf | the fraction of reach length such that all point downstream of x = rfL undergo downward faulting | - |
Δη | the height of faulting | L |
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 | - |
Notes
The sediment is assumed to be uniform with size D. All sediment transport is assumed to occur in a specified fraction of time during which the river is in flood, specified by an intermittency. A Manning-Strickler formulation is used for bed resistance. A generic relation of the general form of that due to Meyer-Peter and Muller is used for sediment transport. The flow is computed using the normal flow approximation.
If the channel slope is negative and the water depth is not a number, “nan”, check the time step and the spatial step length. In particular, the time step may be too large or equivalently the spatial step length may be too small. Change these values and run the model again.
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)
References
- Paola, C., Heller, P. L. & Angevine, C. L. 1992 The large-scale dynamics of grain-size variation in alluvial basins. I: Theory. Basin Research, 4, 73-90.
- Meyer-Peter, E., and Müller, R. 1948 Formulas for bed-load transport. Proceedings, 2nd Congress International Association for Hydraulic Research, Rotterdam, the Netherlands, 39-64.