The present work addresses the coupling of a flamelet database, to a low-Mach approximation of the Navier–Stokes equations using scalar controlling variables. The model is characterized by the chemistry tabulation based on laminar premixed flamelets in combination with an optimal choice of the reaction progress variable, which is determined based on the computational singular perturbation (CSP) method. The formulation of the model focuses on turbulent premixed flames taking into account the effect of heat losses, but it is easily extended to partially premixed and non-premixed regimes. The model is designed for applications in both, Reynolds-averaged Navier–Stokes (RANS) as well as large-eddy simulations (LES) and results for the two methods are compared. A priori analysis of the database is presented to demonstrate the influence of the reaction progress definition and the chemistry tabulation is validated against a one-dimensional premixed laminar flame. The validation of the turbulent case is performed using a turbulent premixed confined jet flame subject to strong heat losses, in which the model shows a good overall performance.