O presente trabalho trata de construir um modelo computacional utilizando o método dos volumes finitos para malhas não-estruturadas, a fim de se calcular a carga hidráulica num meio poroso, considerando este meio não homogêneo e anisotrópico. A anisotropia é uma característica de muitos materiais encontrados na natureza e depende da propriedade estudada no meio. Primeiramente apresenta-se a dedução da equação do transporte advectivo dispersivo e a formulação matemática para a equação de Laplace, esta última utilizada para o cálculo da carga hidráulica. Em seguida, apresenta-se o algoritmo de solução de um programa computacional em linguagem C++ que permite calcular a velocidade do fluxo em cada face de um volume de controle. Finalmente são feitos vários testes para validação do código computacional utilizado, o que levou a crer que o método utilizado é eficaz para os tipos de malhas testados, apresentando algumas diferenças quanto ao erro da solução. / The present work build a computational model using the finite volumes method for unstructured meshes, with the purpose of calculating the hydraulic load in a porous medium, considering it material non - homogeneous and anisotropic. The Anisotropy is a characteristic of many materials found in the nature and it depends on the property studied in this material. First, we present the deduction of the equation of advective-dispersive transport and the mathematical formulation for the Laplaces equation, this last one used for the calculation of the hydraulic load. Soon afterwards, we present the solution algorithm of a computational program in the C++ language that allows to calculate the speed of the flow in each face of the control volume. Finally several tests for validation of the code are made, which makes it that the plausible to assume method is effective for the types of meshes tested, presenting some differences for the wrong solution.

No meio científico são bastante utilizados os modelos matemáticos para avaliar as características hidrodinâmicas de reatores, porém a literatura é pobre em informações relativas à aplicação do equacionamento matemático em regime não-permanente. Neste trabalho foi aplicado um modelo matemático simulando a hidrodinâmica de reatores tubulares com dispersão para avaliar o efeito da variação de vazão afluente sobre os parâmetros hidrodinâmicos. As simulações foram efetuadas considerando vazão e volume constantes, vazão e volume variáveis e vazão variável e volume constante. Foi investigada a influência de dois tipos de ensaios estímulo-resposta, pulso e degrau, para aplicação de modelos matemáticos e determinação das curvas de distribuição do tempo de residência (DTR) experimentais. Teoricamente ambos devem fornecer os mesmos resultados embora o ensaio em pulso costume apresentar maior sensibilidade experimental. Conforme esperado, ambos os ensaios apresentaram os mesmos resultados finais. Também foram avaliados os traçadores empregados nos dois tipos de ensaios estímulo-resposta, pulso e degrau. Foram empregados três traçadores diferentes - verde de bromocresol, azul de bromofenol e eosina Y que proporcionaram curvas com diferentes características. Como ferramenta auxiliar da modelação foram estudas a determinação das curvas DTR experimentais com auxílio de duas técnicas distintas. Para calibrar o modelo matemático proposto foram utilizados dados de reatores em escala de bancada com diferentes configurações submetidas a variações de vazões afluentes. Os dados de um reator UASB em escala piloto - submetido a variações cíclicas de vazão afluente de 40 e 60% - foram empregados para calibrar e verificar o modelo matemático proposto. Os resultados encontrados com o modelo matemático proposto nesta pesquisa demonstraram que a variação de vazão afluente não deve ser negligenciada. O modelo utilizado representou adequadamente o reator UASB. Seus resultados, quando comparados aos modelos matemáticos que não consideram a variação de vazão, mostraram que para flutuações de vazão elevadas, vazões com valores até 60% maiores do que a vazão média, os valores dos coeficientes de difusão diferem significativamente em função das hipóteses empregadas no desenvolvimento do modelo matemático. / Mathematical models for hydrodynamic characteristics evaluation of reactors are commonly used however there is a lack of information in the literature concerning the application of mathematical modeling for non-steady state flow. In this thesis, it is presented a non-steady mathematical model to simulate the hydrodynamic behavior of pipe flow reactors with diffusion in order to evaluate the effect of the variation of the influent flowrate in the hydrodynamic parameters. The simulations were performed considering constant flowrate and volume, variable flowrate and volume and, variable flowrate and constant volume. It was investigated the influence of two types of stimulus-response assays, pulse and step function, on the application of mathematical models and the determination of the experimental retention time distribution curves. Theoretically, both tests should provide the same results although the pulse stimulus assay usually presents higher experimental sensitivity. As it was expected, both tests presented the same final results. The tested tracers were also evaluated in relation to the pulse and step stimulus-response tests. Three tracers were tested (bromocresol green, bromophenol blue and eosin Y) and it was shown that they provided different curves. It was studied the construction of the experimental retention time distribution curves using two procedures as an auxiliary tool for modeling. The proposed mathematical model was calibrated with data of bench scale reactors submitted to different cyclical variations of flowrates. Moreover, data of UASB reactor in pilot scale - submitted to 40 and 60% of cyclical variations of flowrates were utilized to calibrate and verify the obtained mathematical model. The results found with the mathematical model proposed in this research showed that the variations in influent flow rate can not be neglected. The model adequately represented a pilot scale UASB reactor. The results showed that the diffusion coefficients differ significantly for high flowrate fluctuations, when compared to other mathematical models that do not incorporate variable flowrate, depending upon the hypothesis used to derive the model.

Diploma thesis is aimed on description of wastewater treatment plant with gravitational inflow. Focus is on WWTP in the Czech Republic, particularly in the South Moravia region - in this aspect the thesis presents relevant statistics like the number of WWTP per region and the total, capacity of WWTP per region and their substance load on inflow /outflow. The recherche is also limited according to the core subject of theoretical part, which is technology of small WWTP in the category up to 2000 equivalent habitants. The experimental part presents comparsion of particular wastewater treatment plant of this category in the South Moravia region: one with gravitational inflow and one with inflow pump station. The examined wastewater treatment plant are compared according to substance load, hydraulic load and technology. Thesis is concluded with a review of the different operation parameters of the both types of WWTP.