The edge-based face element method for 3D-stream function and flux calculations in porous media flow