Simplify both diffusion laws.
Use correct distance weighted harmonic mean for \lambda and only the inside \lambda for fluxes over the boundary. Requires changing some reference solutions where graded meshes are used.
\lambda