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

Getting solution values per iteration

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 4, 2006, 08:16
Default Getting solution values per iteration
  #1
Joao Bettencourt
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
  Reply With Quote

Old   October 4, 2006, 08:52
Default Re: Getting solution values per iteration
  #2
Joe
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.
  Reply With Quote

Old   October 4, 2006, 09:01
Default Re: Getting solution values per iteration
  #3
Joao Bettencourt
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

  Reply With Quote

Old   October 4, 2006, 10:24
Default Re: Getting solution values per iteration
  #4
Robin
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
  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
Transient, initial variables from a previous solution nakor FloEFD, FloWorks & FloTHERM 0 April 22, 2011 04:34
snappyHexMesh aborting Tobi OpenFOAM Native Meshers: snappyHexMesh and Others 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 02:58
Wall functions Abhijit Tilak Main CFD Forum 6 February 5, 1999 02:16


All times are GMT -4. The time now is 17:52.