A large number of South Africa’s open pit mines are located in fractured rock. Analytical solutions to pit flooding mainly make use of the Dupuit-Forchheimer approximation. A typical example of such a solution is the solution proposed by Marinelli and Niccoli (2000). The challenge with most of these solutions is that the pit geometry is generally chosen as symmetrical and that the solution cannot accommodate geological structures. Numerical solutions to partial differential equations are available in the form of finite-difference, finite-volume and finite-element schemes. These solutions are able to solve for complex geometries and can accommodate geological structures, but require appropriate grid refinement to represent geological structures.