CMS-Flow:Subgrid Turbulence Model

Subgrid Turbulence Model

In CMS-Flow eddy viscosity is calculated as $\nu _{t}=(1-\theta _{m})\nu _{tc}+\theta _{m}\nu _{m}$ where $\theta _{m}$ is weighting factor equal to $\theta _{m}=(H_{s}/h)^{3}$ in which $H_{s}$ is the significant wave height and $\nu _{tc}$ and $\nu _{tw}$ are the current- and wave-related eddy viscosity components respectively. The wave contribution is included using the equation of Kraus and Larson (1991) $\nu _{t}w=\Lambda u_{w}h$, where $\Lambda$ is an empirical coefficient (default is 0.5), and $u_{w}$ is the wave bottom orbital velocity. The current-related eddy viscosity is calculated as a function of the flow gradients, and the bottom shear stress

     $\nu _{tc}=\nu _{t0}+{\sqrt {(c_{0}u_{*}h)^{2}+(l_{h}^{2}|S|)^{2}}}$

where $\nu _{t0}$ is a base value approximately equal to the dynamic viscosity, and $c_{0}$ is an empirical coefficient and $l_{h}$ is the subgrid mixing length. The mixing length is calculated here as

     $l_{h}=\mathrm {K} min({\sqrt {\Delta x\Delta y}},c_{sm}h)$ 

where $c_{sm}$ is an empirical coefficient (Smagorinsky coefficient).

---

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