Model help:DredgeSlotBW: Difference between revisions

From CSDMS
No edit summary
No edit summary
Line 17: Line 17:


==Model introduction==
==Model introduction==
This program calculates the 1D bed evolution of a sand-bed river after installation of a dredge slot. The calculation begins with the assumption of a prevailing mobile-bed normal flow equilibrium before installation of the dredge slot. The flow depth H, volume bedload transport rate per unit width q<sub>b</sub> and volume suspended transport rate per unit width q<sub>s</sub> at normal flow are computed based on input values of discharge Q<sub>ww</sub>, channel width B, bed material sizes D<sub>50</sub> and D<sub>90</sub>, sediment submerged specific gravity R<sub>r</sub> and bed slope S.
This program calculates the 1D bed evolution of a sand-bed river after installation of a dredge slot.  


The sediment is assumed to be sufficiently uniform so that D<sub>50</sub> and D<sub>90</sub> are unchanging in space and time. The input parameter Inter specifies the fraction of any year for which flood flow prevails. At other times of the year the river is assumed to be morphologically dormant.
The reach is assumed to have length L. The dredge slot is excavated at time t = 0, and then allowed to fill in time with no subsequent excavation. The depth of initial excavation below the bottom of the bed prevailing at normal equilibrium is an input variable with the name Hslot. The dredge slot extends from an upstream point equal to r<sub>u*L</sub> to a downstream point rd*Hslot, where ru and rd are user-input values.
The porosity lamp of the sediment deposit is a user-input parameter.
The bedload transport relation used in the calculation is that of Ashida and Michiue (1972). The formulation for entrainment of sediment into suspension is that of Wright and Parker (2004). The formulation for flow resistance is that of Wright and Parker (2004). The flow stratification correction of Wright-Parker is not implemented here for simplicity. A quasi-equilibrium formulation is used to computed the transport rate of suspended sediment from the entrainment rate.
A backwater calculation is used to compute the flow. The water surface elevation at the downstream end of the reach is held constant at the value associated with normal flow equilibrium.
Iteration is required to compute: a) the flow depth prevailing at normal flow; b) the friction slope and depth prevailing at normal flow, b) the friction slope and depth associated with skin friction associated with skin friction from any given value of depth, and b) the minimum Shields number below which form drag is taken to vanish.
==Model parameters==
==Model parameters==
= First tab header =
= First tab header =
Line 62: Line 51:


==Main equations==
==Main equations==
<span class="remove_this_tag">A list of the key equations. HTML format is supported; latex format will be supported in the future</span>
* bedload calculation using Ashida-Michiue formula
::::{|
|width=500px|<math>q_{Bj} = 17 \tau _{*j} ^ \left ( {\frac{3}{2}} \right ) \left ( {\frac{1 - \tau _{*cj}}{\tau _{*j}}} \right ) \left ( 1 - \sqrt{ {\frac{\tau _{*cj}}{\tau _{*j}}}} \right ) \ast p_{j} </math>
|width=50px align="right"|(1)
|}
* Entrainment of suspended sediment
::::{|
|width=500px|<math> E_{si} = {\frac{B \left ( \lambda X_{i} \right ) ^5 }{1 + {\frac{B}{0.3}} \left ( \lambda X_{i} \right ) ^ 5}} </math>
|width=50px align="right"|(2)
|}
::::{|
|width=500px|<math> X_{i} = \left ( {\frac{u_{*sk}}{v_{si}}} R_{pi} ^ \left (0.6\right ) \right ) S_{0} ^ \left (0.08\right ) \left ( {\frac{D_{i}}{D_{50}}}\right ) ^ \left ( 0.2 \right ) </math>
|width=50px align="right"|(3)
|}
* Overall suppression of entrainment due to mixture effects
::::{|
|width=500px|<math> \lambda = 1 - 0.28 \delta _{\Phi} </math>
|width=50px align="right"|(4)
|}
* Particle Reynolds number
::::{|
|width=500px|<math> R_{pi} = {\frac{\sqrt{ R g D_{i} }  D_{i}}{\nu}} </math>
|width=50px align="right"|(5)
|}


