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

Running stably with a large Courant number using pimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 23, 2009, 16:14
Default Running stably with a large Courant number using pimpleFoam
  #1
New Member
 
Michael Lawson
Join Date: Apr 2009
Location: NREL - Boulder, CO
Posts: 11
Rep Power: 8
mlawson is on a distinguished road
I am attempting to simulate laminar incompressible pulsatile flow in a complex geometry. To adequately resolve my geometry I need at least 15 million cells. I am driving the flow by setting the pressure and the outlet and oscillation pressure at the inlet.

Currently I am using pimpleFoam with automatic timestep adjustment to maintain a max Courant number of 2. Using a higher Courant number causes the solution to become unstable.

My problem is that at a max Courant number of 2 my timesteps are between 1E-5 and 1E-6 seconds. I need to simulate at least 1 second of time, which would require at least 100,000 timesteps at my current timestep. Given the size of my model, solving this many timesteps is computationally prohibitive.

I'm using fvSolution and fvSchemes dictionaries directly from the pimpleFoam tutorial in OF-1.6 (see attached).

Does anyone have any suggestions on how I could increase my timesteps by a couple orders of magnitude to make this simulation feasible? Should I be using a different solver? Should I be using a different discretization scheme?

Any help would be appreciated,
-Mike
Attached Files
File Type: txt fvSolution.pimpleFoam.txt (2.2 KB, 86 views)
File Type: txt fvSchemes.pimpleFoamTut.txt (1.9 KB, 56 views)
mlawson is offline   Reply With Quote

Old   July 13, 2012, 06:04
Default
  #2
Senior Member
 
Julien
Join Date: Jun 2012
Location: France
Posts: 153
Rep Power: 5
Djub is on a distinguished road
Hi Mike,

Did you solve your problem?
I am trying to run a case with a very large Reynolds number (2.5e8). So I would like also to run with a large Co (Courrent flow) to limit my time calculation.
Does anyone has a solution? Or maybe an explanation that what I am looking for does not exist (and cannot exist)? Maybe the Co cannot be greater thant 1 in an unsteady flow?
In that case, is there in OpenFoam any Explicit algorithm for unsteady, incompressible, turbulent, newtonian flow ?
Djub is offline   Reply With Quote

Old   November 24, 2014, 10:27
Default
  #3
Member
 
Christian
Join Date: Apr 2014
Posts: 59
Rep Power: 3
crixman is on a distinguished road
Dear all, any news on this topic?
I was wondering if there was a way to detach Co from the pimple solver, being able to use fixed (and larger!) time steps, as I think it's done in Fluent or other solvers.
Any ideas or suggestions?
crixman is offline   Reply With Quote

Old   November 24, 2014, 11:25
Default
  #4
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,132
Rep Power: 20
alexeym will become famous soon enoughalexeym will become famous soon enough
Use first-order schemes. With upwind scheme I was able to run a case (quite simple - von Karman vortex street) with Co up to 50.
alexeym is offline   Reply With Quote

Old   November 24, 2014, 13:05
Default
  #5
Member
 
Christian
Join Date: Apr 2014
Posts: 59
Rep Power: 3
crixman is on a distinguished road
Interesting to know, did you use turbulence schemes as well?
I could give it a try to start the simulation then switch to second order schemes, even if I'm not sure that the Co jump will be too big to be handled by the solver then!
crixman is offline   Reply With Quote

Old   November 24, 2014, 15:11
Default
  #6
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,132
Rep Power: 20
alexeym will become famous soon enoughalexeym will become famous soon enough
No, simulation was at rather low Re (100).
alexeym is offline   Reply With Quote

Old   November 24, 2014, 17:32
Default
  #7
Member
 
Christian
Join Date: Apr 2014
Posts: 59
Rep Power: 3
crixman is on a distinguished road
I'm making an extreme test, I'll keep you updated if I could reach an interesting value of Co (turbomachinery with kOmegaSST)!
crixman is offline   Reply With Quote

Old   November 24, 2014, 18:13
Default
  #8
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,100
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi, you should Check the pimple stuft at the Wikipedia Site.
__________________
Best regards,
Tobias Holzmann

Some interesting OpenFOAM tutorials and videos on www.Holzmann-cfd.de
Tobi is offline   Reply With Quote

Old   November 24, 2014, 18:28
Default
  #9
Member
 
Christian
Join Date: Apr 2014
Posts: 59
Rep Power: 3
crixman is on a distinguished road
Hi Tobi, sure, wanted to learn more and try out to see how far I can push it!
crixman is offline   Reply With Quote

Old   November 24, 2014, 18:33
Default
  #10
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,100
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
It Is also in my blog.
And before I ment the open foam wiki
__________________
Best regards,
Tobias Holzmann

Some interesting OpenFOAM tutorials and videos on www.Holzmann-cfd.de
Tobi is offline   Reply With Quote

Old   November 24, 2014, 18:38
Default
  #11
Member
 
Christian
Join Date: Apr 2014
Posts: 59
Rep Power: 3
crixman is on a distinguished road
I will look it through thanks and post some updates as soon as I got time!
crixman is offline   Reply With Quote

Old   November 29, 2014, 11:37
Default
  #12
Member
 
Christian
Join Date: Apr 2014
Posts: 59
Rep Power: 3
crixman is on a distinguished road
Hi Tobi, I read through your post about Pimple and I am now running at Co=12.5 and everythign seems to converge well!
As I am using upwind schemes for U,k and omega, when switching to linear schemes, should I drop my Co below 1?
crixman is offline   Reply With Quote

Reply

Tags
pimplefoam

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Mesh Refinement Luiz Eduardo Bittencourt Sampaio (Sampaio) OpenFOAM Mesh Utilities 41 January 17, 2013 03:43
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Courant number msrinath80 OpenFOAM Running, Solving & CFD 1 May 29, 2006 17:20
ENO scheme and Courant number. J. Main CFD Forum 2 November 8, 2002 18:07


All times are GMT -4. The time now is 04:21.