Second phase particles (SPP) play an essential role in controlling grain size and properties of polycrystalline nickel-base superalloys. The full-field methods are of great interest as they can reproduce microstructure evolutions and heterogeneities. Several methods have been developed to reproduce the Smith-Zener pinning (SZP) mechanism at the polycrystalline scale, considering the presence of SPP in the microstructure. In this context, complex deterministic numerical frameworks have been developed in the last decades, such as the vertex methods, and different front-capturing approaches, including the multi phase-field (MPF) and level-set (LS) frameworks. The LS approach is an advanced formalism where the physics of the interaction between the GB and SPP can be captured in 2D or 3D and for static or dynamic SPP populations. However, when SPP become very small comparatively to the mean grain size, the cost of such calculations can become prohibitive even in 2D. In this context, a new 2D front-tracking methodology was recently developed to model ReX and GG with the same precision that the LS method but with drastically reduced calculation times. This presentation will be dedicated to the enhancement of this front-tracking approach to consider the SZP with the same precision that the LS approach but with a much more appreciable efficiency.