
[Sponsors] 
August 9, 2010, 12:48 
SIMPLE Algorithm  Questions from a Beginner

#1 
New Member
Emre Turkoz
Join Date: Aug 2010
Posts: 12
Rep Power: 14 
Hi everybody,
I’m beginner in CFD and trying to implement SIMPLE algorithm using MATLAB. My guidebooks are “Numerical Heat Transfer and Fluid Flow” by Patankar and “An Introduction to Computational Fluid Dynamics” by Versteeg & Malalasekera. Using both books simultaneously really helped me a lot but still, I have a couple of questions. I have problems on the boundary conditions. I’m at the critical point where momentum and pressure correction equations are explained. As an explanation to the pressure correction equation, as far as I have understood, Patankar wrote that either the real values of normal velocity or the real values of pressure should be specified at boundaries. If we specify pressure values, then the correction pfield at boundaries should be zero. It’s clear so far. I wonder what happens in that case to the momentum equation’s boundary conditions. How do we specify them? Are they the ones from the previous iteration? And additionally: How do we set pressure difference values at these boundaries? Is it okey if I extend the pressure field from the previous iteration and extrapolate mirror nodes to get pressure difference and add it as a source term? Thank you all 

August 9, 2010, 13:06 

#2  
New Member
Akinola
Join Date: Jul 2010
Posts: 25
Rep Power: 14 
Quote:
0= dp/dx+mhu*d2(vx)/dx2+rho*gx (ie vx=vy=vz=0 at the wall. also the shear velocity derivative will be zero cos they are still at the wall) From the equation above, you will be able to get the pressure outside of the domain. Hope this helps and I ve not confused you too much. Take care 

August 9, 2010, 13:45 

#3 
New Member
Emre Turkoz
Join Date: Aug 2010
Posts: 12
Rep Power: 14 
falopsy, thank you very much for your answer. it really clarifies how i should deal with the pressure values at boundaries in umomentum equation.
i'm also confused at the uvalues part of that equation. are the boundary values of corrected velocity field from the previous iteration new boundary conditions for the next iteration? or should i deal with these nodes as i deal with the inner points, specifying no specific value for u? 

August 9, 2010, 15:13 

#4  
New Member
Akinola
Join Date: Jul 2010
Posts: 25
Rep Power: 14 
Quote:


August 9, 2010, 16:42 

#5 
New Member
Emre Turkoz
Join Date: Aug 2010
Posts: 12
Rep Power: 14 
thank you, again. I'm reading further on the topic. I noticed, Ch.9 of the Versteeg book actually explains boundaries.
thanks! 

August 11, 2010, 11:58 
SIMPLE Algorithm  Questions from a Beginner

#6 
New Member
V GURUPRASAD
Join Date: Jul 2010
Posts: 5
Rep Power: 14 
hi,
the initial pressure values in momentum equns shuid b assumed as zero oly..... 

August 11, 2010, 12:29 

#7 
New Member
Akinola
Join Date: Jul 2010
Posts: 25
Rep Power: 14 
Yes, the initial condition doesnt really matter but what how will you treat it after the first iteration when the pressure will then be nonuniform? You still need a way to treat the boundary condition.


August 11, 2010, 14:01 

#8 
New Member
Emre Turkoz
Join Date: Aug 2010
Posts: 12
Rep Power: 14 
I just solve pressure correction for the inner region, add it up to the pressure from the previous iteration and again impose boundary conditions to this corrected pressure. The results seem reasonable then.
I have another concern. I wrote the algorithm using finite differences with the help of the book "Computational Fluid Dynamics" by Hoffmann. I made no matrix calculations with conjugate gradient, gauss seidel or other iterative stuff. So why do we prefer other methods that includes them? Are they more accurate? Or faster? 

August 13, 2010, 00:21 
HI

#9  
New Member
V GURUPRASAD
Join Date: Jul 2010
Posts: 5
Rep Power: 14 
Quote:
hi, the guessed pr (i.e) ur initial pr. is assumed as zero.. then 4 d nxt iteration u will solve for pressure correction eqn to find new pr. correction value. Then with the pressure correction value v can find new pr value... Actually v r not finding the exact pr values.. V r jus finding the pr. gradient values in SIMPLEX method..... Then u ll compute new U & V values using corrected pr.... Then check for convergence criteria..... Incase of turbulence v shud multiply the under relaxation factor(usually 0.2 to 0.8) with the pr. correction term to reduce the pr. correction value, thus geting soln converged... 

Tags 
matlab, patankar, simple 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
questions about SIMPLE algorithm (by Patankar)  Nutty  Main CFD Forum  5  January 22, 2015 12:22 
SIMPLE algorithm in 3D cylindrical coordinates  zouchu  Main CFD Forum  1  January 20, 2014 17:02 
SIMPLE algorithm  jika  Main CFD Forum  2  October 22, 2009 06:46 
About Phase Coupled SIMPLE (PCSIMPLE) algorithm  Yan Kai  FLUENT  0  April 13, 2007 23:17 
SIMPLE algorithm  Jonathan Castro  Main CFD Forum  3  December 10, 1999 04:59 