Examples of numerical simulations of two-dimensional unsaturated flow with VS2DI code using different interblock conductivity averaging schemes


numerical modelling
Richards equation
hydraulic conductivity

How to Cite

Szymkiewicz, A., Tisler, W., & Burzyński, K. (2015). Examples of numerical simulations of two-dimensional unsaturated flow with VS2DI code using different interblock conductivity averaging schemes. Geologos, 21(3), 161–167. https://doi.org/10.1515/logos-2015-0015


Flow in unsaturated porous media is commonly described by the Richards equation. This equation is strongly nonlinear due to interrelationships between water pressure head (negative in unsaturated conditions), water content and hydraulic conductivity. The accuracy of numerical solution of the Richards equation often depends on the method used to estimate average hydraulic conductivity between neighbouring nodes or cells of the numerical grid. The present paper discusses application of the computer simulation code VS2DI to three test problems concerning infiltration into an initially dry medium, using various methods for inter-cell conductivity calculation (arithmetic mean, geometric mean and upstream weighting). It is shown that the influence of the averaging method can be very large for coarse grid, but that it diminishes as cell size decreases. Overall, the arithmetic average produced the most reliable results for coarse grids. Moreover, the difference between results obtained with various methods is a convenient indicator of the adequacy of grid refinement.



Baker, D., 2000. A Darcian integral approximation to interblock hydraulic conductivity means in vertical infiltration. Computers and Geosciences 26(5):581–590.

Belfort, B. & Lehmann, F., 2005. Comparisonof equivalent conductivities for numerical simulation of one-dimensional unsaturated flow. Vadose Zone Journal 4(4):1191–1200.

Brunner, P., Cook, P.G. & Simmons, C.T., 2009. Hydrogeologic controls on disconnection between surface water and groundwater. Water Resources Research 45: W01422.

Fayer, M.J., 2000. UNSAT-H version 3.0: Unsaturated soil water and heat flow model, theory, user manual, and examples. PNNL 13249. Pacific Northwest National Laboratory, Richland, WA.

Flemisch, B., Fritz, J., Helmig, R., Niessner, J. & Wohlmuth, B., 2007. DUMUX: a multi-scale multi-physics toolbox for flow and transport processes in porous media. [In:] ECOMAS Thematic Conference on Multi-Scale Computational Methods for Solids and Fluids, Paris-Cachan, France, pp. 69–74.

Gasto, J., Grifoll, J. & Cohen, Y., 2002. Estimation of internodal permeabilities for numerical simulations of unsaturated flow. Water Resources Research 38: 1326.

Geo-Slope International, 2008. Vadose Zone Modeling with VADOSE/W 2007.

Haverkamp, R. & Vauclin, M., 1979. A note on estimating finite difference interblock conductivity values for transient unsaturated flow. Water Resources Research 15, 181–187.

Healy, R.W., 1990. Simulation of solute transport in variably saturated porous media with supplemental information on modifications to the U.S. Geological Survey’s computer program VS2D. USGS Water-Resources Investigation Report 90-4025. Reston, USA.

Healy, R.W. & Ronan, A.D., 1996. Documentation of computer program VS2DH for simulation of energy transport in variably saturated porous media; modification of the U. S. Geological Survey’s computer program VS2DT. USGS Water-Resources Investigation Report 96-4230. Reston, USA.

Healy, R.W., 2008. Simulating water, solute, and heat transport in the subsurface with the VS2DI software package. Vadose Zone Journal 7, 632–639.

Healy, R.W. & Essaid, H.I., 2012. VS2DI: Model use, calibration and validation. Transactions of the American Society of Agricultural and Biological Engineers (ASABE) 55, 1249–1260.

Hsieh, P.A., Wingle, A.W. & Healy, R.W., 1999. VS2DI: A graphical software package for simulating fluid flow and solute or energy transport in variably saturated porous media. USGS Water-Resources Investigation Report 99-4130, Reston, USA.

Jackson, C., 2005. Modelling leakage from perched rivers using the unsaturated flow model VS2DTI. British Geological Survey, Groundwater Systems & Water Quality Programme, Internal Report IR/05/019.

Kroes, J.G., van Dam, J.C., Groenendijk, P., Hendriks, R.F.A. & Jacobs, C.M.J., 2008. SWAP version 3.2. Theory description and user manual. Alterra-report 1649, 262 pp, Alterra, Research Institute, Wageningen.

Lappala, E.G., Healy, R.W. & Weeks, E.P., 1987. Documentation of computer program VS2D to solve the equations of fluid flow in variably saturated porous media. USGS Water-Resources Investigation Report 83-4099. U.S. Geological Survey, Reston.

Lin, H.-C.J., Richards, D.R., Yeh, G.-T., Cheng, J.-R., Cheng, H.-P. & Jones, N.L., 1997. FEMWATER: A three-dimensional finite element computer model for simulating density-dependent flow and transport in variably saturated media. Technical Report CHL-97-12, U.S. Army Corps of Engineers, Waterways Experiment Station.

Pruess, K., Oldenburg, C. & Moridis, G., 1999. TOUGH2 user’s guide. Version 2.0. Report LBNL-43134, Lawrence Berkeley National Laboratory, Berkeley, USA.

Richards, L.A., 1931. Capillary conduction of liquids through porous mediums. Journal of Applied Physics 1, 318–333.

Šimůnek, J., van Genuchten, M.Th. & Sejna, M., 2008. Development and applications of the HYDRUS and STANMOD software packages and related codes. Va­dose Zone Journal 7, 587–600.

Szymkiewicz, A., 2009. Approximation of internodal conductivities in numerical simulation of 1D infiltration, drainage and capillary rise in unsaturated soils, Water Resources Research 45, W10403.

Szymkiewicz, A., 2013. Modeling water flow in unsaturated porous media: Accounting for nonlinear permeability and material heterogeneity. Springer, 237 pp.

Szymkiewicz, A. & Burzyński, K., 2011. Computing internodal conductivities in numerical modeling of two dimensional unsaturated flow on rectangular grid. Archives of Civil Engineering 57, 215–225.

Szymkiewicz, A. & Helmig, R., 2011. Comparison of conductivity averaging methods for one-dimensional unsaturated flow in layered soils. Advances in Water Resources 34, 1012–1025.

Van Genuchten, M.Th., 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal 44, 892–898.

Voss, C.I. & Provost, A.M., 2002. SUTRA, a model for saturated–unsaturated, variable-density ground-water flow with solute or energy transport. Water-Resources Investigation Report 2002-4231. U.S. Geological Survey, Reston, USA.