<div class="NavFrame collapsed" style="text-align:left">
<div class="NavFrame collapsed" style="text-align:left">
Line 70: Line 82:
!Symbol!!Description!!Unit
!Symbol!!Description!!Unit
|-
|-
| Q<sub>w</sub>
| E<sub>si</sub>
| flow discharge
| dimensionless entrainment rate of sediment in the ith grain-size range (volume entrainment per bed area per unit fractional bed content)
| m<sup>3</sup> / s
| -
|-
| λ
| overall suppression of entrainment due to mixture effects
| -
|-
| σ<sub>Φ</sub>
| standard deviation of the bed sediment on the sedimentological Φ scale
| -
|-
| u<sub>*sk</sub>
| shear velocity due to sin friction
| -
|-
|-
| I<sub>f</sub>
| R<sub>pi</sub>
| flood intermittency
| particle Reynolds number
| -
| -
|-
|-
| B
| B
| channel width
| constant in newly proposed entrainment relation, equals to 7.8*10<sup>-7</sup>
| -
|-
| S<sub>0</sub>
| the ratio of h/D<sub>50</sub>
| -
|-
| ν
| kinematic viscosity
| -
|-
| D<sub>i</sub>
| characteristic diameter of the ith grain-size range
| -
|-
| D<sub>50</sub>
| median grain diameter of the bed material
| -
|-
| τ<sub>*j</sub>
| dimensionless tractive force for the j diameter
| -
|-
| q<sub>Bj</sub>
| the amount of bed load
| -
|-
| τ<sub>*cj</sub>
| critical dimensionless tractive force for the j diameter
| -
|-
| p<sub>j</sub>
| abundance ratio of the sediment of j diameter
| -
|-
| x
| downstream coordinate
| m
|-
| η
| bed elevation
| m
|-
| q<sub>b</sub>
| volume bedload transport per unit width
| tons/annum
|-
| q<sub>s</sub>
| volume bedload transport per unit width due to skin friction
| tons/annum
|-
| H
| water depth
| m
| m
|-
|-
| D<sub>50</sub>
| H<sub>s</sub>
| median grain size (sand)
| water depth due to skin friction
| mm
| m
|-
| k<sub>si</sub>
| water surface elevation
| m
|-
|-
| D<sub>90</sub>
| Q<sub>w</sub>
| grain size such that 90% is finer (sand)
| flow discharge
| mm
| m<sup>3</sup> / s
|-
|-
| R
| I<sub>f</sub>
| submerged specific gravity of sediment
| flood intermittency
| -
| -
|-
| R<sub>r</sub>
| submerged specific gravity
|
|-
|-
| S
| S
Line 99: Line 183:
|-
|-
| L
| L
| Reach length
| reach length
| m
| m
|-
|-
Line 106: Line 190:
| m
| m
|-
|-
| r<sub>u</sub>
| r<sub>U</sub>
| fraction of reach length defining upstream end of dredge slot
| fraction of reach length for upstream end of dredge slot
| cm / s
| -
|-
|-
| r<sub>d</sub>
| r<sub>d</sub>
| fraction of reach length defining downstream end of dredge slot
| fraction of reach length for downstream end of dredge slot
|  
| -
|-
|-
| λ<sub>p</sub>
| λ<sub>p</sub>
| bed porosity
| bed porosity
| -
|-
| a<sub>U</sub>
| upwinding coefficient: a<sub>u</sub> = 1 corresponds to full upwinding, 0.5 for central difference
| -
| -
|-
|-
| M
| M
| number of intervals
| number of spatial steps (<=2000)
| -
| -
|-
|-
| ∆<sub>x</sub>
| t
| spatial step length
| time step
| m
| yr
|-
|-
| ∆<sub>t</sub>
| Δt
| time step
| time step
| year
| -
|-
| Δx
| spatial step length
| -
|-
|-
| Mtoprint
| Mtoprint
| number of time steps to printout
| number of time steps to printout
| -
| -
|-
|-
| Mprint
| Mprint
| number of printouts
| number of printouts
| -
| -
|-
|-
| a<sub>u</sub>
| upwinding coefficient: a<sub>u</sub> = 1 corresponds to full upwinding
| -
|- 
|}
|}


