# Neumann Boundary Condition for Poisson Equation solution in Polar Coordinates

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

 March 29, 2011, 14:12 Neumann Boundary Condition for Poisson Equation solution in Polar Coordinates #1 Senior Member   Prapanch Nair Join Date: Mar 2009 Location: Bangalore, India Posts: 105 Rep Power: 8 I am attempting to solve the poisson equation in polar coordinates (2D). The problem domain is an annulus. At both 'r' boundaries, neumann homogenous boundary conditions. And theta is periodic. I am attempting to solve this using BiCGStab solver, that I have written myself. For the above BCs the matrix is singular and hence my solution isn't converging. To force an additive constant. I made 'nth' row of my matrix the nth row of identity matrix. And changed the correspoding b(n) to 0. Thus I attempted to make '0' the additive constant. When I did that, the solution at the nth position has a sharp peak in the solution. I have used the same strategy with CG solver on a symmetric system and it worked. It is to be noted that this solution with the peak is a converged solution. When i change the additive constant at n to values other than 0, I am getting a scaled shift in my solutions. Also note that A matrix is nonsymmetric for such a problem in polar coordinates. My BiCGSTAB is working fine with dirichlet BCs either wholly or partially on the boundary.

 July 29, 2011, 10:55 #2 New Member   Amir Join Date: Oct 2010 Posts: 16 Rep Power: 6 Hello, I have the same problem, any guidance is appreciated. no reply since 2009!

July 30, 2011, 19:07
#3
Senior Member

Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 369
Rep Power: 10
Quote:
 Originally Posted by Amir_Ghasemi Hello, I have the same problem, any guidance is appreciated. no reply since 2009!
It is solved just like any other Neumann system is solved. Nothing special.

Not more than that could be said about it. Now it is upto you or user to implement things properly.

I have solved pure neumann equation in 3d and 2d in polor coordinates without any problems of any sort. Can use BicCGstab or multigrid, no probs here.

 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 CFD XUE FLUENT 0 July 8, 2010 06:49 bearcharge Main CFD Forum 0 May 14, 2010 11:32 normunds OpenFOAM Running, Solving & CFD 2 June 4, 2007 06:45 Ciro Cannavacciuolo FLUENT 0 February 8, 2007 04:46 Jan Ramboer Main CFD Forum 11 August 16, 1999 08:59

All times are GMT -4. The time now is 13:48.