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

Minimum Time Step with Navier Stokes

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 30, 2013, 00:48
Default Minimum Time Step with Navier Stokes
  #1
New Member
 
Join Date: Nov 2013
Posts: 5
Rep Power: 4
Atreides is on a distinguished road
Hello,

I am a bit of a newbie at CFD. I am trying to determine what the minimum time step would be when solving the N-S equations for a lid-driven cavity flow (2D) in order to achieve stability. The flow is incompressible, viscous, and the results I am interested in are at steady state. I am using a finite difference method with a co-located grid.

I'm really not sure how to even approach this problem and any advice on what I need to do to find the time step would be greatly appreciated.

Thank you very much for any and all replies.
Atreides is offline   Reply With Quote

Old   November 30, 2013, 05:27
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,417
Rep Power: 31
FMDenaro will become famous soon enoughFMDenaro will become famous soon enough
Quote:
Originally Posted by Atreides View Post
Hello,

I am a bit of a newbie at CFD. I am trying to determine what the minimum time step would be when solving the N-S equations for a lid-driven cavity flow (2D) in order to achieve stability. The flow is incompressible, viscous, and the results I am interested in are at steady state. I am using a finite difference method with a co-located grid.

I'm really not sure how to even approach this problem and any advice on what I need to do to find the time step would be greatly appreciated.

Thank you very much for any and all replies.
The numerical stability region (Reh,cfl) is determined by the type of numerical scheme you use, each one producing a different map of stability.
I don't see your type of discretization, Re number and number of cells.
I can suggest a rude approximation for evaluating a practical dt in case of Re sufficiently high:

dt_p*(umax/dx+vmax/dy)<1 --> dt = alpha*dt_p (alpha = 0.1 - 0.5)
FMDenaro is offline   Reply With Quote

Old   November 30, 2013, 11:18
Default
  #3
New Member
 
Join Date: Nov 2013
Posts: 5
Rep Power: 4
Atreides is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
The numerical stability region (Reh,cfl) is determined by the type of numerical scheme you use, each one producing a different map of stability.
I don't see your type of discretization, Re number and number of cells.
I can suggest a rude approximation for evaluating a practical dt in case of Re sufficiently high:

dt_p*(umax/dx+vmax/dy)<1 --> dt = alpha*dt_p (alpha = 0.1 - 0.5)
Yes, yes, my mistake. I am using CDS for all discretizations as per my instructions. My Reynolds number is only 100 (mu is 0.01, rho is 1, velocity is 1).

I'm not sure exactly what you mean by dt_p? And the value of alpha is dependent on Re?
Atreides is offline   Reply With Quote

Old   November 30, 2013, 11:26
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,417
Rep Power: 31
FMDenaro will become famous soon enoughFMDenaro will become famous soon enough
what about the time-integration scheme?
FMDenaro is offline   Reply With Quote

Old   November 30, 2013, 11:40
Default
  #5
New Member
 
Join Date: Nov 2013
Posts: 5
Rep Power: 4
Atreides is on a distinguished road
I am using an Explicit Time Advance scheme.

I calculate the combination of the advective and viscous terms and its divergence from the initial velocity field.

I solve Poisson's equation for pressure and then compute the velocity field at the next time step.

I hope that answers your question.
Atreides is offline   Reply With Quote

Old   November 30, 2013, 12:42
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,417
Rep Power: 31
FMDenaro will become famous soon enoughFMDenaro will become famous soon enough
Quote:
Originally Posted by Atreides View Post
I am using an Explicit Time Advance scheme.

I calculate the combination of the advective and viscous terms and its divergence from the initial velocity field.

I solve Poisson's equation for pressure and then compute the velocity field at the next time step.

I hope that answers your question.
If you are using a first order Explicit time integration, Re=100 will drive towards a diffusive stability constraint such as:

dt_d*(1/dx^2 + 1/dy^2)/Re <1/2.

As your sistem is non linear such value should be furhter reduced
FMDenaro is offline   Reply With Quote

Old   November 30, 2013, 12:52
Default
  #7
New Member
 
Join Date: Nov 2013
Posts: 5
Rep Power: 4
Atreides is on a distinguished road
Well perhaps I have other errors in my code but this produces a far larger time step than I have been using and it is immediately unstable. From where did you derive this equation for stability?

If I have a square that I have divided into a grid of 100 points with delta x = delta y = 0.1 that gives me a time step of 1 second per your equation?
Atreides is offline   Reply With Quote

Old   November 30, 2013, 13:23
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,417
Rep Power: 31
FMDenaro will become famous soon enoughFMDenaro will become famous soon enough
Quote:
Originally Posted by Atreides View Post
Well perhaps I have other errors in my code but this produces a far larger time step than I have been using and it is immediately unstable. From where did you derive this equation for stability?

If I have a square that I have divided into a grid of 100 points with delta x = delta y = 0.1 that gives me a time step of 1 second per your equation?

The constraint comes from the linear stabilty analysis of the parabolic equation dphi/dt=Gamma* Lap Phi, discretized with the FTCS scheme
Assuming a grid formed by 10x10 cells (dx=dy=0.1) you simply have
dt_d < 0.25. Therefore dt =0.1 should work.
If you see instability after few time-steps then there is an error in the code
FMDenaro is offline   Reply With Quote

Old   November 30, 2013, 13:27
Default
  #9
New Member
 
Join Date: Nov 2013
Posts: 5
Rep Power: 4
Atreides is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
The constraint comes from the linear stabilty analysis of the parabolic equation dphi/dt=Gamma* Lap Phi, discretized with the FTCS scheme
Assuming a grid formed by 10x10 cells (dx=dy=0.1) you simply have
dt_d < 0.25. Therefore dt =0.1 should work.
If you see instability after few time-steps then there is an error in the code
I will investigate. Thanks very much for the help.
Atreides is offline   Reply With Quote

Reply

Tags
finite difference method, navier-stokes, stability

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Rapidly decreasing deltaT for interDyMFoam chrisb2244 OpenFOAM Running, Solving & CFD 3 July 1, 2014 16:40
Sudden jump in Courant number NJG OpenFOAM Running, Solving & CFD 7 May 15, 2014 13:52
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 9 January 3, 2012 09:06
Full pipe 3D using icoFoam cyberbrain OpenFOAM 4 March 16, 2011 10:20
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


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