A crescente demanda mundial por energia, o aumento da produção industrial mundial e o efeito estufa têm gerado um grande problema para a atmosfera no que diz respeito à poluição. As necessidades de conhecimento sobre a física da dispersão de poluentes levaram à criação de ferramentas capazes de reproduzir e responder sobre os possíveis efeitos ao meio ambiente. Túneis de vento de camada limite são ferramentas que possuem condições de simular, em escala reduzida, as emissões de fontes poluidoras tais como chaminés de indústrias e emissões de veículos automotores, bem como permitir a obtenção de soluções aos possíveis danos causados. Este trabalho realizou um estudo exploratório em túnel de vento do comportamento de uma pluma de contaminantes emitidos por uma chaminé isolada, utilizando o hélio como gás traçador. Foram realizados ensaios com diferentes tipos de escoamento, velocidades do vento e velocidades de emissão da pluma. As emissões utilizaram hélio puro e uma mistura binária de hélio e ar. Foram feitos levantamentos do perfil de concentrações no eixo da chaminé, a diferentes posições a sotavento. Avaliou-se também a eficiência de um sistema de anemometria acoplado a uma sonda aspirante como instrumento medidor de concentrações de hélio. A análise dos resultados demonstrou a capacidade do túnel de simular o fenômeno da dispersão de poluentes e a influência das características da rugosidade do terreno no processo de diluição da pluma. Concluiu-se que há a influência aerodinâmica da chaminé sobre a pluma, comprovando-se a existência de três regiões características da pluma, cada uma apresentando um determinado comportamento para as flutuações instantâneas da concentração de hélio. Por último, o sistema de anemometria em conjunto com uma sonda aspirante é somente adequado para estudos de efeitos locais de dispersão de poluentes. / The increasing global energy demand, the growth of the world industrial production and the greenhouse effect have created a great problem for the atmosphere regarding pollution. The need of knowledge about the physics of the pollutants dispersion led to the creation of tools able to simulate and to respond on the possible effects to the environment. Boundary layer wind tunnels are tools capable of simulate, in reduced scale, the emissions from pollutant sources such as chimneys from industries and vehicles, as well as to allow possible hazardous problems to be solved. In this work, an exploratory wind tunnel study is made of the plume behaviour originated from an isolated chimney, using helium as tracer gas. The tests were performed for different types of flows, wind velocities and plume emission velocities. The emissions were from pure helium and a binary mixture of helium and air. The concentrations were measured in several positions leeward of the source. It was also evaluated the efficiency of an anemometry system connected to an aspirating probe as measuring instrument of the helium concentration.The analysis of the results demonstrated the capacity of the wind tunnel to simulate the pollutant dispersion phenomenon and the influence of the terrain roughness in the plume dilution process. It is concluded that the chimney exerts an aerodynamic influence over the plume, proving the existence of three characteristic plume regions, each one presenting a specific behaviour for the instantaneous fluctuations of helium concentrations. Regarding the anemometry system connected to an aspirating probe, it is concluded that it is only adequate for the study of local effects in the dispersion of pollutants.

O objetivo deste trabalho é obter uma nova solução analítica para a equação de advecção-difusão. Para tanto, considera-se um problema bidimensional difusivo-advectivo estacionário com coeficiente de difusão turbulenta vertical variável que modela a dispersão de poluentes na atmosfera. São utilizados três coeficientes difusivos válidos na camada limite convectiva e que dependem da altura, da distância da fonte e do perfil de velocidade. A abordagem utilizada para a resolução do problema é a técnica da Transformada Integral Generalizada, na qual a equação transformada do problema difusivo-advectivo é resolvida pela técnica da Transformada de Laplace com inversão analítica. Nenhuma aproximação é feita durante a derivação da solução, sendo assim, esta é exata exceto pelo erro de truncamento. O modelo ´e avaliado em condições moderadamente instáveis usando o experimento de Copenhagen. Apresentam-se os resultados numéricos e estatísticos, comparando os resultados obtidos com dados experimentais e com os resultados da literatura. O modelo proposto mostrou-se satisfatório em relação aos dados dos experimentos difusivos considerados.

