# SIMPLER Algorithm

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

 July 5, 2008, 22:44 SIMPLER Algorithm #1 David Guest   Posts: n/a Hello. I have some questions about the implementation of the SIMPLER algorithm outlined in Patankar's book, Numerical Heat Transfer and Fluid Flow. I am implementing a one dimensional version with a uniform grid. ||-o-|-o-|-o-|| o shows where the pressure is calculated, | the velocity and || represents the given boundary velocities. 1) Do I use TDMA to calculate the pseudo velocities (u^)? Or do I just iterate through the array, calculating the coefficients (ap, aw, ae, b) and the pseudo velocity as I go? 2) When calculating the pressure field from the pseudo velocity values,the coefficients aE and aW include ae and aw. aE = density * (Area / ae) *Dy *Dz But since the grid is staggered ae does not exist at the point P, where the pressure is calculated. Am I meant to use ap from the point ue (the neighbouring velocity)? 3) In the grid setup shown, there is no boundary pressure value. When I start the TDMA process for pressure, there is no aW coefficient (ie c1 = 0 in the TDMA notation used in the book) and P1 is consequently equal to 1. This means that all subsequent values of Pi are also 1. Is this correct? It doesn't seem like it can be as on the far boundary aE does not exist and QN involves a division by zero. 4) Since I am assuming an incompressible fluid, the density never changes and so terms like (d0 - d) become zero. Is this correct? Thankyou to anyone who can help me with this.

 July 13, 2008, 12:18 Re: SIMPLER Algorithm #2 ztdep Guest   Posts: n/a pressure is calculated at the main point. the boundary of pressure is the second kind. i can help you,please further contack. ztdep@yahoo.com.cn

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Bharath Main CFD Forum 2 July 10, 2009 10:53 Mori Main CFD Forum 1 August 24, 2006 18:45 elisabet Main CFD Forum 0 December 7, 2005 05:54 Erik Main CFD Forum 1 May 23, 2004 03:57 H.Setthapong Main CFD Forum 5 October 18, 2000 06:34

All times are GMT -4. The time now is 13:11.