CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   rhoPimpleFoam quasi square block floating point error (http://www.cfd-online.com/Forums/openfoam/99200-rhopimplefoam-quasi-square-block-floating-point-error.html)

 pcaste12 March 28, 2012 09:19

rhoPimpleFoam quasi square block floating point error

4 Attachment(s)
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!

 pcaste12 April 4, 2012 02:35

Any help?

any help on this? :(

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