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

• Author(s): Wendland, Edson Cezar
• Date: 2015
• Persistent ID: http://hdl.handle.net/10183/130092
• Origin: Oasisbr
• Subject(s): Meios porosos; Diferenças finitas; Mecânica dos fluidos
• Document Type: Master thesis
• Language: Portuguese

Description

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.

The purpose of the present work is to determine numerically the water content profile in a porous medium submitted to an infiltration process. The description of the transient water movement in nonswelling soil is given in terms of the non-linear Fokker-Planck equation or Richards equation. The difficulties in solving this equation are due not only the nonlinearities involved but also to time-dependent boundary conditions and to the fitting of the experimental data for the physical parameters. In this work, we consider a three-level finite difference scheme for the simulation of infiltration in an homogeneous semi-infinite column satisfying constant flux condition at the surface. This is an implicit scheme which is second-order accurate in time and space. The convective term is approximated by four different schemes. The influence of the discretization of the boundary derivative is analysed by four schemes. The results of the numerical simulation by the three-level scheme with convective term approximated by forward difference (upwind) discretizated with ficticious and point boundary exhibited condition a closer approximation to the experimental water content pro1ile than the ones obtained by two-level scheme.