|
[Sponsors] |
Floating point exception error with small DeltaT |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 1, 2009, 21:37 |
Floating point exception error with small DeltaT
|
#1 |
Member
Join Date: Mar 2009
Location: Sydney, New South Wales, Australia
Posts: 42
Rep Power: 17 |
Hi Foamers,
I have been going really well with foam over the past couple omf months until now. I have recently imported a 3D model from a gambit .MSH file. I scaled it successfully using the -[scale] command, and initially it runs very well, until it reaches a point where it bombs out with a Floating point exception error. ----- #0 Foam::error: printStack(Foam::Ostream&) in "/home/swervo/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/swervo/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so" #2 Uninterpreted: [0xffffe400] #3 Foam::Time::adjustDeltaT() in "/home/swervo/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so" #4 Foam::Time::setDeltaT(double) in "/home/swervo/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so" #5 main in "/home/swervo/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/interFoam" #6 __libc_start_main in "/lib/libc.so.6" #7 Foam::regIOobject::readIfModified() in "/home/swervo/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/interFoam" Floating point exception ----- Is there anybody that can shed any light on this error? The model is in the micro domain, solving using interFoam. I have had to use a very small timestep (1E-12) in previous models to reach a convergent solution, but this in itself hasn't thrown up this error in previous models. I notice that in lines #3 and #4 of the error message it mentions DeltaT, but doesn't give any further clues. Does anybody have any ideas? Cheers, R |
|
September 1, 2009, 22:58 |
|
#2 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
hi,
i have a worse result than yours ,the steps i used are as follows: 1.fluentMeshToFoam voim.msh -scale 1000 2.gedit constant/polyMesh/boundary & 3.gedit 0/* & 4.gedit constant/MRFZones & 5.gedit system/controlDict & 6.MRFSimpleFoam but,it displayed : cannot open file file: /home/ry/OpenFOAM/ry-1.5/tutorials/MRFSimpleFoam/voim/system/controlDict at line 0. From function regIOobject::readStream(const word&) in file db/regIOobject/regIOobjectRead.C at line 66. i don't know why? can you help me? regards |
|
September 1, 2009, 22:59 |
|
#3 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
before MRFSimpleFoam,i put $wmake in the console
|
|
September 1, 2009, 23:32 |
|
#4 | |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
Quote:
Best regards, Chiven |
||
September 1, 2009, 23:39 |
|
#5 | |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
Quote:
The timestep is not the key point, please check the initialized fields, the discretization and linear solver settings. Best regards, Chiven |
||
September 2, 2009, 04:38 |
|
#6 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
Dear chiven,
the controlDict file: FoamFile { version 2.0; format ascii; class dictionary; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application simpleFoam; startFrom startTime; startTime 0; stopAt endTime; endTime 500; deltaT 1; writeControl timeStep; writeInterval 50; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression uncompressed; timeFormat general; timePrecision 6; runTimeModifiable yes; // ************************************************** *********************** // I don't know where is the problem.how to modify to it? would you like to help me? regards, jennifer |
|
September 2, 2009, 04:49 |
|
#7 |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
||
September 2, 2009, 04:57 |
|
#8 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
Dear chiven,
I try to modify it ,but the result is all the same :cannot open file file: /home/ry/OpenFOAM/ry-1.5/tutorials/MRFSimpleFoam/voim/system/controlDict at line 0. From function regIOobject::readStream(const word&) in file db/regIOobject/regIOobjectRead.C at line 66. |
|
September 2, 2009, 05:05 |
|
#9 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
it puzzled me several days!
|
|
September 2, 2009, 06:19 |
|
#10 |
Senior Member
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 23 |
If the solver can't find the file, then you're not running it from the correct directory. I see the case path is:
/home/ry/OpenFOAM/ry-1.5/tutorials/MRFSimpleFoam/voim/ Are you sure you're running MRFSimpleFoam from that directory?
__________________
Laurence R. McGlashan :: Website |
|
September 2, 2009, 06:42 |
|
#11 |
New Member
Patrick Mueller
Join Date: Jun 2009
Location: Kaiserslautern
Posts: 8
Rep Power: 17 |
I have a similar problem to yours, therefore I'll just post it here instead of starting a whole new thread.
I'm currently trying to get the MRFSimpleFoam example from the FoamWiki running. Since my iterations fail after ~7 equations I tried decreasing the timestep to get a Courant Number <1. Unfortunately, for that I have to choose a Step less than 1.8e-4. When I decrease the timestep to this number I get the following error: #0 Foam::error: printStack(Foam::Ostream&) in "/home/mueller/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/mueller/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #2 ?? in "/lib/libc.so.6" #3 void Foam::multiply<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<Foam::Tensor<d ouble>, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/home/mueller/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so" #4 Foam::tmp<Foam::GeometricField<Foam::Tensor<double >, Foam::fvPatchField, Foam::volMesh> > Foam: operator*<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<Foam::Tensor<double >, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/mueller/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so" #5 Foam::incompressible::RASModels::kEpsilon::divDevR eff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/home/mueller/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so" #6 main in "/home/mueller/OpenFOAM/mueller-1.6/applications/bin/linux64GccDPOpt/MRFSimpleFoam" #7 __libc_start_main in "/lib/libc.so.6" #8 _start at /build/buildd/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116 Floating point exception I get the same error when I reduce the speed of the rotor to 10 ups and increase the timestep to 2.5e-3. Has anybody an idea why I can't decrease the timestep to an usable amount without getting a floating point exception? Regards |
|
September 2, 2009, 07:22 |
|
#12 | |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
Quote:
If you revise the file in windows platform, maybe it can not be opened. Delete the OLD controlDict file and copy another from the example case of MRFSimple case. Good luck. Chiven |
||
September 2, 2009, 07:50 |
|
#13 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
yes, I put the command on the console:gedit system/controlDict &
then,it displayed the information above, I changed the "application simpleFoam"to"application MRFSimpleFoam",saved it ,at last MRFSimpleFoam,but the result also is:cannot open file. |
|
September 2, 2009, 07:52 |
|
#14 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
copy another from the example case of MRFSimple case?where is the MRFSimple example case?
|
|
September 2, 2009, 07:53 |
|
#15 | |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
Quote:
it is interesting. But i have no idea. best regards, Chiven |
||
September 2, 2009, 07:57 |
|
#16 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
thank all the same!
regards, jennifer |
|
September 2, 2009, 08:01 |
|
#17 |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
||
September 2, 2009, 08:15 |
|
#18 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
how to do the calculation of myself?
|
|
September 2, 2009, 08:23 |
|
#19 |
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 17 |
||
September 2, 2009, 08:39 |
|
#20 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
thank you, chiven!
I sent a email to you just now,it detailed all the steps i adopt,hope for your guidance! |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Floating point exception? | Ricky Wong | FLUENT | 19 | March 28, 2020 16:32 |
turbFoam floating point exception and k-epsilon | Hectux | OpenFOAM Running, Solving & CFD | 4 | April 28, 2009 08:10 |
Floating point exception: Invalid number | Jinfeng | CFX | 4 | January 11, 2009 17:44 |
Floating Point Exception in STAR CCM+ | Sheila | Siemens | 0 | March 31, 2008 07:48 |
[Gmsh] Gmsh and samplesurface | touf | OpenFOAM Meshing & Mesh Conversion | 2 | December 10, 2007 03:27 |