CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Neumann Boundary Condition for Poisson Equation solution in Polar Coordinates

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 29, 2011, 15:12
Default 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: 17
prapanj is on a distinguished road
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.
prapanj is offline   Reply With Quote

Old   July 29, 2011, 11:55
Default
  #2
New Member
 
Amir
Join Date: Oct 2010
Posts: 16
Rep Power: 16
Amir_Ghasemi is on a distinguished road
Hello,
I have the same problem, any guidance is appreciated.
no reply since 2009!
Amir_Ghasemi is offline   Reply With Quote

Old   July 30, 2011, 20:07
Default
  #3
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,285
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by Amir_Ghasemi View Post
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.
arjun is offline   Reply With Quote

Reply

Tags
bicgstab, conjugate gradient, poisson equation, polar coordinates, singularity


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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
External Radiation Boundary Condition (Two sided wall), Grid Interface CFD XUE FLUENT 0 July 8, 2010 07:49
vorticity boundary condition bearcharge Main CFD Forum 0 May 14, 2010 12:32
Slip boundary condition what is inside normunds OpenFOAM Running, Solving & CFD 2 June 4, 2007 07:45
Neumann Boundary Condition Ciro Cannavacciuolo FLUENT 0 February 8, 2007 04:46
Boundary Conditions Jan Ramboer Main CFD Forum 11 August 16, 1999 09:59


All times are GMT -4. The time now is 02:57.