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

bounding epsilon blow up

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

LinkBack Thread Tools Display Modes
Old   September 30, 2010, 11:42
Default bounding epsilon blow up
New Member
Jie (Jay) Zhang
Join Date: Sep 2010
Location: FL, U.S.
Posts: 25
Rep Power: 8
jiez is on a distinguished road
Sponsored Links
It is a case using simpleFoam solver with K-epsilon Turbulence model.

the fellowing result will blow up. Is there anyone can help me to find out where the problem comes from? BC setting('epsilo' attached)?

Time = 0.001
DILUPBiCG: Solving for Ux, Initial residual = 0.64936, Final residual = 0.0285001, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.527993, Final residual = 0.0462129, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.815427, Final residual = 0.0168725, No Iterations 2
DICPCG: Solving for p, Initial residual = 1.42064e-15, Final residual = 1.42064e-15, No Iterations 0
time step continuity errors : sum local = 6.86127e+27, global = 7.64441e+11, cumulative = -2.2062e+23
DILUPBiCG: Solving for epsilon, Initial residual = 4.35678e-07, Final residual = 4.35678e-07, No Iterations 0
bounding epsilon, min: 3.3635e-22 max: 1.18516e+87 average: 7.77262e+80
DILUPBiCG: Solving for k, Initial residual = 5.98618e-09, Final residual = 5.98618e-09, No Iterations 0
ExecutionTime = 193.75 s ClockTime = 195 s

BC setting:


type epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.0012;

type fixedValue;
value uniform 0.0012;

type zeroGradient;
jiez is offline   Reply With Quote
Sponsored Links

Old   September 30, 2010, 11:46
New Member
Jie (Jay) Zhang
Join Date: Sep 2010
Location: FL, U.S.
Posts: 25
Rep Power: 8
jiez is on a distinguished road
I'm new to FOAM.
I will be very appreciated for your help!
jiez is offline   Reply With Quote

Old   September 30, 2010, 20:07
Senior Member
santiagomarquezd's Avatar
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 430
Rep Power: 17
santiagomarquezd will become famous soon enough
Maybe a draft of the geometry, details of problem and


dictionaries can give us a better idea how to solve the blowing up.

Research Scientist
Research Center for Computational Methods (CIMEC) - CONICET/UNL
Tel: 54-342-4511594 Int. 7032
Colectora Ruta Nac. 168 / Paraje El Pozo
(3000) Santa Fe - Argentina.
santiagomarquezd is offline   Reply With Quote

Old   October 8, 2010, 17:11
Senior Member
Travis Carrigan
Join Date: Jul 2010
Location: Arlington, TX
Posts: 147
Rep Power: 9
tcarrigan is on a distinguished road
Try setting the k and e schemes to upwind rather than linear in the fvSchemes directory. I've had problems achieving convergence using second order schemes for the turbulence models.
tcarrigan is offline   Reply With Quote

Old   January 16, 2011, 18:14
New Member
Philipp Bachmann
Join Date: May 2010
Location: Esslingen, Germany
Posts: 7
Rep Power: 9
kroetenechse is on a distinguished road
In my case it was a great help to reduce the relaxation factors. I had the same problem. The relaxation factors are responsible for the search of a the start value in the calculation of the parameters depending of the result in the last time step. If you reduce them, the start value is not so far away to the last result but you´ll need more time for the calculation.

Try it. You find the relaxationFactors in fvSolution and then reduce all factors 0.1 down. If it doesn´t help, try 0.2.

It could be the geometry, too. But try it.
Hope, it helps.
greetz phil
kroetenechse 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
Explanation bounding k and bounding epsilon idrama OpenFOAM 42 July 13, 2017 04:05
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 16 March 4, 2017 09:30
Bounding epsilon or bounding omega Stylianos OpenFOAM 7 February 26, 2017 00:17
MRFSimpleFOAM goes divergenced! renyun0511 OpenFOAM Running, Solving & CFD 0 November 19, 2009 03:11
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 21:21

Sponsored Links

All times are GMT -4. The time now is 19:31.