# momentum equation--pressure correction with SIMPLE

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

 October 24, 2009, 22:31 momentum equation--pressure correction with SIMPLE #1 New Member   anonymous Join Date: Oct 2009 Posts: 2 Rep Power: 0 Hi, I am very new to CFD, I have been unable to find a resource that shows how the discretised momentum equation is derived. Every reference starts with the differential momentum equation, then the next line is what seems to be a random equation with a_i_J multiplying the first term, with no explanation of the meaning of this coefficient. The derivations of the pressure correction equation is similar. There is a d in it, which is related to this a_i_J term in some way that I cannot understand. Can anyone point me to a web based reference, or provide an explaination. Thanks

 October 25, 2009, 05:40 #2 New Member   mhm Join Date: Oct 2009 Posts: 8 Rep Power: 7 Hi fredp I suggest you to study chapters 4, 5 and 6 of “An Introduction to Computational Fluid Dynamics - Versteeg”. It is very good.

October 25, 2009, 10:52
#3
New Member

anonymous
Join Date: Oct 2009
Posts: 2
Rep Power: 0
Quote:
 Originally Posted by mhm Hi fredp I suggest you to study chapters 4, 5 and 6 of “An Introduction to Computational Fluid Dynamics - Versteeg”. It is very good.
I have read this--it wasn't helpful. It has a bunch of undefined coefficients. I would love to find an explanation of the meaning of each of the terms in the momentum equation in words.

 October 30, 2009, 10:14 #4 New Member   zoheb Join Date: Mar 2009 Location: india, mumbai Posts: 24 Rep Power: 8 refer to introduction to computational fluid methods by schafer ...its good....neway i found limited info on the staggered grid....try reading anil dates papers on pressure velocity corrections....here has well enuf explained the terms and their relevance.i used the collocated method ...i find that more convinient...all the best

 October 30, 2009, 14:29 Simple algorithm #5 Member   jk Join Date: Jun 2009 Posts: 64 Rep Power: 8 Dear friends, I have written my SIMPLE algorithm on staggered grid. But encountering the following problem : before that my boundary conditions are as follows: Staggered arrangment is as follows: left and right side (basically a duct) ---> pressure inlet and outlet top and bottom ---> pressure and u velocity Left side ---> pressure inlet (32 pa) right side ---> pressure outlet (0 pa) Top and bottom side ---> dp/dy = 0, dpprime/dy = 0, u=0. URF for pprime (or pressure) = 0.1 and for u and v equal to 0.3 and 0.4 respectively. Problems: (i) pprime values are very very small (so pressure correction also would be very small) (ii) only at the exit my u velocities are reasonable I have run my outer iterations = 3000 and my inner iteration (pprime) = 300. Still no change in the pressure value from the inlet to exit. Please help me in this regard. thanks jyo

 October 31, 2009, 05:53 #6 New Member   zoheb Join Date: Mar 2009 Location: india, mumbai Posts: 24 Rep Power: 8 check with your urf ....it might be a bit low....

 October 31, 2009, 09:46 initialisation problem #7 Member   jk Join Date: Jun 2009 Posts: 64 Rep Power: 8 Hi zobekenobe, Thanks for your mail. I tried with a simple couette problem which is given in the CFD book by Anderson. But the formulation is from Malalasekara. I am getting a linear profile as given in the book. But when i change the initial value (it is given as zero in the book especially pressure) even for a small change my solution is blowing. Really don't know what is happening. My urf's are 0.1 for pressure and 0.3,0.4 for u and v velocity. ( I am using staggered grid and explicit method). I tried with both upwind and central difference. There is no improvement among the two. thanks jyothish

 October 31, 2009, 09:58 #8 Member   jk Join Date: Jun 2009 Posts: 64 Rep Power: 8 Hi, Previous to my last message, I want to convey that effect of dp (between inlet and outlet) is reflected only at the outlet and not in the interior of the domain. thanks jyo

 November 4, 2009, 09:22 #9 Member   jk Join Date: Jun 2009 Posts: 64 Rep Power: 8 Hi, Still my SIMPLE formulation is not giving the correct solution. i am using staggered grid arrangment. I am giving pressure inlet and pressure outlet boundary conditions. Here the pressure inlet is the total pressure. Top and bottom i am giving wall boundary conditions. From Po i am calculating Ps (static pressure) using Po = Ps + 0.5*row*u^2 and this Ps will go into the momentum equations (pressure gradient). Is this right ?. My solution is blowing after some thounsand iterations. I am taking the dimension from anderson book (for couette flow). It is as follows : domain size ----> 0.003m (height) and 0.15 m (length) grid size ----> 21 x 11 for pressure ----> 22 x 11 for u velocity ------> 23 x 12 for v velocity I am using explicit scheme (finite volume) with upwind discretization. Please help me in this regard. Thanks jyothishkumar

 November 12, 2009, 08:29 SIMPLE algorithm in colocated grids #10 Member   jk Join Date: Jun 2009 Posts: 64 Rep Power: 8 Hi, Can you please check the following steps to follow to implement SIMPLE algorithm in colocated grids. Also i have added my questions (if there is any) immediately below the steps. please help me in this regard. 1. First initializing values at the nodes (p, control volume center) Up*, Vp* Pp* etc 2. Finding Up and Vp from the general discretized momentum equation from above guessed value. 3. Find Ue at the cell faces (rhie and chow’s intrapolation ) using the adjacent nodal U values For eg. At face ‘e’ à Ue = f(Up, UE, PE, Pp, PEE, dp, dE, etc) 4. For the pressure correction equation we need Ue ,Up, de(?), dp We can get the value of dp from the discretized momentum equation (CV center) How to get this de (should also be through interpolation like de = (dP + dE) / 2 ?) 5. From the above pressure correction formulation we can get the value of Pprime. 6. Then PP = PP* + (alpha * Pprime) 7. Velocity correction Up = Up* + dP(Pprime(i-1,j)-Pprime(I,j)) -------à without correction Up = (alpha * Up (without correction)) + (1- alpha)*Up old (from previous iteration) 8. Above process is repeated for VP. thanks jyo

 November 16, 2009, 10:51 #11 Member   jk Join Date: Jun 2009 Posts: 64 Rep Power: 8 Hi, Can you pls check my following steps for SIMPLE Algorithm in colocated grids. 1. Find Up's at all points from discretized momentum equation (without pressure term) 2. Apply rho and chi's interpolation for finding the face velocity Ue(interpolated from UP and UE (nodal points or cell centered) 3. substitute this face velocities in the correction equation. Pls help me in this regard. thanks jyothish

 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 kk Main CFD Forum 0 July 2, 2006 23:13 kk Main CFD Forum 0 July 2, 2006 23:10 Bharath Somayaji Main CFD Forum 1 March 1, 2006 07:12 George Main CFD Forum 5 March 22, 2005 07:10 Abhijit Tilak Main CFD Forum 9 January 8, 2001 12:04

All times are GMT -4. The time now is 05:36.