Calculate flux over boundary

For now, expression does not support visualization over facets. You would have to project a function into the appropriate space for now (see for instance:How to compute tangent vectors on facets (2D surface parametrization) - #3 by hherlyng)