The linear relation between the mean rate of product creation and the mean scalar dissipation rate, derived in the seminal paper by K.N.C. Bray ['The interaction between turbulence and combustion', Proceedings of the Combustion Institute, Vol. 17 (1979), pp. 223-233], is the cornerstone for models of premixed turbulent combustion that deal with the dissipation rate in order to close the reaction rate. In the present work, this linear relation is straightforwardly validated by analysing data computed earlier in the 3D Direct Numerical Simulation (DNS) of three statistically stationary, 1D, planar turbulent flames associated with the flamelet regime of premixed combustion. Although the linear relation does not hold at the leading and trailing edges of the mean flame brush, such a result is expected within the framework of Bray's theory. However, the present DNS yields substantially larger (smaller) values of an input parameter c m (or K 2 = 1/(2c m − 1)), involved by the studied linear relation, when compared to the commonly used value of c m = 0.7 (or K 2 = 2.5). To gain further insight into the issue and into the eventual dependence of c m on mixture composition, the DNS data are combined with the results of numerical simulations of stationary, 1D, planar laminar methane-air flames with complex chemistry, with the results being reported in terms of differently defined combustion progress variables c, i.e. the normalised temperature, density, or mole fraction of CH 4 , O 2 , CO 2 or H 2 O. Such a study indicates the dependence of c m both on the definition of c and on the equivalence ratio. Nevertheless, K 2 and c m can be estimated by processing the results of simulations of counterpart laminar premixed flames. Similar conclusions were also drawn by skipping the DNS data, but invoking a presumed beta probability density function in order to evaluate c m for the differently defined c's and various equivalence ratios.