Renormalization of the nonlinear sigma model in 2 + epsilon dimensions. Application to the Heisenberg ferromagnets