Hi,
I want to do a boundary integral for an open surface, namely a half-cylinder, in order to implement boundary conditions in the Nitsche fashion. Parts of my code now look like
geo = CSGeometry()
r = 1/pi
cyl = Cylinder(Pnt(0,0,0), Pnt(1,0,0),r)
top = Plane(Pnt(0,0,0),Vec(0,0,-1))
left = Plane(Pnt(0,0,0),Vec(-1,0,0))
right =Plane(Pnt(1,0,0),Vec(1,0,0))
shell = (cyl * left*right*top)
geo.AddSurface(cyl, shell)
################
mesh = Mesh(geo.GenerateMesh(maxh=0.3))
bnd_D = mesh.Boundaries("dir")
in order to define the mesh and the boundary. I then perform the integral using e.g.
a += SymbolicBFI(InnerProduct(p*mu, v.Trace()), BND, element_boundary = True, definedon = bnd_D)
However, it seems that NGSolve skips this integral. Does anyone know what I am doing wrong?
Best regards,
Alvar