Neste trabalho foi avaliado o impacto na simulação dos poluentes fotoquímicos do uso de diferentes parametrizações da Camada Limite Planetária (CLP) para a Região Metropolitana de São Paulo, utilizando o modelo de mesoescala não-hidrostático Weather Research and Forecasting (WRF) com o acoplamento dos módulos químicos (WRF/Chem). Para análise da acurácia dos resultados obtidos com o modelo, foram utilizadas as medições realizadas pelas estações automáticas da CETESB em superfície, e a validação do perfil vertical da concentração de O3 com base nos dados do experimento de lançamento de ozoniosondas realizado dentro do contexto de um projeto de políticas públicas (recursos FAPESP) realizado durante o período de estudo. A motivação para a realização do presente trabalho foi a grande dificuldade de representação da estrutura da CLP em especial durante o período noturno sobre áreas urbanas. O estudo contou com a utilização inédita do modelo on-line para simulação da composição química da atmosfera da RMSP, e por fim analisou o impacto do acoplamento com o modelo urbano Urban Canopy Model (UCM), com o modelo WRF-Chem. A área de estudo é a Região Metropolitana de São Paulo e adjacências, e o período estudado foi de 28/10 a 01/11/2006. As variações anuais da temperatura do ar, velocidade do vento e radiação solar para o ano de 2006, mostraram que o mês de outubro foi marcado por grande disponibilidade de radiação solar, altos valores de temperatura do ar e ventos intensos. As simulações com casos idealizados mostraram que o esquema numérico de CLP, YSU (Yonsei University) representou melhor a variação da Camada Limite na área urbana, o que influenciou diretamente na dispersão dos poluentes. Já para as simulações de casos reais os resultados mostraram uma boa correlação entre simulado e observado para as concentrações superficiais de O3 e CO e uma grande dificuldade do modelo em representar o perfil vertical da concentração de O3, quando utilizada a versão 2.1 do modelo WRF/Chem. Problema este parcialmente resolvido (podendo melhorar ainda mais) quando é utilizada uma versão mais recente do modelo WRF-Chem, como o caso da versão 3, acoplada com o modelo urbano UCM. Os experimentos realizados para estimar o impacto do tipo do uso do solo sobre a variação da altura da CLP e conseqüente relação com a composição química da atmosfera, mostraram que a altura da CLP é mais alta sobre áreas urbanas em comparação com áreas rurais, portanto influenciando de forma efetiva na dispersão dos poluentes. No presente trabalho foi mostrado o transporte de poluentes existente entre as regiões metropolitanas de São Paulo, Campinas e Sorocaba, devido principalmente a topografia que atua canalizando o vento na direção de Campinas e Sorocaba. A proximidade com as fontes de emissão mostrou ser um fator importante para o nível de concentração de O3, pois à medida que se afasta das fontes de emissão tanto vertical quanto horizontalmente a concentração de O3 aumenta. E devido a sua relativa proximidade com o Oceano Atlântico, a RMSP é beneficiada com a formação da Brisa Marítima que auxilia na dispersão dos poluentes emitidos e formados na RMSP. / This work evaluated the parameterizations´ sensibility of the Planetary Boundary Layer (PBL) in the simulation of photochemical pollutants above the Metropolitan Region of São Paulo, by using the Weather Research and Forecasting (WRF) non-hydrostatic mesoscale model, with the coupling of chemical modules (WRF/Chem). In order to analyze the accuracy of the results obtained with the model, there will be used the measurements made by the automatic air quality stations from CETESB on surface, and the validation of the O3 vertical profile, will be based on ozonesoundings experiments conducted during the period of the study. The motivation for this study was the great difficulty in representing the PBLs structure specially during the night, witch is an new utilization of an on-line model to simulate the chemical composition of the atmosphere of the MASP, and finally, to analyze the impact of the coupling with the Urban Canopy Model (UCM). The study area is the Metropolitan Region of São Paulo and surroundings, and the period was between 10/28 11/01/2006. The annual variations of the air temperature, wind speed and solar radiation for 2006, showed that the month of October was marked by wide availability of solar radiation, high values of air temperature and intense winds. The simulations with idealized cases showed that the PBLs order of height based on the numerical scheme is YSU, MYJ, BOULAC and UCM, which influences directly the dispersion of pollutants with lower concentrations of the CO for YSU and higher for UCM respectively, and for the concentrations of the O3, the scheme YSU shows lower concentrations but unlike the CO, the other schemes are equal. However, for the simulations of real cases, the results showed a good correlation between simulated and observed for the surface concentrations of the O3 and CO, and a great difficulty of the model on the vertical profiles representation of the concentration of O3, when the 2.1 version of the WRF/Chem model is used. This problem was partially solved (can be better improved) when used with the latest version, as version 3, and coupled with the Urban Canopy Model. The experiments performed to estimate the impact of the land use type on the variation of the PBL height and consequent relationship with the chemical composition of the atmosphere, showed that the PBL height is higher on the urban areas compared to rural areas, thus influencing effectively in the dispersal of pollutants, and a lower PBL height was associated with higher concentrations due to lower dispersion of pollutants. This work has demonstrated the transport of pollutants between the metropolitan regions of São Paulo - Campinas and São Paulo - Sorocaba, due to topography that acts targeting the wind in the direction of Campinas. The proximity of emission sources has been important for the level of concentration of O3, as it moves away from sources of emissions both vertically and horizontally increases its concentration. Due the proximity with the Atlantic Ocean, the formation of Sea Breeze helps dispersion of pollutants emitted and formed in the MASP.

