Model help:AgDegNormalFault: Difference between revisions
No edit summary |
m (Text replacement - "http://csdms.colorado.edu/wiki/" to "https://csdms.colorado.edu/wiki/") |
||
(8 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 17: | Line 17: | ||
==Model introduction== | ==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. | 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. | ||
<div id=CMT_MODEL_PARAMETERS> | <div id=CMT_MODEL_PARAMETERS> | ||
Line 164: | Line 162: | ||
==Main equations== | ==Main equations== | ||
Use the same equations as AgDegNormal model | |||
<div class="NavFrame collapsed" style="text-align:left"> | <div class="NavFrame collapsed" style="text-align:left"> | ||
Line 269: | 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 321: | 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 357: | Line 240: | ||
| U | | U | ||
| flow velocity | | flow velocity | ||
| | | L / T | ||
|- | |- | ||
| C<sub>f</sub> | | C<sub>f</sub> | ||
Line 365: | Line 248: | ||
| g | | g | ||
| acceleration of gravity | | acceleration of gravity | ||
| | | L / T<sup>2</sup> | ||
|- | |- | ||
| α<sub>r</sub> | | α<sub>r</sub> | ||
Line 373: | Line 256: | ||
| k<sub>s</sub> | | k<sub>s</sub> | ||
| grain roughness | | grain roughness | ||
| | | L | ||
|- | |- | ||
| n<sub>k</sub> | | n<sub>k</sub> | ||
Line 385: | 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 405: | 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 412: | 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 425: | 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 436: | 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 450: | 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 490: | 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 498: | Line 386: | ||
| Δη | | Δη | ||
| the height of faulting | | the height of faulting | ||
| | | L | ||
|- | |- | ||
|} | |} | ||
Line 508: | 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 534: | Line 418: | ||
</div> | </div> | ||
==Notes== | ==Notes== | ||
If the channel slope is negative and the water depth is not a number, “nan”, check the time step and the spatial step length. | 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== | ==Examples== | ||
Line 540: | 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 546: | Line 432: | ||
==Developer(s)== | ==Developer(s)== | ||
[[User:Gparker|Gary Parker]] | |||
==References== | ==References== | ||
Line 554: | 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.