|
[Sponsors] |
November 1, 2017, 07:20 |
Timestep size and Convergence
|
#1 |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
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 |
|
November 1, 2017, 07:35 |
|
#2 |
Senior Member
david
Join Date: Oct 2012
Posts: 142
Rep Power: 14 |
The CFL condition should get you your answer about the stability of your solution and the time step to use
|
|
November 1, 2017, 07:39 |
|
#3 |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
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 |
|
November 1, 2017, 07:45 |
|
#4 |
Senior Member
david
Join Date: Oct 2012
Posts: 142
Rep Power: 14 |
You can actually specify the courant number in fluent with the density based solver
|
|
November 1, 2017, 07:49 |
|
#5 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Quote:
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? |
||
November 1, 2017, 07:52 |
|
#6 |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
Yes sir, I am actually working on pressure based solver. 2D Laminar flow over bluff-bodies.
__________________
Always Thedal |
|
November 1, 2017, 07:57 |
|
#7 |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
2D laminar flow over multiple bluff bodies.
__________________
Always Thedal |
|
November 1, 2017, 08:21 |
|
#8 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
what is exactly your problem? The solution is always diverging or only for a threshold in the value of dt?
|
|
November 1, 2017, 08:45 |
|
#9 |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
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 |
|
November 1, 2017, 08:55 |
|
#10 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
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 |
|
November 1, 2017, 09:02 |
|
#11 |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
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 |
|
November 1, 2017, 09:17 |
|
#12 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Diverge for further iterations at this time or you see divergence for longer time? Your flow will develop a fully unsteady solution
|
|
November 1, 2017, 09:27 |
|
#13 |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
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 |
|
November 1, 2017, 09:45 |
|
#14 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
the flow appears no yet fully developed... it could be due to an excess of numerical viscosity, what about the spatial discretiazion?
|
|
November 2, 2017, 01:50 |
|
#15 |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
Spacial Discretization:
Gradient - Least Square Cell based Pressure - Second order Momentum - Second order upwind
__________________
Always Thedal |
|
November 2, 2017, 02:24 |
|
#16 |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,285
Rep Power: 34 |
||
November 2, 2017, 04:14 |
|
#17 |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
Can you just elaborate me please, how to do it?
__________________
Always Thedal |
|
November 2, 2017, 04:28 |
|
#18 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
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? |
|
November 2, 2017, 04:40 |
|
#19 | |
Member
Thamilmani M
Join Date: Sep 2017
Location: IIT Bombay, Mumbai
Posts: 52
Rep Power: 9 |
Quote:
The mesh is actually symmetrical only about X axis. Why does it impact?
__________________
Always Thedal |
||
November 2, 2017, 04:42 |
|
#20 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Quote:
The cell Re number seems adequate, what about the inflow velocity? |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Convergence | Centurion2011 | FLUENT | 48 | June 15, 2022 00:29 |
RMS Residual value and Timestep size | Binoy | CFX | 2 | August 7, 2017 13:22 |
Domain size convergence | CRAC | CFX | 1 | April 3, 2016 21:25 |
Time step size & convergence absolute criteria | yuitsang | FLUENT | 5 | April 15, 2013 05:27 |
Combustion Convergence problems | Art Stretton | Phoenics | 5 | April 2, 2002 06:59 |