No presente trabalho apresenta-se um estudo da equação de difusão bidimensional transiente com reação química na aplicação de um caso hipotético. As hipóteses consideradas permitem reproduzir superficialmente o campo de concentração de poluentes liberadas por fontes contínuas em uma região quadrada. Este estudo é motivado pela grande atenção que especialistas têm dado ultimamente para a poluição causada pelo rápido desenvolvimento industrial e tecnológico, com sérios impactos ambientais em rios e lagos. O modelo matemático é discretizado com o Método de Diferenças Finitas – MDF, em malha uniforme, com acurácia de segunda ordem nas variáveis espaciais e primeira ordem na variável temporal. O sistema de equações algébricas é resolvido com o método Gauss-Seidel.  Os resultados demonstram boa concordância com a literatura.