CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   kLowReWallFunction in OF 2.2 (http://www.cfd-online.com/Forums/openfoam-programming-development/115862-klowrewallfunction-2-2-a.html)

 laplacian April 8, 2013 13:17

kLowReWallFunction in OF 2.2

Dear Foamers,
in the new release OF 2.2 new wall-functions are available for epsilon and K.
They seem damping-functions for both Low and High Re mesh.
My question are:
1) Which is the theoretical base for the equation of K in viscous sub-layer (i.e. where come from the numerical coefficients)?
PHP Code:

``` (kLowReWallFunctionFvPatchScalarField.C)00208             scalar C = 11.0; 00209             scalar Cf = (1.0/sqr(yPlus + C) + 2.0*yPlus/pow3(C) - 1.0/sqr(C)); 00210             kw[faceI] = 2400.0/sqr(Ceps2_)*Cf;  ```
2) Is it valid for SST k-w model (y+<5)?

 fumiya April 8, 2013 17:57

Hi,

1) You might want to read the paper(eq. (28)):
http://www.os-cfd.ru/cfd_docs/wall_f..._functions.pdf

2) I think this wall function is designed for use with v2f turbulence model.
I don't know whether or not it is valid to use this wall function with other turbulence models.
Is there anyone familiar with this area?

Hope this helps,
Fumiya

 laplacian April 9, 2013 07:32

Quote:
 Originally Posted by fumiya (Post 419208) Hi, 1) You might want to read the paper(eq. (28)): http://www.os-cfd.ru/cfd_docs/wall_f..._functions.pdf 2) I think this wall function is designed for use with v2f turbulence model. I don't know whether or not it is valid to use this wall function with other turbulence models. Is there anyone familiar with this area? Hope this helps, Fumiya
Reading the paper really helped me.

 xgarnaud September 11, 2013 04:44

I am struggling to reproduce reference results for the skin friction on a Naca airfoil (http://turbmodels.larc.nasa.gov/naca0012_val.html), so I would like to understand better the wall functions in OpenFOAM, in particular for well resolved cases.

in kLowReWallFunction, the friction velocity is computed as
scalar uTau = Cmu25*sqrt(k[faceCellI]);
for both the viscous and log layers, while in nutLowReWallFunction it is computed as