The Monte Carlo method is quite useful in the modeling of particulate systems. It is used here to simulate the particle brekage process during grinding that can be represented by a population balance equation. The simulation technique is foe from discretization of time or size. The results of simulation under restricted conditions of grinding compare very well with the available analytical solution of the population balance equation. The procedure is extended to simulate the grinding process in its entirety. This method provides an alternative to the modeling of the grinding process where the governing population balance equation cannot be readily solved.