dolfin team mailing list archive

["Johan Hoffman" <jhoffman@xxxxxxxxxx>]

Garth and Anders,

I am interested in solving a problem where I a priori specifies an inner
surface (aligned with cell faces) where the solution is discontinuous. The
coupling of the solution over the discontinuity is taken cared of by a
surface integral over that inner surface (like a penalty term).

For this I expect to need:

(1) a data structure with "double nodes" at a predefined surface (=
collection of nodes/faces).

(2) for the nodes on the surface; the ability to assemble only
contribution from either side of the surface for the "double nodes"
respectively (where the two sides of the surface may be represented by a
MeshFunction for example).

(3) a surface integral defined over an inner surface.

I suspect that some parts of (1)-(3) is already in place, as part of the
DG work? What remains to be done to handle this problem?

/Johan