The open-source C++ library OpenFOAM version 16.12, in which a new solver based on Eulerian stochastic field (ESF) method coupled with FGM has been implemented, is used. The method (ESF/FGM) in large eddy simulation framework consists first of all of solving the filtered transport equations for the control variables mixture fraction and progress variable together with the classical filtered transport equations for mass density and momentum. This yields to an unclosed reaction source term. The effect of turbulent scales on the chemical reactions at the sub-grid scale must be taken into account to correctly represent the thermochemical state and close the reaction source term. For this purpose, the knowledge of the sub-grid evolution of the controlling variables is needed. Thus, the Eulerian Stochastic Fields (ESF) method is adopted, which is based on the transport equation of the Favre-filtered joint scalar probability density function (FDF) of the controlling variables. In this approach, the FDF is represented as an ensemble of number of Eulerian stochastic fields for each of the table controlling variables, for which additional stochastic differential equations are solved. The simulations were carried on to solve the set of obtained partial differential equations applying the merged PISO- SIMPLE algorithm and low Mach number formulation.