Imposing a discontinuity in the solution, on an internal interface

If you switch to DG, you need to use a DG variational formulation, see
https://bitbucket.org/fenics-project/dolfin/src/master/python/demo/undocumented/dg-poisson/demo_dg-poisson.py
or a slightly outdated version for DOLFINx that describes some of the math

2 Likes