Resumo A simulação numérica do escoamento deágua no soloé de grande importância para entender e melhorar muitos processos de drenagem e infiltração. Neste trabalho, a abordagem numérica para a simulação do escoamentoágua-ar em uma zona subsaturada foi focada na resolução numérica da equação de Richards. Para resolver o sistema não linear decorrente da discretização da equação que modela o problema, foi utilizado um método numérico de busca direta que não utiliza derivadas, denominado DFSANE. Palavras-chave. Equação de Richards, Meios Porosos, DFSANE, Zona Subsaturada. 1 Introdução O escoamento de fluidos em meios porososé um complexo fenômeno que vem sendo analisado há anos [3], sendo de grande importância para diversasáreas da engenharia e ciência, por exemplo, hidrologia deáguas subterrâneas, engenharia de reservatórios, ciência dos solos, mecânica dos solos e engenharia química [1].