Integral over internal boundary using facet normal

See @MiroK’s answer here: