Hi,
I've got a PDE I'm trying to solve, and the examples + manual are
getting me *almost* there, but there are a couple of things I am having
trouble with. I figure someone here should be able to help me pretty
easily as I think it is a fairly straight forward application of fipy.
here is the equation (solving for Psi, as a function of R and z):
R*(cost1)*Psi + 1/R * (const2)*Psi + del.(1/R * del Psi)
with boundary conditions Psi = 0 on all four sides.
So I have two source terms, and one diffusion term, with a variable
diffusion co-efficient, and that's where I'm having trouble. I don't
quite understand how to set up my diffusion coefficient (1/R), or the R
and 1/R dependance in my two source terms. The main two examples I've
been using for reference are
examples/diffusion/steadyState/otherMeshes/grid3Dinput.py (the second
part) for setting up a 2D problem, and
examples/diffusion/variable/input.py for having a variable diffusion
coefficient. Does anyone know of a better example I should be looking
at? If someone could show me how to set up and/or solve this problem up
I would really appreciate it.
Thanks
Bjørn
--
Bjørn Hansen
University of Washington
Redmond Plasma Physics Laboratory
14700 NE 95th Street, Suite 100
Redmond, WA 98052
425-881-7706 (Fax)425-882-9137
