Hello,
And thanks for answering.
I tried running the example script, but I got this error message
./sphere.py
Traceback (most recent call last):
File "./sphere.py", line 147, in <module>
exec(fipy.tests.doctestPlus._getScript())
File "<string>", line 9, in <module>
NameError: name 'Gmsh2DIn3DSpace' is not defined
Am I making any mistake?
Many thanks
Lorenzo
On 11/01/2011 04:18 PM, Daniel Wheeler wrote:
> 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
> <[email protected] <mailto:[email protected]>> 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
> [email protected] <mailto:[email protected]>
> http://www.ctcms.nist.gov/fipy
> [NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy]
>
>
>
>
> --
> Daniel Wheeler
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy]
