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

Resumo

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.