CFD Online URL
[Sponsors]
Home > Wiki > Near-wall treatment for k-omega models

Near-wall treatment for k-omega models

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search
Line 44: Line 44:
   \omega_p = \frac{\sqrt{k_p}}{{C_\mu^{1/4}}\kappa y_p},
   \omega_p = \frac{\sqrt{k_p}}{{C_\mu^{1/4}}\kappa y_p},
</math>
</math>
-
which are the values for <math>k</math> and <math>\omega</math> according to Wilcox2006 asymptotic analysis of log layer. These wall functions for <math>k</math> and <math>\omega</math> are the results of the solution of model equation for the logarithmic layer.
+
which are the values for <math>k</math> and <math>\omega</math> according to Wilcox(1993) asymptotic analysis of log layer. These wall functions for <math>k</math> and <math>\omega</math> are the results of the solution of model equation for the logarithmic layer.
The above methodology is known to produce spurious results in separated flows, where, by definition, <math>u_\tau = 0</math> at the separation and reattachment point. Many extension of this approach has been proposed.
The above methodology is known to produce spurious results in separated flows, where, by definition, <math>u_\tau = 0</math> at the separation and reattachment point. Many extension of this approach has been proposed.
Line 91: Line 91:
== References ==
== References ==
 +
* {{reference-paper|author=Wilcox, D.|year=1993|title=Turbulence Modeling for CFD|rest='DCV Industries, Inc. La Canada, California'}}
* {{reference-paper|author=Bredberg, J.|year=2000|title=On the Wall Boundary Condition for Turbulence Models|rest='Internal Report, Department of Thermo and Fluid Dynamics, Chalmers University of Tecyhnology Gotebord, Sweden'}}
* {{reference-paper|author=Bredberg, J.|year=2000|title=On the Wall Boundary Condition for Turbulence Models|rest='Internal Report, Department of Thermo and Fluid Dynamics, Chalmers University of Tecyhnology Gotebord, Sweden'}}
* {{reference-paper|author=Menter, F., Esch, T.|year=2001|title=Elements of industrial heat transfer predictions|rest='COBEM 2001, 16th Brazilian Congress of Mechanical Engineering.'}}
* {{reference-paper|author=Menter, F., Esch, T.|year=2001|title=Elements of industrial heat transfer predictions|rest='COBEM 2001, 16th Brazilian Congress of Mechanical Engineering.'}}
* {{reference-paper|author=ANSYS|year=2006|title=FLUENT Documentation|rest=''}}
* {{reference-paper|author=ANSYS|year=2006|title=FLUENT Documentation|rest=''}}

Revision as of 13:31, 31 October 2011

Turbulence modeling
Turbulence
RANS-based turbulence models
  1. Linear eddy viscosity models
    1. Algebraic models
      1. Cebeci-Smith model
      2. Baldwin-Lomax model
      3. Johnson-King model
      4. A roughness-dependent model
    2. One equation models
      1. Prandtl's one-equation model
      2. Baldwin-Barth model
      3. Spalart-Allmaras model
    3. Two equation models
      1. k-epsilon models
        1. Standard k-epsilon model
        2. Realisable k-epsilon model
        3. RNG k-epsilon model
        4. Near-wall treatment
      2. k-omega models
        1. Wilcox's k-omega model
        2. Wilcox's modified k-omega model
        3. SST k-omega model
        4. Near-wall treatment
      3. Realisability issues
        1. Kato-Launder modification
        2. Durbin's realizability constraint
        3. Yap correction
        4. Realisability and Schwarz' inequality
  2. Nonlinear eddy viscosity models
    1. Explicit nonlinear constitutive relation
      1. Cubic k-epsilon
      2. EARSM
    2. v2-f models
      1. \overline{\upsilon^2}-f model
      2. \zeta-f model
  3. Reynolds stress model (RSM)
Large eddy simulation (LES)
  1. Smagorinsky-Lilly model
  2. Dynamic subgrid-scale model
  3. RNG-LES model
  4. Wall-adapting local eddy-viscosity (WALE) model
  5. Kinetic energy subgrid-scale model
  6. Near-wall treatment for LES models
Detached eddy simulation (DES)
Direct numerical simulation (DNS)
Turbulence near-wall modeling
Turbulence free-stream boundary conditions
  1. Turbulence intensity
  2. Turbulence length scale

Contents

Standard wall function approach

There are two possible ways of implementing wall functions in a finite volume code:

  • Additional source term in the momentum equations.
  • Modification of turbulent viscosity in cells adjacent to solid walls.

