[em17]

Hi all,

I am using the WALE SGS model and have a query about the model constants. I would like to change the constant Cw, however, when I looked at the source code I found that the model is also dependent on constants Ck and Ce. It seems that Ce is only used in the calculation for epsilon, so I am not too bothered about this.

I am presuming if I want to change the value of Cw, then I must also change Ck? Does anyone know the relationship between these two constants?

Thank you!
Em

[navidmt]

Hi Emily

Have you found any answer for your question? I came up with the same question and I don't know if I need to run my simulation using the default values of OpenFoam or it is essential to change the values of Cw and consequently the Ck.

I would appreciate if you share your possible solution here with me. Thank you

[clapointe]

If you look at the publication described in the .H file, the authors suggest a "universal" Cw of somewhere between 0.55 and 0.6 if I remember correctly. I am not sure where the default coefficient came from, but you can change it by adding to turbulenceProperties :

WALECoeffs
{
    Cw 0.55;
}

As for Ck, they have an expression for nut in the paper that is some function of Cw and the velocity gradient tensor. If you look at the foam implementation and do some math, you'll likely find that nut is computed such that Ck cancels out.

Caelan

[em17]

Hi,

Caelan is correct! I ended up working through the equations in the source code (WALE.C) to figure it out.

Ck (which is actually the Smagorinsky constant) is used to calculate the subgrid-scale turbulent kinetic energy, . This is then used to calculate the turbulent eddy viscosity, , which is written in your time directories. During the calculation of , the Ck terms will cancel out. So if you are only interested in then you only need to update Cw. However if you are interested in the subgrid-scale turbulent kinetic energy , then you'll need to estimate Ck as well.

You'll also see that Ce is simply used to calculate epsilon, so if you are not interested in epsilon then Ce is irrelevant.

Also, based on my personal experience I have found the OpenFOAM 'standard' constants to work well (Cw=0.325). Although if you have the time/resources available then it's always good to test different values to see how sensitive your simulations are.

I hope this helps,
Emily

[navidmt]

Thank you so much I followed as you suggested.