|
[Sponsors] |
Simulation Time Vs Reynolds number for a given mesh |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 6, 2017, 07:51 |
Simulation Time Vs Reynolds number for a given mesh
|
#1 |
Senior Member
vidyadhar
Join Date: Jul 2016
Posts: 138
Rep Power: 9 |
Hello All,
I would like to know if there is any relation between simulation time (time taken by a machine to execute the code) and Reynolds number for a given mesh. It is taking 50 hrs of time to complete my simulation at Re 400 using 20 cores. I would like to know how much it would take for higher and lower Re. Thanks! Vidyadhar |
|
April 7, 2017, 01:17 |
|
#2 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
It depends on your model.
If you use a model which covers turbulency to some extend the solvers need more steps to converge if there is more turbulency in the flow. For LES simulations the time for calculation increases considerably with turbulency and therewith with the Reynold number. If you use a steady state solver the dependency between Re and the calculation time is not as steep. But you may get problems with convergence at very high Re numbers.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
April 8, 2017, 00:13 |
|
#3 |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,676
Rep Power: 66 |
For a fixed grid, there's no direct correlation. The time needed to do 1 iteration is dependent only on number of cells and the number of equations you solve, etc. The solver doesn't care what the flow Reynolds number is, it just cranks the algorithm the same.
For a steady case and fixed grid, the number of iterations you need to converge to the solution varies somewhat with Reynolds number due to non-linearities and it is hard to say whether you need more or less. For transient simulations, you usually decrease your time-step size as velocity increases in order to maintain the same temporal bandwidth. I.e. you double the velocity, you half the time-step size and do twice as many time-steps. In this case it's linearly proportional to Reynolds. |
|
April 8, 2017, 01:50 |
|
#4 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
Dear LuckyTran
> The time needed to do 1 iteration is dependent only on number of cells a If you use a direct solver for the system of equation this is true. But if you use an iterative solver like conjugated gradients the number of iterations depend on the condition of the system. This condition is more worse with much changing in the flow. Some more advanced solvers for the fluid dynamics may be controlle by the Co number and select the time step by theirself. With more turbulence in the flow they choose lower time steps and need more computational time. This is equivalent to your statement > you usually decrease your time-step size as velocity increases in order to maintain the same temporal bandwidth
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
April 8, 2017, 09:09 |
|
#5 | |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,676
Rep Power: 66 |
The automatic Co is just the software picking time-step for the user because the user is too lazy to do it for themself. There are a lot of acceleration schemes, there are a lot of tricks that can likely speed up and speed down the process. My claim is that these are indirect behaviors, and that the only meaningful difference in solution time is the number of time-steps you do in a transient simulation.
Quote:
As for the pre-conditioning matrix, it is also nearly equivalent to saying the solution time now depends on good/bad is your initial guess. So if you change the Reynolds number, you should also change your initial guess so that the iterative solvers do the same number of iterations. Now where in the algorithm (direct/indirect) is there a switch that says, ah-hah for a different Reynolds number I must do something different? My claim, although it is not 100% true, is that there's not really such a switch because the linear system you are trying to solve is the same (the discretized N-S with all the time and space discretizations). You have these little quirks that make a little bit of difference here and there, but most of the time you are just solving a linear system Ax=b |
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
courant number increases to rather large values | 6863523 | OpenFOAM Running, Solving & CFD | 22 | July 5, 2023 23:48 |
[Other] Contribution a new utility: refine wall layer mesh based on yPlus field | lakeat | OpenFOAM Community Contributions | 58 | December 23, 2021 02:36 |
AMI speed performance | danny123 | OpenFOAM | 21 | October 24, 2020 04:13 |
Moving mesh | Niklas Wikstrom (Wikstrom) | OpenFOAM Running, Solving & CFD | 122 | June 15, 2014 06:20 |
IcoFoam parallel woes | msrinath80 | OpenFOAM Running, Solving & CFD | 9 | July 22, 2007 02:58 |