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

bounding epsilon blow up

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 30, 2010, 11:42
Default bounding epsilon blow up
  #1
New Member
 
Jie Zhang
Join Date: Sep 2010
Location: FL, U.S.
Posts: 20
Rep Power: 6
jiez is on a distinguished road
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:

epsilo:

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

inlet
{
type fixedValue;
value uniform 0.0012;
}

outlet
{
type zeroGradient;
}
jiez is offline   Reply With Quote

Old   September 30, 2010, 11:46
Default
  #2
New Member
 
Jie Zhang
Join Date: Sep 2010
Location: FL, U.S.
Posts: 20
Rep Power: 6
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
Default
  #3
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 15
santiagomarquezd will become famous soon enough
Maybe a draft of the geometry, details of problem and

system/controlDict
system/fvSchemes
system/dvSolution
0/*

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

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   October 8, 2010, 17:11
Default
  #4
Senior Member
 
Travis Carrigan
Join Date: Jul 2010
Location: Arlington, TX
Posts: 128
Rep Power: 7
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 online now   Reply With Quote

Old   January 16, 2011, 18:14
Default
  #5
New Member
 
Philipp Bachmann
Join Date: May 2010
Location: Esslingen, Germany
Posts: 7
Rep Power: 7
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

Reply

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 37 February 5, 2014 13:01
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 13 November 4, 2013 15:13
Bounding epsilon or bounding omega Stylianos OpenFOAM 5 June 9, 2011 00:40
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


All times are GMT -4. The time now is 16:30.