Model help:AgDegNormalFault: Difference between revisions

From CSDMS
No edit summary
m (Text replacement - "http://csdms.colorado.edu/wiki/" to "https://csdms.colorado.edu/wiki/")
 
(12 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:
   * http://csdms.colorado.edu/wiki/Model help:<modelname>
   * 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.
 
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.  


<div id=CMT_MODEL_PARAMETERS>
<div id=CMT_MODEL_PARAMETERS>
==Model parameters==
==Model parameters==
= First tab header =
= Input Files and Directories =
{|{{Prettytable}} class = "wikitable unsortable"  cellspacing="0" cellpadding="0" style="margin:0em 0em 0em 0;"
{|{{Prettytable}} class = "wikitable unsortable"  cellspacing="0" cellpadding="0" style="margin:0em 0em 0em 0;"
|-
|-
!Parameter!!Description!!Unit
!Parameter!!Description!!Unit
|-valign="top"
|-valign="top"
|width="20%"|<span class="remove_this_tag">First parameter</span>
|width="20%"|Input directory
|width="60%"|<span class="remove_this_tag">Description parameter</span>
|width="60%"|path to input files
|width="20%"|<span class="remove_this_tag">[Units]</span>
|width="20%"|
|-
|Site prefix
|Site prefix for Input/Output files
|
|-
|Case prefix
|Case prefix for Input/Output files
|
|-
|}
|}


= Second tab header =
= Run Parameters =
{|{{Prettytable}} class = "wikitable unsortable"  cellspacing="0"  cellpadding="0" style="margin:0em 0em 0em 0;"
{|{{Prettytable}} class = "wikitable unsortable"  cellspacing="0"  cellpadding="0" style="margin:0em 0em 0em 0;"
|-
|-
!Parameter!!Description!!Unit
!Parameter!!Description!!Unit
|-valign="top"
|-valign="top"
|width="20%"|<span class="remove_this_tag">First parameter</span>
|width="20%"|Flood discharge
|width="60%"|<span class="remove_this_tag">Description parameter</span>
|width="60%"|
|width="20%"|<span class="remove_this_tag">[Units]</span>
|width="20%"| m<sup>3</sup> / s
|-
|Intermittency
|
| -
|-
|Channel Width
|
| m
|-
|Grain size
|
| mm
|-
|Bed Porosity
|
| -
|-
|Roughness height
|
| mm
|-
|Ambient Bed Slope
|
|
|-
|Imposed Annual Sediment Transfer Rate from Upstream
|
| tons / annum
|-
|Length of reach
|
| m
|-
|Time step
|
| year
|-
|Number of Time Steps per Printout
|
|
|-
|Number of printout
|
|
|-
|intervals
|
|
|-
|Upwinding coefficient (1 = full upwind, 0.5 = central difference)
|
|
|-
|Coefficient in Manning-Strickler Resistance Relation
|
|
|-
|Coefficient in Sediment Transport Relation
|
|
|-
|Exponent in Sediment Transport Relation
|
|
|-
|Critical Shield stress
|
|
|-
|Fraction of bed shear stress that is skin friction
|
|
|-
|Submerged specific gravity of sediment
|
|
|-
|Height of faulting
|
| m
|-
|Fraction of reach length such that all points downstream undergo downward faulting
|
| -
|-
|Time from beginning of calculation at which faulting occurs
|
| yrs
|-
|}
|}


= Etc. tab header =
= About =
{|{{Prettytable}} class = "wikitable unsortable"  cellspacing="0"  cellpadding="0" style="margin:0em 0em 0em 0;"
|-
!Parameter!!Description!!Unit
|-valign="top"
|width="20%"|Model name
|width="60%"|name of the model
|width="20%"| -
|-
|Author name
|name of the model author
| -
|-
|}
<headertabs/>
<headertabs/>
</div>


==Uses ports==
==Uses ports==
Line 54: Line 162:


==Main equations==
==Main equations==
* Manning-Strickler formulation
Use the same equations as AgDegNormal model
::::{|
|width=500px|<math>C_{f}=\alpha _{r}\left ( \frac{H}{K_{c}} \right )^{\frac{1}{6}}</math>
|width=50px align="right"|(1)
|}
* Total bed material load per unit width
::::{|
|width=500px|<math>\frac{q_{t}}{{\sqrt{RgDD}}}=\alpha_{t}\left ( \frac{\varphi_{s}\tau _{b}}{\rho RgD} -\tau_{c}^* \right )^{n_{t}}</math>
|width=50px align=right|(2)
|}


<div class="NavFrame collapsed" style="text-align:left">
<div class="NavFrame collapsed" style="text-align:left">
Line 71: Line 170:
!Symbol!!Description!!Unit
!Symbol!!Description!!Unit
|-
|-
| X
| Q
| Streamwise coordinate
| flood discharge
| m
| L <sup>3</sup> / T
|-
| t
| time step
| T
|-
| B
| river width
| L
|-
| D
| grain size of the bed sediment
| L
|-
| λ<sub>p</sub>
| bed porosity
| -
|-
| q<sub>w</sub>
| water discharge per unit width
| L<sup>2</sup> / T
|-
| k<sub>c</sub>
| composite roughness height
| L
|-
| G
| imposed annual sediment transfer rate from upstream
| M / T
|-
| G<sub>tf</sub>
| upstream sediment feed rate
| -
|-
| ξ<sub>d</sub>
| downstream water surface elevation
| L
|-
|-
| ΔX
| L
| Spatial step length
| length of reach under consideration
| m
| L
|-
|-
| t
| i
| Temporal coordinate
| number of time steps per printout
| seconds
| -
|-
|-
| C<sub>f</sub>
| p
| Non-dimensional friction coefficient
| number of printouts desired
| -
| -
|-
|-
| Q<sub>w</sub>
| M
| Flood discharge
| number of spatial intervals
| m<sup>3</sup>/s
| -
|-
|-
| I<sub>f</sub>
| R
| Flood intermittency
| submerged specific gravity of sediment
| -
| -
|-
|-
| B<sub>c</sub>
| S<sub>f</sub>
| Channel width
| friction slope
| m
| -
|-
|-
| D
| F<sub>r</sub>
| Characteristic grain size
| Froude number
| mm
|-
| λ<sub>p</sub>
| Bed porosity
| -
| -
|-
|-
| k<sub>c</sub>
| U
| Composite roughness height
| flow velocity
| mm
| L / T
|-
|-
| S<sub>I</sub>
| C<sub>f</sub>
| Ambient bed slope
| bed friction coefficient
| -
| -
|-
|-
| G<sub>tf</sub>
| g
| Imposed annual sediment transport rate
| acceleration of gravity
| tons/annum
| L / T<sup>2</sup>
|-
| α<sub>r</sub>
| coefficient in Manning-Stricker, dimensionless coefficient between 8 and 9
| -
|-
|-
| k<sub>s</sub>
| grain roughness
| L
| L
| Length of reach
|- 
| m
| n<sub>k</sub>
| dimensionless coefficient typically between 2 and 5
| -
|- 
| τ<sup>*</sup>
| Shield number
| -
|-
| ρ
| fluid density
| M / L<sup>3</sup>
|-
| ρ<sub>s</sub>
| sediment density
| M / L<sup>3</sup>
|-
| τ<sub>c</sub>
| critical Shields number for the onset of sediment motion
| -
|-
|-
| Δt
| ψ<sub>s</sub>
| Time step
| the fraction of bed shear stress
| year
| -
|-
|-
| N<sub>toprint</sub>
| q<sub>t</sub> <sup>*</sup>
| number of time steps to printout
| Einstein number
| -
| -
|-
|-
| N<sub>print</sub>
| q<sub>t</sub>
| number of printouts
| volume sediment transport rate per unit width
| L<sup>2</sup> / T
|-
| I<sub>f</sub>
| flood intermittency
| -
| -
|-
|-
| M
| t<sub>f</sub>
| Number of spatial intervals
| cumulative time the river has been in flood
| T
|-
| G<sub>t</sub>
| the annual sediment yield
| M / T
|-
| t<sub>a</sub>
| the number of seconds in a year
| -
| -
|-
|-
| a<sub>u</sub>
| Q<sub>f</sub>
| Upwinding coefficient (1 = full upwind, 0.5 = central difference)
| sediment transport rate during flood discharge
| L<sup>2</sup> / T
|-
| α<sub>t</sub>
| dimensionless coefficient in the sediment transport equation, equals to 8
| -
| -
|-
|-
| α<sub>r</sub>
| n<sub>t</sub>
| Coefficient in Manning-Strickler
| exponent in sediment transport relation, equals to 1.5
| -
| -
|-
|-
| α<sub>s</sub>
| τ<sub>c</sub> <sup>*</sup>
| Coefficient in sediment transport relation
| reference Shields number in sediment transport relation, equals to 0.047
|-
| C<sub>Z</sub>
| dimensionless Chezy resistance coefficient.
|-
| S<sub>l</sub>
| initial bed slope of the river
| -
| -
|-
|-
| η<sub>t</sub>
| η<sub>i</sub>
| Exponent in sediment transport relation
| initial bed elevation
| L
|-
| τ
| shear stress on bed surface
| -
| -
|-
|-
| τ<sup>*</sup><sub>c</sub>
| q<sub>b</sub>
| Reference Shields number in sediment transport relation
| bed material load
| M / T
|-
| Δx
| spatial step length, equals to L / M
| L
|-
| Q<sub>w</sub>
| flood discharge
| L<sup>3</sup> / T
|-
| Δt
| time step
| T
|-
| Ntoprint
| number of time steps to printout
| -
| -
|-
|-
| φ<sub>s</sub>
| Nprint
| Fraction of bed shear stress due to skin friction
| number of printouts
| -
| -
|-
|-
| R
| a<sub>U</sub>
| Submerged specific gravity
| upwinding coefficient (1=full upwind, 0.5=central difference)
| -
| -
|-
|-
| Cz
| α<sub>s</sub>
| Non-dimensional Chézy friction coefficient
| coefficient in sediment transport relation
| -
| -
|-
|-
| G<sub>tf</sub>
| u<sub>*</sub>
| imposed annual sediment transport rate
| shear velocity
| L / T
|-
| α<sub>r</sub>
| coefficient in Manning-Strickler resistance relation
| -
| -
|-
|-
| r<sub>f</sub>
| τ<sub>b</sub> <sup>*</sup>
| the fraction of reach length such that all points downstream of x = r<sub>f</sub>L undergo downward faulting
| non-dimensional total shear stress
| -
| -
|-
|-
| t<sub>f</sub>
| r<sub>f</sub>
| the time from beginning of calculation at which faulting occurs
| the fraction of reach length such that all point downstream of x = r<sub>f</sub>L undergo downward faulting
| -
| -
|-
|-
| ∆η
| Δη
| height of faulting
| the height of faulting
| m
| L
|-
|-
|}
|}


