Governing equations for a fully coupled flowing-reaction-deformation behavior with mass transfer in heap leaching are developed. The model equations are solved using an explicit finite difference method under the conditions of invariable application rate and constant hydraulic head. The distribution of the degree of the saturation, as well as the distributions of the concentration of the reagent and the solute is given. A cubic relationship between the mineral recovery and the leaching duration is obtained based on the numerical results. The relationship can be used to predict the recovery percentage of the valuable metal.