Questões ambientais estão no centro das discussões nas últimas décadas. A poluição atmosférica, causada pela expansão pós-revolução industrial fez surgir a necessidade de aprender a descrever, usando modelos matemáticos, esse fenômeno. Com esse conhecimento pode-se propor soluções que mitiguem a poluição e os danos colaterais causados ao ambiente. A dispersão de poluentes modelada por soluções analíticas, a partir das equações de advecção-difusão oferecem um conhecimento sobre cada componente que constrói a equação, característica inexistente em outras abordagens, como a numérica. Entretanto ela era incapaz de descrever propriedades que se referem à turbulência, as estruturas coerentes, causadas por componentes não-lineares suprimidas por construção das equações governantes do modelo. Este trabalho estudou uma forma de recuperar características associadas à turbulência através de uma componente fundamental em estruturas coerentes, a fase. Essa é incluída no modelo que passa a descrever manifestações da turbulência em processos de dispersão através de flutuações de pequena escala na concentração da solução do modelo sesquilinear, que é determinístico-estocástico. No decorrer do trabalho há um estudo através de variações de parâmetros para compreender os efeitos da fase no modelo. Ele também foi aplicado ao experimento de Copenhagen e a dois cenários reais com a intenção de compreender o modelo frente à variáveis micrometeorológicas assim como aprimorá-lo para simular a dispersão de poluentes oriundos de fontes de forma realística. / Environmental issues have been at the center of discussions in the last few decades. Atmospheric pollution, caused by post-industrial revolution, has increased the necessity to describe, using mathematical models, this phenomenon. With this knowledge is possible to propose solutions mitigating the pollution and collateral damages caused in the environment. The pollutant dispersion modeled by analytical solutions, from advection-diffusion equations, offers a knowledge about each component that constructs the equation, a characteristic that does not exist in other approaches, such as numerical. However it was unable to describe properties that refer to turbulence, coherent structures, caused by nonlinear components suppressed by constructing the model governing equations. This work studied a way to recover characteristics associated with turbulence through a fundamental component in coherent structures, the phase. This is included in the model which describes manifestations of turbulence in the dispersion process through the presence of small-scale concentration fluctuations in the sesquilinear model, which is deterministicstochastic. In the course of this work there is a study through variations of parameters to understand the phase effects in the model. It was also applied to Copenhagen experiment and to two real scenarios with the intention of understanding the model regarding micrometeorological variables as well as improving it to simulate the pollutant dispersion from sources in a realistic way.

