CFD Online Logo CFD Online URL
Home > Wiki > Yap correction

Yap correction

From CFD-Wiki

Jump to: navigation, search
Turbulence modeling
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

The Yap correction [Yap. C. J. (1987)] consists of a modification of the epsilon equation in the form of an extra source term, S_\epsilon, added to the right hand side of the epsilon equation. The source term can be written as:

\rho S_\epsilon \equiv 0.83 \, \rho \, \frac{\epsilon^2}{k} \, \left(\frac{k^{1.5}}{\epsilon \, l_e} - 1 \right) \, \left(\frac{k^{1.5}}{\epsilon \, l_e} \right)^2


l_e \equiv c_\mu^{-0.75} \, \kappa \, y_n

y_n is the normal distance to the nearest wall.

This source term should be added to the epsilon equation in the following way:

\frac{\partial}{\partial t} \left( \rho \epsilon \right) +
\frac{\partial}{\partial x_j} 
 \rho \epsilon u_j - \left( \mu + \frac{\mu_t}{\sigma_\epsilon} \right) 
 \frac{\partial \epsilon}{\partial x_j}
\left( C_{\epsilon_1} f_1 P - C_{\epsilon_2} f_2 \rho \epsilon \right)
+ \rho E
+ \rho S_\epsilon

Where the epsilon equation has been written in the same way as is in the CFD-Wiki article on low-Re k-epsilon models.

The Yap correction is active in nonequilibrium flows and tends to reduce the departure of the turbulence length scale from its local equilibrium level. It is an ad-hoc fix which seldom causes any problems and often improves the predictions.

Yap showed strongly improved results with the k-epsilon model in separated flows when using this extra source term. The Yap correction has also been shown to improve results in a stagnation region. Launder [Launder, B. E. (1993)] recommends that the Yap correction should always be used with the epsilon equation.

Implementation issues

The Yap source term contains the explicit distance to the nearest wall, y_n. This distance is sometimes difficult to efficiently calculate in complex geometries. In structured grids, the coordinate distance to the nearest wall can be used as an approximation. Otherwise, a brute force calculation must be used which greatly benefits from a multi grid approach. In topologies with domain boundaries that are not walls the problem becomes more complex, because the non-wall boundaries will block the direct path to the wall boundaries. A simple loop over length must now be accompanied by topological path checking. This makes the Yap correction most suitable for use in a structured code where some normal wall distance is readily available. There are several alternative formulations that can be used instead though (anyone have the references??).

When implementing the Yap correction it is common to use it only if the source term is positive. Hence:

\rho S_\epsilon^{implemented} = max(\rho S_\epsilon, 0)


Launder, B. E. (1993), "Modelling Convective Heat Transfer in Complex Turbulent Flows", Engineering Turbulence Modeling and Experiments 2, Proceedings of the Second International Symposium, Florence, Italy, 31 May - 2 June 1993, Edited by W. Rodi and F. Martelli, Elsevier, 1993, ISBN 0444898026.

Yap, C. J. (1987), Turbulent Heat and Momentum Transfer in Recirculating and Impinging Flows, PhD Thesis, Faculty of Technology, University of Manchester, United Kingdom.

My wiki