Just to clarify: as far as FiPy is concerned, "interiorFaces" have a cell on either side of them, whereas "exteriorFaces" have a cell on only one side; the other side is "outside" the mesh. It doesn't matter whether the mesh has simple topological connectivity (square, sphere, etc.) or complex topological connectivity (torus). "interiorFaces" separate "mesh" from "more mesh"; "exteriorFaces" separate "mesh" from "not mesh".
On Nov 19, 2014, at 10:53 AM, Daniel Wheeler <[email protected]> wrote: > Hi Kyle, > > #Boundary Conditions > phi.constrain(X,mesh.exteriorFaces) > phi.constrain(X,mesh.interiorFaces) > > The above is the problem. You are constraining the internal faces, which > makes no sense in FiPy. I am not even sure how FiPy behaves when that > constraint is added. However, I assume that is not what you want to do. If > you remove that constraint the result seems to look nice. BTW, > "exteriorFaces" refers to both the inner and outer external faces of the > annulus. Perhaps you thought "interiorFaces" referred to the interior annulus > faces that are external to the mesh. > > Also, you don't need to import the results into Mathematica. You can just use > "viewer = fipy.Viewer(phi); viewer.plot()" to look at the results. > > Cheers, > > Daniel > > -- > Daniel Wheeler > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