O presente trabalho propõe um novo modelo para dispersão de poluentes na atmosfera, tal modelo foi idealizado no trabalho de dissertação da autora e continuou-se seu desenvolvimento nesta pesquisa. O modelo é baseado na solução semi-analítica da equação de advecção-difusão para emissão contínua, com resolução através do método de separação de variáveis e da transformada de Fourier. As condições de contorno são tratadas como infinitas reflexões do poluente no solo e no topo da camada limite atmosférica. Adiante, estas reflexões são utilizadas de modo parcial, na tentativa de considerar fenômenos da dispersão que não podem ser explicitados no modelo determinístico, de forma que os contornos podem ser entendidos como estocásticos, ou seja, pode-se interpretar os contornos como uma amostragem de uma distribuição. Além disso, é realizada uma otimização nos contornos parcialmente reflexivos, com o objetivo de desenvolver uma metodologia de otimização e determinar os valores ótimos para a reflexão parcial. Os resultados obtidos foram, primeiramente, comparados com os experimentos de Copenhagen e Hanford. Posteriormente, comparou-se o modelo com dados de concentração coletados em uma fábrica de celulose, a CMPC Celulose Riograndense. Simulou-se, também, a dispersão de poluentes emitidos por uma usina termelétrica no Brasil, que faz parte do programa de pesquisa e desenvolvimento tecnológico do setor de energia elétrica da Agência Nacional de Energia Elétrica (ANEEL). / The present work proposes a new model for pollutant dispersion in the atmosphere, this model was idealized in the dissertation work of the author and continued its development in this research. The model is based on the semi-analytic solution of the advectiondiffusion equation for continuous emission, with resolution through the method of separation of variables and the Fourier transform. The boundary conditions are treated as infinite reflections of the pollutant in the soil and at the top of the atmospheric boundary layer. These reflections are used in a partial way in the attempt to consider phenomena of dispersion that can not be explained in the deterministic model, so that the boundaries can be understood as stochastic, that is, one can interpret the boundaries as a sampling of a distribution. In addition, an optimization is performed in the partially reflective boundaries, with the purpose of developing an optimization methodology and determining the optimal values for the partial reflection. The results obtained were firstly compared with the experiments of Copenhagen and Hanford. Subsequently, the model was compared with concentration data collected at a cellulose production plant. The dispersion of pollutants emitted by a thermoelectric plant in Brazil was also simulated, which is part of the research and technological development program of the electric energy sector of the National Electric Energy Agency (ANEEL).

O presente estudo visa à obtenção de um modelo utilizando o acoplamento dos códigos WRF e CALPUFF com o objetivo de obter as características do campo meteorológico e da dispersão dos poluentes NO2 e SO2 de uma Usina Termelétrica localizada em Linhares. Foi utilizada uma grade de 100×100 células, com resolução de 1 km durante 90 h. Para vias de comparação, outra simulação foi feita sem a entrada do modelo WRF no CALMET, utilizando uma grade 15×15 sem modificar nenhum outro parâmetro. Os resultados obtidos foram comparados com dados do aeroporto de Vitória e com a legislação ambiental vigente. Para a camada limite planetária, o resultado obtido da simulação WRF/CALMET se encontra dentro da faixa encontrada na literatura. Os demais resultados, média de velocidade e direção do vento, para as duas simulações diferem entre si e dos dados do aeroporto, fato que já era esperado devido à distância e orografia entre a Usina e os dados da estação utilizados como entrada nos modelos. A concentração máxima obtida para os poluentes estudados se encontram dentro dos padrões de qualidade do ar. Devido à falta de dados observacionais, não foi possível uma genuína validação dos resultados, mas, sabendo a localização dos picos, foi possível sugerir locais de amostragem para futura comprovação dos resultados. / This study aims to create a model using the coupling of the WRF and CALPUFF codes in order to obtain characteristics of the meteorological field and the dispersion of pollutants NO2 and SO2 of a power plant located at the city of Linhares. The field consists of a grid of 100×100 cells resolution of 1 km for 90 h. Another simulation was made without using WRF as an input into CALMET, in the purpose of comparison, using a 15×15 grid and no change of other parameters. The results were compared to data from the airport of Vitória and against environmental legislation. For the planetary boundary layer the results of WRF/CALMET simulation are within the range found in the literature. The results of average wind speed and direction obtained by both simulations are different from each other and from the data of the airport,such as expected due to the distance and orography of the power plant and station data used as input into the models. The maximum concentrations of the pollutants are within air quality standards. Due to lack of observational data, genuine validation of the results is not feasible, but knowing the location of the concentration peak, it was possible to propose suitable sampling sites for future verification.

