Um Modelo Matemático para a Simulação do Transporte de Solutos em Solos Não-Saturados

Resumo

Uma formulação de elementos finitos semi-discretos é proposta para aproximar o sistema acoplado não linear de Equações Diferenciais Parciais, composto da Equação de Richards para o fluxo e da equação do transporte de Advecção-Dispersão-Sorção para a concentração, que modela o transporte de solutos no solo. A equação de Richards é apresentada em três formas [1]: forma h, devida a carga de pressão; forma θ, baseado na umidade e forma mista, que é a combinação das duas variáveis. Neste trabalho, a equação de Richards é resolvida na forma mista que tem vantagens em relação à forma h, pois garante as propriedades de conservação de massa. Em ambas equações a variável espacial é aproximada pelo método de elementos finitos e a integração do tempo é resolvido usando a famılia do Método Trapezoidal Generalizado.