CFD Online Discussion Forums

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)
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!

pcaste12 April 4, 2012 02:35

Any help?
 
any help on this? :(


All times are GMT -4. The time now is 06:37.