Hello,

I'm new to both FiPy and Gmsh.
after days of struggling with solving Laplace eq. for a 3D space and a
sphere,
with robin boundary condition on the sphere's surface, and fixed value b.c.
near "infinity",
I think it is time I'd ask for advice:

given a sphere with radius r=R,
and the B.Cs:
1.    -\alpha c +  \frac{\partial c}{\partial r}= 0 , where r=R
2.    c=0                                                    , where r>>R

I want to solve Laplace equation: {\nabla^2 c = 0}
outside the sphere.

now, what would be the better approach? and how can I overcome the following
problems?

1.generate a sphere using gmsh.
   define the robin BC using getExteriorFaces() on the sphere's surface.
   the problem: how can I define BC near infinity? (r>>R),  e.g. for R=1,
r=100 is infinite enough.

2. generate a sphere with spherical hole in the middle.
    define the BC in the "infinity" using getExteriorFaces()
    the problem: how can i define the robin BC on the inner hole-sphere
surface?

or is there a better way?

I would really appreciate your help,

thank you,
Mor
BGU University, physics dep.
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
  [NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy]

Reply via email to