Muitos dispositivos na indústria aeroespacial trabalham com combustão turbulenta nos seus sistemas de propulsão. O conhecimento sobre este fenômeno é muito importante para explicar seu comportamento físico e químico. Neste contexto, o estudo de modelos que incluam transporte de escalares na zona de reação é importante para compreender e desenvolver máquinas mais eficientes. Ao utilizar dinâmica de fluidos computacional (CFD) em conjunto com os métodos experimentais é possível reduzir o custo e tempo de desenvolvimento destas máquinas. O objetivo principal deste trabalho é adaptar um modelo de combustão turbulenta para chamas difusivas não pré-misturadas, que utiliza a metodologia das equações médias de Reynolds, para um modelo que representa as grandezas turbulentas utilizando simulação de grandes escalas. O modelo foi implementado no código aberto OpenFoam sob o nome de reactingFoamLES. A estrutura numérica foi desenvolvida a partir de um modelo adaptado para o tratamento de combustão para chamas não pré-misturadas com química detalhada.