CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Wiki > Yap correction

Yap correction

From CFD-Wiki

Revision as of 16:21, 18 October 2006 by Jola (Talk | contribs)
Jump to: navigation, search

The Yap correction [Yap. C. J. (1987)] consists of a modifification 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

Where

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} 
\left[
 \rho \epsilon u_j - \left( \mu + \frac{\mu_t}{\sigma_\epsilon} \right) 
 \frac{\partial \epsilon}{\partial x_j}
\right]
=
\left( C_{\epsilon_1} f_1 P - C_{\epsilon_2} f_2 \rho \epsilon \right)
\frac{\epsilon}{k}
+ \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. In an unstructured 3D solver this distance is usually not available and it can be ambiguous how to compute it in more complex topologies. This makes the Yap correction most suitable for use in a structured code where the 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)

References

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