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

role of the time step in convergence criteria !

Register Blogs Community New Posts Updated Threads Search

Like Tree8Likes
  • 1 Post By medaouarwalid
  • 1 Post By FMDenaro
  • 1 Post By medaouarwalid
  • 2 Post By piu58
  • 1 Post By medaouarwalid
  • 1 Post By juliom
  • 1 Post By medaouarwalid

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 11, 2017, 22:37
Wink role of the time step in convergence criteria !
  #1
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
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 !
Светлана likes this.
medaouarwalid is offline   Reply With Quote

Old   June 12, 2017, 05:36
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,769
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
The time step has a great relevance on the solution, both in terms of accuracy and in the numerical stability.
Светлана likes this.
FMDenaro is offline   Reply With Quote

Old   June 12, 2017, 10:12
Default
  #3
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
The time step has a great relevance on the solution, both in terms of accuracy and in the numerical stability.
please more details , How do I know in advance choose the right time step
Светлана likes this.
medaouarwalid is offline   Reply With Quote

Old   June 12, 2017, 11:59
Default
  #4
Senior Member
 
piu58's Avatar
 
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15
piu58 is on a distinguished road
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)
piu58 is offline   Reply With Quote

Old   June 12, 2017, 12:04
Default
  #5
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by piu58 View Post
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.
Thank you, what i am learning from this forum is larger than what i learnd in my entire years of study . What a great forum
Светлана likes this.
medaouarwalid is offline   Reply With Quote

Old   June 12, 2017, 15:03
Default
  #6
Senior Member
 
Julio Mendez
Join Date: Apr 2009
Location: Fairburn, GA. USA
Posts: 290
Rep Power: 18
juliom is on a distinguished road
Send a message via Skype™ to juliom
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
Светлана likes this.
juliom is offline   Reply With Quote

Old   June 12, 2017, 16:45
Default
  #7
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by juliom View Post
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
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 !?!?
Светлана likes this.

Last edited by medaouarwalid; June 12, 2017 at 18:09.
medaouarwalid is offline   Reply With Quote

Old   June 13, 2017, 00:19
Default
  #8
Senior Member
 
piu58's Avatar
 
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15
piu58 is on a distinguished road
> 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)
piu58 is offline   Reply With Quote

Old   June 13, 2017, 05:55
Default
  #9
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by piu58 View Post
> 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.
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??
medaouarwalid is offline   Reply With Quote

Old   June 13, 2017, 07:10
Default
  #10
Senior Member
 
piu58's Avatar
 
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15
piu58 is on a distinguished road
> 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)
piu58 is offline   Reply With Quote

Old   June 13, 2017, 07:15
Default
  #11
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by piu58 View Post
> 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.
Balance between cell zize and velocity ? I didnt hear this befor, can you explain to me ?
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 .
medaouarwalid is offline   Reply With Quote

Old   June 13, 2017, 07:48
Default
  #12
Senior Member
 
piu58's Avatar
 
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15
piu58 is on a distinguished road
Quote:
Originally Posted by medaouarwalid View Post
Balance between cell zize and velocity ? I didnt hear this befor, can you explain to me ?
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 .

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)
piu58 is offline   Reply With Quote

Old   June 13, 2017, 08:01
Default
  #13
Senior Member
 
dilaw meda
Join Date: Jun 2017
Location: algeria
Posts: 145
Rep Power: 8
medaouarwalid is on a distinguished road
Quote:
Originally Posted by piu58 View Post
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.
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
medaouarwalid is offline   Reply With Quote

Reply


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
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


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