
[Sponsors] 
August 12, 2015, 10:19 
SIMPLE vs PIMPLE

#1 
Member
Nicole Andrew
Join Date: Sep 2014
Location: Pretoria, South Africa
Posts: 58
Rep Power: 8 
Hello Foamers!
I am trying to solve for the steady state flow of slag in a furnace. I was having some issues so I returned to a simple case of a rectangular box, with one wall at 1700 degC and the opposite wall at 1500 degC, with a starting temperature of 1600 degC throughout. I used buoyantBoussinesqSimpleFoam but it seems to generate very high velocities (1e5+) before starting to converge to a reasonable solution and even after 1000s of iterations it is still not converging well. On the other hand, I tried buoyantBoussinesqPimpleFoam and it seems to be converging to a reasonable solution much much faster. Has anyone had any similar experiences with the SIMPLE algorithm seeming to converge better/faster than PIMPLE? I have not had a chance to study Pimple much yet. I can post an example case, but thought I would just check if anyone has some quick ideas for me, or some literature I could look at. Thanks! 

August 12, 2015, 13:13 

#2 
Senior Member
anonymous
Join Date: Aug 2014
Posts: 205
Rep Power: 9 
SIMPLE solvers and steady solvers, and thus the time you are seeing in the output does not have a physical meaning
PIMPLE solvers are unsteady solvers. Your case has a big difference in temperature and it might be easier to obtain a "steady" solution using PIMPLE, because it helps to create a diagonal dominant matrix and because your case might be highly unsteady There are also other factors that might influence the case such as mesh, schemes and solution configuration etc 

August 13, 2015, 04:12 

#3 
Member
Nicole Andrew
Join Date: Sep 2014
Location: Pretoria, South Africa
Posts: 58
Rep Power: 8 
Hi,
Thanks so much. The 'time' I was referring to is the time it took my computer to solve. I think you are probably right about my case being highly unsteady, I am glad to hear that this is something that other people may have experienced before. That's interesting about PIMPLE generating more diagonally dominant matrices, I will look into that a bit more. I still need to investigate various schemes etc. But thanks so much for the advice 

August 13, 2015, 04:30 

#4 
Senior Member
anonymous
Join Date: Aug 2014
Posts: 205
Rep Power: 9 
Dear Nicole,
PIMPLE helps to generate a diagonally dominant matrix because of the implicit term that the time derivate operator adds to the diagonal of the matrix of the equation discretized with the finite volume method. I'm glad that I could help you 

August 17, 2015, 12:06 

#5 
New Member
Harshad Joshi
Join Date: Nov 2010
Location: India
Posts: 9
Rep Power: 12 
Hi Nicole,
In addition to above, I believe you are trying the wrong solver. buoyantBoussinesqSimpleFoam is the incompressible solver with effects of buoyancy added in. Rule of thumb for this solver recommended by OpenFOAM trainers is a deltaT of 40C. In the description, the deltaT is clearly 200C. Why not try switching to the buoyantSimpleFoam solver? It is the compressible solver with buoyancy. H. 

August 17, 2015, 18:29 

#6 
Senior Member
anonymous
Join Date: Aug 2014
Posts: 205
Rep Power: 9 
Well I wouldn't really say the same that the OpenFOAM, trainers you are referring to.
I will be quick on it. The Boussinesq approximation can be found in the rhok term which reads rhok=1beta*(TTref) The second term of rhok should be <<1 so that the boussinesq approximation maintains true. The order of magnitude of the second term is beta*deltaT and thus it depends on the properties of the flow. Moreover sometimes even though beta*deltaT is "big" one could stay using the boussinesq approximation to obtain some first results. Although it is true that if the second term becomes big the convergence of the solver becomes difficult 

August 20, 2015, 10:55 

#7 
Member
Nicole Andrew
Join Date: Sep 2014
Location: Pretoria, South Africa
Posts: 58
Rep Power: 8 
Hi Joshi,
Thanks for your inputs. ssss has pretty much answered how I would have. Except to add that the temperature differences in the actual fluid are much smaller than the boundary temperature differences due to the good mixing and heat transfer. So I think my simulation falls in the 40 degree deltaT range in any case. I am not using the full buoyant solver because of issues that I presented here: http://www.cfdonline.com/Forums/ope...edvolume.html 

Tags 
boussinesq, buoyantboussinesqpimple, buoyantboussinesqsimple, simple algorithm 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Continuity divergence in simple gas flow  Okka  FLUENT  7  January 11, 2019 02:01 
What is PIMPLE? Can it be used for supersonic case?  shenzhou1987  OpenFOAM Running, Solving & CFD  3  January 26, 2015 02:48 
questions about pimple !!  1988  OpenFOAM PreProcessing  4  July 1, 2014 03:57 
Pimple?  n.makhtoomi  OpenFOAM Programming & Development  5  June 5, 2014 11:58 
SIMPLE algorithm in 3D cylindrical coordinates  zouchu  Main CFD Forum  1  January 20, 2014 18:02 