The source term in the first approach is simply the difference between logarithmic and linear interpolation of velocity gradient multiplied by viscosity (the difference between shear stresses). The second approach does not attempt to reproduce the correct velocity gradient. Instead, turbulent viscosity is modified in such a way as to guarantee the correct shear stress.

Using the compact version of log-law


\frac{U}{u_\tau} = \frac{1}{\kappa} \ln E y^{+}
(1)

where E=9.8 is equivalent to additive constants, and and using \tau_w = \rho u_\tau^2 we obtain:


  \tau_w = \frac{\rho u_\tau \kappa U}{\ln Ey^{+}},

On the other hand, the linear interpolation for shear stress, rembembering that U|_{y=0}=0, is:


  \tau_w = (\nu_t + \nu)\frac{U_p}{y_p}.

Comparing the above equations we obtain an expression for turbulent viscosity can be obtained:


  \nu_t = \nu\left( \frac{y^{+}\kappa}{\ln Ey^{+}} - 1\right).
(2)

Note that u_\tau has been been incorporated in y^+. The latter remains the only unknown in the equation and has to be estimated for the current velocity field. In the standard approach this cannot be done explicitly and instead an implicit way of obtaining y^{+} has to be employed.

After multiplying log law (1) by y_p/\nu and after reorganising some terms we get:


  \frac{\kappa U_p y_p}{\nu} = y^{+}\ln{Ey^{+}}.
(3)

This equation can be solved numerically with respect to y^+ for example via root searching algorithms e.g. Newton method for specified U_p, y_p and \nu. One iteration in a Newton method for (3) is


  y^{+}_{n+1} = \frac{\frac{\kappa U_p y_p}{\nu} + y^{+}_{n}}{1 + \ln E y_n^{+}}.

Thus obtained y^{+} is then substituted to (2) Eventually the estimated u_\tau serves also to define the values of turbulent quantities in the cell adjacent to the wall:


  k_p = \frac{u_\tau^2}{\sqrt{C_\mu}}

  \omega_p = \frac{\sqrt{k_p}}{{C_\mu^{1/4}}\kappa y_p},

which are the values for k and \omega according to Wilcox(1993) asymptotic analysis of log layer. These wall functions for k and \omega are the results of the solution of model equation for the logarithmic layer.

The above methodology is known to produce spurious results in separated flows, where, by definition, u_\tau = 0 at the separation and reattachment point. Many extension of this approach has been proposed.


Automatic wall treatments

Menter suggested a mechanism that switches automatically between HRN and LRN treatments.

The full description to appear soon. The idea is based on blending:


  \omega_\text{vis} = \frac{6\nu}{\beta y^2}

  \omega_\text{log} = \frac{u_\tau}{C_\mu^{1/4} \kappa y}

  \omega_p = \sqrt{\omega_{\text{vis}}^2 + \omega_{\text{log}}^2},



  u_\tau = \sqrt[4]{(u_\tau^{\text{vis}})^4 + (u_\tau^{\text{log}})^4},

FLUENT

Both k- omega models (std and sst) are available as low-Reynolds-number models as well as high-Reynolds-number models.

The wall boundary conditions for the k equation in the k- omega models are treated in the same way as the k equation is treated when enhanced wall treatments are used with the k- epsilon models.

This means that all boundary conditions for
- wall-function meshes will correspond to the wall function approach, while for the
- fine meshes, the appropriate low-Reynolds-number boundary conditions will be applied.

In Fluent, that means:

If the Transitional Flows option is enabled in the Viscous Model panel, low-Reynolds-number variants will be used, and, in that case, mesh guidelines should be the same as for the enhanced wall treatment
(y+ at the wall-adjacent cell should be on the order of y+ = 1. However, a higher y+ is acceptable as long as it is well inside the viscous sublayer (y+ < 4 to 5).)

If Transitional Flows option is not active, then the mesh guidelines should be the same as for the wall functions.
(For [...] wall functions, each wall-adjacent cell's centroid should be located within the log-law layer, 30 < y+ < 300. A y+ value close to the lower bound y+ = 30 is most desirable.)


References

  • Wilcox, D. (1993), "Turbulence Modeling for CFD", 'DCV Industries, Inc. La Canada, California'.
  • Bredberg, J. (2000), "On the Wall Boundary Condition for Turbulence Models", 'Internal Report, Department of Thermo and Fluid Dynamics, Chalmers University of Tecyhnology Gotebord, Sweden'.
  • Menter, F., Esch, T. (2001), "Elements of industrial heat transfer predictions", 'COBEM 2001, 16th Brazilian Congress of Mechanical Engineering.'.
  • ANSYS (2006), "FLUENT Documentation", .
My wiki