On Fri, Jun 14, 2013 at 7:30 PM, Richard Garner <[email protected] > wrote:
> Ideally, one > could make the 2d region's equation a boundary condition for the 3d (in > the correct cell faces of the 3d region of course), but I don't know if > one could specify a constraint that has a 2D transvers Laplacian of the > dependent variable. Hi Richard, It is possible to do this with FiPy. You will probably need to build up source terms in the cells along the adjoining boundaries and identify the boundary regions. Alternatively, I thought that maybe I could have two > coupled equations, one for the 2d region and one for the 3d region, but > then one would have a 3d mesh and the other a 2d mesh. Is that possible? > Yes. You can have multiple meshes and interpolate between them with relative ease. > Finally, perhaps one could make the thin region be just one cell thick in > the thin dimension, and simply treat it as a 3d region. Doesn't seem like > it would add a great deal of numerical burden with the one extra layer. > Obviously not. You are assuming a linear gradient in the z-direction with this approximation. > Anyway, I am curious about how the "gurus" would approach this. As I > mentioned, I do this in Comsol often, but it's not quite clear what they > are doing "under the hood" in this instance. That is one reason why I > like Fipy and am starting to get into it. > Good to hear that you are looking into FiPy. I hope we can support your efforts in the future. I think all the above methods should, hopefully, give the same outcome. Probably the third method requires the least effort and should be the first one to try. -- Daniel Wheeler
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
