Model help:AgDegNormGravMixSubPW

From CSDMS
Revision as of 18:13, 12 April 2011 by Happyfei (talk | contribs)
The CSDMS Help System

AgDegNormGravMixSubPW

It is the calculator for evolution of upward-concave bed profiles in rivers carrying sediment mixtures in subsiding basins.

Model introduction

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.

Model parameters

Parameter Description Unit
First parameter Description parameter [Units]
Parameter Description Unit
First parameter Description parameter [Units]

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

A list of the key equations. HTML format is supported; latex format will be supported in the future

Notes

The river is assumed to be morphologically active for If fraction of time, during which the flow is approximated as constant. Otherwise, the river is assumed to be morphologically dead.

The river flows into a basin that is subsiding with rate s. The basin has constant width; the ratio of basin width to river width is rB. The river has sinuosity W. For each unit of bedload deposited, L units of washload (typically sand transported in suspension) is deposited across the depositional basin.

In particular, the program computes the time evolution of the spatial profiles of bed elevation, bed slope, total bedload transport rate and grain size distribution of the surface (active) layer of the bed.

If run for a sufficient length of time, the river profile approaches a steady-state balance between subsidence. At this steady state the profile displays both an upward-concave elevation profile and downstream fining of the surface material.

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.

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 ηd.

Gravel bedload transport of mixtures is computed with a user-specified selection of the Parker (1990), or Wilcock-Crowe (2003) surface-based formulations for gravel transport.Sand and finer material must first be excluded from the grain size distributions, which then must be renormalized for gravel content only, in the case of the Parker (1990) relation. In the case of the Wilcock-Crowe (2003) relation, the sand is retained in the computation.

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. These constraints are easy to relax.

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.

The flow is calculated using the normal flow (local equilibrium) approximation.

In performing the calculation, the following control parameters must be specified: M = number of spatial intervals, so that the spatial step length = L/M; dt = time step length; Ntoprint = number of time steps to a printout; Nprint = number of printouts in the calculation.

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

Key papers

Links