
[Sponsors] 
March 22, 2005, 04:21 
pressure poisson equation on nonstaggered grids

#1 
Guest
Posts: n/a

Hi everyone,
I am a newbie in pressurecorrection method (solution of Poisson type equation for pressure) and trying to develop a new CFD code in cylindrical coordinates and on nonstaggered grids based on this method. I had coded the algorithm, everything looks fine to me but there must be something wrong since the code explodes. Most of the books I come up with discusses the issue in cartesian coordinated and on staggered grids due to wellknown reasons. My question is, can anyone suggest any sources especially discussing the issue in terms of cylindrical coordinates(books etc.) which I can start with? Or are there any source codes which are publicly available that will serve for my purpose? I know it's a bit general but what might be the critical issues which I should pay special attention? Thanx everyone in advance, abilge 

March 22, 2005, 18:59 
Re: pressure poisson equation on nonstaggered gri

#2 
Guest
Posts: n/a

Check out Suhas Patankar's book.. It is one of the classical texts on this.


March 23, 2005, 01:07 
Re: pressure poisson equation on nonstaggered gri

#3 
Guest
Posts: n/a

What solver are you using for solving the poisson equation? Have you tried FISHPACK? Ususally for nonstaggered grid, the RhieChow interpolation is used to avoid pressurevelocity decoupling.


March 23, 2005, 03:43 
Re: pressure poisson equation on nonstaggered gri

#4 
Guest
Posts: n/a

Hi Sachin,
Thanx for your fast response first of all... Correct me if I'm wrong but as far as I remember that book fully describes SIMPLE method for incompressible flows involving heat transfer and again in cartesian coordinates.I guess you are right, it is a classical text on the subject matter since every paper I come across cites that book. Perhaps I was not clear with my question. My goal is to solve the pressure Poisson equation in cylindrical coordinates to obtain the pressure field. And while doing this I do not want to deal with staggered grids. So I believe what I want to do is different from SIMPLE algorithm. I know that this is doable by referring to a paper by Abdallah (1987). But in that paper things are still in cartesian coordinates. Do I ask for too much Thanx for your time once again. abilge 

March 23, 2005, 03:55 
Re: pressure poisson equation on nonstaggered gri

#5 
Guest
Posts: n/a

Hi agg,
I want to thank you too for your fast response. Let me describe the situation once more for clarity: I have obtained the pressure Poisson equation in cylindrical coordinates by taking the divergence of momentum equations and invoking equation of continuity for zero divergence condition. And I intend to use Neuman type boundary conditions except for the exit (we are talking about a laminar pipe flow where the exit is open to atmosphere so the exit boundary condition is Driechlecht). For the solution of the Pressure Poisson equation I use successive over relaxation (SOR) with a relaxation coefficient of 1.7. Coming to your comment; No I have not heard about FISHPACK. What is it used for? And you've talked about "RhieChow interpolation" to avoid velocitypressure decoupling? Can you suggest any references for that? I believe this is what I'm looking for since I heard the most troublesome aspect of the subject at hand occurs while implementing the Neumann type boundary conditions. Perhaps what you have suggested might help.. Than you once again for your time.. abilge 

March 23, 2005, 10:31 
Re: pressure poisson equation on nonstaggered gri

#6 
Guest
Posts: n/a

Hi,
Have you tried underrelaxing? say alfa_p' = 0.6 (alfa_p<1). "SIMPLE" like algorithms usually require the pressure corrections to be underrelaxed. However,Energy equation may zip through with overrelaxing. Hope this helps Abhijit 

March 23, 2005, 11:42 
Re: pressure poisson equation on nonstaggered gri

#7 
Guest
Posts: n/a

Hi Abhijit,
Well as I mentioned previously, I tried overrelaxation upto now. Perhaps underrelaxation might help but since my eventual goal is to couple energy equation, this might cause problem as you told. However, it is worth trying I believe since energy equation is not present in the solution at the moment so I'll give it a try and let you know about the outcome. Thanx a lot for your time... abilge 

March 23, 2005, 12:12 
Re: pressure poisson equation on nonstaggered gri

#8 
Guest
Posts: n/a

