O presente trabalho consistiu no desenvolvimento de rotinas computacionais para a previsão de emissões em motores de combustão interna a partir do acoplamento da metodologia do equilíbrio químico com a cinética química. A modelagem matemática da metodologia do equilíbrio químico foi baseada em balanços de massas e em relações de equilíbrio gerando um sistema não-linear com vinte e duas incógnitas e vinte e duas equações resolvidas pelo método de Newton-Raphson. Foi possível calcular as concentrações de vinte e uma espécies químicas geradas na combustão de diversos combustíveis, misturas de combustíveis, como etanol e gasolina, para uma determinada pressão, temperatura e razão combustível/ar. Após o desenvolvimento da rotina de equilíbrio químico realizou-se um acoplamento com uma rotina de cinética química para a previsão de espécies que possuem uma cinética lenta, como os compostos nitrogenados e o monóxido de carbono. Utilizou-se seis equações para a previsão de espécies pela cinética química sendo três delas pertencentes ao mecanismo de Zeldovich.