|
[Sponsors] |
October 4, 2006, 09:16 |
Getting solution values per iteration
|
#1 |
Guest
Posts: n/a
|
Hi,
I want to use the value of CL and/or Cd as criteria for convergence. I'm thinking of calculating something like (CL(i)-CL(i-1))/CL(i-1), where i is the current iteration and i-1 is the previous. How can i get the values of variables from the previous and current iteration and do such a calculation? Thanks Joćo |
|
October 4, 2006, 09:52 |
Re: Getting solution values per iteration
|
#2 |
Guest
Posts: n/a
|
AFAIK you will have to use MMS to create user variables on the stack user /USER_DATA to store the previous iterations values and then use a user CEL function to read them e.g.:
(CL(i)-CL(i-1))/CL(i-1) CL(i) will be a standard CEL expression defined in CFXpre. CL(i-1) will be a call to a user CEL routine e.g. getPreviousCL(args) You had better have hairy balls because MMS isnt for the feint of heart. |
|
October 4, 2006, 10:01 |
Re: Getting solution values per iteration
|
#3 |
Guest
Posts: n/a
|
Joe,
Thanks. I was gessing it wouldn't be easy...Just what is MMS, so that I can get this started? Joao |
|
October 4, 2006, 11:24 |
Re: Getting solution values per iteration
|
#4 |
Guest
Posts: n/a
|
It is not recommended to converge based on a global parameter like this. It is entirely possible that the lift coefficient could stabilize for a while then change. For instance, the boundary layer may continue to develop slowly and have little effect on lift, then separate, which will have a dramatic effect on lift. In such a case, your convergence criteria would completely miss the separation.
The best criteria for convergence are the MAX and RMS residuals. It would be better to determine what MAX or RMS residual level is needed to get your lift coefficient within a required tolerance and stick with this value. Keep in mind that simply observing that a value doesn't change does not mean you have a solution. Regards, Robin |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Transient, initial variables from a previous solution | nakor | FloEFD, FloWorks & FloTHERM | 0 | April 22, 2011 05:34 |
[snappyHexMesh] snappyHexMesh aborting | Tobi | OpenFOAM Meshing & Mesh Conversion | 0 | November 10, 2010 04:23 |
strange node values @ solid/fluid interface - help | JB | FLUENT | 2 | November 1, 2008 13:04 |
IcoFoam parallel woes | msrinath80 | OpenFOAM Running, Solving & CFD | 9 | July 22, 2007 03:58 |
Wall functions | Abhijit Tilak | Main CFD Forum | 6 | February 5, 1999 02:16 |