Esta tese tem o objetivo de apresentar um modelo matemático, para simular a dispersão de poluentes na atmosfera, que considera a variação temporal do campo de vento e dos coeficientes de difusão turbulenta, além disso, representar uma fonte móvel através de fontes pontuais. Sendo assim, usa-se a ideia do método da decomposição de Adomian e a técnica GILTT (Generalized Integral Laplace Transform Technique) no intuito de resolver a equação de advecção difusão, a qual descreve o fenômeno citado. Ainda, implementa-se o modelo proposto com o conjunto de dados do experimento de OLAD (Over Land Alongwind Dispersion) e, por fim, comparam-se os resultados obtidos e os dados de concentração coletados no experimento mencionado. / This thesis aims to present a mathematical model to simulate the dispersion of pollutants in the atmosphere, which considers the temporal variation of the wind field and the eddy diffusivity. Moreover, it represents a moving source through point sources. To reach this goal, it uses the idea of the Adomian decomposition method together with the GILTT technique (Generalized Integral Laplace Transform Technique) in order to solve the advection-diffusion equation, which describes the phenomenon. It further implements the model proposed with the dataset of OLAD (Dispersion Over Land Alongwind) experiment and finally the results obtained and the concentration observed are compared.

Atualmente, com a melhoria das técnicas de paralelização, cada vez mais vem sendo empregado o uso de clusteres de computadores para programas exigentes de alta performance. Porém, para um programa rodar adequadamente em paralelo, o mesmo deve ser escrito de maneira paralelizável. A proposta deste trabalho é o desenvolvimento de uma plataforma computacional que facilite a escrita de programas orientados à análise de dados climáticos e simulações de fenômenos atmosféricos escritos em paralelo. A plataforma foi pensada de maneira a ser o mais exível possível, permitindo que novas funcionalidades sejam adicionadas através de plugins. Em um futuro próximo, pretende-se disponibilizá-la para uso público com licença GNU General Public License (GNU GPL) e código aberto. Para demonstrar seu potencial de uso, foi realizado um estudo exploratório com dados provenientes do projeto de reanálise Twentieth Century Reanalysis version 2 (20CRv2). Tal estudo visa obter informações sobre que parâmetros e escalas de tempo são importantes para a descrição do fenômeno da zona de convergência do Atlântico sul (ZCAS). / Since parallel computing technologies are improving very fast, computer clusteres are getting more and more employed for highly demanding computational tasks. Nevertheless, for a software to run in a cluster, it needs to be written in a parallel way, which is not necessarily a simple task. Therefore, the aim of this work is to develop a computer platform capable of manipulating, analysing and evaluating climate observational data and controlling simulation in a parallel manner. Features can be added to the platform through plugins, making it very exible and extensible for a very large range of tasks. A simple application is purposed to demonstrate the platform's potential uses. The application consists of a exploratory study of data from the 20th century reanalysis project concerning the South American Convergence Zone (SACZ), a very important phenomena in the south hemisphere.

