CFD Online Discussion Forums - rhoPimpleFoam quasi square block floating point error

I will try to follow the recommendation made in the thread
http://http://www.cfd-online.com/Forums/openfoam/98988-how-give-enough-info-get-help.html

So
===0===
I will not create a dropbox account but add the image direclty as far as my case is really simple
===1===
I am working with rhoPimple even if for the time bieng it is not transcient for the purpose of test
===2===
The geometry is very simple: two boxes attached with a very small difference in the size in the middle (see the box.png)
Very important if I make a perfectly square box then it works!
Also find the checkMesh log file in attachement (chM.txt).
I am using the version 2.1 of OpenFoam.
===3===
Now the part where I am not totally sure, I've used the example of the Tjunction so you will find the fvSolution and the fvShemes alsmost identical as the Tjunction one (fv.zip)
===4===
I will also add the 0 files, because I am still new (I think it will be very hard to be old for an old man like me) in openFoam so I could have made mistake (0.zip)

===5===
Finnaly the error:

Code:

PIMPLE: iteration 50

DILUPBiCG:  Solving for Ux, Initial residual = 0.0273749, Final residual = 2.68484e-06, No Iterations 1

DILUPBiCG:  Solving for Uy, Initial residual = 0.0250746, Final residual = 3.50406e-06, No Iterations 2

DILUPBiCG:  Solving for Uz, Initial residual = 0.0273749, Final residual = 2.68484e-06, No Iterations 1

DILUPBiCG:  Solving for h, Initial residual = 0.00403133, Final residual = 9.32419e-06, No Iterations 2

DICPCG:  Solving for p, Initial residual = 0.122236, Final residual = 6.46543e-07, No Iterations 93

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0

time step continuity errors : sum local = 2.02885e-07, global = -5.91572e-09, cumulative = 0.0011558

rho max/min : 2 0.5

DILUPBiCG:  Solving for epsilon, Initial residual = 0.541168, Final residual = 5.44544e-06, No Iterations 3

DILUPBiCG:  Solving for k, Initial residual = 0.640926, Final residual = 4.25578e-06, No Iterations 3

PIMPLE: not converged within 50 iterations

ExecutionTime = 22.07 s  ClockTime = 22 s



 Averages of T :  outlet = 298.251  inlet = 298

 Averages of p :  outlet = 101300  inlet = 106128

 Averages of rho :  outlet = 1.18057  inlet = 1.19037

 Averages of U :  outlet = (-4.19546e-05 21.4418 -4.21684e-05)  inlet = (0 27.5132 0)

 MassFlows:   outlet = 0.00102701  inlet = -0.001

Courant Number mean: 2.83429 max: 53.8633

deltaT = 4.00441e-06

Time = 0.00115827



diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0

PIMPLE: iteration 1

DILUPBiCG:  Solving for Ux, Initial residual = 0.888993, Final residual = 0.0180781, No Iterations 1

DILUPBiCG:  Solving for Uy, Initial residual = 0.378975, Final residual = 0.00749124, No Iterations 1

DILUPBiCG:  Solving for Uz, Initial residual = 0.888993, Final residual = 0.0180781, No Iterations 1

DILUPBiCG:  Solving for h, Initial residual = 0.842746, Final residual = 0.0214127, No Iterations 1

#0  Foam::error::printStack(Foam::Ostream&) in "/media/data/software/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#1  Foam::sigFpe::sigHandler(int) in "/media/data/software/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

#2   in "/lib/libc.so.6"

#3  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/media/data/software/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"

#4  Foam::hPsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/media/data/software/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"

#5  

 in "/media/data/software/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"

#6  __libc_start_main in "/lib/libc.so.6"

#7  

 in "/media/data/software/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"

Exception en point flottant

Of course this test is for my quest of solving a bigger problem, I am still on the same problem since 9 month (I should propably change my job :p ) and even if I made a lots of home work, it is still quite dark inside:confused:.

I ve made this test with rhoSimpleFoam and rhoSimplecFoam and it looks good. The purpose of using rhoPimple is because I need a quite high pressure (5bar~500000pa) and i would like to reach this pressure gradualy.

As usual many thank for your help!