Difference between revisions of "SED TEST1 CASE"

From WikiROMS
Jump to navigationJump to search
(Replaced content with '<div class="title">Suspended Sediment Test in Channel</div> The previous contents of this page have been move to their proper location: TEST_CHAN_CASE.')   (change visibility)
 
(16 intermediate revisions by one other user not shown)
Line 1: Line 1:
<div class="title">Suspended Sediment Test in Channel</div>
<div class="title">Suspended Sediment Test in Channel</div>


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.
The previous contents of this page have been move to their proper location: [[TEST_CHAN_CASE]].
 
[[Image:Test_case_1.gif|center]]
 
==Domain==
 
The model domain is a long, narrow rectangular box, with a flat bottom.
 
{| border="1"
!Model Parameter
!Variable
!Value
|-
|Length
|<math>l</math>
|10000 m
|-
|Width
|<math>w</math>
|1000 m
|-
|Temperature
|<math>T</math>
|<math>20^{o} C</math>
|}
 
==Bottom Sediment==
 
{| border="1"
!Model Parameter
! Variable
! Value
|-
|Size
|<math>D_{50}</math>
|0.15 mm
|-
|Density
|<math>\rho_{s}</math>
|<math>2650 \,kg/m^{3}</math>
|-
|Settling Velocity
|<math>w_{s}</math>
|0.001 m/s
|-
|Critical Shear Stress
|<math>\tau_{c}</math>
|<math>0.05 \,N/m^{2}</math>
|-
|Porosity
|<math>\phi</math>
|0.90
|}
 
An infinite supply of sediment (not erosion, no armoring) is used.
 
==Forcing==
 
Inflow boundary condition only <br>
No rotation <math>(f=0)</math><br>
No wind<br>
No heating/cooling<br>
 
==Boundary Conditions==
 
Inflow maintained as steady flow, depth-mean flow <math><u> = 1 \,m/s</math> <br>
Outflow = open <br>
Sides = walls with no fluxes, no friction <br>
Bottom roughness <math>Z_{0} = 0.005</math> <br>
Sediment flux calculated by model <br>
Surface = free surface, no fluxes <br>
 
==Analysis Data==
 
Velocity profile at x = 8000 m <br>
Suspended sediment profile at x = 8000 m <br>
Time series of net rate of sediment flux from bed and/or concentration at specified levels <br>
 
==Physical Constants==
 
Gravitational acceleration <math>g \,= \,9.81 m/s^{2}</math><br>
Von Karman's constant <math>\kappa \,= \,0.41</math><br>
Water density <math>\rho_{w} \,= \,998 kg/m^{3}</math><br>
Dynamic viscosity (and minimum diffusivity)<math>\nu \,= \,1e-6 \,m^{2}/s</math><br>
 
