|
[Sponsors] |
when we can say an unsteady problem has converged? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 25, 2013, 10:16 |
when we can say an unsteady problem has converged?
|
#1 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
how much should Ux,Uy,global or local continuity residual reach that we can say convergency has occured?
I have attached my solution's figure. any opinion is appreciated. whats the condition of convergency in unsteady compressible problems? thanks. |
|
January 26, 2013, 09:03 |
|
#2 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
is it true that at each time step at a 2-D unsteady problem initial residual of all of p,Ux,Uy should be less than 10^-3?than final residuals are not important?
|
|
January 26, 2013, 10:03 |
|
#3 |
Senior Member
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 200
Rep Power: 18 |
Hi Ehsan,
A reduction of the residuals by three to five orders of magnitude is recommended in book of Ferziguer and Peric. I suggest that you make some test varying the convergence condition (tolerances) and evaluate how they impact the final results. Best regards, Paulo |
|
January 26, 2013, 13:49 |
|
#4 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
hi dear paulo
thanks.you mean this reduction should occur in initial residauls? I don't grasp convergency subject obviously in unsteady problems such that in steady flows is clear .shall not each time reach to convergency? Can't we suppose at general each time iteration at a unsteady problem like a small duration steady problem alone?then why we compare time steps with eachother.i'm totatally mixed up with this comparison.please tell me some information. Thanks. |
|
January 27, 2013, 06:59 |
|
#5 |
Senior Member
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 200
Rep Power: 18 |
Hi Ehsan,
Sometimes we do a transient analysis but we are only interested in the final steady state results. In this case the transient solution is used because it improves convergence or because the transient algorithm is faster than a steady state algorithm. There are other situations in which we really need to know how the flow evolutes with time, I suppose this is your case. Then you have to be sure that the convergence condition doesn’t cause errors that add up each time step causing a big error in the solution. So, what I was suggesting is that you run a simulation, and then run another dividing all tolerances by ten. Then you compare the results (Ux,Uy etc) at some important time instants, for instance t=0.1, t=0.2 etc. The tolerances must not affect the results, if they do, you must lower them until they don’t affect the results anymore. Another important variable is the time step, if the time step is not small enough, it may cause errors that affect the solution. It would be also interesting to do some comparisons, halving the time step and comparing the results to evaluate it’s influence. Best regards, Paulo |
|
January 27, 2013, 13:28 |
|
#6 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
thank you dear Paulo
Do you mean residuals(thats easier tracking by pyFoam) by results or p,Ux,... Themselves? I've divided tolerances by ten as you said to me.now residuals of variables are same(trend and value) exactly.but global countinuity is lower at second one from -9 order to -10.has the problem converged in your opinion and your experiment? Thank you. |
|
January 27, 2013, 21:16 |
|
#7 | |
Member
Eric Robertson
Join Date: Jul 2012
Posts: 95
Rep Power: 15 |
Personally, I prefer:
Ux, Ux, Uz = 1e-5 or 1e-6 if possible p = 1e-3 or 1e-4 (it is usually quite difficult to get it lower than this) Turbulent variables (k, omega, epsilon etc) 1e-5 or 1e-6 I think you're also asking where this is seen, so for example, you would look for your initial residuals as such: Quote:
You may also track convergence if you see your residuals don't get any smaller or larger over a large period of time, though it is best to try and get the actual residual figure down as low is possible. |
||
January 28, 2013, 15:21 |
|
#8 |
Senior Member
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 200
Rep Power: 18 |
Hi Ehsan,
About your question, I was suggesting that you compare the results of p,Ux, Uy etc, not the residuals. |
|
January 28, 2013, 17:22 |
|
#9 | |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
Quote:
thank you again.i should compare results in folders?or another easier way is possible? I use adjustableRunTime with automatic Co.is it better to change it to constant time step for more stability? Last edited by immortality; January 29, 2013 at 07:22. |
||
January 29, 2013, 07:19 |
|
#10 |
Senior Member
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 200
Rep Power: 18 |
Hi Ehsan,
To compare all results would be very difficult. I suggest that you choose some important locations in the flow and compare the results in the same locations. You could also draw some streamlines with paraview and compare their position using one tolerance or another. About the time step I think it would be better to use a constant time step for making comparisons, but if the results are similar the time steps would be equal in both simulations. I didn't understand what you mean by more instability, normally we want to avoid instabilities. |
|
January 29, 2013, 09:14 |
|
#11 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
thanks.sorry for the mistake!I corrected it.so I don't know how to realize a number in the column list at each folder is related to which point.should i select randomly a number for instance 20th number at the list?
|
|
January 30, 2013, 08:01 |
|
#12 |
Senior Member
Paulo Vatavuk
Join Date: Mar 2009
Location: Campinas, Brasil
Posts: 200
Rep Power: 18 |
Hi Ehsan,
If it's difficult to choose a few points, you could draw some velocity profiles using paraview and make comparisons between the profiles obtained with different tolerances. |
|
February 2, 2013, 04:21 |
|
#13 |
Member
Anant Diwakar
Join Date: Jan 2013
Posts: 68
Rep Power: 13 |
Hi
I think for unsteady flow, the best way of monitoring convergence is to monitor the variation of time-averaged value of some parameter (eg. Cl_avg, Cd_avg etc.). Regards Anant |
|
February 2, 2013, 04:23 |
|
#14 |
Member
Anant Diwakar
Join Date: Jan 2013
Posts: 68
Rep Power: 13 |
Hi,
I think the best way of monitoring convergence in an unsteady flow is to monitor the variation of time-averaged quantity of some parameter, like Cl_avg, Cd_avg etc. Regards Anant |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[GAMBIT] HELP!I have a problem when starting the unsteady iteration | usst139 | ANSYS Meshing & Geometry | 3 | January 21, 2010 05:04 |
SOS: About VOF-Dynamic mesh unsteady problem | whpsoft | FLUENT | 4 | May 9, 2006 20:26 |
Newbie Unsteady Problem | Anthony | FLUENT | 7 | March 18, 2006 13:39 |
Steady problem and unsteady mode | mini | FLUENT | 6 | April 3, 2005 21:28 |
Question concerning about unsteady problem.... | ghlee | Main CFD Forum | 6 | October 2, 1998 06:36 |