CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Bounding epsilon and convergence

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

Like Tree1Likes
  • 1 Post By shyam

LinkBack Thread Tools Display Modes
Old   February 13, 2013, 06:06
Question Bounding epsilon and convergence
New Member
Join Date: Jan 2013
Posts: 4
Rep Power: 4
kirli is on a distinguished road
Hello Foamers.

I am kind of new in OpenFoam..

I am trying to run SimpleFoam turbulence simulation on a 3D geometry.
It is a kind of pipe with flow splitter. The fluid is water (one phase).

I got this kind of error:

DILUPBiCG: Solving for Ux, Initial residual = 1.96286e-14, Final residual = 1.96286e-14, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 3.56233e-14, Final residual = 3.56233e-14, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 1.38946e-13, Final residual = 1.38946e-13, No Iterations 0
DICPCG: Solving for p, Initial residual = 1, Final residual = 2.63715e+06, No Iterations 1001
time step continuity errors : sum local = 7.66448e+73, global = 1.60172e+66, cumulative = 1.60172e+66
[3] #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/"
[3] #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/"
[3] #2 in "/lib/"
[3] #3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/"

What is the reason for the negative K?
kirli is offline   Reply With Quote

Old   February 13, 2013, 08:39
Default Additional information
New Member
Join Date: Jan 2013
Posts: 4
Rep Power: 4
kirli is on a distinguished road
I believe that the problem is in the mesh.

I have made a CAD model (solidWorks) and then mesh it in salome. Import it in UNV format to OpenFoam
Attached Images
File Type: jpg Capture2.JPG (32.8 KB, 17 views)
File Type: jpg Capture1.JPG (78.0 KB, 18 views)
kirli is offline   Reply With Quote

Old   June 5, 2013, 07:29
New Member
shyam prasad
Join Date: Mar 2009
Posts: 25
Rep Power: 8
shyam is on a distinguished road
Hi Kirli,
Most of us go through this problem of bounding epsilon or omega. Some of things which helped me are the following, which are taken from various threads of CFD-online.

check your mesh using the command checkMesh and look for non-orthogonality, if its between
a. 0 to 50 full corrected scheme is applicable,
b. 50 to 70 limited correction is required,
c. 70 to 80 stability possible, accuracy compromised,
d. above 80 stability very difficult to attain

For case b use the following settings in fvSchemes

default faceLimited leastSquares 0.5;
default Gauss linear limited 0.33;
default limited 0.33;

for case c and d I would suggest to revisit your mesh.
Romarius likes this.
shyam is offline   Reply With Quote


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
groovyBC issue - k and epsilon sagnikmazumdar OpenFOAM Running, Solving & CFD 24 March 1, 2015 08:16
problem with Min/max rho tH3f0rC3 OpenFOAM 7 February 23, 2013 06:37
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24
Compressible epsilon blows up swahono OpenFOAM 10 November 26, 2010 06:38
Epsilon Convergence Trouble Carlos FLUENT 4 August 27, 2007 11:22

All times are GMT -4. The time now is 18:33.