'''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==
 
===Channel Initial Conditions===
 
Channel initial conditions:
The test channel was modeled by establishing a grid parameterized with dx = 100 m , dy=100 m, f0 = 0, and h = 10 m (flat bottom). Initial conditions set a vertical logarithmic velocity profile for u (not required but provided reasonable starting values), v = 0, zeta (water surface height) = 0, SSC in the water column = 0, and bed thickness = 1 m (to provide unlimited supply). The model was forced with 2 methods :
 
—Simulation 1: Imposing a constant flow of 10 m³/s/m of width. This simulation allowed that water surface elevation to vary. Radiation boundary conditions were imposed for the water level along with a constant flow imposed by a depth averaged velocity ubar = 1.0 m/s at the upstream and downstream boundaries.
 
—Simulation 2: Imposing a constant bed slope and water surface slope of 4x10-5 m/m. This simulation forced the water surface elevation and hence the bottom stress. Radiation boundary conditions were imposed for the depth-averaged velocity along with a clamped water surface condition at each boundary. The bed slope of 4x10-5 was selected to produce a depth-averaged velocity of 1m/s (similar to simulation 1) with a Z0 = 0.005.
 
For both simulations, vertical mixing was parameterized using six different closure schemes: MY25, KKLpara, KKLmin, KE, KW88, and ANA
MY25 = Mellor Yamada Level 2.5 Closure, parabolic wall proximity function
KKLpara = Generic Length Scale parameterized as Mellor Yamada Level 2.5 with a parabolic wall proximity function
KKLmin = Generic Length Scale parameterized as Mellor Yamada Level 2.5 with a linear wall proximity function (minimum distance to each surface)
KE = Generic Length Scale parameterized as k-epsilon closure
KW88 = Generic Length Scale parameterized as k-omega closure
ANA = analytical expression of a parabolic vertical eddy diffusivity and viscosity profile:
 
Kz = k u* z (1 - z/D)
where Kz = vertical eddy viscosity, m²/s
z = height above the bottom, m
D = depth of flow, m, = 10 m
and u* = friction velocity which was calculated according to the logarithmic profile.
 
[[Image:Equation_ustar.gif|center]]
 
where k = 0.41 and u = 1.0 m/s.
 
Only suspended-sediment transport was modeled. At the center of each bottom level horizontal grid cell, the bed shear stress was estimated and used to calculate sediment resuspension using the relation:
 
[[Image:Jw_eq_1_E.gif|center]]
 
where E = erosive flux, kg/m²/s
Eros_rate = erosion rate, kg/ m²/s, set at 5.0 e - 4
por_1 = 1-porosity = 0.6
τw = shear stress exerted by the water, = ρ Cd V |V|
ρ= fluid density, kg/m³, = 998 kg/m³
V = magnitude of velocity = sqrt(u² + v²), m/s
τc = critical shear stress for erosion, = 0.05 N/m²
The erosive flux was multiplied the bed area and by dt for each grid to provide a mass of sediment (kg) that would be exchanged between the bed and the water column. Deposition was calculated from the bottom boundary condition of the settling velocity algorithm. The depositional flux was divided by the grid horizontal area to yield a mass flux in kg. The erosive flux and depositional fluxes were added at each time step to obtain a net mass transfer that was added/subtracted from the bed and added/subtracted to the bottom cell in the water column. Bed erosion was not limited because the initial depth was set at 1.0 m and the model was not run long enough to erode to this level.
 
===Results for Simulation 1: Constant flow of <math>10 \,m^{3}/s/m</math> of width.===
 
[[Image:Sed_test1.gif|center|frame|Figure 1. Vertical profiles of velocity, eddy viscosity, eddy diffusivity, turbulent kinetic energy, dissipation, and suspended-sediment concentration (SSC) for simulation 1.]]
 
For all of the vertical profiles, the solutions for MY25 are consistent with the solutions of GLS as KKLpara. This is to be expected because there are no buoyancy terms (see test case 2 - estuary case, for differences that arise due to the length scale limitations imposed by the 2 methods). Typically, the vertical profiles from the KW88 method agree most closely with the analytical solution. The turbulent kinetic energy for the MY25 (and KKLpara) method has a lower value near the bed that is compensated for by a larger water surface slope (not shown).
 
===Results for Simulation 2: Constant bed slope and water surface slope of<math>4x10^{-5} \,m/m</math>.===
 
[[Image:Test_chan.gif|center|frame|Figure 2. Vertical profiles of velocity, eddy viscosity, eddy diffusivity, turbulent kinetic energy, dissipation, and suspended-sediment concentration (SSC) for simulation 2.]]
 
For all of the vertical profiles, the solutions for MY25 are consistent with the solutions of GLS as KKLpara (as for simulation 1). Typically, the vertical profiles from the KW88 method agree most closely with the analytical solution. The turbulent kinetic energy for all simulations are now consistent because the stress is the same for each closure.
 
===Results for open-channel flow test case.===
 
{|border="1"
!
!<math>\partial\eta\backslash\partial x</math>
!<math> u^{*} (m^{2}\backslash s^{2})</math>
|-
|calculated
|4.00e-5
|0.0625
|-
|ANA
|4.21e-5
|0.0643
|-
|<math>\kappa-\varepsilon</math>
|3.98e-5
|0.0626
|-
|MY25
|3.00e-5
|0.0544
|}

Latest revision as of 16:04, 21 September 2010

Suspended Sediment Test in Channel

The previous contents of this page have been move to their proper location: TEST_CHAN_CASE.