Este trabalho descreve os modelos matemáticos e computacionais utilizados para a simulação numérica do comportamento de motores alternativos de combustão interna. São comparados quatro métodos de diferenças finitas, Filtering Remedy and Methodology (FRAM), MacCormack, Lax-Wendroff de um passo e Lax-Wendroff de dois passos, para a solução das equações do escoamento nos dutos dos sistemas de admissão e exaustão. É apresentada uma metodologia desenvolvida para a obtenção de estabilidade na solução da fronteira dos dutos, em face às variações adversas das propriedades do elemento ao qual o duto está acoplado. Os modelos para plenos, silenciadores, junções e turboalimentadores também são descritos. São mostrados resultados do programa de simulação computacional desenvolvido, para diversas configurações de motores e escoamentos, indicando sua precisão e sua capacidade de tratar as mais diversas configurações de motores.