Localizados na interface entre a terra e o mar, sistemas estuarinos estão propensos a grandes influências de substâncias que são introduzidas por diversas fontes, incluindo rios, deposição atmosférica, efluentes e escoamento de água superficial. Estas áreas, em geral, possuem um forte gradiente de concentração de substâncias químicas, onde essas são transformadas, retidas ou eliminadas. O Canal de Piaçaguera está inserido no Sistema Estuarino de Santos e São Vicente (SESS), localizado na Bacia Hidrográfica da Baixada Santista, situada na região Região Metropolitana da Baixada Santista (RMBS), que compreende os municípios de Cubatão, Santos, São Vicente, Guarujá e Praia Grande. No SESS foram implantadas, na década de 50, diversas indústrias de base (siderurgia, petroquímica e de fertilizantes), que assim como o Porto de Santos, maior da América Latina, detém alto potencial poluidor, o que fez desse estuário um grande receptor de resíduos tóxicos e efluentes líquidos contaminados. A região do Canal de Piaçaguera, área considerada mais crítica em termos de poluição desse estuário, recebe efluentes de diversas indústrias e tem como outras importantes fontes poluidoras as atividades de dragagem, deposição atmosférica, águas provenientes de drenagem superficial, esgotos clandestinos e o Porto de Santos. O presente trabalho buscou ampliar o conhecimento da poluição das águas superficiais do Canal de Piaçaguera e verificar a influência da hidrodinâmica na contaminação dessa área e em outras áreas do estuário. Primeiramente foi realizada uma avaliação dos dados de qualidade das águas superficiais da Companhia Ambiental do Estado de São Paulo (Cetesb) medidos em pontos dos rios Mogi, Piaçaguera e Cubatão, os quais despejam suas águas no canal, dos anos de 2000 a 2010. Para estimar as concentrações e área atingidas foram realizadas simulações da dispersão de poluentes em pontos dentro e fora do canal, a partir da modelagem numérica acoplada, de hidrodinâmica e dispersão. O diagnóstico de qualidade da água indica comprometimento da qualidade ambiental dos corpos dágua estudados, visto que apresentam teores acentuados de poluentes diversos e de fontes difusas, destacando-se as industriais e dos terminais portuários. O rio Piaçaguera contribuiu com o maior número de poluentes e mais elevada concentração deles. A modelagem demonstrou que o Canal de Piaçaguera apresenta circulação com energia bem menor que outras áreas do estuário, sendo notório o menor espalhamento de poluentes nessa área, de modo que a baixa velocidade das correntes no Canal é associada a menor dispersão de partículas. Considerando o lançamento contínuo de poluentes e fatores agravantes como, baixa hidrodinâmica, a bioacumulação e resuspensão dos sedimentos contaminados, é evidente que essa área necessita de maior controle e rigidez na disposição dos materiais com alto potencial poluidor. / Located at the interface between land and sea, estuarine systems are prone to large influences of substances that are introduced by various sources, including rivers, atmospheric deposition, wastewater and surface water runoff. These areas generally have a strong concentration gradient of chemicals where these are transformed, retained or eliminated. The Piaçaguera\'s channel is inserted at the estuarine system of Santos and São Vicente (SESS), located in the Santos Basin, located in the metropolitan area of Santos (RMBS), which includes the cities of Cubatao, Santos, Sao Vicente, Guaruja and Praia Grande. At the SESS were implanted, in the 1950s, several basic industries (steel, petrochemicals and fertilizers), that as the Port of Santos, Latin America\'s largest, has high pollution potential as well as the Port of Santos, what made this a great receiver estuary of toxic waste and contaminated wastewater. The region Piaçagueras channel, an area considered most critical in terms of pollution of the estuary receives effluents from various industries and its other major pollution sources of dredging activities, atmospheric deposition, water from surface drainage, underground running down and the Port of Santos . This study aimed to expand knowledge of the pollution of surface waters of the Channel Piaçaguera and the influence of hydrodynamics on the contamination of the area. We first carried out an evaluation of the data quality of surface waters of the State Environmental Company of São Paulo (Cetesb) at the points measured in the river Mogi, Piaçaguera and Cubatão, which discharge their waters into the canal, of the years 2000 to 2010. To estimate the concentrations and area affected were carried out simulations of pollutant dispersion in points inside and outside the channel from the coupled numerical modeling of hydrodynamics and dispersion. The diagnosis of water quality indicates impairment of environmental quality of the water bodies studied, since they exhibit high values concentration of various pollutants and diffuse sources, especially the industrial and port terminals. The river Piaçaguera contributed the largest number of pollutants and the higher their concentration theirs. The modeling demonstrated that the Piaçagueras channel has circulation with far less energy than other areas of the estuary, being remarkable the lowest spread of pollutants in this area, so that the low kinetic energy of currents in the channel is associated with a lower dispersion of particles. Considering the continuous discharge of pollutants and aggravating factors such as low hydrodynamics, bioaccumulation and resuspension of contaminated sediments, it is clear that this area needs greater rigidity in the control and disposal of materials with high pollution potential.

