Sediment Transport 1: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 112: Line 112:
\frac {\partial}{\partial t}\left (\frac {hC_{tk}}{\beta_{tk}}    \right) + \frac {\partial(hV_j C_{tk})} {\partial x_j} = \frac {\partial}{\partial x_j} \left[ v_s h \frac {\partial (r_{sk}C_{tk})} {\partial x_j}  \right] + \alpha_t \omega_{sk} (C_{t*k} - C_{tk})  (2-47)
\frac {\partial}{\partial t}\left (\frac {hC_{tk}}{\beta_{tk}}    \right) + \frac {\partial(hV_j C_{tk})} {\partial x_j} = \frac {\partial}{\partial x_j} \left[ v_s h \frac {\partial (r_{sk}C_{tk})} {\partial x_j}  \right] + \alpha_t \omega_{sk} (C_{t*k} - C_{tk})  (2-47)
</math>
</math>
for j=1,2; k=1,2,......N , where N is the number of sediment size classes and
::t = time [s]
::h = water depth [m]
::<math>x_j</math> = Cartesian coordinate in the j<sup>th</sup> direction [m]
::<math>c{tk}</math> = depth-averaged total-load sediment mass concentration for size class k defined as <math>c_{tk} = q_{tk} /(Uh)</math>  in which <math>q_{tk}</math>  is the total-load mass transport [kg/m<sup>3</sup>]
::<math>c_{tk*}</math> = depth-averaged total-load sediment mass concentration for size class k and described in the Equilibrium Concentration and Transport Rates section [kg/m<sup>3</sup>]
::<math>\beta_{tk}</math>  = total-load correction factor described in the Total-Load Cor-rection Factor section [-]
::<math>r_{sk}</math>  = fraction of suspended load in total load for size class k and is described in Fraction of Suspended Sediments section [-]
::<math>v_s</math>  = horizontal sediment mixing coefficient described in the Hori-zontal Sediment Mixing Coefficient section [m<sup>2</sup>/s]
::<math>\alpha_t</math>  = total-load adaptation coefficient described in the Adaptation Coefficient section [-]
::<math>\omega_{sk}</math>  = sediment fall velocity [m/s]

Revision as of 13:16, 23 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]

The entrainment and deposition functions are calculated as

(2-44)

where

z = vertical coordinate from the bed [m]
= vertical sediment mixing coefficient [m2/s]
c = local sediment concentration [kg/m3]
= sediment fall velocity [m/s]
= calculated sediment concentration at an elevation a above the bed [kg/m3]
= equilibrium (capacity) sediment concentration at an elevation a above the bed [kg/m3]

Bed Change Equation

The bed change is calculated as

(2-46)

where

= 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) [-]
is the bed load mass transport rate [kg/m/s]

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).

Nonequilibrium Total-Load Transport Model

Total-load Transport Equation

The single-sized sediment transport model described in Sánchez and Wu (2011a) was extended to multiple-sized sediments within CMS by Sánchez and Wu (2011b). In this model, the sediment transport is separated into current- and wave-related transports. The transport due to currents includes the stirring effect of waves; and the wave-related transport includes the transport due to asymmetric oscillatory wave motion as well as steady contributions by Stokes drift, surface roller, and undertow. The current-related bed and suspended transports are combined into a single total-load transport equation, thus reducing the computational costs and simplifying the bed change computation. The two-dimensional horizontal (2DH) transport equation for the current-related total load is

for j=1,2; k=1,2,......N , where N is the number of sediment size classes and

t = time [s]
h = water depth [m]
= Cartesian coordinate in the jth direction [m]
= depth-averaged total-load sediment mass concentration for size class k defined as in which is the total-load mass transport [kg/m3]
= depth-averaged total-load sediment mass concentration for size class k and described in the Equilibrium Concentration and Transport Rates section [kg/m3]
= total-load correction factor described in the Total-Load Cor-rection Factor section [-]
= fraction of suspended load in total load for size class k and is described in Fraction of Suspended Sediments section [-]
= horizontal sediment mixing coefficient described in the Hori-zontal Sediment Mixing Coefficient section [m2/s]
= total-load adaptation coefficient described in the Adaptation Coefficient section [-]
= sediment fall velocity [m/s]