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

turbFoam floating point exception and k-epsilon

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 28, 2009, 03:43
Default turbFoam floating point exception and k-epsilon
  #1
New Member
 
Join Date: Mar 2009
Posts: 28
Rep Power: 8
Hectux is on a distinguished road
Hi,

I got two questions:
Does anyone have an idea what the reason for this floating point exception could be:

Courant Number mean: 5.18024e+43 max: 4.56829e+47
DILUPBiCG: Solving for Ux, Initial residual = 0.877633, Final residual = 9.78282e-09, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.938388, Final residual = 9.62189e-09, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.952206, Final residual = 2.34812e-08, No Iterations 2
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00298175, No Iterations 1001
time step continuity errors : sum local = 1.19866e+81, global = -3.05383e+62, cumulative = -3.05383e+62
DICPCG: Solving for p, Initial residual = 9.50744e-07, Final residual = 9.50744e-07, No Iterations 0
time step continuity errors : sum local = 2.17993e+170, global = -1.77391e+136, cumulative = -1.77391e+136
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam:ILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam:ILUPreconditioner:ILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::Istream&) in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::lduMatrix:reconditioner::addasymMatrixCons tructorToTable<Foam:ILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::Istream&) in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::Istream&) in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#8 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libfiniteVolume.so"
#9 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressibleTurbulenceModels.so"
#10 Foam::turbulenceModels::kEpsilon::correct() in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libincompressibleTurbulenceModels.so"
#11 main in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/turbFoam"
#12 __libc_start_main in "/lib/i686/libc.so.6"
#13 Foam::regIOobject::readIfModified() in "/home/adam/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/turbFoam"
Floating point exception

And my second question:

What is influenced by the choice of InternalField in the file for k und epsilon?
At the moment I had small numbers like 1e-11 for both.
Does someone have literature for this case?


Thanks and best regards.
Hectux is offline   Reply With Quote

Old   April 28, 2009, 04:10
Default
  #2
Senior Member
 
Rishi .
Join Date: Mar 2009
Posts: 138
Rep Power: 8
hellorishi is on a distinguished road
Hi,

My guess is that your deltaT is very large !

Courant Number mean: 5.18024e+43 max: 4.56829e+47
and
time step continuity errors : sum local = 2.17993e+170, global = -1.77391e+136, cumulative = -1.77391e+136

I would suggest to check your log file and check the Courant number from first time step. If its increasing exponentially( my guess), then use 1/10th of current delta t till you get Courant number near 1 or less. Courant number should be below 1.

Hope this helps
Rishi
hellorishi is offline   Reply With Quote

Old   April 28, 2009, 04:20
Default
  #3
New Member
 
Join Date: Mar 2009
Posts: 28
Rep Power: 8
Hectux is on a distinguished road
Thanks for the really quick answer.

I will decrease my timesteps...and look up Courant number
Hectux is offline   Reply With Quote

Old   April 28, 2009, 05:10
Default
  #4
New Member
 
Join Date: Mar 2009
Posts: 28
Rep Power: 8
Hectux is on a distinguished road
Still the same problem.

I am trying the solve it with adjustable deltaT as described by niklas in this forum
Hectux is offline   Reply With Quote

Old   April 28, 2009, 07:10
Default
  #5
Senior Member
 
Gavin Tabor
Join Date: Mar 2009
Posts: 181
Rep Power: 8
grtabor is on a distinguished road
Looks to me like you have a divergence problem with your solution; you are unlikely to solve this by adjusting the timestep (particularly by the time you get to a courant number e+47! Go back through the log file and see if you can identify which equation is having the problem at the start (you might notice that one of the equations is not converging, or stops converging, or suddenly requires 500 iterations to converge where it used to need 50). Is the mesh strongly non-orthogonal? If so, generate a better one or increase the non-orthogonality corrections. Try running with upwind differencing on the div term for a bit, or use simpleFoam to start the solution off.

Gavin
grtabor 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
ERROR: Floating point exception: Overflow Elisabetta CFX 19 March 30, 2015 06:34
Error : floating point exception; overflow. maryam CFX 3 May 18, 2010 07:44
Floating point exception: Invalid number Jinfeng CFX 4 January 11, 2009 17:44
floating point exception mojtaba CFX 2 November 15, 2008 07:58
floating point exception Mojtaba CFX 0 November 1, 2008 02:54


All times are GMT -4. The time now is 14:43.