An integro-differential model for non-Fickian tracer transport in porous media: validation and numerical simulation. (English) Zbl 1398.76216

Summary: Diffusion processes have traditionally been modeled using the classical parabolic advection-diffusion equation. However, as in the case of tracer transport in porous media, significant discrepancies between experimental results and numerical simulations have been reported in the literature. Therefore, in order to describe such anomalous behavior, known as non-Fickian diffusion, some authors have replaced the parabolic model with the continuous time random walk model, which has been very effective. Integro-differential models (IDMs) have been also proposed to describe non-Fickian diffusion in porous media. In this paper, we introduce and test a particular type of IDM by fitting breakthrough curves resulting from laboratory tracer transport. Comparisons with the traditional advection-diffusion equation and the continuous time random walk are also presented. Moreover, we propose and numerically analyze a stable and accurate numerical procedure for the two-dimensional IDM composed by a integro-differential equation for the concentration and Darcy’s law for flow. In space, it is based on the combination of mixed finite element and finite volume methods over an unstructured triangular mesh.


76S05 Flows in porous media; filtration; seepage
76M10 Finite element methods applied to problems in fluid mechanics
65R20 Numerical methods for integral equations
45K05 Integro-partial differential equations
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
65N30 Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs


Full Text: DOI Link


[1] Douglas, Superconvergence in the pressure in the simulation of miscible displacement, SIAM Journal on Numerical Analysis 22 pp 962– (1985) · Zbl 0624.65124
[2] Ewing, Galerkin method for miscible displacement problems in porous media, SIAM Journal on Numerical Analysis 17 pp 351– (1980) · Zbl 0458.76092
[3] Peaceman, Improved treatment of dispersion in numerical calculation of multidimensional miscible displacement, Society of Petroleum Engineers Journal 6 pp 213– (1996)
[4] Ewing, The Mathematics of Reservoir Simulation (1984) · Zbl 0533.00031
[5] Settari, Development and application of variational methods for simulation of miscible displacement in porous media, Society of Petroleum Engineers Journal 17 pp 228– (1977)
[6] Chen, Analysis of non-Fickian diffusion problems in a composite medium, Computer Physics Communications 150 pp 31– (2003) · Zbl 1196.65161
[7] Hassanizadeh, On the transient non-Fickian dispersion theory, Transport in Porous Media 23 pp 107– (1996)
[8] Neuman, Perspective on theories of non-Fickian transport in heterogeneous media, Advances in Water Resources 32 pp 670– (2009)
[9] Bromly, Non-Fickian transport in homogeneous unsaturated repacked sand, Advances in Water Resources 40 pp W07402– (2004)
[10] Cortis, Anomalous transport in ’classical’ soil and sand columns, Soil Science Society of America Journal 68 pp 1539– (2004)
[11] Cortis, Quantitative characterization of pore-scale disorder effects on transport in ’homogeneous’ granular media, Physical Review E 70 pp 041108– (2004)
[12] Gelhar, A critical review of data on field-scale dispersion in aquifers, Water Resources Research 28 pp 1955– (1992)
[13] Huang, Experimental investigation of solute transport in large, homogeneous and heterogeneous, saturated soil columns, Transport in Porous Media 18 pp 283– (1995)
[14] Levy, Measurement and analysis of non-Fickian dispersion in heterogeneous porous media, Journal of Contaminant Hydrology 64 pp 203– (2003)
[15] Mackay, A natural gradient experiment on solute transport in a sand aquifer: 1. approach and overview of plume movement, Water Resources Research 22 pp 2017– (1986)
[16] Sidle, Spatially varying hydraulic and solute transport characteristics of a fractured till determined by field tracer tests, Funen, Denmark, Water Resources Research 34 pp 2515– (1998)
[17] Silliman, Laboratory investigation of longitudinal dispersion in anisotropic porous media, Water Resources Research 23 pp 2145– (1987)
[18] Sternberg, Laboratory observation of nonlocal dispersion, Transport in Porous Media 23 pp 135– (1996)
[19] Strack, A mathematical model for dispersion with a moving front in groundwater, Water Resources Research 28 pp 2973– (1992)
[20] Barbeiro, H1-second order convergent estimates for non-Fickian models, Applied Numerical Mathematics 61 pp 201– (2011) · Zbl 1209.65141
[21] Cannon, A priori L2 error estimates for finite-element methods for nonlinear diffusion equations with memory, SIAM Journal on Numerical Analysis 27 pp 595– (1990) · Zbl 0709.65122
[22] Ewing, Finite volume element approximations of nonlocal in time one-dimensional flows in porous media, Computing 64 pp 157– (2000) · Zbl 0969.76052
[23] Ewing, Sharp L2-error estimates and superconvergence of mixed finite element methods for non-Fickian flows in porous media, SIAM Journal on Numerical Analysis 40 pp 1538– (2002) · Zbl 1044.76032
[24] Ferreira, A second order approximation for quasilinear non-Fickian models, Computer Methods in Applied Mechanics 13 pp 471– (2013)
[25] Ferreira, Supraconvergence and supercloseness in Volterra equations, Applied Numerical Mathematics 62 pp 1718– (2012) · Zbl 1254.65131
[26] Sinha, Some new error estimates of a semidiscrete finite volume element method for a parabolic integro-differential equation with nonsmooth initial data, SIAM Journal on Numerical Analysis 43 pp 2320– (2006) · Zbl 1106.65121
[27] Ferreira JA Pinto L Non-Fickian tracer transport in porous media Proceedings of the 14th International Conference on Computational and Mathematical Methods in Science and Engineering Rota (Cádiz), Spain 2014
[28] Dentz, Time behavior of solute transport in heterogeneous media: transition from anomalous to normal transport, Advances in Water Resources 27 pp 155– (2004)
[29] Gao, Comparison of alternative models for simulating anomalous solute transport in a large heterogeneous soil column, Journal of Hydrology 377 pp 391– (2009)
[30] Cortis, Computing ’anomalous’ contaminant transport in porous media: the CTRW MATLAB toolbox, Ground Water 43 pp 947– (2005)
[31] Scheidegger AE An evaluation of the accuracy of the diffusivity equation for describing miscible displacement in porous media Proceedings of Theory of Fluid Flow in Porous Media Norman, Oklahoma 1959
[32] Huang, Evidence of one-dimensional scale-dependent fractional advection-dispersion, Journal of Contaminant Hydrology 85 pp 53– (2006)
[33] Logan, Solute transport in porous media with scale-dependent dispersion and periodic boundary conditions, Journal of Hydrology 184 pp 261– (1996)
[34] Mishra, Analysis of solute transport with a hyperbolic scale-dependent dispersion model, Hydrological Processes 4 pp 45– (1990)
[35] Zoua, Analytical solutions to non-Fickian subsurface dispersion in uniform groundwater flow, Journal of Hydrology 179 pp 237– (1996)
[36] Arbogast, Computational methods for multiphase flow and reactive transport problems arising in subsurface contaminant remediation, Journal of Computational and Applied Mathematics 74 pp 19– (1996) · Zbl 0866.76041
[37] Arbogast, A fully mass and volume conserving implementation of a characteristic method for transport problems, SIAM Journal of Scientific Computing 28 pp 2001– (2066) · Zbl 1293.76097
[38] Bergamaschi, A mixed finite element-finite volume formulation of the black-oil model, SIAM Journal of Scientific Computing 20 pp 970– (1998) · Zbl 0959.76039
[39] Lipnikov, Local flux mimetic finite difference methods, Numerische Mathematik 112 pp 115– (2009) · Zbl 1165.65063
[40] Brezzi, Mixed and Hybrid Finite Element Methods (1991) · Zbl 0788.73002
[41] Raviart, A mixed finite element method for 2nd order elliptic problems, Mathematical Aspects of Finite Element Methods, Lecture Notes in Mathematics 606 pp 292– (1977) · Zbl 0362.65089
[42] Crouzeix, Une méthode multipas implicite-explicite pour l’approximation des équations d’évolution paraboliques, Numerische Mathematik 35 pp 257– (1980) · Zbl 0419.65057
[43] Barth, The design and application of upwind schemes on unstructured meshes, AIAA Paper 93-0668 (1993)
[44] Dawson, Godunov mixed methods for immiscible displacement, International Journal for Numerical Methods in Fluids 11 pp 835– (1990) · Zbl 0704.76059
[45] Hundsdorfer, Monotonicity-preserving linear multistep methods, SIAM Journal on Numerical Analysis 41 pp 605– (2003) · Zbl 1050.65070
[46] Ewing R Lazarov R Pasciak J Vassilev A Mathematical modeling, numerical techniques, and computer simulation of flows and transport in porous media Proceedings of Computational Techniques and Applications Melbourne 1995
[47] Bramble, Analysis of the inexact Uzawa algorithm for saddle point problems, SIAM Journal on Numerical Analysis 34 pp 1072– (1997) · Zbl 0873.65031
[48] Brenner, A multigrid algorithm for the lowest-order Raviart-Thomas mixed triangular finite element method, SIAM Journal on Numerical Analysis 29 pp 647– (1992) · Zbl 0759.65080
[49] Arbogast, Enhanced cell-centered finite differences for elliptic equations on general geometry, SIAM Journal of Scientific Computing 19 pp 404– (1998) · Zbl 0947.65114
[50] Arbogast, Efficient mixed methods for groundwater flow on triangular or tetrahedral meshes, Water Science and Technology Library 1 pp 3– (1994)
[51] Acosta, The maximum angle condition for mixed and nonconforming elements: application to the Stokes equations, SIAM Journal on Numerical Analysis 37 pp 18– (2000) · Zbl 0948.65115
[52] Douglas, Global estimates for mixed methods for second order elliptic equations, Mathematics of Computation 44 pp 39– (1985) · Zbl 0624.65109
[53] Gastaldi, Optimal L-error estimates for nonconforming and mixed finite element methods of lowest order, Numerische Mathematik 50 pp 587– (1987) · Zbl 0597.65080
[54] Hubbard, Multidimensional slope limiters for MUSCL-type finite volume schemes on unstructured grids, Journal of Computational Physics 155 pp 54– (1999) · Zbl 0934.65109
[55] Buffard, Monoslope and multislope MUSCL methods for unstructured meshes, Journal of Computational Physics 229 pp 3745– (2010) · Zbl 1189.65204
[56] Park, Multi-dimensional limiting process for hyperbolic conservation laws on unstructured grids, Journal of Computational Physics 229 pp 788– (2010) · Zbl 1185.65150
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.