CMS-Flow:Subgrid Turbulence Model

Subgrid Turbulence Model

In CMS-Flow eddy viscosity is calculated as the sum of the kinematic viscosity ${\nu }_0$, the current-related eddy viscosity ${\nu }_c$ and the wave-related eddy viscosity ${\nu }_w$

     ${\nu }_t={\nu }_0+{\nu }_c+{\nu }_w$  

The wave contribution is included using the equation of Kraus and Larson (1991)

     ${\nu }_w=\mathrm{\Lambda }{u}_{w}h$

where $\mathrm{\Lambda }$ is an empirical coefficient (default is 0.5), and $u_w$ is the wave bottom orbital velocity and $h$ is the water depth.

The current-related eddy viscosity is calculated as a function of the flow gradients, and the bottom shear stress

     ${\nu }_c={\nu }_0+\sqrt{(c_0{u}_{*}h{)}^2+(c_1\mathrm{\Delta }|S|{)}^2}$

where ${\nu }_0$ is a base value approximately equal to the kinematic viscosity, $c_0$ is an empirical coefficient, $c_{sm}$ is an empirical coefficient (Smagorinsky coefficient), $\mathrm{\Delta }$ is the local cell area, and $|S|$ is equal to

     $|S|=\sqrt{(\frac{\partial U}{\partial x}{)}^2+(\frac{\partial V}{\partial y}{)}^2+\frac{1}{2}(\frac{\partial U}{\partial y}+\frac{\partial V}{\partial x}{)}^2}$

---

References

LARSON, M.; HANSON, H., and KRAUS, N. C., 2003. Numerical modeling of beach topography change. Advances in Coastal Modeling, V.C. Lakhan (eds.), Elsevier Oceanography Series, 67, Amsterdam, The Netherlands, 337-365.

---

CMS-Flow