An efficient method for simulating steady unsaturated flow in random porous media: Using an analytical perturbation solution as initial guess to a numerical model

Thomas Harter, T. ‐C Jim Yeh

Research output: Contribution to journalArticle

16 Scopus citations

Abstract

Numerical simulation of flow through multidimensional heterogeneous soils under unsaturated conditions is a computationally intensive task. The governing unsaturated flow equation is nonlinear. The degree of nonlinearity depends on the unsaturated hydraulic properties of the soil and the degree of heterogeneity. In this paper, we present a highly efficient method to simulate unsaturated steady state flow through random porous media. Our method utilizes the analytical (approximate) solution derived by a perturbation‐spectral method as an initial guess solution for a numerical model to simulate two‐dimensional vertical infiltration problems. It is found that this approach, which we call “ASIGNing,” reduces the required CPU time by one to two orders of magnitude. ASIGNing is demonstrated to operate successfully under a wide variety of boundary conditions which may substantially deviate from those imposed on the initial guess solution, A large range of mean and variances in the independent variables ln Ks and α or alternatively In α, has been tested and it is shown that the method works well for variances of the unsaturated hydraulic conductivity σln2K′ ≤ 5 and average 〈α〉 ≤ 0.1 [cm−1].

Original languageEnglish (US)
Pages (from-to)4139-4149
Number of pages11
JournalWater Resources Research
Volume29
Issue number12
DOIs
StatePublished - Dec 1993
Externally publishedYes

ASJC Scopus subject areas

  • Water Science and Technology

Fingerprint Dive into the research topics of 'An efficient method for simulating steady unsaturated flow in random porous media: Using an analytical perturbation solution as initial guess to a numerical model'. Together they form a unique fingerprint.

  • Cite this