A general numerical model is described for the dissolution kinetics of spherical particles in binary systems for any combination of first order reactions at the particle-matrix interface and long distance diffusion within the matrix. The model is applicable to both finite and infinite media and handles both complete and partial particle dissolution. It is shown that interfacial reactions can have a strong effect on the dissolution kinetics, the solute concentration at the particle-matrix interface and the solute concentration profile in the matrix.