|
[Sponsors] |
April 18, 2011, 07:15 |
K Epsilon convergance issue
|
#1 |
New Member
Join Date: Jan 2011
Posts: 6
Rep Power: 15 |
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 |
|
April 18, 2011, 09:11 |
|
#2 | |
Senior Member
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 23 |
Quote:
Code:
outlet U - zeroGradient or inletOutlet p - fixedVlaue 0 k - zerogradient epsilon - zerogradinet nut - calculated (value 0) |
||
April 18, 2011, 09:28 |
|
#3 |
New Member
Join Date: Jan 2011
Posts: 6
Rep Power: 15 |
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 |
|
Tags |
k-epsilon, les. rans |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
SimpleFoam k and epsilon bounded | nedved | OpenFOAM Running, Solving & CFD | 16 | March 4, 2017 09:30 |
How to get reference to k and epsilon in the epsEqn and kEqn | cfd_explorer | OpenFOAM Programming & Development | 0 | March 10, 2011 11:16 |
How to get reference to k and epsilon in the epsEqn and kEqn | cfd_explorer | OpenFOAM | 0 | March 10, 2011 10:58 |
Compressible epsilon blows up | swahono | OpenFOAM | 10 | November 26, 2010 06:38 |
SimpleFoam k and epsilon bounded | nedved | OpenFOAM Running, Solving & CFD | 1 | November 25, 2008 21:21 |