[gfilip]

pisoFoam crashes when I try to use the locDynOneEqEddy LES model.

libOpenFOAM.so and libincompressibleRASModels.so both show up in the attached log.

When I run the case with say dynSmagorinsky, everything works fine.


in my constant/LESProperties:

***************************************
LESModel locDynOneEqEddy
delta cubeRootVol;

locDynOneEqEddyCoeffs
{
ce 1.05;
filter simple;
}
****************************************

I have no issues with this turb model in a channel flow problem with both pisoFoam and channelFoam. I have checked the initial conditions several times but I cannot find anything that could be responsible for crash in this case.

Any hints are appreciated.

[Tarak]

Hii,

You can try to solve this in 2 ways:

1. In the wall boundary condition for turbulent kinetic energy, give a very small value e.g 1e-20 instead of 0.

OR

2. Start solution using oneEqEddy or dynOneEqEddy model (or any other simpler model) and after some hundreds of iterations switch to locDynOneEqEddy model.

Thanks,
Tarak

[vut]

Hi Tarak,

Your response for gfilip's question interested me alot.

May you give me more details of switching to locDynOneEqEddy model after hundreds of iterations?

More generally, I am searching how to initialize a LES simulation from a known solution?

Thank in advance,

VUT