Accurate numerical prediction of surface heat transfer in the presence of film cooling within aero-engine has long been a goal of the aero-engine industry. It requires accurate simulation of the turbulent mixing and reaction processes between freestream and the cooling flow. In this study, the Stress Blended Eddy Simulation (SBES) turbulence model is used together with the Flamelet Generated Manifold (FGM) combustion model to calculate the surface heat flux upstream and downstream of an effusion cooling hole. All simulations are run using ANSYS Fluent. The test case corresponds to an experimental rig run at MIT. Calculations are performed for both reacting and non-reacting jet cooling cases across a range of jet-to-stream blowing ratios, and compared with the experimental data. The effects of mesh resolution are also investigated. Calculations are also performed across a range of Damkohler number (i.e. flow to chemical time ratio) from zero to 30, with unity blowing ratio, and the differences in the maximum surface heat flux magnitude in the reacting and non-reacting cases at a specific location downstream of the hole are investigated. Results from these analyses show good correlation with the experimental heat flux data upstream and downstream of the cooling hole, including the heat flux augmentation due to local reaction. Results from the Damkohler number sweep also show a good match with the experimental data across the range investigated.