I would like to use Hedge to simulate the following case, but I am not
sure whether it is possible...
What I would like to do is using a particular numerical flux on some
specific elements interfaces, for example to simulate a crack in a
I defined an internal boundary (straight line) in my .msh file and
Hedge is treating it correctly when I use
However, I can't figure out how to define a flux *between the two
edges of the boundary*. I looked at file hedge/mesh/__init__.py to
understand how Hedge deals with internal boundaries and it seems that
the adjacency information on the concerned elements is lost, so that
those elements behave as it there had no neighbor to share fluxes with
on this interface (hedge/mesh/__init__.py:277+).
I imagined two approaches to implement this feature:
* restoring the interface information, by adding
"interfaces.append(els_faces)" between lines 289 and 290 of
hedge/mesh/__init__.py but I am not sure about the implications of
this on the rest of the program
* defining a periodicity condition between the two edges of the
boundary, but I can't make Gmsh mesh this and I'm afraid this will not
Do you think using Hedge for such case is possible? Do you have any
hint to help me implement this?
Thanks in advance