Investigating pore network models as numerical tools to solve steady saturated groundwater flow