
[Sponsors] 
March 2, 2013, 12:30 
Time step size & convergence absolute criteria

#1 
New Member
yui tsang
Join Date: Mar 2013
Posts: 1
Rep Power: 0 
Hi,
I am a student currently working on a simulation of a confined room with a steady inlet using fluent. My model is struggling to get converged (10^4) with a steady approach and therefore I am trying to use a transient approach instead. My question is how I should determine the time step size and convergence absolute criteria. I have read from previous relevant post in which time step size was said to be equal to delta x/velocity, where delta x is size of the smallest cell. however, I wonder how I could determine the velocity for that? On top of my question about the convergence absolute criteria, I also noticed that residuals is not the only indicator for convergence. It was suggested in another post we could check if the solution is converged by implementing some control/monitor points. I wonder what this means explicitly. Background of case: Room with demensions of 4.6x3.36x3m with Displacement Ventilation system. Inlet and outlet are defined as massflow inlet and pressureoutlet. the only heat source in the room is human  defined as heat flux of 160W/m2 on the head Thank you very much for your time on my post! Matthew Last edited by yuitsang; March 2, 2013 at 14:02. 

March 4, 2013, 04:42 

#2 
Senior Member
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,280
Rep Power: 18 
Good questions!
1. Your "nonconverged" solution, how do the results look like? Completely wrong or not? If not, you could use the velocity you get there to estimate your velocity for the transient case. The "dx / dv" comes from the Courant criterion, or CourantFriedrichLewy = CFL criterion (http://en.wikipedia.org/wiki/Courant...Lewy_condition). This is the maximum time step you can use, when you use a explicit time integration. It means, that a fluid particle does not move more than one cell during a single time step and is needed for stability in explicit time integration. As I understand it, it is not relevant for implicit timeintegration, which you should use anyway most of the times. Main advantage of explicit time integration is, that it is easier to implement. Since you bought Fluent, this is fortunately not relevant to you... Now, what I suggest is to take one characteristical lenght scale "L" of your system (such as the room size) and have  let's say  20 time steps for the fluid to travel through the domain, so that dt = L / (20*v). You can increase that number (20) later, if you don't get the results you are looking for or have numerical problems. 2. Small residuals can however mean, that your solution is still not final, thus it changes from iteration to iteration (but slowly). Sometimes you can monitor some characteristical value inside you domain, such as the (area integrated) heatflux through the outlet, or some force on a body and plot it each iteration. At best, this is your value of interest and some (area or volume) integral value. In Fluent you can go to "Monitors > Surface Monitors" or "Volume Monitors" and create one. If this value stops changing when iterating, your solution also converges. This is a good additional (to the residuals) tool. 3. Can you post a picture of you setup and of your mesh? You should be able to get convergency, also for the steady state...
__________________
The skeleton ran out of shampoo in the shower. 

March 6, 2013, 11:02 

#3 
Senior Member
OJ
Join Date: Apr 2012
Location: United Kindom
Posts: 475
Rep Power: 11 
You may need to do the timestep study to understand what timestep is just enough to bring down your residuals to the convergence criteria you are happy with (say 1e5 etc). Or, easier way is to use adaptive timestepping to let FLUENT determine its ideal timestep.
This way we ensure that every timestep is converging. Regarding the convergence criteria, for transient analysis there are different scenarios. If the case can finally come to steady, you should see the monitors straight. If the case is periodic, you will keep seeing the periodic oscillations. In this case you will need to monitor the periodicity in the monitors, and then stop when have simulated enough long. If it is chaotic, the monitor oscillations can be indeterministic and you may need to monitor certain volumetric parameters to make sure the simulation has run long enough. For your case, you may also consider plotting volumetric mass imbalance and total volumetric heat generated in the system as well. There may be a point where total heat in the room will stay constant with increase in time, this may be yoru convergence. OJ 

March 21, 2013, 01:32 
transient convergence enquiry

#4 
Senior Member
Join Date: Nov 2009
Posts: 125
Rep Power: 8 
Dear all,
with reference to the attachment, what are the oscillations mean please? i'm performing a fluidthermal coupling 3D transient analysis for a total time of 2400s. as seen on the attachment, less than 100s, i'm getting oscillations in the TEnergy RMS residual plot on all components. on the fluid side, the turbulence RMS residuals and others have reached 1e4 criteria, but not for heat transfer, and this calculation is still ongoing. Please help! Thanks.
__________________
Thank you for your kind attention. Kind regards, mactech001 Currently using: ANSYS v13 

March 21, 2013, 05:53 

#5 
Senior Member
OJ
Join Date: Apr 2012
Location: United Kindom
Posts: 475
Rep Power: 11 
This is a FLUENT forum and from your image, it seems that you are using CFX. If may be appropriate to post this question on CFX forum.
http://www.cfdonline.com/Forums/cfx/ Have you read all the documentation on CFD FAQ here: http://www.cfdonline.com/Wiki/Ansys_FAQ The residuals don't seem right. In transient simulation, you should make sure that your imbalances are reduced to within 1% for every timestep. I can't see the plot of imbalances in your image. OJ 

April 15, 2013, 04:27 

#6 
Senior Member
Join Date: Nov 2009
Posts: 125
Rep Power: 8 
Dear OJ,
thanks for the advice. what if the imbalances are not within 1% between every timestep? what does that indicate i should check please?
__________________
Thank you for your kind attention. Kind regards, mactech001 Currently using: ANSYS v13 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Superlinear speedup in OpenFOAM 13  msrinath80  OpenFOAM Running, Solving & CFD  18  March 3, 2015 06:36 
How to write k and epsilon before the abnormal end  xiuying  OpenFOAM Running, Solving & CFD  8  August 27, 2013 15:33 
OF 1.6  Ubuntu 9.10 (64bit)  GLIBCXX_3.4.11 not found  piprus  OpenFOAM Installation  22  February 25, 2010 14:43 
NavierStokes time step size  Martin  Main CFD Forum  2  June 6, 2008 03:38 
Phase locked average in run time  panara  OpenFOAM  2  February 20, 2008 15:37 