K-epsilon Model

From Thermal-FluidsPedia

(Difference between revisions)
Jump to: navigation, search
(Created page with 'In the turbulent models presented in the preceding section, the eddy diffusivity, <math>\varepsilon _{M}</math>, was obtained by algebraic expressions. While they are very easy t…')
Line 18: Line 18:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(2)}}
|}
|}
where <math>C_{\mu }</math> is a dimensionless empirical constant and L plays the same role as the mean free path in the kinematic viscosity of the gas. Therefore, both K and L must be determined in order to use eq. (4.436) to evaluate the eddy diffusivity.  
where <math>C_{\mu }</math> is a dimensionless empirical constant and L plays the same role as the mean free path in the kinematic viscosity of the gas. Therefore, both K and L must be determined in order to use eq. (4.436) to evaluate the eddy diffusivity.  
Line 33: Line 33:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(3)}}
|}
|}
For a two-dimensional steady state flow, the convection of K is
For a two-dimensional steady state flow, the convection of K is
Line 43: Line 43:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(4)}}
|}
|}
Since eddy diffusion in the x-direction is negligible for a boundary layer type flow, only diffusion in the y-direction needs to be considered:
Since eddy diffusion in the x-direction is negligible for a boundary layer type flow, only diffusion in the y-direction needs to be considered:
Line 53: Line 53:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(5)}}
|}
|}
where <math>\sigma _{K}</math> is a dimensionless empirical constant. The rate of K production can be obtained by the eddy shear stress and the time averaged velocity gradient
where <math>\sigma _{K}</math> is a dimensionless empirical constant. The rate of K production can be obtained by the eddy shear stress and the time averaged velocity gradient
Line 63: Line 63:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(6)}}
|}
|}
where eq. (4.436) was used to evaluate the eddy diffusivity. The rate of K destruction, <math>\varepsilon </math>, can be evaluated by analyzing a fluid packet with diameter L, oscillating in the turbulent flow field. If the oscillating velocity is K1/2, the drag force acting on the fluid packet will be <math>C_{D}\rho L^{2}(K^{1/2})^{2}</math>, where CD is the drag coefficient that is approximately equal to 1. The mechanical power dissipated per unit mass is
where eq. (4.436) was used to evaluate the eddy diffusivity. The rate of K destruction, <math>\varepsilon </math>, can be evaluated by analyzing a fluid packet with diameter L, oscillating in the turbulent flow field. If the oscillating velocity is K1/2, the drag force acting on the fluid packet will be <math>C_{D}\rho L^{2}(K^{1/2})^{2}</math>, where CD is the drag coefficient that is approximately equal to 1. The mechanical power dissipated per unit mass is
Line 73: Line 73:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(7)}}
|}
|}
The K-equation for a boundary layer type flow can be obtained by substituting eqs. (4.438) – (4.441) into eq. (4.437)
The K-equation for a boundary layer type flow can be obtained by substituting eqs. (4.438) – (4.441) into eq. (4.437)
Line 83: Line 83:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(8)}}
|}
|}
By following a similar procedure, the dissipation equation can be obtained as the following:
By following a similar procedure, the dissipation equation can be obtained as the following:
Line 93: Line 93:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(9)}}
|}
|}
where <math>\sigma _{\varepsilon },\text{ }C_{1},\text{ and }C_{2}</math>are additional dimensional empirical constants. After K and ε are obtained, the eddy diffusivity can be obtained by eliminating L between eqs. (4.436) and (4.441), i.e.
where <math>\sigma _{\varepsilon },\text{ }C_{1},\text{ and }C_{2}</math>are additional dimensional empirical constants. After K and ε are obtained, the eddy diffusivity can be obtained by eliminating L between eqs. (4.436) and (4.441), i.e.
Line 103: Line 103:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(10)}}
|}
|}
where CD in eq. (4.441) has been set to 1. Equations (4.442) – (4.444) become the equations for the K-ε model. It is different from the algebraic equations presented in the preceding subsection because additional partial differential equations must be solved. Similar to eq. (4.407), both eqs. (4.442) and (4.443) are parabolic because diffusion in the x-direction is neglected. Jones and Launder (1972) suggested the following values for the K-ε model:
where CD in eq. (4.441) has been set to 1. Equations (4.442) – (4.444) become the equations for the K-ε model. It is different from the algebraic equations presented in the preceding subsection because additional partial differential equations must be solved. Similar to eq. (4.407), both eqs. (4.442) and (4.443) are parabolic because diffusion in the x-direction is neglected. Jones and Launder (1972) suggested the following values for the K-ε model:
Line 117: Line 117:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(11)}}
|}
|}
The generation of K represented by the second term on the right-hand side of eq. (4.442) becomes
The generation of K represented by the second term on the right-hand side of eq. (4.442) becomes
Line 135: Line 135:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(12)}}
|}
|}
If the length scale is taken as
If the length scale is taken as
Line 149: Line 149:
</center>
</center>
-
|{{EquationRef|(1)}}
+
|{{EquationRef|(13)}}
|}
|}
where the left and right hand sides, respectively, represent the rates of production and dissipation of K, i.e. generation of K is balanced by dissipation of K in the viscous sublayer. Therefore, the mixing length theory can be used in the sublayer region and the fully turbulent region can be described using the K-ε model.
where the left and right hand sides, respectively, represent the rates of production and dissipation of K, i.e. generation of K is balanced by dissipation of K in the viscous sublayer. Therefore, the mixing length theory can be used in the sublayer region and the fully turbulent region can be described using the K-ε model.

