Modelling variably saturated flow in porous media for heap leach analysis

A comprehensive model of solution flow through heap leach systems requires a numerical solution with the ability to predict variably saturated-unsaturated flow through porous media domains that contain materials with spatially variable properties. The types of flow associated with stockpile leach processes lead to flow problems, such as infiltration into dry soil, drainage, perched water tables and flow through heterogeneous materials. Computational methods for modelling this class of flow behaviour are conventionally based on the classical Richards' equation. However, the governing equations are highly non-linear, difficult to solve, and require iterative numerical solution methods. The pressure-based form of the Richards' equation suffers from poor mass balance, while the mixed form can possess convergence difficulties. An adaptive transformed mixed algorithm is described, which reduces the non-linearity of the problem, optimizes the time step size, and provides a fast, numerically robust scheme that significantly reduces computation (or CPU) time. This method is shown to give fast, accurate solutions on a number of complex flow problems. The utility of this algorithm is illustrated through its implementation within the PHYSICA computational modelling software, which incidentally, provides a framework for a comprehensive heap leach model.