Hi,
I forgot to mention, that one set of relaxation parameters (for u,v,w,p' & Temperature/Energy) in general will not work for all problems, you need to play with those. Relaxation parameters are very much problem dependent. If you use SIMPLE you may need to underrelax momentum equations as well (again its problem dependant). This is usually the case unless you use sophisticated algorithms (MG,Preconditioners....) Problems with large gradients in any of the primitive variables usually require underrelaxation. If underrelaxing all the variables does'nt work, then it could be bugs in the code. To begin with you can try a very simple case of 2D flow between fixed parallel plates. Abhijit 

March 23, 2005, 13:40 
Re: pressure poisson equation on nonstaggered gri

#9 
Guest
Posts: n/a

Hi Abhijit ,
That seems like just the kind of information I'm looking for. As you might appreciate you may have to go through thousands of references just to come across such an information. By the way I guess I am about the discover the bug you are talking about on a 2D laminar pipe problem. If I handle it then I may proceed with adding energy equation. Then I'll decide whether I should move on to more sophisticated algorithms such as multigrid methods or so because it seems like I'm gonna suffer from the slow rate of convergence. I'll let you know.. Thanx once again for your time.. abilge 

March 24, 2005, 15:32 
Re: pressure poisson equation on nonstaggered gri

#10 
Guest
Posts: n/a

For FISHPACK (uses cyclic reduction technique): http://www.scd.ucar.edu/css/software/fishpack/
For MUDPACK (uses multigrid techniques): http://www.scd.ucar.edu/css/software/mudpack/ You can use the above solvers to compare results with the method you use now, so you know your implementation is correct. I have used Fishpack, it uses secondorder discretization and is fast. You can also try Mudpack. The RhieChow interpolation is not just for the boundary, but in general to have the discretization of Pressure Poisson equation consistent with that of divergence and pressure gradient operator. You can find the refence in any standard CFD text book or search on google. 

March 25, 2005, 11:20 
Re: pressure poisson equation on nonstaggered gri

#11 
Guest
Posts: n/a

Hi agg,
Well I still cannot believe you sent me the links to source files I hope I can manage to incorporate these subprograms in my code. If I can it would be wonderful. By the way I solved my previous problem (solution of cylindrical pressure Poisson equation on nonstaggered grids). However, as you might guess, the code suffers from low rate of convergence since I do not use multigrid methods or so..Perhaps the source that you have pointed out may help me handle this issue. In the mean time I looked for documents regarding to RhieChow interpolation. Now I'm going over those.. Thanx for everything especially for your time.. abilge 

March 29, 2005, 09:33 
Re: pressure poisson equation on nonstaggered gri

#12 
Guest
Posts: n/a

Hi agg,
I just want to ask you something: Suppose that I want to use the nonstaggered version of FISHPACK to solve the Poisson equation for pressure and supplied the righthandside function to FISHPACK an d it produced a pressure field (and we don't know for sure whether this pressure field is the desired field or not). My question is: Do I have to compare the results produced by the consecutive calls of the FISHPACK and carry out an error analysis i.e. the global error is less than a certain value? This might sound like a very nerdy question to you (I want to apologize in advance if it is) since I could not go into FISHPACK (I was busy fixing the bugs with my existing code, thanx GOD I've fixed it) but I could not find any sort of error definitions in the arguement list of the package. Your answer might help me plan things in advance before I jump into things. I would like to thank you for everything once again. You were like a candle in my darkness by showing me an alternative way (seriously!). abilge 

March 29, 2005, 12:34 
Re: pressure poisson equation on nonstaggered gri

#13 
Guest
Posts: n/a

No you do not have to check consequtive calls. What you can do is for the same right hand side, you can call Fishpack once and also your solver once and then check the two solutions. If there is a significant difference, it could be because of boundary conditions (if you use iterative technique and Fishpack uses direct methods) OR you have implemented Fishpack incorrectly in which case, you will have to use a right hand side for which you already know a solution as a check OR the implementation of your pressure solver is incorrect (but since you say that you have fixed the bugs, unlikely to be the case, but no harm in checking again right!)
OR You can skip Fishpack altogether and just use your solver for a standard Laplace problem for which analytical solution exists. You will find such a problem in "Computational fluid dynamics" by T.J. Chung and more in "Conduction of heat in solids" by H.S. Carslaw and J.C. Jaeger 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Pressure work term in turbulent K.E. equation  lost.identity  Main CFD Forum  0  March 8, 2011 13:21 
Does star cd takes reference pressure?  monica  Siemens  1  April 19, 2007 11:26 
FEM pressure poisson equation (Implicit)  cfdbeginner  Main CFD Forum  0  August 9, 2005 13:32 
Gas pressure question  Dan Moskal  Main CFD Forum  0  October 24, 2002 22:02 
pressure poisson equation!!  frederic felten  Main CFD Forum  1  November 17, 2000 04:17 