SED TEST1 CASE

Suspended Sediment Test in Channel

This case provides a fundamental check of the ability of a model to 1) represent a simple flow, 2) flux material from the bed, and 3) develop a suspended-sediment profile.

Domain

The model domain is a long, narrow rectangular box, with a flat bottom.

Model Parameter	Variable	Value
Length	$l$	10000 m
Width	$w$	1000 m
Temperature	$T$	$20^{o}C$

Bottom Sediment

Model Parameter	Variable	Value
Size	$D_{50}$	0.15 mm
Density	$\rho _{s}$	$2650\,kg/m^{3}$
Settling Velocity	$w_{s}$	0.001 m/s
Critical Shear Stress	$\tau _{c}$	$0.05\,N/m^{2}$
Porosity	$\phi$	0.90

An infinite supply of sediment (not erosion, no armoring) is used.

Forcing

Inflow boundary condition only
No rotation $(f=0)$
No wind
No heating/cooling

Boundary Conditions

Inflow maintained as steady flow, depth-mean flow $<u>=1\,m/s$
Outflow = open
Sides = walls with no fluxes, no friction
Bottom roughness $Z_{0}=0.005$
Sediment flux calculated by model
Surface = free surface, no fluxes

Analysis Data

Velocity profile at x = 8000 m
Suspended sediment profile at x = 8000 m
Time series of net rate of sediment flux from bed and/or concentration at specified levels

Physical Constants

Gravitational acceleration $g\,=\,9.81m/s^{2}$
Von Karman's constant $\kappa \,=\,0.41$
Water density $\rho _{w}\,=\,998kg/m^{3}$
Dynamic viscosity (and minimum diffusivity)$\nu \,=\,1e-6\,m^{2}/s$

Note: If a model incorporates physical constants that differ from these, and/or automatically calculates some values specified here, please specify the values used.

Results

	$\partial \eta \backslash \partial x$	$u^{*}(m^{2}\backslash s^{2})$
calculated	4.00e-5	0.0625
ANA	4.21e-5	0.0643
$\kappa -\varepsilon$	3.98e-5	0.0626
MY25	3.00e-5	0.0544