'''Output'''
'''Output'''
Line 195: Line 395:
|-
|-
| η
| η
| Bed surface elevation
| bed surface elevatioon
| m
| L
|-
|-
| S
| H
| Bed slope
| water depth
| -
| L
|-
|-
| H
| ξ
| Water depth
| water surface elevation
| m
| L
|-
|-
| τ<sub>b</sub>
| τ<sub>b</sub>
| Total (skin friction + form drag) Shields number
| bed shear stress
| M / (T<sup>2</sup> L)
|-
| S
| bed slope
| -
| -
|-
| q<sub>t</sub>
| total bed material load
| m<sup>2</sup>/s
|-
|-
|}
|}
Line 218: 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. 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
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.  


The water depth is calculated using a Chézy formulation, when only the Chézy coefficient is specified in the input text file. The Manning-Strickler formulation is implemented, when only the roughness height, k<sub>c</sub>, and the coefficient α<sub>r</sub> are given in the input file. When all the three parameters are present, the program will ask the user which formulation they would like to use.
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 226: 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: http://csdms.colorado.edu/wiki/Special:Upload</span>
* <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 232: Line 432:


==Developer(s)==
==Developer(s)==
<span class="remove_this_tag">Name of the module developer(s)</span>
[[User:Gparker|Gary Parker]]


