Dear list, The constraints system is extremely elegant way of applying boundary conditions to a problem. I was just wondering if there is any background information (notes, a paper or thesis) that covers how this is implemented?
I have been solving advection-diffusion-reaction type problems using the finite volume method (see my online notes here, http://danieljfarrell.github.io/FVM) and have found a elegant way of adding general boundary conditions to this equation. See here, http://danieljfarrell.github.io/FVM/advection_diffusion_boundary_conditions.html#general-matrix-form But this only really works when there is a transient term. As soon as the equation changes form then a I need to rethink how to apply the boundary conditions. The Fipy system seems to get round this problem I just wondered how it is working. I have read the source code and noticed the use of "masks" but I couldn't really glean much more than that. If possible I wish to implement a Fipy like constraints system in my code, but I want to stay in semi-discrete form so that I can used the method of lines. Can the constraints system be applied to the semi-discrete form? Best wishes, Dan
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
