Dear All, Many months ago I asked some questions about how to solve Laplace equations in 3D on a "swiss cheese" domain. See http://bit.ly/rX3KfN . I did little progress, mainly because I could not work much on this problem. I think I was a little to ambitious to start with that. I now would like to start tackling the 2D ultrasimplified version of that problem: you have a 10 by 10 square, whose side is [-5,5] in some coordinates. Then in the origin (0,0) you have a hole, i.e. a circle of radius 1. You would like to solve
D\nabla^2\rho=0 (where \rho=\rho(x,y) is scalar you can think of a a density of some kind), with boundary conditions rho=1 along the sides of the square and rho=0 along the circle (I can set D=1 if I want to) [Laplace equation, nothing more and nothing less]. I can equally think, if I want to, that the whole is filled with some material with D=0 and use anyway the condition rho=0 at the interface between the two materials. Once I solve Laplace equation, my next step is to calculate the incoming flux to the hole. Any suggestions about the practical implementation of this problem? Is there already any script similar to what I have in mind? Best Regards Lorenzo _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy]