==References==
==References==
<span class="remove_this_tag">Key papers</span>
* 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.


==Links==
==Links==
<span class="remove_this_tag">Any link, eg. to the model questionnaire, etc.</span>
* [[Model:AgDegNormalFault]]
* [[Model_help:AgDegNormal]]


[[Category:Modules]] [[Category:Utility components]]
[[Category:Utility components]]

Latest revision as of 17:17, 19 February 2018

The CSDMS Help System

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

Parameter Description Unit
Input directory path to input files
Site prefix Site prefix for Input/Output files
Case prefix Case prefix for Input/Output files
Parameter Description Unit
Flood discharge m3 / s
Intermittency -
Channel Width m
Grain size mm
Bed Porosity -
Roughness height mm
Ambient Bed Slope
Imposed Annual Sediment Transfer Rate from Upstream tons / annum
Length of reach m
Time step year
Number of Time Steps per Printout
Number of printout
intervals
Upwinding coefficient (1 = full upwind, 0.5 = central difference)
Coefficient in Manning-Strickler Resistance Relation
Coefficient in Sediment Transport Relation
Exponent in Sediment Transport Relation
Critical Shield stress
Fraction of bed shear stress that is skin friction
Submerged specific gravity of sediment
Height of faulting m
Fraction of reach length such that all points downstream undergo downward faulting -
Time from beginning of calculation at which faulting occurs yrs
Parameter Description Unit
Model name name of the model -
Author name name of the model author -

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

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:

See also: Help:Images or Help:Movies

Developer(s)

Gary Parker

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.

Links

Retrieved from "https://csdms.colorado.edu/csdms_wiki/index.php?title=Model_help:AgDegNormalFault&oldid=216486"