CMS-Flow:Subgrid Turbulence Model

From CIRPwiki
Jump to navigation Jump to search

Subgrid Turbulence Model

In CMS-Flow eddy viscosity is calculated as

     νt=(1θm)νc+θmνw  

where θm is weighting factor equal to θm=(Hs/h)3 in which Hs is the significant wave height and νc and νw are the current- and wave-related eddy viscosity components respectively. The wave contribution is included using the equation of Kraus and Larson (1991)

     νw=Λuwh

where Λ is an empirical coefficient (default is 0.5), and uw 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

     νc=ν0+(c0u*h)2+(csmΔ|S|)2

where ν0 is a base value approximately equal to the kinematic viscosity, c0 is an empirical coefficient, csm is an empirical coefficient (Smagorinsky coefficient), Δ is the local cell area, and |S| is equal to

     |S|=2(Ux)2+2(Vy)2+(Uy+Vx)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