On Tue, Oct 12, 2010 at 5:26 PM, BIN ZHANG <zhn...@gmail.com> wrote: > Dear Daniel: > > Thanks a lot for your suggestions. Now I'm actually able to solve the > problem with normal boundary conditions. But I still have one extra > question, is it possible for me to use complicated boundaries like the one > shown in the attachment. For the white circle on the left, I would have > phi=0, and for the white circle on the right, I would have phi=1. > > I tried to define these boundaries using mesh.getInteriorFaces(), but I got > an error about: "Face list has interior faces" > > Is there a way to enforce these complicated boundaries in fipy? >
You can't have interior boundary conditions in FiPy. Cast the boundary condition to a volume integral and apply a source term based on location or use Benny's suggestion, which can more accurately capture the boundary. -- Daniel Wheeler
