DG on nonconforming mesh

I am trying to implement a DG method on nonconforming meshes. See the picture below for an example of the mesh. Are there any ideas on how to implement the jumps and averages in NGSolve for the interface where the meshes on both sides do not match well? Thank you very much.


Have a look at


best Christopher

1 Like

Hi Christopher,

Thank you very much for the instructive example. I managed to set up the contribution over the interface in the bilinear form. I still have one further question, how do I set up an integral only over S_1 (see the attach picture)?

I know I can approximate this by setting up an integral over the whole bottom interface of the upper rectangle and multiply an indicator function. But it would be nice if this can also be done by ContactBoundary.

Thanks agian for your help.