On Jan 23, 2013, at 3:39 PM, Laurence wrote: > I'm happy with the mask and how it ascertains to the mesh I've already > created. > > My issue is that I am creating a diffusion field on the mesh, and > at one place in the mesh, a localised reaction is taking place > depending on the concentration at that point. This will deplete > the concentration of the solution at that specific location, > and so cause further diffusion. I have looked into using the levelset > functions, but the surfactant equations seem not to specifically describe > the equations I require. > > I thought if I created a subset of the mesh, I could write PDEs that > only affected those cells. > > Is there a way I can write my PDEs such that the area the equations > effect is limited to certain areas.
You can just use the mask to limit the coefficients to the area of interest, so you could write a diffusion equation with a source that applies everywhere and have the source be zero everywhere but where the reaction is taking place. This just came up a couple of days ago: http://thread.gmane.org/gmane.comp.python.fipy/2878/focus=2885 _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]