
Unsaturated Flow Through Spherical Inclusions with Contrasting Sorptive Numbers
Author(s) -
Furman Alex,
Warrick A. W.
Publication year - 2005
Publication title -
vadose zone journal
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.036
H-Index - 81
ISSN - 1539-1663
DOI - 10.2136/vzj2004.0076
Subject(s) - hydraulic conductivity , inclusion (mineral) , nonlinear system , flow (mathematics) , vadose zone , mechanics , mathematics , richards equation , boundary value problem , thermodynamics , materials science , mathematical analysis , geometry , geotechnical engineering , geology , water content , physics , groundwater , soil science , soil water , quantum mechanics
The analytic element method (AEM) is used to model unsaturated flow through a spherical inclusion of contrasting hydraulic properties. The steady state Richards' equation is combined with the Gardner model for unsaturated hydraulic conductivity to form the Helmholtz equation. The later is solved by means of the AEM. The background and inclusion materials are assumed to have different saturated hydraulic conductivities and different sorptive numbers; hence, the conditions are more general than treatments of spherical inclusions. Continuity of the interfacial head boundary condition leads to a nonlinear system of equations, whose solution requires an iterative solution. Analysis includes the effect of the hydraulic properties and of the background flux and evaluation of computational efficiency for contrasting hydraulic properties. To examine water contents, a methodology is presented for matching Gardner and van Genuchten parameters. The new solution is more realistic than the previous solution for a spherical inclusion with a sorptive number the same as the background but is computationally more tedious.