Revision as of 04:27, 9 April 2010

In the turbulent models presented in the preceding section, the eddy diffusivity, \varepsilon _{M}, was obtained by algebraic expressions. While they are very easy to apply, their shortcoming is their lack of universal applicability because all of them work well only near the wall. As evidenced by Figs. 4.35 and 4.36, the agreement between the wall function and the experimental results is not very good for large y+. When the above algebraic models are applied to the internal flow, they will break down near the center of the tube. Thus, it is imperative to develop a turbulent model that is applicable in the region that is sufficiently far away from the wall. An ideal turbulent model should be applicable to any region in the turbulent flow although it may not be able to be represented by algebraic equations. Significant advancement of high-speed computers makes developments and applications of non-algebraic turbulent models possible. Among different non-algebraic turbulent models, the K-ε model is one of the most widely used models and it will be discussed below. The K-ε model is based on the analogy between the motion of the fluid packet in the turbulent flow and the random motion of an ideal gas. The kinetic theory of gases stated that the kinematic viscosity of gas can be obtained by ν = cλ / 3 where c is the mean speed of the molecules and λ is the mean free path. For turbulent flow, the kinetic energy due to the fluid packet random motion is

K=\frac{1}{2}(\overline{{u}'^{2}}+\overline{{v}'^{2}}+\overline{{w}'^{2}})

(1)

and the mean speed of the fluid packet is K1/2. The eddy diffusivity for the momentum can be defined in a similar way to the kinematic viscosity of gas (Kolmogorov, 1942), i.e.,

\varepsilon _{M}=C_{\mu }K^{1/2}L

(2)

where Cμ is a dimensionless empirical constant and L plays the same role as the mean free path in the kinematic viscosity of the gas. Therefore, both K and L must be determined in order to use eq. (4.436) to evaluate the eddy diffusivity. For three-dimensional flow, the equation for kinetic energy associated with the velocity fluctuation was derived in Section 3.5.2. While we can obtain the equation for K by simplifying eq. (2.528) for two-dimensional turbulent boundary flow (Oosthuizen and Naylor, 1999), a more direct approach suggested by Bejan (2004) will be introduced here. If we consider a control volume in a horizontal slender flow region (turbulent boundary layer), balance of the kinetic energy for velocity fluctuation is



Convection of K = Eddy diffusion of K + Rate of K generation − Rate of K destruction

(3)

For a two-dimensional steady state flow, the convection of K is

\text{Convection of }K\text{ = }\bar{u}\frac{\partial K}{\partial x}+\bar{v}\frac{\partial K}{\partial y}

(4)

Since eddy diffusion in the x-direction is negligible for a boundary layer type flow, only diffusion in the y-direction needs to be considered:

\text{ Eddy diffusion of }K=\frac{\partial }{\partial y}\left( \frac{\varepsilon _{M}}{\sigma _{k}}\frac{\partial K}{\partial y} \right)

(5)

where σK is a dimensionless empirical constant. The rate of K production can be obtained by the eddy shear stress and the time averaged velocity gradient

\text{Rate of }K\text{ generation }=\left( \varepsilon _{M}\frac{\partial \bar{u}}{\partial y} \right)\cdot \left( \frac{\partial \bar{u}}{\partial y} \right)=C_{\mu }K^{1/2}L\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}

(6)

where eq. (4.436) was used to evaluate the eddy diffusivity. The rate of K destruction, \varepsilon , can be evaluated by analyzing a fluid packet with diameter L, oscillating in the turbulent flow field. If the oscillating velocity is K1/2, the drag force acting on the fluid packet will be CDρL2(K1 / 2)2, where CD is the drag coefficient that is approximately equal to 1. The mechanical power dissipated per unit mass is

\varepsilon =\frac{[C_{D}\rho L^{2}(K^{1/2})^{2}]K^{1/2}}{\rho L^{3}}=C_{D}\frac{K^{3/2}}{L}

(7)

The K-equation for a boundary layer type flow can be obtained by substituting eqs. (4.438) – (4.441) into eq. (4.437)

\text{ }\bar{u}\frac{\partial K}{\partial x}+\bar{v}\frac{\partial K}{\partial y}\text{ }=\frac{\partial }{\partial y}\left( \frac{\varepsilon _{M}}{\sigma _{k}}\frac{\partial K}{\partial y} \right)+\varepsilon _{M}\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}-\varepsilon

