[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: polar coordinate singularity




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