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

divergence problem with turbulence

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 6, 2011, 08:15
Default divergence problem with turbulence
  #1
Senior Member
 
Join Date: Mar 2011
Posts: 158
Rep Power: 15
tH3f0rC3 is on a distinguished road
Hi,

when I solve my case with f.rx. simpleFoam with turbulence switched off there is no problem.
But as soon as I turn on turbulence there are divergence problems, propably due to a negative epsilon or k value because of bad itteration.
I also know that this problem can be solved by using methods with higher arrangement (order,... don't know how it is called. ax + c - ax˛ +bx +c) or by limiting the gradient.

How can I limit the gradient? -> In fvSolutions?
How can use methods with higher arrangement?

Does someone know?
Thanks for the answer!

Best Regards,
tH3f0rC3
tH3f0rC3 is offline   Reply With Quote

Old   April 6, 2011, 08:24
Default
  #2
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 23
maddalena will become famous soon enough
Quote:
Originally Posted by tH3f0rC3 View Post
Hi,

when I solve my case with f.rx. simpleFoam with turbulence switched off there is no problem.
But as soon as I turn on turbulence there are divergence problems, propably due to a negative epsilon or k value because of bad itteration.
I also know that this problem can be solved by using methods with higher arrangement (order,... don't know how it is called. ax + c - ax˛ +bx +c) or by limiting the gradient.

How can I limit the gradient? -> In fvSolutions?
How can use methods with higher arrangement?

Does someone know?
Thanks for the answer!

Best Regards,
tH3f0rC3
Hi,
maybe this is a good thread where to start: http://www.cfd-online.com/Forums/ope...ime-steps.html


mad
maddalena is offline   Reply With Quote

Old   April 6, 2011, 09:29
Default
  #3
Senior Member
 
Join Date: Mar 2011
Posts: 158
Rep Power: 15
tH3f0rC3 is on a distinguished road
Hi maddalena,

"try lowering the tolerances of the linear solvers. Something in the order of 1e-12 might be appropriate."
Can you tell me where i can lower the tolerance? Is it in the fvSolution?

When I change in fvSolution the solver crashes (at time step 20) even earlier as before (time step 27).
Here is the error message:
Code:
Time = 18
DILUPBiCG:  Solving for Ux, Initial residual = 0.311333, Final residual = 0.0143 016, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.184059, Final residual = 0.0120 39, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.282154, Final residual = 0.0103 257, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.699041, Final residual = 0.0292449,  No Iterations 10
time step continuity errors : sum local = 6.17123e+12, global = -4.29599e+11, cu mulative = -1.38856e+12
DILUPBiCG:  Solving for epsilon, Initial residual = 0.0355234, Final residual = 0.00232541, No Iterations 1
bounding epsilon, min: -1.20957e+34 max: 2.38968e+38 average: 9.01526e+33
DILUPBiCG:  Solving for k, Initial residual = 0.249567, Final residual = 0.01152 8, No Iterations 1
bounding k, min: -9.14125e+25 max: 6.68445e+29 average: 6.57551e+26
ExecutionTime = 83 s  ClockTime = 84 s
Time = 19
DILUPBiCG:  Solving for Ux, Initial residual = 0.956852, Final residual = 0.0501 376, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.99477, Final residual = 0.01063 42, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.957249, Final residual = 0.0912 169, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.965656, Final residual = 0.0432162,  No Iterations 2
time step continuity errors : sum local = 5.79332e+53, global = 2.72507e+50, cum ulative = 2.72507e+50
DILUPBiCG:  Solving for epsilon, Initial residual = 1, Final residual = 0.005137 37, No Iterations 46
DILUPBiCG:  Solving for k, Initial residual = 2.50403e-13, Final residual = 2.50 403e-13, No Iterations 0
ExecutionTime = 87.22 s  ClockTime = 89 s
Time = 20
DILUPBiCG:  Solving for Ux, Initial residual = 0.605269, Final residual = 0.0264 482, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.592722, Final residual = 0.0560 965, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.598655, Final residual = 0.0273 214, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.136164, Final residual = 0.00529025 , No Iterations 2
time step continuity errors : sum local = 7.53995e+60, global = -4.89601e+50, cu mulative = -2.17093e+50
#0  Foam::error::printStack(Foam::Ostream&) in "/data/caehgb.za/studienarbeit-di rk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOA M.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/data/caehgb.za/studienarbeit-dirk/work /OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  __restore_rt in "/lib64/tls/libc.so.6"
#3  Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigne d char) const in "/data/caehgb.za/studienarbeit-dirk/work/OpenFoam_ParaView/Open FOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/data/caehgb.za/s tudienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccD POpt/libfiniteVolume.so"
#5  Foam::incompressible::RASModels::kEpsilon::correct() in "/data/caehgb.za/stu dienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPO pt/libincompressibleRASModels.so"
#6  main in "/data/caehgb.za/studienarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/ OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/simpleFoam"
#7  __libc_start_main in "/lib64/tls/libc.so.6"
#8  Foam::regIOobject::writeObject(Foam::IOstream::streamFormat, Foam::IOstream: :versionNumber, Foam::IOstream::compressionType) const in "/data/caehgb.za/studi enarbeit-dirk/work/OpenFoam_ParaView/OpenFOAM/OpenFOAM-1.7.1/applications/bin/li nux64GccDPOpt/simpleFoam"
Floating point exception
Best Regards,
tH3f0rC3

Last edited by tH3f0rC3; April 6, 2011 at 09:45.
tH3f0rC3 is offline   Reply With Quote

Old   April 6, 2011, 11:00
Default
  #4
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 23
maddalena will become famous soon enough
Quote:
Originally Posted by tH3f0rC3 View Post
time step continuity errors : sum local = 6.17123e+12, global = -4.29599e+11, cu mulative = -1.38856e+12

If you have such a big continuity error at time step 20, then your problem is on the BC probably...
check them once again.

mad
maddalena is offline   Reply With Quote

Old   April 6, 2011, 11:14
Default
  #5
Senior Member
 
Join Date: Mar 2011
Posts: 158
Rep Power: 15
tH3f0rC3 is on a distinguished road
Ok, I will check this.

Best Regards,
tH3f0rC3
tH3f0rC3 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Nozzle Flow - Divergence Problem Ijaz FLUENT 9 January 11, 2014 05:36
Turbulence model for mixing problem nileshjrane OpenFOAM Running, Solving & CFD 1 September 7, 2010 18:48
High Turbulence Intensity Problem bwg FLUENT 1 January 13, 2010 14:09
divergence problem brusly FLUENT 0 October 10, 2006 01:25
a problem about turbulence modeling Chen Yi Main CFD Forum 8 May 30, 2003 03:11


All times are GMT -4. The time now is 08:23.