The present paper describes the application of local perturbations in Large Eddy Simulation to enhance transition to turbulence applied on the heated wall of a natural convection boundary layer. Simulations are carried out with a Prandtl number Pr=0.71 and a maximum Grashof number of Gr=4.1×1011. The purpose of this work is to accelerate the transition process in the numerical domain, reducing the needed simulation effort to reach the turbulent regime in a natural convection boundary layer flow. A parametric study has been performed to identify the fastest transition to turbulence in the proposed technique. The obtained results are compared with experimental, empirical and other numerical reference data to assess the method’s performance. The analysis of the transition method contemplates its effect in the mean flow field, as well as in the resolved Reynolds stresses and turbulent heat fluxes.