Este trabalho utiliza um procedimento alternativo para lidar com uma classe de sistemas hiperbólicos não lineares, empregando um modelo preliminar para descrever o movimento de poluentes em uma atmosfera, representado por um sistema hiperbólico não linear de m+2 equações diferenciais parciais, a saber: conservação de massa e momentum para uma mistura com múltiplos componentes de gases e ar e m equações de balanço de massa para os gases. A simulação deste sistema segue um procedimento sistemático, tratando um problema simultâneo como um problema seqüencial, através da fatoração do operador em uma parte ordinária não homogênea (dependente do tempo) e o operador hiperbólico homogêneo associado. Este último é simulado por um esquema de Glimm de evolução no tempo, especificamente indicado para a simulação de problemas hiperbólicos, apresentando precisão assegurada matematicamente, baseado em uma teoria com sólido embasamento termodinâmico. Este esquema é o que melhor preserva a identidade do choque, isto é, a magnitude e a posição da onda de choque. Entretanto, além de sua limitação a problemas unidimensionais, sua implementação requer a solução completa, ou aproximada, do problema de Riemann associado. Neste trabalho emprega-se um aproximante de Riemann – buscando evitar a exigência de uma solução completa do problema de Riemann associado – para cada dois passos consecutivos, aproximando, assim a solução do problema de Riemann associado por funções constantes por partes, sempre satisfazendo a condição de salto sem, necessariamente, satisfazer a condição de entropia. Comparações entre os resultados obtidos através da solução exata do problema de Riemann associado com aqueles obtidos pelo aproximante de Riemann mostram o bom desempenho desta metodologia.