|
[Sponsors] |
pimpleFoam vs simpleFoam vs pisoFoam vs icoFoam? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 4, 2009, 17:40 |
pimpleFoam vs simpleFoam vs pisoFoam vs icoFoam?
|
#1 |
Senior Member
Pei-Ying Hsieh
Join Date: Mar 2009
Posts: 334
Rep Power: 18 |
Hi,
Can someone explain the difference among pimpleFoam, simpleFoam, pisoFoam, and icoFoam? When to select which solver? Thanks! Pei |
|
September 5, 2009, 05:01 |
|
#2 |
Senior Member
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 323
Rep Power: 18 |
Hi
icoFoam is transient solver for incompressible, laminar flow of Newtonian fluids. pimpleFoam is large time-step transient solver for incompressible, flow using the PIMPLE(merged PISO-SIMPLE) algorithm. pisoFoam is transient solver for incompressible flow. Turbulence modelling is generic, i.e. laminar, RAS or LES may be selected. simpleFoam is steady-state solver for compressible, turbulent flow Regards Ata Last edited by ata; September 18, 2009 at 05:35. |
|
September 21, 2009, 04:53 |
|
#3 |
Member
David GISEN
Join Date: Jul 2009
Location: Germany
Posts: 70
Rep Power: 17 |
simpleFoam is for incompressible, turbulent flow.
|
|
September 21, 2009, 05:51 |
pimpleFoam vs simpleFoam vs pisoFoam vs icoFoam? Reply to Thread
|
#4 |
Senior Member
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 323
Rep Power: 18 |
Hi David
Thanks it's true. Regards Ata |
|
February 18, 2010, 11:02 |
Using LES with simplefoam
|
#5 |
Senior Member
Join Date: Mar 2009
Location: Norway
Posts: 138
Rep Power: 17 |
May LES be used with simpleFoam?
|
|
February 18, 2010, 11:20 |
|
#6 |
Senior Member
Florian Krause
Join Date: Mar 2009
Location: Munich
Posts: 103
Rep Power: 17 |
Hello,
simpleFoam is a steady-state solver, as ata already mentioned..... Thus, No you cannot use simpleFoam for LES. Use pisoFoam instead! Best, Florian |
|
September 10, 2010, 06:36 |
|
#7 |
Senior Member
Nilesh Rane
Join Date: Apr 2010
Posts: 122
Rep Power: 16 |
guys,
I am getting different results for same case using icoFoam and simpleFoam. icoFoam seems to give better results. What could be the reason?? Isn't the two suppose to give similar results??? BTW i am using the default solvers coming with OF170. Only change is i changed fvscheme for U to upwind (which is the scheme used in simpleFoam anyway).
__________________
Imagination is more important than knowledge..
|
|
September 5, 2012, 12:14 |
|
#8 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 14 |
Hi!
I have the same kind of question: what is the interest of pisoFoam compared to pimpleFoam -and vice-versa. I am trying to simulate the interaction between wind and buildings: my simulation has to be unsteady, incompresible, with large to very large Reynolds number and complex geometry. Should I use pisoFoam or pimpleFoam? And why? In your case, Nilesh, icoFoam is unsteady laminar, whereas simpleFoam is stead turbulent. Equations are not similar! I know people who have choosen to prefer pimpleFoam over pisoFoam (http://www.symscape.com/node/948). But they are running RANS, I would like to perform LES. Any advice or explaination ? |
|
September 5, 2012, 17:41 |
|
#9 |
Senior Member
|
Using LES or RANS doesn't matter. People who normally prefer pimpleFOAM say that pimple Algorithm is more robust and efficient. They also say that it can be used for larger (compared to piso, where courant number of > 1 could result in diverged simulation) time steps. If larger time step is of interest (as in some sudo transient cases ) then surely pimple is an attractive algorithm.
You said that you want to do LES. In LES you might also need to think about time scales and you must keep your time step smaller than the time scales of the large eddies which you are actually simulating. What I have observed is that people normally keep the courant number less that 0.5-0.4 . Pimple offers you more control by providing nOuterIterations parameter. Where as if nOuterIterations is = 1, you pimple is just simply a piso Algorithm. So definitely some addition cost on computation with more outer iterations could improve your results, but it is sometimes an important decision whether you want to improve the results further or not. regards, Awais |
|
September 6, 2012, 05:21 |
|
#10 | |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 14 |
Hi Awais. Thanks a lot for your advices. But I did not catch everything
"as in some sudo transient cases" . What means "sudo" ? My goal is to model very large Reynolds problem: about 1e5 or 1e6. Thus, keeping a little Co is out of my CPU capabilities. Today I am running with Co close to 1. TJunction tutorial works with a Co of 5, and only 1 nOuterIterations . So this tuto is actually a piso algorithm? I think about using a Co of 5 (same as the tuto), and 1 nOuterIterations (just to be not only piso but real pimple). Do you think it'd be OK? Another aspect: I am not familiar with the concept of "final" solvers. For the momment I used the same as "normal" ones: Quote:
Last edited by Djub; September 6, 2012 at 05:58. Reason: adding a question |
||
September 6, 2012, 10:12 |
|
#11 |
Senior Member
|
Hi Julien,
By Sudo Transient I mean a simulation which is physically not transient, but you use a transient flow solver like pimpleFoam. To see this in action you can also have a look at the angleDuct tutorial in rhoPimpleFoam/ras. There the courant number is 10. But if you are using LES, you would need time Averaging of the solution. Plus time scale of the structures in flow is also important. If you are doing RANS simulation then may be the setup you have just said (i.e. Co 5 and 1 nouterIteration) might work. But again in transient simulation you will be interested in averaged values of the flow, so just think whether you will be able to capture all the flow fluctuations with that time step or not. My understanding is that you should keep Co < 1 (specially) when your flow is highly transient like flow over a cylinder produces vortex street and if the time step is high we might miss many modes of the fluctuations and this could be reflected in our averaged values (or may be it could also take more time to get a time averaged solution). Regarding pFinal and UFinal I think these are somehow the stringent criteria for last Iteration on the equation in one pimple Loop so as to get the best possible solution in the numerical sense. Regards, Awais |
|
September 10, 2012, 04:03 |
|
#12 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 |
Small side remark - he probably meant to write 'pseudo'
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. |
|
September 10, 2012, 11:25 |
|
#13 |
Senior Member
|
Thanks for pointing that out Anton.
Regards, Awais |
|
September 12, 2012, 04:59 |
|
#14 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 14 |
Hi everybody,
I ran some cases with Pimple. It rocks! From a working PISO case, the only change is about nOuterIterations (in fvsolution ). In my case, I use 2 nOuterItreations. You have also to choose for final solvers, I suppose for the last iteration. But PIMPLE permits the use of the entry adjustTimeStep yes in the controlDict. In this case, you control directly the CFL (entry maxCo). I am now working with a maxCo of 10 ! And it is still stable and realist . And of course much faster ! In this case, your time step is not constant so you may prefer to use writeControl adjustableRunTime (instead of timestep). In this case, writeInterval is in seconds. Djub Last edited by Djub; September 12, 2012 at 11:35. |
|
March 20, 2013, 13:39 |
Courant number
|
#15 |
Member
Suranga Dharmarathne
Join Date: Jan 2011
Location: TX, USA
Posts: 39
Rep Power: 15 |
Isn't it a problem if you increased Courant number>1
Br, Suranga. |
|
March 21, 2013, 05:53 |
|
#16 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 14 |
Pimple being an Implicit numerical method (cf Wikipedia), you can use a large Co, greater than 1. Co needs to be smaller than 1 for explicit method (see here).
By the way, as owayz said, I took a great care to define correctly the time scales of the large eddies I am simulating. In my case, there was no problem at Co=10. Nevertheless, Co=100 crashed, as well as Co=50. At Co=30, the results were not satisfactory. So the limit (in my case, which was a simple rectangular rod in a laminar flow, for vortex shedding estimation), the limit is between Co=10 (OK) and Co=30 (KO) |
|
March 27, 2013, 14:09 |
physical example
|
#17 |
Senior Member
Join Date: Nov 2012
Location: Bavaria
Posts: 145
Rep Power: 14 |
Hi,
most stupid question comes probably at last. I've read the discussion more then once but I am still not clear about the usage of these different solvers. Could anybody try to explain the advantages/disadvantages on the basis of a practical situation. Supposed there is a duct flow (L > 5000mm, W =250mm, H = 20mm, fluid = water, point of origin is located in the center of the duct) with... a) small Reynoldsnumber (100 - 200) at inlet, no obstacle (e.g. rod or cube) in the middle of the duct, b) small Reynoldsnumber (100 -200) at inlet, obstacle (cube with side length a=20mm) positioned in the center of the duct, c) high Reynoldsnumber (> 3000) at inlet, no obstacle located in the duct, d) high Reynoldsnumber (> 3000) at inlet, obstacle (cube a=20mm) positioned in the center of the duct. For what case would you use which solver and why? You could also change the physical setting of one of the cases if that helps to show the advantage of a specific solver. Thanks a lot in advance! Aylalisa |
|
March 28, 2013, 03:59 |
|
#18 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 |
If there is no motion of the obstacle, then it doesn't influence your choice of solver. I believe you can use pimpleFoam for all cases. You might be able to use even larger timesteps (thus saving some computational time) by using simpleFoam for the high Reynolds cases.
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. |
|
April 3, 2013, 12:03 |
modify the situation
|
#19 |
Senior Member
Join Date: Nov 2012
Location: Bavaria
Posts: 145
Rep Power: 14 |
What kind of flow phenomenon makes you use pisoFoam and icoFoam respectively?
|
|
April 3, 2013, 12:26 |
|
#20 |
Senior Member
Julien
Join Date: Jun 2012
Location: France
Posts: 152
Rep Power: 14 |
icoFoam is laminar, whereas pisoFoam can deal with turbulence.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Laminar simpleFoam and inviscid simpleFoam | herenger | OpenFOAM Running, Solving & CFD | 7 | July 11, 2013 07:27 |
SimpleFoam as Newtonian laminar flow solver | titio | OpenFOAM Running, Solving & CFD | 2 | March 8, 2013 05:44 |
Density in icoFoam Densidad en icoFoam | manuel | OpenFOAM Running, Solving & CFD | 8 | September 22, 2010 05:10 |
Error running simpleFoam in parallel | skabilan | OpenFOAM Running, Solving & CFD | 2 | August 29, 2008 10:42 |
Kubuntu uses dash breaks All scripts in tutorials | platopus | OpenFOAM Bugs | 8 | April 15, 2008 08:52 |