
[Sponsors] 
December 6, 2010, 02:09 
calculation of time step

#1 
Member
Aamer Shahzad
Join Date: Mar 2010
Posts: 58
Rep Power: 7 
Hello all.....
i am trying to calculate coefficient of lift and drag on corrugated airfoil of insect at Reynolds number of 34000. the case is inherently unsteady. i am trying 2 turbulence models at present. i.e SA and KW SST transitional. i have following questions. 1. How can i calculate the correct time step size, no of time steps and iterations per time step for my case. (any reference, tutorial or book etc ??? or any advice on this) 2. What is the residual criteria that has to be set for unsteady simulation.. is it absolute or relative, and what value approximately? 3. What is the role of Courent number in determining the correct time step size.? 

December 6, 2010, 07:18 

#2 
Member
Dan
Join Date: Oct 2010
Location: UK
Posts: 41
Rep Power: 6 
Below is my derrivation for the stability requirements of for direct solution of the navier stokes equations on a cartesian grid,
The dimensionalised navier stoke equation for the conservation of momentum, is given by for stability it is required that at all points, sub, where is some suitible factor of saftey, then. assuming each term has equal weight then, rearranging gives four requirements for Viscous term requirement the discretization of $\partial ^2 u$ is given by, which will always have magnitude, substituting gives, which simplifies to, Diffusion term requirement the discretization of is given by, which will always have magnitude; substituting gives, which simplifies to, Pressure term Requirment can be discretized as, therfore substitution gives, External Force term Requirment therfore if the bellow inequality is satisfied then the discretized equations will be stable, where is the maximum velocity found in the system, A slighly more relaxed stability requirement for the 2 dimensional case was given in "Griebel et al: Numerical Simulation in Fluid Dynamics: SIAM 1998" but there was no acompanying derrivation. Hope this helps 

December 6, 2010, 08:53 

#3 
Member
Aamer Shahzad
Join Date: Mar 2010
Posts: 58
Rep Power: 7 
Dear dandalf....
thanks for the quick reply..... but ur reply gives mathematical solution.... how to apply this to a practical problem in fluent. how can i calculate the time steps size, no of time steps and no of iterations per time step, while looking at ur derivations??? 

December 6, 2010, 09:37 

#4 
Member
Dan
Join Date: Oct 2010
Location: UK
Posts: 41
Rep Power: 6 
Presumably you know the reuired spatial resolution to capture the details, of your flow. So you have some
which is your smallest spacial discretization on your mesh, if not you should be able to find a suitabble spatial resolution using the relation, you can then simply substitute this into the abouve equations to get a stable which is your time step. as for the number of itterations required, after each itteration you need to check some convergence criteria, if for instance you are solving an incompressible flow using the partial step method, you can check for conservation of mass after each itteration of the pressure poisson, and keep going unill it is below some acceptable error. As for total number of time steps, in your case with an insect wing, you have a known time period for 1 beat cycle, you have calculated a suitable time step. it should be fairly trivial to calculate a number of time steps to cover a sufficient number of beat cycles to convince yourself that you have captured all unsteady beat interaction effects. If this hasn't helped You might find some good sugestions of how to proceed in, Rajat Mittal, Computational Modeling in BioHydrodymanics: Trends, Chalenges and Recent Advances,IEEE JOURNAL OF OCIANIC ENGINEERING, 29 , 595604, 2004. 

December 7, 2010, 05:11 

#5 
New Member
Join Date: Aug 2010
Posts: 4
Rep Power: 6 
i am sorry but this explanation ist also so difficult.
__________________
Mechanical Man 

December 7, 2010, 07:06 
Hi

#6 
Member
Dynampally Pavitran
Join Date: Mar 2010
Location: India
Posts: 73
Rep Power: 7 
1. How can i calculate the correct time step size, no of time steps and iterations per time step for my case. (any reference, tutorial or book etc ??? or any advice on this)
1. Time step size should be decided on how much temporal accuracy is enough, usually you need to do a timestep dependency study. 2. No. of time steps decides your total duration( usually 3 to 4 flow through times is good enough). 3. No. of iterations per time step (In each time step your residuals should be less than 1E4), usually this is achieved in less than 10 iterations, actually it depends on the problem also. 2. What is the role of Courent number in determining the correct time step size. 1. Courant number which is also called CFL number is stability criteria for unsteady simulations. it depends on timestep size & cellsize. 2. for explicit solvers you should see that CFL < 1 and for implicit solvers it can be higher. Information regarding this can be found in any CFD text book (eg: CFD by John D Anderson) 

December 8, 2010, 03:27 
Finding smallest division in mesh

#7 
New Member
Irshad
Join Date: Dec 2009
Location: Bangalore
Posts: 20
Rep Power: 7 
Hi guys,
I too have the same doubt in finding the smallest element size "dx" in the mesh. Once we find this, we can easily calculate the time step "dt" from CFL criterion as (dt*U)/dx = C Where C courant number But the problem is finding smallest division in mesh. Could anybody help in that? 

December 8, 2010, 03:55 
Hi

#8 
Member
Dynampally Pavitran
Join Date: Mar 2010
Location: India
Posts: 73
Rep Power: 7 
Usually minimum delx are found near the boundaries, eg: for airfoils you can find min delx near leading edge or trailing edge.


December 8, 2010, 04:26 

#9 
New Member
Muhammad Nafees
Join Date: Aug 2010
Posts: 2
Rep Power: 0 
Dear Pavitran.....
thanks for enlightening us on the subject...... your explanation will hopefully help... 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Time step size and max iterations per time step  pUl  FLUENT  27  May 12, 2015 20:04 
How to write k and epsilon before the abnormal end  xiuying  OpenFOAM Running, Solving & CFD  8  August 27, 2013 15:33 
SLTS+rhoPisoFoam: what is rDeltaT???  nileshjrane  OpenFOAM Running, Solving & CFD  4  February 25, 2013 05:13 
calculation diverge after continue to run  zhajingjing  OpenFOAM  0  April 28, 2010 04:35 
Is it to possible to change time step DeltaT during a calculation  dbxmcf  OpenFOAM Running, Solving & CFD  4  November 4, 2007 15:21 