|
[Sponsors] |
role of the time step in convergence criteria ! |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 11, 2017, 22:37 |
role of the time step in convergence criteria !
|
#1 |
Senior Member
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8 |
hi, excuse me for my questions, i am new in simulations, i'am using fluent to simulate an air jet, i use transient time study and K-epsilon model, when i set the time step to 0.0001s the solution converged rapidally , and every 20 iterations i get the message "solution converged" , with 10^-4 for continuity , 10^-6 for velocities and 10^-8 for energie , but when i increased the time step to 0.01 it stills not converged until now. could any one explain to me please !
|
|
June 12, 2017, 05:36 |
|
#2 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,769
Rep Power: 71 |
The time step has a great relevance on the solution, both in terms of accuracy and in the numerical stability.
|
|
June 12, 2017, 10:12 |
|
#3 |
Senior Member
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8 |
||
June 12, 2017, 11:59 |
|
#4 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
You should understand the concept of the Courant number. Some of the solvers control the time step by this number, which is a fine idea.
Than you have to think about the relation between Co number and convergence All in all it can be said: The finer the mesh the lower the max Co number (or the time step) need to be. A smaller time step may take less wall time for calculating a solution.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
June 12, 2017, 12:04 |
|
#5 | |
Senior Member
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8 |
Quote:
|
||
June 12, 2017, 15:03 |
|
#6 |
Senior Member
|
Besides what professor Denaro said, the time step is critique in explicit schemes but not in implicit scheme. The latter is stable at high time step. However, is not as accurate as the explicit one
|
|
June 12, 2017, 16:45 |
|
#7 |
Senior Member
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8 |
I am using a bounded secpnd order implicit fomulation + i found that best time step need to be less that ( min size/ velocity) , in my case it means less than 0.0001 !!! So to see the result of 1 hour i need 36 milion iterations !?!?
Last edited by medaouarwalid; June 12, 2017 at 18:09. |
|
June 13, 2017, 00:19 |
|
#8 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
> So to see the result of 1 hour i need 36 milion iterations !?!?
Look at your mesh and at your boundary conditions. Something is wrong here.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
June 13, 2017, 05:55 |
|
#9 |
Senior Member
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8 |
Why? I have the smalest size of the mesh is 0,000875m, velocity is 7m/s so the recomanded time step will be less than 0.0001 . This time step will incrementes every time the solution converge ( in my case it was every 9 iterations when time step was 10^-5) so to reach 1 minutes i need 6 milions time step, true or not ?? Or i dont need to wait until 1 minute to see the result of the velocities??
|
|
June 13, 2017, 07:10 |
|
#10 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
> Why?
Because it is unnormal that a well formulated case converges such slowly. There may be skewed cells or an inappropriate balance between cell size and velocity. Nonphysical boundary conditions lead often to convergence problems. That are the directions I go first.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
June 13, 2017, 07:15 |
|
#11 | |
Senior Member
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8 |
Quote:
Boundry conditions are simple in my case , i have a room of dimensions 2 meter heigh, 4meter length and 4 meter width, and an air jet blowing vertically through a nozzle placed in the ceiling, so the ceiling i set it wall, the inlet is in the liddle of the ceiling , walls of the room as pressure outlet because the room is big , floor also as pressure outlet . |
||
June 13, 2017, 07:48 |
|
#12 | |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
Quote:
What means pressure outlet? Constant pressure? Waht velocity b.c. did you set? Do you have a compressible medium? If not it may be you move in volume without any way to let it out.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
||
June 13, 2017, 08:01 |
|
#13 |
Senior Member
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8 |
I did pressure outlet at the walls of the room because the velocity near to it was close to zero, also the pressure was constant in the room (atmospherique pressure) , the inlet like i said is a nozzle placed in the middle of the ceiling, i set it velocity inlet and give it 7m/s intensity 4% and hydraulique diameter 0.004m
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Transient simulation not converging | skabilan | OpenFOAM Running, Solving & CFD | 14 | December 16, 2019 23:12 |
p_rgh initial residual no change with different settings | manuc | OpenFOAM Running, Solving & CFD | 3 | June 26, 2018 15:53 |
High Courant Number @ icoFoam | Artex85 | OpenFOAM Running, Solving & CFD | 11 | February 16, 2017 13:40 |
Floating point exception error | lpz_michele | OpenFOAM Running, Solving & CFD | 53 | October 19, 2015 02:50 |
plot over time | fferroni | OpenFOAM Post-Processing | 7 | June 8, 2012 07:56 |