[romant]

Hej,

I looked through the different wall function for k-epsilon turbulence models, however I feel I am missing something. In Wilcox (Turbulence modeling for CFD 2006) as well as many other publications, in wall functions the velocity seems to be set in the first grid point. In OpenFOAM, the following properties are changed with wall functions in the example of the k-epsilon turbulence model,

1. epsilon: based on the logarithmic law of the wall, epsilon is set in the first cell node (epsilonWallFunction)
2. G: the turbulence production term is determined from the law of the wall and also set in the first node from the wall (epsilonWallFunction)
3. k: is only set to zeroGradient (kqRWallFunction)
4. velocity is not touched in any of these functions?
5. the turbulent viscosity (nut, mut) are set based on the law of the wall (mutUWallFunction)

The question is, why?

[romant]

Hej,

I think I found the solution to this. I found it in Pope, Turbulent Flow, 2000 under the way that wall functions are created. The original question was why the velocity is not set by any of the wall function. The answer is, that wall functions are apparently more robust by setting the shear stress and not the velocity, which makes sense when looking at separation points and other nasty flow phenomena, in which the velocity can become essentially 0. The way this is for the k-epsilon model done is

nutkWallFunction

1. determine the nominal value for the shear velocity, u_shear
2. find y+ using the shear velocity
3. find the nominal values for the velocity in the first grid point (based on the log law in this case), U_p,nominal
4. set the shear stresses -<u_i u_j> (not the normal stresses in the Reynolds stresses) to u_shear * U_p,real / U_p,nominal , In OF this happens through the turbulent viscosity (nu_t)

[cfdonline2mohsen]

Thanks roman
it was very helpful.

[immortality]

hi
Whould you explain more and Why stage 4 is done?

[romant]

Originally Posted by immortality

hi
Whould you explain more and Why stage 4 is done?

Step 4 is introduced in order to have a more stable solver. Instead of setting the velocity in the first point away from the wall, a force (the wall shear force) is introduced in the equation. This way, the solver won't have any 0 velocities introduced by accident (think of stagnation points), this could destabilize the solver as you would be introducing 0's into the matrix.

The easiest way to introduce this source term into the equation is by modifying the turbulent viscosity, see Pope, Turbulent Flows, 2000, pp 442-443 and Versteeg and Malalasekera, An introduction to computational fluid dynamics: The finite volume method, pp. 276-277

[huangxianbei]

Originally Posted by romant

Hej,

1. set the shear stresses -<u_i u_j> (not the normal stresses in the Reynolds stresses) to u_shear * U_p,real / U_p,nominal , In OF this happens through the turbulent viscosity (nu_t)

Hi,romant:
Thank you for your explain. I have a question about the wall function for non-linear models. As the non-linear model can't be depicted by a simple turbulent viscosity, how can we set the shear stress?
say we obtained the shear stress as in your form: u_shear * U_p,real / U_p,nominal, so we should set the Reynolds shear stress by u_shear * U_p,real / U_p,nominal - nu*dU/dy, here, y represents the wall-normal direction. However, the non-linear Reynolds stress is not symmetric, it has 6 shear stress components, which one should be modified? I'm confused by this problem for some days. I'm looking forward for your reply.

Thank you

Xianbei

[zoptirik]

Originally Posted by huangxianbei

Hi,romant:
Thank you for your explain. I have a question about the wall function for non-linear models. As the non-linear model can't be depicted by a simple turbulent viscosity, how can we set the shear stress?
say we obtained the shear stress as in your form: u_shear * U_p,real / U_p,nominal, so we should set the Reynolds shear stress by u_shear * U_p,real / U_p,nominal - nu*dU/dy, here, y represents the wall-normal direction. However, the non-linear Reynolds stress is not symmetric, it has 6 shear stress components, which one should be modified? I'm confused by this problem for some days. I'm looking forward for your reply.

Thank you

Xianbei

Did you solve your wall function problem. I have tried to implement a non-linear model and I have stuck the same problem. What kind of wall function should I use.

thank you