|
[Sponsors] |
May 11, 2015, 09:46 |
Time step too small with maxCo=1
|
#1 |
Member
daniel fernex
Join Date: Oct 2014
Location: Braunschweig, Germany
Posts: 36
Rep Power: 11 |
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 |
|
May 12, 2015, 06:21 |
|
#2 |
Senior Member
|
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 |
|
May 14, 2015, 05:10 |
|
#3 |
Member
daniel fernex
Join Date: Oct 2014
Location: Braunschweig, Germany
Posts: 36
Rep Power: 11 |
Hi Tom,
thanks for the answer. I've already looked at this awesome tutorial, but it didn't solve my problem.
Regards, Daniel |
|
May 15, 2015, 07:37 |
|
#4 |
Senior Member
|
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 |
|
May 15, 2015, 09:41 |
|
#5 |
Member
daniel fernex
Join Date: Oct 2014
Location: Braunschweig, Germany
Posts: 36
Rep Power: 11 |
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 |
|
May 15, 2015, 10:05 |
|
#6 |
Member
daniel fernex
Join Date: Oct 2014
Location: Braunschweig, Germany
Posts: 36
Rep Power: 11 |
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 |
|
|
|
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 08:35 |
How to write k and epsilon before the abnormal end | xiuying | OpenFOAM Running, Solving & CFD | 8 | August 27, 2013 15:33 |
dynamic Mesh is faster than MRF???? | sharonyue | OpenFOAM Running, Solving & CFD | 14 | August 26, 2013 07:47 |
same geometry,structured and unstructured mesh,different behaviour. | sharonyue | OpenFOAM Running, Solving & CFD | 13 | January 2, 2013 22:40 |
pisoFoam with k-epsilon turb blows up - Some questions | Heroic | OpenFOAM Running, Solving & CFD | 26 | December 17, 2012 03:34 |