On Mon, Mar 7, 2011 at 9:28 AM, Jonathan Guyer <gu...@nist.gov> wrote:
>> (2) I actually have D=0 inside the 3 spheres (hence nothing diffuses >> inside them) with the boundary condition \rho=0 on their surface. > > In principle, either is doable. > > It's the \rho=0 boundary condition that could prove to be a problem with (2). > FiPy doesn't generally allow for internal boundary conditions, although I > think the new pho.constrain() model should work on internal faces, but we > haven't tested it at all in that application. I don't think (2) will be a problem, just set a very large source wherever D=0. -- Daniel Wheeler
