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

SIMPLE algorithm for collocated storage

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 6, 2008, 16:05
Default SIMPLE algorithm for collocated storage
  #1
alex
Guest
 
Posts: n/a
Hi,

I'm trying to write a 2d FV solver but I can't get the SIMPLE algorithm to work. I get the basics of how it should work but it just keeps adding to the pressure at each iteration until pressure=NaN.

Do you know of a good reference that describes how to code a SIMPLE solver for a collocated, structured grid arrangement. (has to be collocated as I can get it to work with a staggered grid). Preferably written in as simple english as possible!

Thanks

Alex
  Reply With Quote

Old   April 6, 2008, 23:28
Default Re: SIMPLE algorithm for collocated storage
  #2
hadian
Guest
 
Posts: n/a
the difference between staggered and collocated is that you should have correction term for calculating fluxes on cell faces. have you read the paper of Rhie and Chow?

Rhie, C.M. and Chow, W.L. (1983). "Numerical Study of the Turbulent Flow Past an Airfoil Trailing Edge Separation", AIAA Journal 21(11), 1525-1532.
  Reply With Quote

Old   April 7, 2008, 02:20
Default Re: SIMPLE algorithm for collocated storage
  #3
sudhakar
Guest
 
Posts: n/a
hi, A very good reference for SIMPLE in colocated grid is the book 'computational fluid dynamics' by ANIL DATE. If you cant get this, you can refer to some of the papers published by the author.
  Reply With Quote

Old   April 8, 2008, 15:36
Default Re: SIMPLE algorithm for collocated storage
  #4
alex
Guest
 
Posts: n/a
Hadian, Sudhakar,

Thank you both for your help. The book you recommended by Anil Date was indeed a very good reference. Highly recommended to anyone else interested in SIMPLE for colocated grids.

There is one thing I'm a little unsure of though. In implementing the dP'/dn=0 boundary condition, (on the west boundary for example) I set the west coefficient of the pressure correction to zero. However I still need a value of AW for AP(=AW+AE+AN+AS) (where AW & AP, etc. are coefs. of the pressure correction equation). AW involves ap (ap = coef. of momentum equation) at the west face, which is obtained by linear interpolation between nodal values for internal CVS. How do I obtain ap at a boundary?

Thanks

  Reply With Quote

Old   April 8, 2008, 23:48
Default Re: SIMPLE algorithm for collocated storage
  #5
sudhakar
Guest
 
Posts: n/a
hi alex, I think you are little confused. The issue is straightforward. Near the west boundary, aW=0; substitute this in your calculation of aP, then it becomes aP=aE+aS+aN. Once you implement this, you will experience that your iterative solver takes slightly longer time to converge. hope it helps.
  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
SIMPLE algorithm in 3D cylindrical coordinates zouchu Main CFD Forum 1 January 20, 2014 18:02
The correction on pressure equation of SIMPLE algorithm in MRFSimpleFOAM solver renyun0511 OpenFOAM Running, Solving & CFD 0 November 10, 2010 02:47
About Phase Coupled SIMPLE (PC-SIMPLE) algorithm Yan Kai Main CFD Forum 0 April 18, 2007 03:48
About Phase Coupled SIMPLE (PC-SIMPLE) algorithm Yan Kai FLUENT 0 April 13, 2007 23:17
SIMPLE algorithm Jonathan Castro Main CFD Forum 3 December 10, 1999 05:59


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