The accuracy of large-eddy simulation (LES) of a turbulent premixed Bunsen flame is investigated in this paper. To distinguish between discretization and modeling errors, multiple LES, using different grid sizes h but the same filterwidth , are compared with the direct numerical simulation (DNS). In addition, LES using various values of but the same ratio / h are compared. The chemistry in the LES and DNS is parametrized with the standard steady premixed flamelet for stochiometric methane-air combustion. The subgrid terms are closed with an eddyviscosity or eddy-diffusivity approach, with an exception of the dominant subgrid term, which is the subgrid part of the chemical source term. The latter subgrid contribution is modeled by a similarity model based upon 2 , which is found to be superior to such a model based upon . Using the 2 similarity model for the subgrid chemistry the LES produces good results, certainly in view of the fact that the LES is completely wrong if the subgrid chemistry model is omitted. The grid refinements of the LES show that the results for = h do depend on the numerical scheme, much more than for h = /2 and h = /4. Nevertheless, modeling errors and discretization error may partially cancel each other; occasionally the = h results were more accurate than the h ≤ results. Finally, for this flame LES results obtained with the present similarity model are shown to be slightly better than those obtained with standard β-pdf closure for the subgrid chemistry.