|March 28, 2012, 09:19||
rhoPimpleFoam quasi square block floating point error
Join Date: Apr 2011
Posts: 38Rep Power: 7
I will try to follow the recommendation made in the thread
I will not create a dropbox account but add the image direclty as far as my case is really simple
I am working with rhoPimple even if for the time bieng it is not transcient for the purpose of test
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.
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)
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)
Finnaly the error:
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
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!
|Thread||Thread Starter||Forum||Replies||Last Post|
|Problem with UDF compiling for kTkLW model||Wantami||FLUENT||0||July 18, 2011 05:11|
|How to install CGNS under windows xp?||lzgwhy||Main CFD Forum||1||January 11, 2011 19:44|
|Gmsh and samplesurface||touf||Open Source Meshers: Gmsh, Netgen, CGNS, ...||2||December 10, 2007 03:27|
|user defined function||cfduser||CFX||0||April 29, 2006 10:58|
|error while compiling the USER Sub routine||CFD user||CFX||3||November 25, 2002 16:16|