This way we can formulate the flux dependent on both the matrix as
well as the facet permeability. This significantly reduces the condition
number of the resulting matrix and thus improves newton convergence.