This dissertation presents a computational model for the solution of the equations for the turbulent combustion of a pre-mixed gas. The method presented is a hybrid random choice method With the assumption of incompressibility, the Navier-Stokes equations governing the fluid motion may be uncoupled from the equations governing the chemistry. The solution is modelled using a modification of the random vortex method The reaction-diffusion equations governing time and spatial dependence of temperature and species mass-fractions are solved numerically using the random choice method with allowance for the reaction-diffusion terms. This technique is tested on two problems with analytic solutions. The combustion model is applied first to a mixing layer flow and then to a flow past a back-facing step. The method demonstrates an ability to resolve turbulent effects on flame front propagation