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

Time step size????

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 23, 2009, 00:02
Default Time step size????
  #1
Member
 
MrFluent
Join Date: Mar 2009
Posts: 33
Rep Power: 8
mr_fluent is on a distinguished road
I am wondering about one question.
Lets say i have a cartesian mesh domain, whose smallest cell size is 1E-4meters.

And if i have maximum velocity in domain aroun say 100m/s.

And in momentum equations if i am treating convection terms as explicit but diffusion terms fully implicit. How do i decide what maximum time step size i could use.

Is there any pointers where i could find a definitive formula to use in this case.
mr_fluent is offline   Reply With Quote

Old   March 23, 2009, 00:25
Default
  #2
New Member
 
Prakash Ayappan
Join Date: Mar 2009
Posts: 25
Rep Power: 8
technocrat.prakash is on a distinguished road
Quote:
Originally Posted by mr_fluent View Post
I am wondering about one question.
Lets say i have a cartesian mesh domain, whose smallest cell size is 1E-4meters.

And if i have maximum velocity in domain aroun say 100m/s.

And in momentum equations if i am treating convection terms as explicit but diffusion terms fully implicit. How do i decide what maximum time step size i could use.

Is there any pointers where i could find a definitive formula to use in this case.
For that you have to integrate the momentum eqn. And make sure the constant co-efficient term to be positve. A rough estimate will be (for 2D case energy eqn)
time step = (del X)^2/(Alpha). where alpha is the diffusivity.
technocrat.prakash is offline   Reply With Quote

Old   March 23, 2009, 00:51
Default
  #3
Member
 
MrFluent
Join Date: Mar 2009
Posts: 33
Rep Power: 8
mr_fluent is on a distinguished road
Quote:
Originally Posted by technocrat.prakash View Post
For that you have to integrate the momentum eqn. And make sure the constant co-efficient term to be positve. A rough estimate will be (for 2D case energy eqn)
time step = (del X)^2/(Alpha). where alpha is the diffusivity.

Thanks for replying. The reason for my question was that is there any other criteria that also should be taken care of.

My momentum term diagonal matrix is all positive for two reasons. By treating convection term explicit, i only have diffusion terms in matrix, they produce positive diagonal plus its a transient case of there is contribution from 1/delt_t too.
mr_fluent is offline   Reply With Quote

Old   March 23, 2009, 03:55
Default
  #4
Senior Member
 
ztdep's Avatar
 
p ding
Join Date: Mar 2009
Posts: 287
Rep Power: 9
ztdep is on a distinguished road
Send a message via Yahoo to ztdep Send a message via Skype™ to ztdep
i prefer a full implicit method, so we do not have the problem about time step size
ztdep is offline   Reply With Quote

Old   March 23, 2009, 04:05
Default Courant criteria
  #5
New Member
 
Join Date: Mar 2009
Posts: 4
Rep Power: 8
Hochola is on a distinguished road
Do you have to consider the Courant criteria for this case?
Hochola is offline   Reply With Quote

Old   March 23, 2009, 04:05
Default
  #6
Member
 
MrFluent
Join Date: Mar 2009
Posts: 33
Rep Power: 8
mr_fluent is on a distinguished road
Quote:
Originally Posted by ztdep View Post
i prefer a full implicit method, so we do not have the problem about time step size

Actually i also prefer but since i have to use a very small time step, irrespective of method i use, i was thinking if i could ditch convection term then i am all left with diffusional terms. Now if i only take diffusional terms, my matrix would be more of less similar for u,v,w equations. If this is the case i could reduce operation count greatly.
So whole idea was if i am forced to use small step could i make some shortcut.



But it seems i might have to go with implicit, because in my case del_x/umax is too small, (smaller than time step i have to use).

thanks for replying though.
mr_fluent is offline   Reply With Quote

Old   March 23, 2009, 06:06
Default
  #7
Member
 
MrFluent
Join Date: Mar 2009
Posts: 33
Rep Power: 8
mr_fluent is on a distinguished road
thanks guys, i finally ended up making it fully implicit and in almost similar calculation efficiency if it were only diffusional matrix.
I checked it and its working fine. Thanks again.
mr_fluent is offline   Reply With Quote

Old   March 23, 2009, 06:25
Default
  #8
New Member
 
Prakash Ayappan
Join Date: Mar 2009
Posts: 25
Rep Power: 8
technocrat.prakash is on a distinguished road
Good. Keep Going.
technocrat.prakash is offline   Reply With Quote

Old   March 23, 2009, 08:20
Default
  #9
Senior Member
 
ztdep's Avatar
 
p ding
Join Date: Mar 2009
Posts: 287
Rep Power: 9
ztdep is on a distinguished road
Send a message via Yahoo to ztdep Send a message via Skype™ to ztdep
Quote:
Originally Posted by mr_fluent View Post
Actually i also prefer but since i have to use a very small time step, irrespective of method i use, i was thinking if i could ditch convection term then i am all left with diffusional terms. Now if i only take diffusional terms, my matrix would be more of less similar for u,v,w equations. If this is the case i could reduce operation count greatly.
So whole idea was if i am forced to use small step could i make some shortcut.



But it seems i might have to go with implicit, because in my case del_x/umax is too small, (smaller than time step i have to use).

thanks for replying though.
i do not think you can reduce the time greatly, since the explicit method has the limitation in the time step.
ztdep is offline   Reply With Quote

Old   March 23, 2009, 08:48
Default
  #10
Member
 
MrFluent
Join Date: Mar 2009
Posts: 33
Rep Power: 8
mr_fluent is on a distinguished road
Quote:
Originally Posted by ztdep View Post
i do not think you can reduce the time greatly, since the explicit method has the limitation in the time step.
Actually i am already supposed to use time step sizes smaller than 1e-4 or so.

Anyway i was trying to implement approaximate factorization as kim and moin used in their paper (1984 dns paper). For that i was supposed to solve them by tdma. But since i am solving by tdma, keeping convection terms in solver matrix does not really add anything to calculation time. So i did that. And now it is fully implicit and working well.
mr_fluent is offline   Reply With Quote

Reply

Tags
courant number, time step size

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
Computational time sunnysun OpenFOAM Running, Solving & CFD 5 March 16, 2009 04:32
Modeling in micron scale using icoFoam m9819348 OpenFOAM Running, Solving & CFD 7 October 27, 2007 00:36
Transient simulation not converging skabilan OpenFOAM Running, Solving & CFD 12 September 17, 2007 17:48
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 10:46.