A spatial stabilization via bubble functions of linear finite element methods for nonlinear evolutionary convection-diffusion equations is discussed. The method of lines with SUPG discretization in space leads to numerical schemes that are not only difficult to implement, when considering nonlinear evolutionary equations, but also do not produce satisfactory results. The method we propose can be seen as an alternative to this kind of methods. Once the numerical approximation belonging to a linear finite element space enriched with bubble functions is obtained, the bubble part is discarded. The linear part is shown to give a stabilized approximation to the solution being approached. The bubble functions are deduced using a linear steady convection-diffusion model in such a way that the linear part of the approximation to the linear steady model (after static condensation of bubbles) reproduces the SUPG method in the convection-dominated regime. However, for the nonlinear evolutionary equations we consider in the paper the method we propose is not equivalent to the SUPG method. Some numerical experiments are provided in the paper to show the efficiency of the procedure.