CMS-Flow:Subgrid Turbulence Model: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
m (Reverted edits by Alex (Talk) to last revision by Mitchell brown)
No edit summary
Line 2: Line 2:
== Subgrid Turbulence Model ==
== Subgrid Turbulence Model ==


In CMS-Flow eddy viscosity is calculated as <math> \nu_t = (1-\theta_m)\nu_{tc} + \theta_m \nu_m </math>  where <math>\theta_m</math>  is weighting factor equal to <math>\theta_m = (H_s/h)^3 </math>  in which <math>H_s</math>  is the significant wave height and <math>\nu_{tc}</math>  and <math>\nu_{tw}</math>  are the current- and wave-related eddy viscosity components respectively. The wave contribution is included using the equation of Kraus and Larson (1991)  <math> \nu_tw = \Lambda u_w h </math>, where  <math>\Lambda</math> is an empirical coefficient (default is 0.5), and  <math>u_w</math> 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
In CMS-Flow eddy viscosity is calculated as  


       <math> \nu_{tc} = \nu_{t0} + \sqrt{ (c_0 u_* h)^2 + (\sqrt{c_{sm} \Delta x \Delta y}  |S| )^2 } </math>
       <math> \nu_t = \nu_c + \nu_w </math>


where <math>\nu_{t0}</math> is a base value approximately equal to the dynamic viscosity, and <math>c_0</math> is an empirical coefficient and <math>l_h</math> is the subgrid mixing length. The mixing length is calculated here as 
where <math>\nu_c</math> and <math>\nu_w</math> are the current- and wave-related eddy viscosity components respectively. The wave contribution is included using the equation of Kraus and Larson (1991)


       <math>l_h = c_{sm} \Delta x \Delta y </math>  
       <math> \nu_w = \Lambda u_w h </math>
 
where  <math>\Lambda</math> is an empirical coefficient (default is 0.5), and  <math>u_w</math> is the wave bottom orbital velocity and <math>h</math> is the water depth. The current-related eddy viscosity is calculated as a function of the flow gradients, and the bottom shear stress
 
      <math> \nu_{c} = \nu_{0} + \sqrt{ (c_0 u_* h)^2 + (c_{sm}\Delta |S| )^2 } </math>
 
where  <math>\nu_{0}</math> is a base value approximately equal to the kinematic viscosity, <math>c_0</math> is an empirical coefficient, <math>c_{sm}</math>  is an empirical coefficient (Smagorinsky coefficient), <math> \Delta </math> is the local cell area, and <math>|S|</math> is equal to
 
      <math> |S| = \sqrt{ \biggl( \frac{ \partial U}{\partial x} \biggr) ^2 +  \biggl( \frac{ \partial V}{\partial y} \biggr) ^2 + \frac{1}{2} \biggl( \frac{ \partial U}{\partial y} + \frac{ \partial V}{\partial x}  \biggr) ^2 } </math>


where <math>c_{sm}</math>  is an empirical coefficient (Smagorinsky coefficient).


----
----

Revision as of 18:39, 5 May 2010

Subgrid Turbulence Model

In CMS-Flow eddy viscosity is calculated as

       

where and are the current- and wave-related eddy viscosity components respectively. The wave contribution is included using the equation of Kraus and Larson (1991)

     

where is an empirical coefficient (default is 0.5), and is the wave bottom orbital velocity and is the water depth. The current-related eddy viscosity is calculated as a function of the flow gradients, and the bottom shear stress

     

where is a base value approximately equal to the kinematic viscosity, is an empirical coefficient, is an empirical coefficient (Smagorinsky coefficient), is the local cell area, and is equal to

     



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