Standard k-epsilon model
From CFD-Wiki
(Difference between revisions)
m (fixed case on headings) |
|||
| Line 1: | Line 1: | ||
| - | == Transport | + | == Transport equations for standard k-epsilon model == |
For k <br> | For k <br> | ||
| Line 31: | Line 31: | ||
</math> | </math> | ||
| - | == Effect of | + | == Effect of buoyancy == |
:<math> | :<math> | ||
| Line 46: | Line 46: | ||
</math> | </math> | ||
| - | == Model | + | == Model constants == |
:<math> | :<math> | ||
C_{1 \epsilon} = 1.44, \;\; C_{2 \epsilon} = 1.92, \;\; C_{\mu} = 0.09, \;\; \sigma_k = 1.0, \;\; \sigma_{\epsilon} = 1.3 | C_{1 \epsilon} = 1.44, \;\; C_{2 \epsilon} = 1.92, \;\; C_{\mu} = 0.09, \;\; \sigma_k = 1.0, \;\; \sigma_{\epsilon} = 1.3 | ||
</math> | </math> | ||
Revision as of 15:49, 25 May 2006
Contents |
Transport equations for standard k-epsilon model
For k
![\frac{\partial}{\partial t} (\rho k) + \frac{\partial}{\partial x_i} (\rho k u_i) = \frac{\partial}{\partial x_j} \left[ \left(\mu + \frac{\mu_t}{\sigma_k} \right) \frac{\partial k}{\partial x_j}\right] + P_k + P_b - \rho \epsilon - Y_M + S_k](/W/images/math/0/3/3/03313e8802538459d0a202c34efc1274.png)
For dissipation 
![\frac{\partial}{\partial t} (\rho \epsilon) + \frac{\partial}{\partial x_i} (\rho \epsilon u_i) = \frac{\partial}{\partial x_j} \left[\left(\mu + \frac{\mu_t}{\sigma_{\epsilon}} \right) \frac{\partial \epsilon}{\partial x_j} \right] + C_{1 \epsilon}\frac{\epsilon}{k} \left( P_k + C_{3 \epsilon} P_b \right) - C_{2 \epsilon} \rho \frac{\epsilon^2}{k} + S_{\epsilon}](/W/images/math/0/6/c/06ca9efd2f29b3816707de0452572c77.png)
Modeling turbulent viscosity
Turbulent viscosity is modelled as:
Production of k
Where 
is the modulus of the mean rate-of-strain tensor, defined as :
Effect of buoyancy
where Prt is the turbulent Prandtl number for energy and gi is the component of the gravitational vector in the ith direction. For the standard and realizable - models, the default value of Prt is 0.85.
The coefficient of thermal expansion, 
, is defined as
