The tabulated premixed conditional moment closure (T-PCMC) method has been shown to provide the capability to model turbulent, premixed methane flames with detailed chemistry and reasonable runtimes in Reynolds-averaged Navier–Stokes (RANS) environment by Martin et al. (2013, “Modeling an Enclosed, Turbulent Reacting Methane Jet With the Premixed Conditional Moment Closure Method,” ASME Paper No. GT2013-95092). Here, the premixed conditional moment closure (PCMC) method is extended to large eddy simulation (LES). The new model is validated with the turbulent, enclosed reacting methane backward facing step data from El Banhawy et al. (1983, “Premixed, Turbulent Combustion of a Sudden-Expansion Flow,” Combust. Flame, 50, pp. 153–165). The experimental data have a rectangular test section at atmospheric pressure and temperature with an inlet velocity of 10.5 m/s and an equivalence ratio of 0.9 for two different step heights. Contours of major species, velocity, and temperature are provided. The T-PCMC model falls into the class of table lookup turbulent combustion models in which the combustion model is solved offline over a range of conditions and stored in a table that is accessed by the computational fluid dynamic (CFD) code using three controlling variables: the reaction progress variable (RPV), variance, and local scalar dissipation rate. The local scalar dissipation rate is used to account for the affects of the small-scale mixing on the reaction rates. A presumed shape beta function probability density function (PDF) is used to account for the effects of subgrid scale (SGS) turbulence on the reactions. SGS models are incorporated for the scalar dissipation and variance. The open source CFD code OpenFOAM is used with the compressible Smagorinsky LES model. Velocity, temperature, and major species are compared to the experimental data. Once validated, this low “runtime” CFD turbulent combustion model will have great utility for designing the next generation of lean premixed (LPM) gas turbine combustors.