Channel flow - Rouse profile

In this scenario, the water depth and the surface slope have been chosen identical to those in the purely physical channel scenario introduced in section 12.1.2.

Under certain conditions, the suspended matter equation

$\displaystyle \partial_t C + \partial_z\left(w_c C - \nu_t \partial_z C\right) = 0,$ (273)

has an analytical solution:

Let the eddy diffusivity profile be parabolic with

$\displaystyle \nu_t = \kappa u_* (-z)\frac{D+z_0+z}{D+z_0}$ (274)

with the depth $ D$, the bottom roughness length $ z_0$, the van Karman number $ \kappa$ and the bottom friction velocity $ u_*$ and the vertical coordinate $ z$. Then the analytical solution of (273) is

$\displaystyle \displaystyle \frac{C}{C_0} = \left(\frac{-z}{D+z_0+z} \right)^{-w_c/(\kappa u_*)},$ (275)

where $ C_0$ is the suspended matter concentration at $ z=-(D+z_0)/2$i and depends on the initial conditions for $ C$. The Rouse number is then defined as:

$\displaystyle R =\frac{-w_c}{u_*}.$ (276)

When running the rouse scenario with a two-equation turbulence closure model, then the analytical solution for the Rouse profile is only approximated, since the eddy diffusivity deviates from (275). In order to be closer to the analytical solution, it is ncessary to chose in gotmturb.nml the analytical parabolic profile for eddy diffusivity, which is done by the follwing settings in gotmturb.nml:

  turb_method=     2,
  tke_method=      1,
  stab_method=     1

This Rouse scenario may be calculated with Eulerian concentrations or with Lagrangian particles, depending on the setting of bio_eulerian in bio.nml. When using the Lagrangian particle method, it is advisable to average the concentration over all time steps which belong to one output time step, by setting bio_lagrange_mean=.true.

It is also possible to include a sink term at the bed (for simulating the effect of grazing by benthic filter feeder), the seetings for this have to be made in mussels.nml.

Karsten Bolding 2012-01-24