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

Time step too small with maxCo=1

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

Like Tree4Likes
  • 1 Post By tomf
  • 1 Post By fernexda
  • 2 Post By fernexda

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 11, 2015, 10:46
Default Time step too small with maxCo=1
  #1
Member
 
daniel fernex
Join Date: Oct 2014
Location: Braunschweig, Germany
Posts: 36
Rep Power: 11
fernexda is on a distinguished road
Dear all,
I'm running unsteady simulation of vertical axis wind turbine with sliding meshes and using the pimpleFoam solver with OF 2.1.

I'm trying to reproduce a simulation from a paper, using the same conditions (mesh size, inlet velocity, ...).

My problem is that when I try to refine the mesh in the boundary layer, the time step decreases like hell (down to 1e-9 when I'm solving the boundary layer to y+=1), and the simulation time prediction is many centuries... This seems normal because I'm fixing the Courant number to one. However, from what I've seen in the literature, such simulations can be achieved much faster with the same boundary layer resolution with other solvers such as fluent for instance.

So the question is: is this small time step OpenFOAM specific ? Does anyone know why such simulations can be achieved with fluent for example ? I'm stuck because of this, and forced not to resolve the boundary layer and to use a wall function which impacts the results...

Thanks in advance for your answers.

Regards,
Daniel
fernexda is offline   Reply With Quote

Old   May 12, 2015, 07:21
Default
  #2
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 634
Rep Power: 32
tomf will become famous soon enoughtomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi Daniel,

Have a look at the links below and start playing around with maxCo with different settings for nOuterCorrectors and tolerance levels and different values for the relaxation factors. Just make sure you do get the "Pimple loop converged in XX iterations" comment.

http://www.cfd-online.com/Forums/blo...hm-part-i.html
http://www.cfd-online.com/Forums/blo...m-part-ii.html

Regards,
Tom
tomf is offline   Reply With Quote

Old   May 14, 2015, 06:10
Default
  #3
Member
 
daniel fernex
Join Date: Oct 2014
Location: Braunschweig, Germany
Posts: 36
Rep Power: 11
fernexda is on a distinguished road
Hi Tom,
thanks for the answer. I've already looked at this awesome tutorial, but it didn't solve my problem.
  • If I use PISO, I keep maxCo to 1 which takes too long because of the small time steps.
  • If I use SIMPLE, I keep maxCo to 5, what increases the time step. This is good, but the inner iterations take more time until the 'pimple loop converges'. And if I wand the convergence to be reached faster I need to increase the residual tolerance and therefore decrease the accuracy. So in the end, it doesn't change a lot. And I fear that if I increase the Courant number too much, it will not converge. Do you know what's the limit ? From what I've read, it can be dangerous to set Co>5. Are there any other tricks to get a reasonable simulation time ?
Thanks !


Regards,

Daniel
fernexda is offline   Reply With Quote

Old   May 15, 2015, 08:37
Default
  #4
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 634
Rep Power: 32
tomf will become famous soon enoughtomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi Daniel,

It is pretty much case-specific on how high you can go. I have ran simulations with maxCo at 250 and still got accurate results, but on other projects I had to go back to somewhere around 10. The best way is running some tests to find the cut-off. You have some parameters to play around with (relaxation factors, tolerance levels, number of correctors/nonOrthogonalcorrerctors, etc.). My advice would be to just take a simplified version of your problem and do a sensitivity study, I expect that you will learn quite a lot.

Regards,
Tom
vivek05 likes this.
tomf is offline   Reply With Quote

Old   May 15, 2015, 10:41
Default
  #5
Member
 
daniel fernex
Join Date: Oct 2014
Location: Braunschweig, Germany
Posts: 36
Rep Power: 11
fernexda is on a distinguished road
Hi Tom,

that sounds like a good idea. I didn't know so high courant number could be used, depending on the case. I'll try to play around with the different parameters and see what's the maximum Co.

Thanks a lot,
Daniel
vivek05 likes this.
fernexda is offline   Reply With Quote

Old   May 15, 2015, 11:05
Default
  #6
Member
 
daniel fernex
Join Date: Oct 2014
Location: Braunschweig, Germany
Posts: 36
Rep Power: 11
fernexda is on a distinguished road
Hi,

I've juste found out that using maxCo=300 still gives me accurate results ! I've reduced my simulation time from 2 weeks to 8 hours !

I thought the long simulation times were OpenFOAM specific, but the solution is just to find the right parameters. That's great news, because I can stick with OpenFOAM for the future !

Regards,
Daniel
wind_ and vivek05 like this.
fernexda is offline   Reply With Quote

Reply

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
pimpleFoam: turbulence->correct(); is not executed when using residualControl hfs OpenFOAM Running, Solving & CFD 3 October 29, 2013 09:35
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 16:33
dynamic Mesh is faster than MRF???? sharonyue OpenFOAM Running, Solving & CFD 14 August 26, 2013 08:47
same geometry,structured and unstructured mesh,different behaviour. sharonyue OpenFOAM Running, Solving & CFD 13 January 2, 2013 23:40
pisoFoam with k-epsilon turb blows up - Some questions Heroic OpenFOAM Running, Solving & CFD 26 December 17, 2012 04:34


All times are GMT -4. The time now is 18:13.