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

Timestep size and Convergence

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 1, 2017, 06:20
Lightbulb Timestep size and Convergence
  #1
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
Does time-step size affect the convergence? Like, if a del t is beyond some limit, will it diverge?
If so, How to calculate the optimum time-step size for an unsteady flow problem?
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 1, 2017, 06:35
Default
  #2
Senior Member
 
david
Join Date: Oct 2012
Posts: 142
Rep Power: 13
davidwilcox is on a distinguished road
The CFL condition should get you your answer about the stability of your solution and the time step to use
davidwilcox is offline   Reply With Quote

Old   November 1, 2017, 06:39
Default
  #3
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
Sorry, I am working with ANSYS, where Courant number can't be found before going about the solution. Also, My discretization schemes are Second-order implicit. Hence, Will Courant number affect?
And, Thank you so much for your reply, sir.
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 1, 2017, 06:45
Default
  #4
Senior Member
 
david
Join Date: Oct 2012
Posts: 142
Rep Power: 13
davidwilcox is on a distinguished road
You can actually specify the courant number in fluent with the density based solver
davidwilcox is offline   Reply With Quote

Old   November 1, 2017, 06:49
Default
  #5
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by thedal View Post
Sorry, I am working with ANSYS, where Courant number can't be found before going about the solution. Also, My discretization schemes are Second-order implicit. Hence, Will Courant number affect?
And, Thank you so much for your reply, sir.

The time step enters into the accuracy and stability of a numerical scheme.
If a scheme is consistent and stable, you have convergence.
Explicit schemes are subject to stability constraints that depend both on the convective part (CFL) and the diffusive part of the equations. Implicit schemes are generally unconditionally stable and the time step affects the accuracy.
What kind of flow problem are you simulating?
FMDenaro is offline   Reply With Quote

Old   November 1, 2017, 06:52
Default
  #6
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
Yes sir, I am actually working on pressure based solver. 2D Laminar flow over bluff-bodies.
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 1, 2017, 06:57
Default
  #7
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
2D laminar flow over multiple bluff bodies.
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 1, 2017, 07:21
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
what is exactly your problem? The solution is always diverging or only for a threshold in the value of dt?
FMDenaro is offline   Reply With Quote

Old   November 1, 2017, 07:45
Default
  #9
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
Mine, Flow over multiple square cylinders at low Reynolds number cases.
For del t 0.01, I have run it for more than 20000 timesteps. The solution has not converged. So, I tried to run at del t 0.1, But even then the problems dont converge if sbyD i.e., the distance between cylinders/D is 3 or lesser than that. I have tried to give for that del t 0.1 and some 30000 timesteps. Yet I find the residuals actually diverge after 20000 time steps. And It never look like converging.
Hence, My doubt. the del t and convergence has any relation?
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 1, 2017, 07:55
Default
  #10
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
But do you mean convergence towards a steady state? And how you are sure that your flow problem has a physical steady state at that Re number?
A single bluff body has onset of unsteady vortex shedding for Re as low as 40-50
FMDenaro is offline   Reply With Quote

Old   November 1, 2017, 08:02
Default
  #11
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
Sir, I am working at Re 100.
Yes. Convergence towards a steady state.
And I am trying to validate a simulation and experimental work already published. Hence, I think it should give proper vortex shedding at that Re.
This is the residuals for del t = 0.1 for Re 100 flow over two cylinders separated by 3D. It actually start diverging lately. That I don't know whether it will converge eventually or not (i.e., vortex will start shedding steadily).

https://www.dropbox.com/s/ow6ul3wikg...duals.png?dl=0
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 1, 2017, 08:17
Default
  #12
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Diverge for further iterations at this time or you see divergence for longer time? Your flow will develop a fully unsteady solution
FMDenaro is offline   Reply With Quote

Old   November 1, 2017, 08:27
Default
  #13
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
For a longer time.
It looks like the vortex never sheds also.
I have attached the vorticity contours also at the end of 2300 s.

https://www.dropbox.com/s/p4bc7krs50...02cyl.jpg?dl=0
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 1, 2017, 08:45
Default
  #14
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
the flow appears no yet fully developed... it could be due to an excess of numerical viscosity, what about the spatial discretiazion?
FMDenaro is offline   Reply With Quote

Old   November 2, 2017, 00:50
Default
  #15
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
Spacial Discretization:

Gradient - Least Square Cell based
Pressure - Second order
Momentum - Second order upwind
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 2, 2017, 01:24
Default
  #16
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,274
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by thedal View Post
Spacial Discretization:

Gradient - Least Square Cell based
Pressure - Second order
Momentum - Second order upwind

Shift your mesh in any direction by 1.0E-10 m and try again.
arjun is offline   Reply With Quote

Old   November 2, 2017, 03:14
Default
  #17
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
Quote:
Originally Posted by arjun View Post
Shift your mesh in any direction by 1.0E-10 m and try again.
Can you just elaborate me please, how to do it?
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 2, 2017, 03:28
Default
  #18
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by thedal View Post
Can you just elaborate me please, how to do it?

I suppose that he wants you to "destroy" the grid symmetry with respect to the body.
However, how about the cell Re number? It should be O(1) everywhere to ensure a very low numerical diffusion.
What about the inflow profile?
arjun likes this.
FMDenaro is offline   Reply With Quote

Old   November 2, 2017, 03:40
Default
  #19
Member
 
thedal's Avatar
 
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 8
thedal is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I suppose that he wants you to "destroy" the grid symmetry with respect to the body.
However, how about the cell Re number? It should be O(1) everywhere to ensure a very low numerical diffusion.
What about the inflow profile?
https://www.dropbox.com/s/5kg6b25vr6...%20Re.jpg?dl=0

The mesh is actually symmetrical only about X axis. Why does it impact?
__________________
Always
Thedal
thedal is offline   Reply With Quote

Old   November 2, 2017, 03:42
Default
  #20
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,777
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by thedal View Post
https://www.dropbox.com/s/5kg6b25vr6...%20Re.jpg?dl=0

The mesh is actually symmetrical only about X axis. Why does it impact?

The cell Re number seems adequate, what about the inflow velocity?
FMDenaro 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
Convergence Centurion2011 FLUENT 48 June 14, 2022 23:29
RMS Residual value and Timestep size Binoy CFX 2 August 7, 2017 12:22
Domain size convergence CRAC CFX 1 April 3, 2016 20:25
Time step size & convergence absolute criteria yuitsang FLUENT 5 April 15, 2013 04:27
Combustion Convergence problems Art Stretton Phoenics 5 April 2, 2002 05:59


All times are GMT -4. The time now is 20:05.