Line 150: Line 238:
In the calculation of River Bed Elevation Variation with a Dredge Slot: the river is assumed to be sand-bed. The calculation proceeds using a backwater formulation. Flow resistance is computed using the Wright-Parker (2004) formulation. The bedload transport rate is computed using the Ashida-Michiue (1972) formulation. The rate of entrainment into suspension is computed using the Wright-Parker formulation without the stratification correction.
In the calculation of River Bed Elevation Variation with a Dredge Slot: the river is assumed to be sand-bed. The calculation proceeds using a backwater formulation. Flow resistance is computed using the Wright-Parker (2004) formulation. The bedload transport rate is computed using the Ashida-Michiue (1972) formulation. The rate of entrainment into suspension is computed using the Wright-Parker formulation without the stratification correction.


The calculation begins with the assumption of a prevailing mobile-bed normal flow equilibrium before installation of the dredge slot. The flow depth H, volume bedload transport rate per unit width q<sub>b</sub> and volume suspended transport rate per unit width q<sub>s</sub> at normal flow are computed based on input values of discharge Q<sub>ww</sub>, channel width B, bed material sizes D<sub>50</sub> and D<sub>90</sub>, sediment submerged specific gravity R<sub>r</sub> and bed slope S.
The sediment is assumed to be sufficiently uniform so that D<sub>50</sub> and D<sub>90</sub> are unchanging in space and time. The input parameter Inter specifies the fraction of any year for which flood flow prevails. At other times of the year the river is assumed to be morphologically dormant.
The reach is assumed to have length L. The dredge slot is excavated at time t = 0, and then allowed to fill in time with no subsequent excavation. The depth of initial excavation below the bottom of the bed prevailing at normal equilibrium is an input variable with the name Hslot. The dredge slot extends from an upstream point equal to r<sub>u*L</sub> to a downstream point rd*Hslot, where ru and rd are user-input values.
The porosity lamp of the sediment deposit is a user-input parameter.
The bedload transport relation used in the calculation is that of Ashida and Michiue (1972). The formulation for entrainment of sediment into suspension is that of Wright and Parker (2004). The formulation for flow resistance is that of Wright and Parker (2004). The flow stratification correction of Wright-Parker is not implemented here for simplicity. A quasi-equilibrium formulation is used to computed the transport rate of suspended sediment from the entrainment rate.
A backwater calculation is used to compute the flow. The water surface elevation at the downstream end of the reach is held constant at the value associated with normal flow equilibrium.
Iteration is required to compute: a) the flow depth prevailing at normal flow; b) the friction slope and depth prevailing at normal flow, b) the friction slope and depth associated with skin friction associated with skin friction from any given value of depth, and b) the minimum Shields number below which form drag is taken to vanish.
* Note on model running
The model is allowed up to 2000 spatial steps, any more than that will cause a memory overwrite and the data cannot be trusted
The Normal flow data is included in the initial output, at the bottom of the file.
The sediment is assumed to be sufficiently uniform such that D<sub>50</sub> and D<sub>90</sub> do not change in space or time.
The bedload transport calculations use an Ashida-Michiue formulation (1972), the entrainment of suspended sediment uses that of Wright and Parker (2004), as does the formulation for flow resistance.


==Examples==
==Examples==
Line 164: Line 274:


==References==
==References==
<span class="remove_this_tag">Key papers</span>
* Ashida, K. and M. Michiue, 1972, Study on hydraulic resistance and bedload transport rate in alluvial streams, Transactions, Japan Society of Civil Engineering, 206: 59~69 (in Japanese)
* Wright, S. and G. Parker, 2004, Flow resistance and suspended load in sand-bed rivers: simplified stratification model, Journal of Hydraulic Engineering, 130(8): 796~805.


==Links==
==Links==

Revision as of 12:36, 27 April 2011

The CSDMS Help System

DredgeSlotBW

This model is a calculator for aggradation and degradation of sediment mixtures in gravel-bed streams subject to cyclic hydrographs.

Model introduction

This program calculates the 1D bed evolution of a sand-bed river after installation of a dredge slot.

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

  • bedload calculation using Ashida-Michiue formula
