Simulação de infiltração em meios porosos por diferença finita

• Autor(a) principal: Wendland, Edson Cezar
• Data de Publicação: 1991
• Tipo de documento: Dissertação
• Idioma: por
• Título da fonte: Biblioteca Digital de Teses e Dissertações da UFRGS
• Texto Completo: http://hdl.handle.net/10183/130092
• topic: Meios porosos; Diferenças finitas; Mecânica dos fluidos

O objetivo do presente trabalho é determinar numericamente o perfil de umedecimento em um meio poroso submetido a um processo de infiltração. A descrição do movimento transiente de água em solo não-expansivo é dada em termos da equação não-linear de Fokker-Planck ou da equação de Richards. As dificuldades na solução da equação são não somente devido aos termos não-lineares envolvidos mas também às condições de fronteira dependentes do tempo e ao ajuste dos dados experimentais para os parâmetros fisicos. Neste trabalho é considerado um esquema de diferença finita de três niveis para a simulação de infiltração em uma coluna homogênea, semi-infinita satisfazendo a condição de fluxo constante na superficie. O esquema proposto é implicito com erro de truncamento de segunda ordem no tempo e espaço. O termo convectivo é aproximado por quatro diferentes esquemas. A influência da discretização da derivada na fronteira também é analisada através de quatro esquemas. Os resultados obtidos na simulação numérica pelo esquema de três niveis com termo convectivo aproximado por diferença ascendente (upwind) e condição de fronteira discretizada com ponto ficticio apresentam melhor aproximação aos perfis de umedecimento experimentais que os obtidos por esquemas de dois niveis.