Controle Adaptativo do Erro Numérico das Estimativas de Contaminantes em Aquíferos

Abstract

Vários modelos computacionais que implementam o transporte de soluto em meio poroso saturado surgem constantemente em publicações científicas devido à importância na compreensão e previsão do transporte de constituintes dissolvidos em água subterrânea. As soluções numéricas obtidas por esquemas computacionais não estão imunes aos erros de discretização. No entanto, a confiabilidade nos resultados obtidos das complexas operações provenientes da dinâmica de fluidos computacional pode ser aumentada através de estimadores de erro a posteriori que indicam a precisão da solução numérica de um modelo matemático para simular o fenômeno físico de interesse. Neste trabalho é avaliado um estimador residual para a equação parabólica que descreve os fenômenos de advecção-dispersão-reação (ADR) em meio poroso saturado, considerando o transporte em regime de pequena advecção. Para a solução numérica da equação ADR é implementada um código computacional, em linguagem PYTHON, utilizando a discretização espaço-temporal que emprega o esquema -estável. Em situações de aquíferos isotrópicos e homogêneos, a boa concordância entre a solução numérica e a solução analítica disponível na literatura, é identificada também pelo estimador de erro residual.