I am writing code to understand the pressure correction algorithm in a cylindrical coordinate system and using axi-symmetric approximation to solve the flow through corrugated pipes for highly viscous fluids. I have solved it in cartesian domain using the immersed boundary approach and it shows a very good match with the results. But the same doesn't seem to work with CCS.
I am still figuring out the problem. First of all there is this "1/r" term which I think would lead to discontinuity at the center, r=0. Secondly while setting up the pressure poisson equation I need to add those extra terms that contain "1/r" as the coefficients. There is one little astonishment I get when I remove this term where my solution for velocity field almost matches with the correct solution. But the pressure terms never agree and go on increasing.
If anyone has any idea please suggest it to me.