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

Best approach to study temporal and spatial convergence

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 15, 2013, 15:20
Default Best approach to study temporal and spatial convergence
  #1
Member
 
navid
Join Date: Jan 2010
Posts: 89
Rep Power: 7
ndabir is on a distinguished road
Hi all,

I am doing an unsteady CFD simulation and I am totally confused of what approach to take to study both temporal and spatial convergence. here is the options that I think might be helpful:

1) First refine the grid with a specific time step until the spatial error goes to zero. Then take that grid and use different time steps to study temporal convergence.

2) Other approach can be the opposite of the above. Which means, first study the time step effect on a specific grid and then after reducing time step error to almost zero, then fix the refined time step and reduce the mesh size.

Which one is better? Is there any better approach?
My simulations take relatively long time to give results (2 days) so I need a good approach to avoid long simulations.
ndabir is offline   Reply With Quote

Old   July 15, 2013, 18:01
Default
  #2
Member
 
Reza
Join Date: Mar 2009
Location: Appleton, WI
Posts: 99
Rep Power: 8
triple_r is on a distinguished road
I don't have much experience in this kind of problems, but here are my two cents:

Time step might be dependent on the spacial resolution, specially if you are using explicit methods, so the second method is not something that I would try. When you increase the grid resolution while keeping the time step the same, you are essentially increasing the CFL number, and that might cause instabilities in the code.

Also, try using an adaptive time-stepping method: select a time step and march one step, then use two steps two march the same time and compare your result of interest, if the difference is smaller than a threshold, use the larger step for the next time step. If the difference is too large, use four steps (or three) and repeat this until you get a difference small enough for your case.

This might seem very time consuming, but if you are planning on running the whole simulation with different time steps to check the time-step independency, then you are already solving with the smaller and larger steps, this might even help you save solution time.
triple_r is offline   Reply With Quote

Old   July 15, 2013, 19:25
Default
  #3
Member
 
navid
Join Date: Jan 2010
Posts: 89
Rep Power: 7
ndabir is on a distinguished road
Thanks for your response. My simulation is a bubble collapse near a wall. The bubble collapses because of the pressure difference between inside bubble (low pressure) and outside bubble inside water higher pressure. I use VOF (Volume of Fluid) to capture the interface. I use SIMPLE algorithm in FLUENT which is basically an implicit solver. Previously I used explicit method to only discretize VOF equation but because of stability issues, I changed it to implicit. Now I am not facing any divergence problems but temporal and spatial convergence is an issue.
The problem with my simulation is the process is very slow at the beginning of simulation but most of the important physics happen only in the last 20% time of the simulation. So I definitely need variable time step during one simulation. This makes it harder to do the temporal convergence.
ndabir is offline   Reply With Quote

Old   July 16, 2013, 03:47
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,640
Rep Power: 23
FMDenaro will become famous soon enough
First I strongly suggest to chech for a simple analytical solution!

Then:

Spatial accuracy: use a fixed dt, taken as small as possible
Temporal accuracy: use a fixes space step taken as small as possible

Remember that the local truncation error is a function of L = L(dt,dx,dy,dz) therefore to chech for a single discretization parameter you must ensure that the other do not enter into the convergence (you see that when the convergence is not reached)

In general the accuracy analysis is performed by fixing dt/h = constant
FMDenaro is offline   Reply With Quote

Old   July 16, 2013, 09:43
Default
  #5
Member
 
navid
Join Date: Jan 2010
Posts: 89
Rep Power: 7
ndabir is on a distinguished road
how to determine the value of dt/h=constant? I mean what should the value of this "constant" be?
ndabir is offline   Reply With Quote

Old   July 16, 2013, 10:32
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,640
Rep Power: 23
FMDenaro will become famous soon enough
the value is constrained by the stability requirement of the scheme
FMDenaro is offline   Reply With Quote

Old   July 16, 2013, 10:46
Default
  #7
Member
 
Reza
Join Date: Mar 2009
Location: Appleton, WI
Posts: 99
Rep Power: 8
triple_r is on a distinguished road
And if the time-advancing scheme is unconditionally stable, then the time step size comes from the fastest physical phenomenon happening during the simulation. If this fastest event has a time scale of Dt, then your time step, dt, should be smaller to make sure you are capturing that event.

BTW, I think I remember you being able to specify a CFL number in FLUENT instead of specifying a constant time step. You might be better off with that since, as you said, in the early stages it is quite calm so time step can be large for a given CFL, but when the bubble starts imploding, things move faster and for CFL to remain constant the software has to bring the time step down. Then, for time-step independency study, you can run your simulation for different upper limits for CFL. I hope it makes sense.
triple_r is offline   Reply With Quote

Reply

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
Temporal and Spatial averages shahzad_munir Main CFD Forum 0 September 18, 2012 06:09
Comparison of different spatial and temporal discretization schemes harly OpenFOAM Running, Solving & CFD 2 September 24, 2011 12:15
temporal and spatial discretization. pity Main CFD Forum 0 June 10, 2007 22:53
Spatial and temporal discretization Giuseppe FLUENT 2 January 23, 2006 06:18
1st order temporal & 2nd order spatial Prateep Chatterjee FLUENT 0 January 19, 2003 01:31


All times are GMT -4. The time now is 20:35.