This paper presents a level set framework for the modelling of doping effect during surface diffusion phenomena in a granular packing. The molecular flux of the doped compound is related to the chemical potentials of all the diffusion species. The evolution of the grain compact is simulated in three dimensions, based on the resulting kinetic law relating the surface diffusion velocity to the thermodynamic driving force. An anisotropic adaptive mesh, based on the level set function properties is used to refine the mesh in the surroundings of the grain surface. The simulations have been perfomed by using parallel computing strategy.