The fracture network (left) and calculated flow characteristics (right).  Constant heads are assigned along the left (20 m) and right (0 m) boundaries.  Notice that small head gradients coincide with areas with large fracture density, while large head gradients fall in areas of non-conducting matrix. The modeling results suggest that fracture media may not be treated as an "equivalent porous media" with uniform conductivity and head distribution.