calculation of time step

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

 December 6, 2010, 02:09 calculation of time step #1 Member   Aamer Shahzad Join Date: Mar 2010 Posts: 58 Rep Power: 8 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: 7 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 lth and rajann_786 like this.

 December 6, 2010, 08:53 #3 Member   Aamer Shahzad Join Date: Mar 2010 Posts: 58 Rep Power: 8 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: 7 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 Bio-Hydrodymanics: Trends, Chalenges and Recent Advances,IEEE JOURNAL OF OCIANIC ENGINEERING, 29 , 595-604, 2004. o_mars_2010 likes this.

 December 7, 2010, 05:11 #5 New Member   Join Date: Aug 2010 Posts: 4 Rep Power: 8 i am sorry but this explanation ist also so difficult. rajann_786 likes this. __________________ Mechanical Man

 December 7, 2010, 07:06 Hi #6 Member   Dynampally Pavitran Join Date: Mar 2010 Location: India Posts: 74 Rep Power: 8 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 time-step 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 1E-4), 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 time-step size & cell-size. 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) lth and rajann_786 like this.

 December 8, 2010, 03:27 Finding smallest division in mesh #7 New Member   Irshad Join Date: Dec 2009 Location: Bangalore Posts: 19 Rep Power: 8 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: 74 Rep Power: 8 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 Linear Mode

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

 Similar Threads Thread Thread Starter Forum Replies Last Post pUl| FLUENT 31 August 21, 2015 04:46 xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33 nileshjrane OpenFOAM Running, Solving & CFD 4 February 25, 2013 05:13 zhajingjing OpenFOAM 0 April 28, 2010 04:35 dbxmcf OpenFOAM Running, Solving & CFD 4 November 4, 2007 15:21

All times are GMT -4. The time now is 15:38.