We investigate how the flamelet combustion model coupled with the collision distance algorithm reproduces turbulent reacting flows. The wrinkling length scale of the flame is modelled as function of the turbulent fluctuating velocity due to shear stress and wake effects calculated by a porous computational mesh. The porous media is based on the Minkowski difference to check the collision among convex objects. The combustion and porosity models are implemented in the framework of an in house Fortran code (Shock Towards Kinetic Explosion Simulator) that solves the full set of Navier-Stokes equations for reacting flows. Results are presented for large Reynolds numbers in relatively small combustion chambers with Karlovitz number around unity. Comparison with experimental data leads to good agreement. The novel approach for calculation of porosities is fully described as it was coded. (C) 2019 Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved.