Recently a new numerical model devoted to the full field modeling of microstructural evolutions at the polycrystal scale has been proposed and validated [1]. The latter is based on a level set description of interfaces in a finite element framework. Firstly introduced to model 2D and 3D primary recrystallization with nucleaction [1] and [2], it has then been extended to consider the grain growth stage [3] and [4]. The ability of this approach to model the Zener pinning phenomenon without any assumption concerning the shape of second phase particles was also demonstrated [5]. This model has nevertheless an elevated computational cost and requires many numerical parameters whose calibration is not straightforward. In the present paper, some major improvements of the model which address these two points are discussed. A comparative study is also provided in order to illustrate the gains achieved in terms of computational efficiency and robustness.