Abstract Unsaturated flow of liquid in a bed of uniform and spherical ore particles is studied numerically and experimentally. An unsteady and two-dimensional model is developed based on the mass conservation equations of liquid phase in the bed and in the particles. The model equations are solved using a fully implicit finite difference method giving the distribution of the degree of saturation in the particles and in the bed and the vertical velocity of flow in the bed, as well as, the effect of periodic infiltration on the above distributions. To calibrate the computational model, several column tests are performed using periodic infiltration of water on 40 cm high columns composed of ore having particles smaller than 25 mm. The numerical analysis shows that (a) the results obtained from numerical modelling under the same operating conditions as used for column tests, are in good agreement with those from experimental procedure, (b) the degree of saturation of the bed and the time required to reach steady state conditions depend on the inflow of water and intrinsic permeability of the bed and (c) the velocity fluctuations and the fluctuations of the degree of saturation in the bed depend on the inflow of water, period of infiltration, height and intrinsic permeability of the bed.