|
[Sponsors] |
Neumann Boundary Condition for Poisson Equation solution in Polar Coordinates |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 29, 2011, 15: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: 17 |
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, 11:55 |
|
#2 |
New Member
Amir
Join Date: Oct 2010
Posts: 16
Rep Power: 16 |
Hello,
I have the same problem, any guidance is appreciated. no reply since 2009! |
|
July 30, 2011, 20:07 |
|
#3 | |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,286
Rep Power: 34 |
Quote:
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. |
||
Tags |
bicgstab, conjugate gradient, poisson equation, polar coordinates, singularity |
|
|
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 |