Poisson solver

Steve

Hi, I am writing a code for solving viscous incompressible flows using Navier Stokes equations in (u,v,p). Using normal cartesian non-staggered grid.

I am using the fractional step method, so I will be solving for intermediate velocities and then a pressure Poisson equation to get the new velocities.

Presently I am using the point by point SOR method for my Poisson solver. I would like to ask what Poisson solution method would be better, if I am still required to solve on point-by-point basis. I am thinking of using conjugate gradient, and for my code, multigrid method is not advisable.

Another thing is, for conjugate gradient, do I need to include my boundary equations in Ax=b ??

Thanks a lot!!