Modelo hidrodinâmico e de transporte bidimensional de grade não estruturada para lagos rasos

Abstract

Devido à proximidade de grandes cidades, ecossistemas aquáticos são alvos de despejos industriais, esgotos domésticos e resíduos sólidos. Modelos matemáticos são largamente utilizados para entender os padrões de fluxo e o transporte de substâncias nestes ecossistemas. Neste trabalho, desenvolveu-se um modelo hidrodinâmico e de transporte de poluentes bidimensional para corpos d’água rasos que utiliza grades não estruturadas. O modelo foi desenvolvido em dois módulos: hidrodinâmico e de transporte de massa. O módulo hidrodinâmico utiliza as equações completas de Saint Venant para quantificação da circulação em corpos d’águas. No módulo hidrodinâmico aqui desenvolvido é feito um aperfeiçoamento de um esquema semi-implícito de diferenças finitas utilizado no modelo IPH-ECO para um esquema semi-implícito de volumes finitos. O módulo de transporte de massa, por sua vez, representa o transporte de um escalar conservativo sob influência dos processos físicos do meio. Dois esquemas foram empregados e comparados para a solução numérica da equação de transporte: um esquema de diferenças centrais e um esquema com limitador de fluxo. O modelo foi aplicado em um estudo de caso simplificado formado por uma lagoa circular e no rio Guaíba. O módulo hidrodinâmico foi avaliado através de testes de conservação de volume feitos comparando com os resultados encontrados nno modelo IPHECO na aplicação realizada na lagoa circular e na solução analítica no rio Guaíba. O módulo hidrodinâmico do modelo de grades não estruturadas (IPH-UnTRIM2D) mostrou-se mais conservativo que o modelo IPH-ECO em todos os cenários simulados na lagoa circular. A solução hidrodinâmica e de transporte do modelo de grades não estruturadas no estudo de caso idealizado no rio Guaíba representou meandros, sulcos, estreitamentos e alargamentos de canais com comprometimento com uma situação real. Quanto à simulação de transporte de escalares e poluentes, foi comprovada a eficiência do esquema com limitador de fluxo sobre os esquemas de diferenças centrais. Os testes de conservação de massa para diferentes cenários no rio Guaíba indicaram erros pouco significativos em ambos os esquemas quando comparados com a ordem de grandeza do ecossistema.

Aquatics ecossistems around large towns are always exposed waste domestic and industrial and solid residues. Mathematical models are widely used to understand flow patterns and transport of mass in these ecosystems. In this work a model was developed on unstructured two-dimensional complex dynamic model for shallow lakes. The model was divided in two modules: (a) hydrodynamic module, describing quantitative flows and water level; (b) a mass transport module. The hydrodynamic module solves Saint Venant equations to quantify water circulation in a water body. Futhermore, the hydrodynamic module was improved to semiimplicit finite volume scheme. The mass transport module represents transport of a conservative scalar under the physical processes influence. Two schemes were used and compared to the transport equation numerical solution : a central difference scheme and a flux limiter scheme. The model was applied in the circular lake and river Guaíba. The hydrodynamic module was validated by evaluating tests of volume conservation made comparing with results of the model IPH-ECO and analytical solution. It proved to be more conservative than the model IPH-ECO in all scenarios simulated in the circular lake. The hydrodynamic and transport solution in river Guaíba study case represented bends, grooves, and channels widening and narrowing. In the transport simulation, the flux limiter scheme was more efficient than central differences scheme. Futhermore, mass conservation tests for different scenarios in Guaíba indicated minor errors in both schemes when compared with ecosystem magnitude.