Messed up the reply again. ---------- Forwarded message ---------- From: Daniel Wheeler <daniel.wheel...@gmail.com> Date: Tue, Nov 1, 2011 at 11:18 AM Subject: Re: Simple 2D Problem To: Lorenzo Isella <lorenzo.ise...@gmail.com>
I think Gmsh should be able to do this. Try it. Tale a look at < http://matforge.org/fipy/browser/trunk/examples/cahnHilliard/sphere.py> to get you started with the meshing. When you run the problem, be sure to use "DiffusionTermNoCorrection" rather than the regular "DiffusionTerm". On Mon, Oct 31, 2011 at 3:41 PM, Lorenzo Isella <lorenzo.ise...@gmail.com>wrote: > 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 > fipy@nist.gov > http://www.ctcms.nist.gov/fipy > [NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy] > -- Daniel Wheeler -- Daniel Wheeler
_______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy]
