Re: polar coordinate singularity

[Sanjay Govindjee <sanjay@xxxxxxxxxxxxxxx>]

There should be nothing special needed to do this in general for finite 
difference (or
finite element) methods
other than the usual need for one sided derivatives at the boundary of 
the domain.

-sg

Barry Smith wrote:
>   A simple check in google of finite differences polar singularity
> gave lots of possible references. One near the top that may be
> relevent was 
> http://enstrophy.colorado.edu/~mohseni/PSpdf/MyPapers/JCP2000.pdf
>
>   PETSc itself doesn't address the nitty-gritty details of 
> differencing schemes
> so I don't think has any particular tool to help manage the differencing.
>
>    Barry
>
>
> On Tue, 9 May 2006, Sean Dettrick wrote:
>
> > Hi,
> > I am solving Poisson's equation
> > Div. Grad(Phi)=rho
> > in the (R,Theta) plane with KSP, and it is working well if I avoid 
> > the axis. But I would like to put a point at R=0.   Is there a 
> > recommended way to do that?
> > Thanks,
> > Sean