CMS-Flow:Subgrid Turbulence Model: Difference between revisions
(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…') |
(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 | 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 | ||
(A1) | (A1) | ||
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 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. | ||
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 | 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 22:59, 23 October 2009
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 (A1) 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.
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.