Sediment Transport 1: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 53: Line 53:


::<math>\frac {\partial (hC)}{\partial t} + \frac {\partial (hV_j C)}{\partial x_j} = \frac{\partial}{\partial x_j} \left(v_s h \frac{\partial C}{\partial x_j}\right) + E_b - D_b </math>  (2-43)
::<math>\frac {\partial (hC)}{\partial t} + \frac {\partial (hV_j C)}{\partial x_j} = \frac{\partial}{\partial x_j} \left(v_s h \frac{\partial C}{\partial x_j}\right) + E_b - D_b </math>  (2-43)
where
::t = time[s]
::h = water depth [m]
::<math>x_j</math> = Cartesian coordinate in the
::<math>E_b</math> entrainment or pick-up function [kg/m<sup>2</sup>/s]
::<math>D_b</math> deposition or settling function [kg/m<sup>2</sup>/s]

Revision as of 19:53, 22 July 2014

Sediment Transport

Overview

For sand transport, the wash-load (i.e. sediment transport which does not contribute to the bed-material) can be assumed to be zero, and therefore, the total-load transport is equal to the sum of the bed- and the suspended-load transports: .

There are currently three sediment transport models available in CMS:

(1) Equilibrium total load

(2) Equilibrium bed load plus non-equilibrium suspended load, and

(3) Non-equilibrium total-load.

The first two models are single-size sediment transport models and are only available with the explicit time-stepping schemes. The third is multiple-sized sediment transport model and is available with both the explicit and implicit time-stepping schemes.

Equilibrium Total-load Transport Model

In this model, both the bed load and suspended load are assumed to be in equilibrium. The bed change is solved using a simple mass balance equation known as the Exner equation.

(2-42)

for , where N is the number of sediment size classes and

t = time [s]


h = total water depth [m]
=Cartesian coordinate in the jth direction [m]
= equilibrium total-load transport rate [kg/m/s]
= bed elevation with respect to the vertical datum [m]
= bed porosity [-]
= morphologic acceleration factor [-]
= sediment density [~2650 kg/m3 for quartz sediment]
= empirical bed-slope coefficient (constant) [-]

Because the model assumes that both the sediment transport is equilibri-um, it only recommended for coarse grids with resolutions larger than 50-100 m where the assumption of equilibrium sediment transport is more appropriate. As mentioned above the equilibrium total-load sediment transport model is a single-size sediment transport model and is only available with the explicit time-stepping scheme. For more information on the equilibrium sediment transport model, the reader is referred to Buttolph et al. (2007).

Equilibrium Bed-load plus Nonequilibrium Suspended Load Transport Model

Calculations of suspended load and bed load are conducted separately. The bed load is assumed to be in equilibrium and is included in the bed change equation while the suspended load is solved through the solution of an advection-diffusion equation. Actually the advection diffusion equation is a non-equilibrium formulation, but because the bed load is assumed to be in equilibrium, this model is referred to the "Equilibrium A-D" model.

Suspended-load Transport Equation

The transport equation for the suspended load is given by

(2-43)

where

t = time[s]
h = water depth [m]
= Cartesian coordinate in the
entrainment or pick-up function [kg/m2/s]
deposition or settling function [kg/m2/s]