Hello,
I'm trying to use fipy.Gmsh2d to define the following mes[image: grid.png]h.
I could not figure out the simple way to define the Dirichlet boundary
condition on the red edges (Lines 4,5,6) and Neumann conditions elsewhere
(black edges).
thanks in advance,
Alex
geo = '''
Point(1) = {160, 112, 0, 10};
Point(2) = {0, 112, 0, 10};
Point(3) = {0, 0, 0, 10};
Point(4) = {47, 0, 0, 10};
Point(5) = {47, 58, 0, 10};
Point(6) = {94, 58, 0, 10};
Point(7) = {94, 0, 0, 10};
Point(8) = {160, 0, 0, 10};
Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 5};
Line(5) = {5, 6};
Line(6) = {6, 7};
Line(7) = {7, 8};
Line(8) = {8, 1};
// define the boundary
Line Loop(10) = {1, 2, 3, 4, 5, 6, 7, 8};
// define the domain
Plane Surface(1) = {10};
Physical Surface("front") = {1};
'''
mesh = Gmsh2D(geo)
x, y = mesh.cellCenters
import matplotlib.pyplot as plt
plt.figure()
plt.scatter(x,y)
--
To unsubscribe from this group, send email to [email protected]
View this message at https://list.nist.gov/fipy
---
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
