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

How to use PIMPLE properly?

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree10Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 15, 2019, 05:34
Default
  #21
Senior Member
 
floquation's Avatar
 
Kevin van As
Join Date: Sep 2014
Location: TU Delft, The Netherlands
Posts: 252
Rep Power: 17
floquation will become famous soon enough
Quote:
Originally Posted by Santiago View Post
You use the SIMPLE loops (or as you call them, outer) in an attempt to stabilize (make the simulation more robust) the run by introducing artificial viscosity via a relaxation flux. Is important to recall that some of the time schemes in OF add some sort of numerical flux in the convective term. Thus, repetitive SIMPLE iterations will add to this spurious flux. Note that This "error" is inversely proportional to delta t. In an orthogonal grid, such flux may become a 'false friend' if your Co becomes too small. So, if your grid is not bad, Why do you need more than one SIMPLE loop?
I can see how such an artificial flux may result in outer/SIMPLE iterations not converging. When the relaxation is turned off for the final timestep (here: iteration 200), the residuals do again converge, because the artificial flux is "turned off"?
What remains is the question "why use more than one iteration": That's to make sure that your system of equations are solved properly. Within one outer/SIMPLE iteration, the alpha/VOF equation would only be solved once using the old pressure and velocity fields. (Evidently, pressure and velocity are already being iterated using the inner/PISO iterations.) To correctly solve for the entire system of equations, you'd thus need to iterate, don't you?
If I combine what you said with what I think to be true, then one should use little outer/SIMPLE iterations together with no relaxation for the best result? (Unless relaxation is needed, because the mesh is poor / system is ill-posed / etc., for which the system needs stabilisation.)
floquation is offline   Reply With Quote

Old   February 27, 2019, 19:03
Default
  #22
New Member
 
Join Date: Aug 2018
Posts: 9
Rep Power: 5
mszeto715 is on a distinguished road
Thank you for the response Joaran. How have people learned to use PIMPLE correctly before Tobias Holzmann's publication? It's surprising to me that the usage explained by Dr. Holzmann is not a template in any of the OpenFOAM tutorials. Perhaps this is an artifact of OpenFOAM being open-source code. I have even found that it is not applicable to some of the PIMPLE-related solvers like heatTransfer/chtMultiRegionFoam (though I guess they are just outdated). I would like to better understand the context behind this correct usage (i.e., using residual controls, setting the number of outer loops to some high value). Also, are there other sources that describe strategies for configuring the settings? Many thanks for your attention. -Mimi
mszeto715 is offline   Reply With Quote

Old   February 28, 2019, 06:44
Default
  #23
New Member
 
Joaquín Aranciaga
Join Date: Oct 2018
Posts: 20
Rep Power: 4
joaran is on a distinguished road
Quote:
Originally Posted by mszeto715 View Post
Thank you for the response Joaran. How have people learned to use PIMPLE correctly before Tobias Holzmann's publication? It's surprising to me that the usage explained by Dr. Holzmann is not a template in any of the OpenFOAM tutorials. Perhaps this is an artifact of OpenFOAM being open-source code. I have even found that it is not applicable to some of the PIMPLE-related solvers like heatTransfer/chtMultiRegionFoam (though I guess they are just outdated). I would like to better understand the context behind this correct usage (i.e., using residual controls, setting the number of outer loops to some high value). Also, are there other sources that describe strategies for configuring the settings? Many thanks for your attention. -Mimi
Hi mszeto715,



I actually haven't found any better information about the configuration of fvSolution than Holzmann's book. I think that the best way to learn to use OpenFOAM is playing with it, making mistakes, and search in this forum and/or googling your errors to find an answer. I suppose though it'd be better to take a course, if you can afford it. And ultimately all the information is contained in the installation files, be it a comment, or a program line, so I'd recommend you try to learn how to read the files and make an effort to understand what the main lines do. For that, it's a really good idea to spend some time taking a tutorial on C++.
Sorry for not being such a help, this was just my (tiny) experience using OpenFOAM.


Joa
joaran is offline   Reply With Quote

Old   January 15, 2020, 15:51
Default Slow simulation
  #24
New Member
 
Aditya Srivastava
Join Date: Sep 2019
Posts: 1
Rep Power: 0
tatadocomo is on a distinguished road
I am also simulting rising bubble in interfoam. The mesh is very refined as of size D/40 where D is diameter of bubble. The p_rgh is taking 1000 iteration to converge which is slowing my simulation. Kindly suggest what to do
tatadocomo is offline   Reply With Quote

Reply

Tags
pimple

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
PIMPLE – the value of the final under-relaxation factor Zbynek OpenFOAM 7 January 20, 2020 10:38
pimpleDyMFoam computation randomly stops babapeti OpenFOAM Running, Solving & CFD 5 January 24, 2018 05:28
error while running modified pimple solver R_21 OpenFOAM Programming & Development 0 May 28, 2015 06:59
A question on the PIMPLE algorithm GerhardHolzinger OpenFOAM Running, Solving & CFD 4 February 13, 2015 06:49
Help for the small implementation in turbulence model shipman OpenFOAM Programming & Development 25 March 19, 2014 10:08


All times are GMT -4. The time now is 23:45.