CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

K Epsilon convergance issue

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   April 18, 2011, 06:15
Default K Epsilon convergance issue
  #1
New Member
 
Join Date: Jan 2011
Posts: 6
Rep Power: 5
Ollie is on a distinguished road
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
Ollie is offline   Reply With Quote

Old   April 18, 2011, 08:11
Default
  #2
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 11
maddalena is on a distinguished road
Quote:
Originally Posted by Ollie View Post

outlet
U - fixedValue 1.461, 0, 0)
p - fixedVlaue 0
k - zerogradient
epsilon - zerogradinet
nut - calculated (value 0)
Not sure if that can be the answer, however: you should not fix both pressure and velocity at the outlet! What happens if you set the following?
Code:
outlet
U - zeroGradient or inletOutlet
p - fixedVlaue 0
k - zerogradient 
epsilon - zerogradinet
nut - calculated (value 0)
mad
maddalena is offline   Reply With Quote

Old   April 18, 2011, 08:28
Default
  #3
New Member
 
Join Date: Jan 2011
Posts: 6
Rep Power: 5
Ollie is on a distinguished road
Hi maddalena,

Firstly thanks for your quick reply. I have tried what you have suggested already, I have tried changing the boundary conditions around for a while, and the result is still the same.

Could it be the solver pisoFoam, its looks like pisoFoam is creating a solution which looks like a statistical average, instead of a transient solution ?


regards

Ollie
Ollie is offline   Reply With Quote

Reply

Tags
k-epsilon, les. rans

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 13 November 4, 2013 14:13
How to get reference to k and epsilon in the epsEqn and kEqn cfd_explorer OpenFOAM Programming & Development 0 March 10, 2011 10:16
How to get reference to k and epsilon in the epsEqn and kEqn cfd_explorer OpenFOAM 0 March 10, 2011 09:58
Compressible epsilon blows up swahono OpenFOAM 10 November 26, 2010 05:38
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 20:21


All times are GMT -4. The time now is 17:41.