SUMMARYWithin the context of the sintering process simulation, this paper proposes a numerical strategy for the direct simulation of the matter transport by surface diffusion, in two and three dimensions. The level set formulation of the surface diffusion problem is first established. The resulting equations are solved by using a finite element method. A stabilization technique is then introduced, in order to avoid the spurious oscillations of the grain boundary that are a consequence of the dependence of the surface velocity on the fourth-order derivative of the level set function. The convergence and the accuracy of this approach are proved by investigating the change in an elliptic interface under surface diffusion. Cases in direct relation with the sintering process are analyzed besides: sintering between two grains of the same size or of two different sizes. Finally, 3D simulations involving a small number of particles show the ability of the proposed strategy to deal with strong deformations of the grain surface (formation of necks) and to access directly important parameters such as the closed porosity rate.