Computing internodal conductivities in numerical modeling of two dimensional unsaturated flow on rectangular grid

This paper compares numerical solutions of transient two-dimensional unsaturated flow equation by using different averaging schemes for internodal conductivities. Averaging methods such as arithmetic mean, geometric mean, upstream weighting, and integrated mean are taken into account, as well as a recent approach based on steady-state approximation. The latter method proved the most flexible, producing relatively accurate solutions for both downward and upward flow cases.

W artykule porównano rozwiązania numeryczne dwuwymiarowego równania przepływu nienasyconego z użyciem różnych metod uśredniania współczynnika przewodności hydraulicznej między sąsiednimi węzłami. Uwzględniono średnią arytmetyczną, geometryczną, 'pod prąd' oraz całkową, a także niedawno zaproponowany sposób uśredniania oparty na analizie przepływu ustalonego. Ten ostatni sposób okazał się najbardziej uniwersalny, gdyż umożliwiał otrzymanie stosunkowo dokładnych wyników dla przepływu wody w strefie nienasyconej zarówno w kierunku w dół jak i w górę.