Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions. (English) Zbl 1352.76110

Summary: The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.


76S05 Flows in porous media; filtration; seepage
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
74L10 Soil and rock mechanics


Full Text: DOI Link


[1] Szymkiewicz, Unified macroscopic model for unsaturated water flow in soils of bimodal porosity, Hydrological Sciences Journal 51 (6) pp 1106– (2006)
[2] Lunati, Multiscale finite-volume method for compressible multiphase flow in porous media, Journal of Computational Physics 216 pp 616– (2006) · Zbl 1220.76049
[3] Hou, A multiscale finite element method for elliptic problems in composite materials and porous media, Journal of Computational Physics 134 pp 169– (1997) · Zbl 0880.73065
[4] Ryan, A hybrid micro-scale model for transport in connected macro-pores in porous media, Journal of Contaminant Hydrology 126 pp 61– (2011)
[5] Auriault, Heterogeneous medium: is an equivalent macroscopic description possible?, International Journal of Engineering Science 29 (7) pp 785– (1991) · Zbl 0749.73003
[6] Lewandowska, Homogenization modeling and parametric study of moisture transfer in an unsaturated heterogeneous porous medium, Transport in Porous Media 45 pp 321– (2001)
[7] Arbogast, Gravitational forces in dual-porosity systems: I. model derivation by homogenization, Transport in Porous Media 13 pp 179– (1993)
[8] Arbogast, Derivation of the double porosity model of single phase flow via homogenization theory, SIAM Journal on Mathematical Analysis 21 (4) pp 823– (1990) · Zbl 0698.76106
[9] Clark, Fluid flow in a layered medium, Quarterly of Applied Mathematics 52 (4) pp 777– (1994) · Zbl 0815.35038
[10] Cook, Microstructure diffusion models with secondary flux, Journal of Mathematical Analysis and Applications 189 pp 731– (1995) · Zbl 0831.35138
[11] Showalter, Diffusion of fluid in a fissured medium with microstructure, SIAM Journal on Mathematical Analysis 22 (6) pp 1702– (1991) · Zbl 0764.35053
[12] Pesyńska, Multiscale elliptic-parabolic systems for flow and transport, Electronic Journal of Differential Equations 2007 (147) pp 1– (2007)
[13] Lewandowska, Modeling of unsaturated water flow in double-porosity soils by the homogenization approach, Advances in Water Resources 27 pp 283– (2004)
[14] Szymkiewicz, Micromechanical approach to unsaturated water flow in structured geomaterials by two-scale computations, Acta Geotechnica 3 pp 37– (2008)
[15] Szymkiewicz, Two-scale modeling of unsaturated water flow in a double-porosity medium under axisymmetric conditions, Canadian Geotechnical Journal 45 pp 238– (2008)
[16] Szymkiewicz A Modeling of unsaturated water flow in highly heterogeneous soils 2004
[17] Pruess K Narasimhan TN A practical method for modeling fluid and heat flow in fractured porous media Sixth SPE Symposium on Reservoir Simulation 1982 1 37
[18] Peter, Homogenization of secondary-flux models of partially fissured media, International Journal of Numerical Analysis and Modeling 5 pp 150– (2008) · Zbl 1165.76051
[19] Auriault, Homogenization of Coupled Phenomena in Heterogenous Media (2009)
[20] Szymkiewicz, Calculating effective conductivity of heterogeneous soils by homogenization, Archives of Hydro-Engineering and Enviromental Mechanics 52 (2) pp 111– (2005)
[21] Genuchten, A closed-form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Science Society of America Journal 44 (5) pp 892– (1980)
[22] Mualem, A new model for predicting the hydraulic conductivity of unsaturated porous media, Water Resources Research 12 (3) pp 513– (1976)
[23] Durner, Proceedings of the International Workshop on Characterization and Measurement of the Hydraulic Properties of Unsaturated Porous Media pp 817– (1997)
[24] Dhatt, Finite Element Method (2012)
[25] Hochbruck, Exponential Rosenbrock-type methods, SIAM Journal on Numerical Analysis 47 (1) pp 786– (2009) · Zbl 1193.65119
[26] Carr, Efficient simulation of unsaturated flow using exponential time integration, Applied Mathematics and Computation 217 (14) pp 6587– (2011) · Zbl 1419.76494
[27] Carr, A dual-scale modelling approach for drying hygroscopic porous media, SIAM Multiscale Modeling and Simulation 11 (1) pp 362– (2013) · Zbl 1267.76106
[28] Carr, A variable-stepsize Jacobian-free exponential integrator for simulating transport in heterogeneous porous media: application to wood drying, Journal of Computational Physics 233 pp 66– (2013) · Zbl 1286.76110
[29] Kelley, Solving Nonlinear Equations with Newton’s Method (2003) · Zbl 1031.65069
[30] Guezaine, Gmsh: A 3-d finite element mesh generator with built-in pre- and post-processing facilities, International Journal for Numerical Methods in Engineering 79 (11) pp 1309– (2009) · Zbl 1176.74181
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.