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

Floating point exception error with small DeltaT

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 1, 2009, 20:37
Default Floating point exception error with small DeltaT
  #1
Member
 
Join Date: Mar 2009
Location: Sydney, New South Wales, Australia
Posts: 42
Rep Power: 8
rassilon is on a distinguished road
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
rassilon is offline   Reply With Quote

Old   September 1, 2009, 21:58
Default
  #2
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
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
renyun0511 is offline   Reply With Quote

Old   September 1, 2009, 21:59
Default
  #3
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
before MRFSimpleFoam,i put $wmake in the console
renyun0511 is offline   Reply With Quote

Old   September 1, 2009, 22:32
Default
  #4
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by renyun0511 View Post
hi,
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
Please check the controlDict file, if there is no problem, delete it, then copy another from the corresponding tutorial cases and change it for your case.

Best regards,
Chiven
chiven is offline   Reply With Quote

Old   September 1, 2009, 22:39
Default
  #5
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by rassilon View Post
Hi Foamers,


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.


R

The timestep is not the key point, please check the initialized fields, the discretization and linear solver settings.

Best regards,
Chiven
chiven is offline   Reply With Quote

Old   September 2, 2009, 03:38
Default
  #6
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
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
renyun0511 is offline   Reply With Quote

Old   September 2, 2009, 03:49
Default
  #7
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by renyun0511 View Post
Dear chiven,
application simpleFoam;

jennifer
try it

application MRFSimpleFoam;

and you had better reference to example case of MRFSimpleFoam.

best regards,
Chiven
chiven is offline   Reply With Quote

Old   September 2, 2009, 03:57
Default
  #8
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
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.
renyun0511 is offline   Reply With Quote

Old   September 2, 2009, 04:05
Default
  #9
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
it puzzled me several days!
renyun0511 is offline   Reply With Quote

Old   September 2, 2009, 05:19
Default
  #10
Senior Member
 
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 14
l_r_mcglashan will become famous soon enough
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
l_r_mcglashan is offline   Reply With Quote

Old   September 2, 2009, 05:42
Default
  #11
New Member
 
Patrick Mueller
Join Date: Jun 2009
Location: Kaiserslautern
Posts: 8
Rep Power: 8
PaMuelle is on a distinguished road
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
PaMuelle is offline   Reply With Quote

Old   September 2, 2009, 06:22
Default
  #12
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by renyun0511 View Post
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.
How do you revise the controlDict? Do you use the editor of linux, like vi?

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
chiven is offline   Reply With Quote

Old   September 2, 2009, 06:50
Default
  #13
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
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.
renyun0511 is offline   Reply With Quote

Old   September 2, 2009, 06:52
Default
  #14
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
copy another from the example case of MRFSimple case?where is the MRFSimple example case?
renyun0511 is offline   Reply With Quote

Old   September 2, 2009, 06:53
Default
  #15
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by renyun0511 View Post
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.

it is interesting. But i have no idea.

best regards,
Chiven
chiven is offline   Reply With Quote

Old   September 2, 2009, 06:57
Default
  #16
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
thank all the same!

regards,
jennifer
renyun0511 is offline   Reply With Quote

Old   September 2, 2009, 07:01
Default
  #17
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by renyun0511 View Post
copy another from the example case of MRFSimple case?where is the MRFSimple example case?

you can try to run the case

OpenFOAM/OpenFOAM-1.5/tutorials/MRFSimpleFoam/mixerVessel2D

and then do the calculation of yourself.

best regards,
Chiven
chiven is offline   Reply With Quote

Old   September 2, 2009, 07:15
Default
  #18
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
how to do the calculation of myself?
renyun0511 is offline   Reply With Quote

Old   September 2, 2009, 07:23
Default
  #19
Senior Member
 
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 8
chiven is on a distinguished road
Quote:
Originally Posted by renyun0511 View Post
how to do the calculation of myself?

reference to the user guide

http://www.opencfd.co.uk/openfoam/doc/user.html

enjoy,
chiven
chiven is offline   Reply With Quote

Old   September 2, 2009, 07:39
Default
  #20
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
thank you, chiven!
I sent a email to you just now,it detailed all the steps i adopt,hope for your guidance!
renyun0511 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
Floating point exception? Ricky Wong FLUENT 9 February 15, 2014 23:49
turbFoam floating point exception and k-epsilon Hectux OpenFOAM Running, Solving & CFD 4 April 28, 2009 07:10
Floating point exception: Invalid number Jinfeng CFX 4 January 11, 2009 17:44
Floating Point Exception in STAR CCM+ Sheila CD-adapco 0 March 31, 2008 06:48
Gmsh and samplesurface touf Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 December 10, 2007 03:27


All times are GMT -4. The time now is 00:04.