CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

rhoPimpleFoam quasi square block floating point error

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 28, 2012, 10:19
Default rhoPimpleFoam quasi square block floating point error
  #1
Member
 
Pierre Castellani
Join Date: Apr 2011
Location: Paris
Posts: 38
Rep Power: 15
pcaste12 is on a distinguished road
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 ) and even if I made a lots of home work, it is still quite dark inside.

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!
Attached Images
File Type: png box.png (6.3 KB, 14 views)
Attached Files
File Type: txt chM.txt (2.8 KB, 7 views)
File Type: zip fv.zip (1.4 KB, 9 views)
File Type: zip 0.zip (3.7 KB, 7 views)
pcaste12 is offline   Reply With Quote

Old   April 4, 2012, 03:35
Default Any help?
  #2
Member
 
Pierre Castellani
Join Date: Apr 2011
Location: Paris
Posts: 38
Rep Power: 15
pcaste12 is on a distinguished road
any help on this?
pcaste12 is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with UDF compiling for kTkLW model Wantami FLUENT 0 July 18, 2011 06:11
How to install CGNS under windows xp? lzgwhy Main CFD Forum 1 January 11, 2011 19:44
[Gmsh] Gmsh and samplesurface touf OpenFOAM Meshing & Mesh Conversion 2 December 10, 2007 03:27
user defined function cfduser CFX 0 April 29, 2006 11:58
error while compiling the USER Sub routine CFD user CFX 3 November 25, 2002 16:16


All times are GMT -4. The time now is 00:13.