[pnm] BoundaryFlux class only works for Dirichlet Boundaries
The BoundaryFlux class for evaluating fluxes entering or leaving the network only works if Dirichlet conditions are set on the relevant boundaries.
This is because the class evaluates the local residual at the respective DOFs. For Dirichlet BCs, the value is not zero but equals the boundary flux.
One could maybe check if a Dirichlet BC is set and proceed as before and just evaluate the source
term instead otherwise.
Note: source
corresponds to a neumann
BC at the boundary. The PNM does not support neumann
but requires a source
instead because there are no boundary scvfs
.