[robingilbert]

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/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  ?? in "/lib/libc.so.6"
#3  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#4  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#5  main in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/applications/bin/linux64GccDPOpt/rhoPimpleFoam"
#6  __libc_start_main in "/lib/libc.so.6"
#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]

Hi,

It§s 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
matej

[robingilbert]

thank you,

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

[trex930]

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?

T