CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

Navier-Stokes time step size

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

LinkBack Thread Tools Search this Thread Display Modes
Old   June 4, 2008, 12:14
Default Navier-Stokes time step size
Posts: n/a

I'm implementing a 3D finite volume solver for the Navier-Stokes (NS) equations. I'm using a 4-stage RK for the time integration and I would like the time step to be adaptive and automatic, i.e. determined run time per time step.

However, I don't know how to choose an appropriate time step for NS. I'm not comfortable inferring a time step from some combination of, say, the time step limits for the 1D wave and heat equations using various schemes. For one thing I could arrive at a condition more restrictive than necessary. For example, I found that Euler's scheme for the wave equation is unconditionally unstable whereas it is conditionally stable for advection-diffusion. Or even worse, I could make a choice that is not restrictive enough.

The easiest approach would be to take the time step as the minimum of the convective and diffusive limits. But is this good enough ? And how is the diffusive time step limit really computed in 3D?

Thanks in advance
  Reply With Quote

Old   June 5, 2008, 13:27
Default Re: Navier-Stokes time step size
Posts: n/a
Dear Martin,

The time step in 3D is often chosen as the minimum of timesteps determined from convective and diffusive limits. There is no major difference between 2d and 3d definitions, only the characteristic length scale (when considering arbitrary ploygonal meshes) are a little different. Yo could compute dt_vis from the concept that VNN= nu*dt_vis/ char_len^2, where nu is the kinematic viscosity and char_len is a characteristic length for a volume, typically chosen as sqrt(cell_area) of cubeorot(cell_volume) depending on 2d/3d. Typically you could set VNN=CFL and get dt_vis for every cell. You could then set dt=min(dt_inv_min,dt_vis_min), which is global time stepping strategy or use the smallest of dt_inv and dt_vis in each cell as its own dt; which is local time stepping. For NS solvers, local time stepping relieves much of the grid--induced stiffness and results in a more rapid convergence. Of course the CFL in your case would be bounded from above by 1. In our codes however, we do not invoke the diffusive time limit even for a NS solver, we generally employ implicit time integration with a convective time limit coupled with local time stepping to achieve rapid steady state viscous solutions.

Hope this helps.



  Reply With Quote

Old   June 6, 2008, 03:38
Default Re: Navier-Stokes time step size
Posts: n/a
Hello Ganesh,

Thank you very much for your reply. It confirms my beliefs on how this is done. For now I have included only explicit time stepping but implicit methods are certainly something to consider once I get my basic solver validated.

Regards, Martin
  Reply With Quote


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
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 05:36
Time step, Number of time step, Maxximum Iterations per time step sandisk FLUENT 0 July 18, 2011 02:57
critical error during installation of openfoam Fabio88 OpenFOAM Installation 21 June 2, 2010 03:01
natural convection - time step size co2 FLUENT 7 June 3, 2004 21:56
how to increase time step size? co2 FLUENT 6 May 17, 2004 07:25

All times are GMT -4. The time now is 19:36.