|
[Sponsors] |
January 16, 2017, 15:39 |
Solving Stokes Flows in OpenFOAM
|
#1 |
Senior Member
Join Date: Jul 2013
Posts: 124
Rep Power: 12 |
Hi all,
I am interested in solving Stokes flows (zero Reynolds number) using OpenFOAM. It's my understanding that there are several different ways to do this, for example you can just delete for example fvm::div(phi, U) from the UEqn in icoFoam and recompile and have a Stokes flow solver. However, I've noticed that doing so, the amount of time to converge is still a function of viscosity nu, which I would not expect. Are there other changes that need to be made to for example icoFoam in order to accurately solve zero Reynolds number flows, or better yet, what would be the optimum approach for solving steady state Stokes flows using OpenFOAM, i.e. time resolution is not important. I would appreciate any help. |
|
January 16, 2017, 21:29 |
|
#2 | |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34 |
Quote:
This greatly affect your simulation time. You have to use urf = 1 and depend on explicit urf. This comes with stability problems but a lower value of explicit urf would help. |
||
January 16, 2017, 21:40 |
|
#3 |
Senior Member
Join Date: Jul 2013
Posts: 124
Rep Power: 12 |
Hi arjun, thanks for the message. I'm not sure I fully understand, I did not know there were implicit and explicit under relaxation factors. I have only used under relaxation with simpleFoam. Do you mean that the best solution is to modify simpleFoam and eliminate the convection terms there? Could you please explain the difference between explicit and implicit urf?
Thanks again. |
|
January 16, 2017, 21:58 |
|
#4 | |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34 |
Quote:
I know it is implicit because all these solvers apply it in implicit form. |
||
January 16, 2017, 22:07 |
|
#5 |
Senior Member
Join Date: Jul 2013
Posts: 124
Rep Power: 12 |
Thanks, so if you wanted to solve a steady Stokes flow problem would you use simpleFoam?
|
|
January 16, 2017, 22:39 |
|
#6 | |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34 |
Quote:
You can solve stokes problem with simpleFOAM as the convective part will be calculated to be insignificant. Whole thing will be difficult to converge. |
||
January 16, 2017, 22:57 |
|
#7 |
Senior Member
Join Date: Jul 2013
Posts: 124
Rep Power: 12 |
That makes sense. So which solver would you use?
|
|
January 16, 2017, 23:33 |
|
#8 |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,273
Rep Power: 34 |
||
March 12, 2018, 10:04 |
|
#9 |
Member
Pablo Alarcón
Join Date: Mar 2018
Location: Liège
Posts: 59
Rep Power: 8 |
I'm also interested in this problem, but I want to construct a solver using (as a base) only the stokes flow, is there any possibility of having a solver considering only the stokes flow? or is it necessary to develop it from scratch?, something like taking the icofoam solver and adapt it with the SIMPLE algorithm?
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
problem with Min/max rho | tH3f0rC3 | OpenFOAM | 8 | July 31, 2019 09:48 |
settlingFoam unstable? | bendel_boy | OpenFOAM Running, Solving & CFD | 38 | July 8, 2016 05:07 |
Micro Scale Pore, icoFoam | gooya_kabir | OpenFOAM Running, Solving & CFD | 2 | November 2, 2013 13:58 |
calculation stops after few time steps | sivakumar | OpenFOAM Running, Solving & CFD | 7 | March 17, 2013 06:37 |
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 | bookie56 | OpenFOAM Installation | 8 | August 13, 2011 04:03 |