|
[Sponsors] |
turbFoam floating point exception and k-epsilon |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 28, 2009, 03:43 |
turbFoam floating point exception and k-epsilon
|
#1 |
New Member
Join Date: Mar 2009
Posts: 28
Rep Power: 17 |
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. |
|
April 28, 2009, 04:10 |
|
#2 |
Senior Member
Rishi .
Join Date: Mar 2009
Posts: 149
Rep Power: 17 |
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 |
|
April 28, 2009, 04:20 |
|
#3 |
New Member
Join Date: Mar 2009
Posts: 28
Rep Power: 17 |
Thanks for the really quick answer.
I will decrease my timesteps...and look up Courant number |
|
April 28, 2009, 05:10 |
|
#4 |
New Member
Join Date: Mar 2009
Posts: 28
Rep Power: 17 |
Still the same problem.
I am trying the solve it with adjustable deltaT as described by niklas in this forum |
|
April 28, 2009, 07:10 |
|
#5 |
Senior Member
Gavin Tabor
Join Date: Mar 2009
Posts: 181
Rep Power: 17 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
ERROR: Floating point exception: Overflow | Elisabetta | CFX | 23 | July 14, 2018 06:57 |
Error : floating point exception; overflow. | maryam | CFX | 3 | May 18, 2010 07:44 |
Floating point exception: Invalid number | Jinfeng | CFX | 4 | January 11, 2009 16:44 |
floating point exception | mojtaba | CFX | 2 | November 15, 2008 06:58 |
floating point exception | Mojtaba | CFX | 0 | November 1, 2008 01:54 |