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

Adaptive time step

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 27, 2005, 23:26
Default Adaptive time step
  #1
styoung317
Guest
 
Posts: n/a
Hi,everyone! The following question confused me: In a dynamic process, we may set the time step according to Courant's condition to assure stability , |v*dt/dx|<1, in which v is the largest characteristic velocity; as the time ellapses, the velocity is changing ,too, so can we change the time step accordingly? Or some implicit principle exists? Thanks in advance. Young
  Reply With Quote

Old   June 28, 2005, 03:16
Default Re: Adaptive time step
  #2
shekharc
Guest
 
Posts: n/a
off course u can change time-step everytime.
  Reply With Quote

Old   June 28, 2005, 08:41
Default Re: Adaptive time step
  #3
styoung317
Guest
 
Posts: n/a
Thank you for response! Actually,I found that it was not feasible. For example,if a time step of 1e-5 is used, after an interval of 0.2s or so,the value of velocity will be one tenth of the initial one,but a corresponding time step of 1e-4 results in divergence.

My strategy of pressure correction may account for that,I think it somewhat affects the stability of calculation.
  Reply With Quote

Old   June 28, 2005, 10:46
Default Re: Adaptive time step
  #4
shekharc
Guest
 
Posts: n/a
u keep n eye on each component of the velocity throughout all the grids, evaluate u/dx, v/dy, and w/dz for all the grids. Now, say c = max(|u|/dx, |v|/dy, |w|/dz). So, u use dt = CFL/c.

For viscous flow calculations, there is more constrain u shud take into account...the Grid-Fourior number. dt < [dx*dx*dy*dy*dz*dz]/[2*KIN_VISC*(dx*dx + dy*dy + dz*dz)].

At last, your required dt would be the minimum of the two.

  Reply With Quote

Old   June 28, 2005, 16:13
Default Re: Adaptive time step
  #5
ganesh
Guest
 
Posts: n/a
Dear styoung317,

One way of solving dynamic problems or unsteady problems is to use the concept of pseudo-time stepping or Dual time stepping by which you introduce a dual step in addition to each physical step. At each physical time step, you will be trying to drive the pseudo-time derivative to zero, thereby solving a steady problem and then for this case you can use implicit time stepping and accelerate convergence. Of course in sucha case the physical time step is generally held constant to a vlue based on some physical considerations the problem possesses. Trying a fully explicit scheme is not wrong, but as you would have already seen, it would go too slow, for your time step will be restricted by physical considerations. Even for a steady case we compute the dt as you mentioned by the CFL criterion, using the largest eigen value, which again is a function of the solution and will change with iteration. Hope this helps.

Regards,

Ganesh
  Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
directMapped problem panda60 OpenFOAM Bugs 4 July 8, 2010 11:23
Time step in transient simulation shib FLUENT 0 June 17, 2010 14:07
DPM UDF particle position using the macro P_POS(p)[i] dm2747 FLUENT 0 April 17, 2009 02:29
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 12:16
natural convection - time step size co2 FLUENT 7 June 3, 2004 22:56


All times are GMT -4. The time now is 06:54.