Submitted by Jean Medeiros (jeanletras@uepb.edu.br) on 2016-08-18T11:50:05Z No. of bitstreams: 1 PDF - André Luiz da Silva.pdf: 2669742 bytes, checksum: 796dc6a0d580f489f471969a901b6a51 (MD5) / Made available in DSpace on 2016-08-18T11:50:05Z (GMT). No. of bitstreams: 1 PDF - André Luiz da Silva.pdf: 2669742 bytes, checksum: 796dc6a0d580f489f471969a901b6a51 (MD5) Previous issue date: 2016-04-28 / Currently there is a great concern for groundwater contamination due to environmental pollution. Because of the need to deal with problems of this nature, it is necessary to develop research to enable understanding of the spread of pollutants in porous media motivating the need to develop mathematical models and computer codes to simulate reactive transport multispecies pollutants. This work aimed to propose a mathematical model for the numerical simulation of the spread of pollutants in the saturated region of the soil, considering the advection mechanisms, diffusion, sorption and chemical reaction. The mathematical model obtained was discretized using the method of finite volume. The results show that the method of Finite Volume, presented a very satisfactory performance in the discretization of pollutant reactive two-dimensional transport problem of the conservation equations in groundwater, the numerical method was compared to a semi-infinite medium getting good results, simulated experimental design guided the work results. It also highlighted the good performance of the computer language as a support tool for the evaluation of environmental impacts caused by pollution sources, and you can view the extent and behavior of a plume of contamination and carry out numerous transactions in short time at low cost and producing results very close to reality. / Atualmente existe uma grande preocupação com a contaminação das águas subterrâneas devido a poluição ambiental. Devido a necessidade de se lidar com problemas desta natureza, se faz necessário o desenvolvimento de pesquisas que possibilitem o entendimento da propagação de poluentes em meios porosos motivando a necessidade do desenvolvimento de modelos matemáticos e códigos computacionais para simularem o transporte reativo de múltiplas espécies poluentes. Este trabalho teve como objetivo propor um modelo matemático para a simulação numérica da propagação de poluentes na região saturada do solo, contemplando os mecanismos de advecção, difusão, sorção e reação química. O modelo matemático obtido foi discretizado utilizando o método dos Volumes Finitos. Os resultados obtidos demonstram que o método dos Volumes Finitos, apresentou um desempenho bastante satisfatório na Discretização das equações de conservação do problema do transporte bidimensional reativo de poluentes em águas subterrâneas, o método numérico foi comparado para um meio semi-infinito obtendo bons resultados, o planejamento experimental simulado norteou os resultados do trabalho. Ficou evidenciado também a boa performance da linguagem computacional como ferramenta de apoio para a avaliação de impactos ambientais causados por fontes poluidoras, sendo possível visualizar a extensão e o comportamento de uma pluma de contaminação e realizar inúmeras operações em curto período de tempo, a baixo custo e produzindo resultados bem próximos da realidade.