[Ollie]

Hi OpenFoamers,

I am having some trouble with K Epsilon when using PisoFOAM, I am initialising in SimpleFOAM and taking the final time step obtained from SimpleFOAM to be used as the 0 file for PisoFOAM.

The flow type is high Reynolds number flow over a cylinder within a control volume.

The flow properly initialises in simpleFOAM and converges as expcted. However after 10,000 iterations in pisoFoam there is little difference if any in the flow structure, a steady-state results appears to be gained. I am using an Euler scheme and have ran for 10 seconds. There is sort of instantaneous convergance.

The expected flow structure should be completely turbulent, with eddies being shed alternately from the cylinder, however a pointed "tail" remains in the cylinder wake.

I have ran a similar model in LES with no complications, thus justifing the mesh is of a high standard.

Has anyone else had this problem, if so, how did you resolve it?


My initial/boundary conditions are below.

internal Field
U - (1.461, 0, 0)
p - 0
epsilon - 0.024
k - 0.020491
nut - 0

inlet
U - fixed value (1.461, 0, 0)
p - zerogradient
k - fixed value 0.020491
epsilon - fixed value 0.024
nut - calculated 0

Cylinder
U - fixed value - (0,0,0)
p - zerogradient
k - kqRWallFunction (value 0.020491)
epsilon - epsilonWallFunction (value 0.024)
nut- nutWallFunction (value 0)

outlet
U - fixedValue 1.461, 0, 0)
p - fixedVlaue 0
k - zerogradient
epsilon - zerogradinet
nut - calculated (value 0)

Cheers
Ollie