[math]\displaystyle{ q_{Bj} = 17 \tau _{*j} ^ \left ( {\frac{3}{2}} \right ) \left ( {\frac{1 - \tau _{*cj}}{\tau _{*j}}} \right ) \left ( 1 - \sqrt{ {\frac{\tau _{*cj}}{\tau _{*j}}}} \right ) \ast p_{j} }[/math] (1)
  • Entrainment of suspended sediment
[math]\displaystyle{ E_{si} = {\frac{B \left ( \lambda X_{i} \right ) ^5 }{1 + {\frac{B}{0.3}} \left ( \lambda X_{i} \right ) ^ 5}} }[/math] (2)
[math]\displaystyle{ X_{i} = \left ( {\frac{u_{*sk}}{v_{si}}} R_{pi} ^ \left (0.6\right ) \right ) S_{0} ^ \left (0.08\right ) \left ( {\frac{D_{i}}{D_{50}}}\right ) ^ \left ( 0.2 \right ) }[/math] (3)
  • Overall suppression of entrainment due to mixture effects
[math]\displaystyle{ \lambda = 1 - 0.28 \delta _{\Phi} }[/math] (4)
  • Particle Reynolds number
[math]\displaystyle{ R_{pi} = {\frac{\sqrt{ R g D_{i} } D_{i}}{\nu}} }[/math] (5)

Notes

In the calculation of River Bed Elevation Variation with a Dredge Slot: the river is assumed to be sand-bed. The calculation proceeds using a backwater formulation. Flow resistance is computed using the Wright-Parker (2004) formulation. The bedload transport rate is computed using the Ashida-Michiue (1972) formulation. The rate of entrainment into suspension is computed using the Wright-Parker formulation without the stratification correction.

The calculation begins with the assumption of a prevailing mobile-bed normal flow equilibrium before installation of the dredge slot. The flow depth H, volume bedload transport rate per unit width qb and volume suspended transport rate per unit width qs at normal flow are computed based on input values of discharge Qww, channel width B, bed material sizes D50 and D90, sediment submerged specific gravity Rr and bed slope S.

The sediment is assumed to be sufficiently uniform so that D50 and D90 are unchanging in space and time. The input parameter Inter specifies the fraction of any year for which flood flow prevails. At other times of the year the river is assumed to be morphologically dormant.

The reach is assumed to have length L. The dredge slot is excavated at time t = 0, and then allowed to fill in time with no subsequent excavation. The depth of initial excavation below the bottom of the bed prevailing at normal equilibrium is an input variable with the name Hslot. The dredge slot extends from an upstream point equal to ru*L to a downstream point rd*Hslot, where ru and rd are user-input values.

The porosity lamp of the sediment deposit is a user-input parameter.

The bedload transport relation used in the calculation is that of Ashida and Michiue (1972). The formulation for entrainment of sediment into suspension is that of Wright and Parker (2004). The formulation for flow resistance is that of Wright and Parker (2004). The flow stratification correction of Wright-Parker is not implemented here for simplicity. A quasi-equilibrium formulation is used to computed the transport rate of suspended sediment from the entrainment rate.

A backwater calculation is used to compute the flow. The water surface elevation at the downstream end of the reach is held constant at the value associated with normal flow equilibrium.

Iteration is required to compute: a) the flow depth prevailing at normal flow; b) the friction slope and depth prevailing at normal flow, b) the friction slope and depth associated with skin friction associated with skin friction from any given value of depth, and b) the minimum Shields number below which form drag is taken to vanish.

  • Note on model running

The model is allowed up to 2000 spatial steps, any more than that will cause a memory overwrite and the data cannot be trusted

The Normal flow data is included in the initial output, at the bottom of the file.

The sediment is assumed to be sufficiently uniform such that D50 and D90 do not change in space or time.

The bedload transport calculations use an Ashida-Michiue formulation (1972), the entrainment of suspended sediment uses that of Wright and Parker (2004), as does the formulation for flow resistance.

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

  • Ashida, K. and M. Michiue, 1972, Study on hydraulic resistance and bedload transport rate in alluvial streams, Transactions, Japan Society of Civil Engineering, 206: 59~69 (in Japanese)
  • Wright, S. and G. Parker, 2004, Flow resistance and suspended load in sand-bed rivers: simplified stratification model, Journal of Hydraulic Engineering, 130(8): 796~805.

Links

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