CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Floating point exception (

robingilbert April 16, 2010 15:34

Floating point exception
Hello guys,

I am trying to run a simulation of a heatsource placed inside an airconditioned room. i used the rhoPimpleFoam solver. my mesh quality is good ( hexahedral meshes) and checkMesh showed no problems at all. but on running the simulation i am getting the following error:


rho max/min : 3 0.3
bounding p, min: -1e+06 max: 3e+06 average: 1e+05
DILUPBiCG:  Solving for epsilon, Initial residual = 0.9, Final residual = 9e-06, No Iterations 645
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 1e-06, No Iterations 6
DILUPBiCG:  Solving for Ux, Initial residual = 0.9, Final residual = 0.04, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.02, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.9, Final residual = 0.02, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.1, Final residual = 0.003, No Iterations 2
#0  Foam::error::printStack(Foam::Ostream&) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/"
#1  Foam::sigFpe::sigFpeHandler(int) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/"
#2  ?? in "/lib/"
#3  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/"
#4  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/"
#5  main in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/applications/bin/linux64GccDPOpt/rhoPimpleFoam"
#6  __libc_start_main in "/lib/"
#7  _start at /build/buildd/eglibc-2.10.1/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception

please help me solve this problem

matejfor April 19, 2010 02:56


Its pretty easy. The SutherlandTransport routine in function calculate and correct is dividing by zero, which is not allowed in this solar system. To avoid legal consequences, the program crashed. I guess you are running with pressure or temperature to nonphysical values either caused by wrong BCs or initials. I would set all the thermo properties constant to see what is going on and turn the sutherland back later.

more info about the crash info is here.

good luck

robingilbert April 21, 2010 20:55

thank you,

I figured out the problem.... i had given a non-physical value which caused the problem... thanks for helping me out.

trex930 June 24, 2010 03:52

Hey robingilbert,

How did you manage to fix it? I have been playing around with my boundary conditions, and this thing keeeeeps on popping up. Which condition did you change that made it work? was it in the Pressure file?


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