CMS-Flow:Subgrid Turbulence Model: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
(Created page with ' == Subgrid Turbulence Model == The 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 facto…')
Line 2: Line 2:
== Subgrid Turbulence Model ==
== Subgrid Turbulence Model ==


The 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 mathH_s/math  is the significant wave height and   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 (set to 0.5 here), and   is the wave bottom orbital velocity. The current-related eddy viscosity is calculated using a subgrid turbulence in which the turbulent eddy viscosity is a function of the flow gradients and the grid size model of the form  
The 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 = \del u_w h </math>, where <math>del</math> is an empirical coefficient (set to 0.5 here), and <math>u_w</math> is the wave bottom orbital velocity. The current-related eddy viscosity is calculated using a subgrid turbulence in which the turbulent eddy viscosity is a function of the flow gradients and the grid size model of the form  
(A1)
      <math> \nu_{tc} = \nu_{t0} + sqrt{(c_0 u_{*} h)^2 + (l_h^2 abs{S})^2} </math>
where  is a base value approximately equal to the dynamic viscosity, and  is an empirical coefficient and  is the subgrid mixing length. The mixing length is calculated here as  where  is an empirical coefficient (Smagorinsky coefficient). Note that if  and  are set to zero and  is set to 0.578 , where  is the bed friction coefficient, (A1) reduces to the FALCONER (1980) equation originally used in CMS-Flow.


 
where  <math>\nu_{t0}</math> is a base value approximately equal to the dynamic viscosity, and <math>c_0</mathis an empirical coefficient and  <math>l_h</math> is the subgrid mixing length. The mixing length is calculated here as  <math>l_h = \kappa min(sqrt{\del x \del y}, c_{sm} h) </math> where <math>c_{sm}</mathis an empirical coefficient (Smagorinsky coefficient).
        math \frac{\partial}{\partial t} \biggl( \frac{ h C_t }{\beta _t} \biggr) + \frac{\partial (U_j h C_t)}{\partial x_j} = \frac{\partial }{\partial x_j} \biggl[ \nu _s h \frac{\partial (r_s C_t)}{\partial x_j} \biggr] + \alpha _t \omega _s (C_{t*} - C_t) /math
 
where mathh/math is the total water depth (math h = \zeta + \eta /math), mathC_t/math is the total load concentration, mathC_{t*} /math is the sediment transport capacity, math\beta _t/math is the total load correction factor, math \nu _s /math is the diffusion coefficient, mathr_s/math is the fraction of suspended sediments, math\alpha_t/math is the total load adaptation coefficient, and math\omega_s/math is the sediment fall velocity. The concentration capacity may be calculated with either the Lund-CIRP (Carmenen and Larson 2007), the van Rijn (1985), or the Watanabe (1987) transport equations. The advantage of a bed-material approach is that the suspended- and bed-load transport equations are combined into a single equation and there is one less empirical parameter to estimate.

Revision as of 23:07, 23 October 2009

Subgrid Turbulence Model

The eddy viscosity is calculated as where is weighting factor equal to in which is the significant wave height and and are the current- and wave-related eddy viscosity components respectively. The wave contribution is included using the equation of KRAUS and LARSON (1991) Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \nu_tw = \del u_w h } , where is an empirical coefficient (set to 0.5 here), and is the wave bottom orbital velocity. The current-related eddy viscosity is calculated using a subgrid turbulence in which the turbulent eddy viscosity is a function of the flow gradients and the grid size model of the form

     

where is a base value approximately equal to the dynamic viscosity, and is an empirical coefficient and is the subgrid mixing length. The mixing length is calculated here as Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle l_h = \kappa min(sqrt{\del x \del y}, c_{sm} h) } where is an empirical coefficient (Smagorinsky coefficient).