(8)

By following a similar procedure, the dissipation equation can be obtained as the following:

\text{ }\bar{u}\frac{\partial \varepsilon }{\partial x}+\bar{v}\frac{\partial \varepsilon }{\partial y}\text{ }=\frac{\partial }{\partial y}\left( \frac{\varepsilon _{M}}{\sigma _{\varepsilon }}\frac{\partial \varepsilon }{\partial y} \right)+C_{1}\varepsilon _{M}\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}\frac{\varepsilon }{K}-C_{2}\frac{\varepsilon ^{2}}{K}

(9)

where \sigma _{\varepsilon },\text{ }C_{1},\text{ and }C_{2}are additional dimensional empirical constants. After K and ε are obtained, the eddy diffusivity can be obtained by eliminating L between eqs. (4.436) and (4.441), i.e.

\varepsilon _{M}=C_{\mu }\frac{K^{2}}{\varepsilon }

(10)

where CD in eq. (4.441) has been set to 1. Equations (4.442) – (4.444) become the equations for the K-ε model. It is different from the algebraic equations presented in the preceding subsection because additional partial differential equations must be solved. Similar to eq. (4.407), both eqs. (4.442) and (4.443) are parabolic because diffusion in the x-direction is neglected. Jones and Launder (1972) suggested the following values for the K-ε model:

C_{\mu }=0.09,\text{ }C_{1}=1.44,\text{ }C_{2}=1.92,\text{ }\sigma _{K}=1,\text{ and }\sigma _{\varepsilon }=1.3

Since the effect of molecular viscosity, ν, was neglected in eqs. (4.442) and (4.443), the K-ε model cannot be applied in the viscous sublayer where molecular viscosity dominates the velocity profiles. In the sublayer region, Prandtl’s mixing length theory is valid and eq. (4.413) can be rewritten as

\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}=\frac{\varepsilon _{M}^{2}}{l^{4}}

(11)

The generation of K represented by the second term on the right-hand side of eq. (4.442) becomes

\varepsilon _{M}\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}=\frac{\varepsilon _{M}^{3}}{l^{4}}

Substituting eq. (4.436) into the right-hand side of the above equation, one obtains:

\varepsilon _{M}\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}=\frac{C_{\mu }^{3}K^{3/2}L^{3}}{l^{4}}

Considering eq. (4.441), the above equation can be rewritten as

\varepsilon _{M}\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}=\frac{C_{\mu }^{3}}{C_{D}}\left( \frac{L}{l} \right)^{4}\varepsilon

(12)

If the length scale is taken as

L=l\left( \frac{C_{D}}{C_{\mu }^{3}} \right)^{1/4}

eq. (4.446) becomes

\varepsilon _{M}\left( \frac{\partial \bar{u}}{\partial y} \right)^{2}=\varepsilon

(13)

where the left and right hand sides, respectively, represent the rates of production and dissipation of K, i.e. generation of K is balanced by dissipation of K in the viscous sublayer. Therefore, the mixing length theory can be used in the sublayer region and the fully turbulent region can be described using the K-ε model.