zbMATH — the first resource for mathematics

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport. (English) Zbl 1453.86051
Summary: A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., \(^{36}\mathrm{Cl}^-)\) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., \(^{90}\mathrm{Sr}^{++})\) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.
86A60 Geological problems
76S05 Flows in porous media; filtration; seepage
86-08 Computational methods for problems pertaining to geophysics
CrunchFlow; PETSc
Full Text: DOI
[1] Grisak, G.; Pickens, J., An analytical solution for solute transport through fractured media with matrix diffusion, J. Hydrol., 52, 1-2, 47-57 (1981)
[2] Tang, D.; Frind, E.; Sudicky, EA, Contaminant transport in fractured porous media: analytical solution for a single fracture, Water Resour. Res., 17, 3, 555-564 (1981)
[3] Feenstra, S.; Cherry, J.; Sudicky, E.; Haq, Z., Matrix diffusion effects on contaminant migration from an injection well in fractured sandstone, Groundwater, 22, 3, 307-316 (1984)
[4] Nilson, R.; Lie, K., Double-porosity modelling of oscillatory gas motion and contaminant transport in a fractured porous medium, Int. J. Numer. Anal. Methods Geomech., 14, 8, 565-585 (1990)
[5] Blessent, D.; Jørgensen, PR; Therrien, R., Comparing discrete fracture and continuum models to predict contaminant transport in fractured porous media, Groundwater, 52, 1, 84-95 (2014)
[6] Robinson, NI; Sharp, JM Jr; Kreisel, I., Contaminant transport in sets of parallel finite fractures with fracture skins, J. Contam. Hydrol., 31, 1-2, 83-109 (1998)
[7] Steefel, CI; Lichtner, PC, Multicomponent reactive transport in discrete fractures: I. Controls on reaction front geometry, J. Hydrol., 209, 1-4, 186-199 (1998)
[8] Yoshida, H.; Yamamoto, K.; Yogo, S.; Murakami, Y., An analogue of matrix diffusion enhanced by biogenic redox reaction in fractured sedimentary rock, J. Geochem. Explor., 90, 1-2, 134-142 (2006)
[9] Weatherill, D.; Graf, T.; Simmons, CT; Cook, PG; Therrien, R.; Reynolds, DA, Discretizing the fracture-matrix interface to simulate solute transport, Groundwater, 46, 4, 606-615 (2008)
[10] Fairley, J.P.: Fracture/matrix interaction in a fracture of finite extent. Water Resour. Res. 46(8), (2010)
[11] Chambon, JC; Broholm, MM; Binning, PJ; Bjerg, PL, Modeling multi-component transport and enhanced anaerobic dechlorination processes in a single fracture-clay matrix system, J. Contam. Hydrol., 112, 1-4, 77-90 (2010)
[12] Chen, F.; Falta, RW; Murdoch, LC, Numerical analysis of contaminant removal from fractured rock during boiling, J. Contam. Hydrol., 134, 12-21 (2012)
[13] Liu, C-T; Yeh, H-D, Development of approximate solutions for contaminant transport through fractured media, Appl. Math. Model., 39, 2, 438-448 (2015) · Zbl 1432.74064
[14] Joshi, N.; Ojha, C.; Sharma, P.; Madramootoo, CA, Application of nonequilibrium fracture matrix model in simulating reactive contaminant transport through fractured porous media, Water Resour. Res., 51, 1, 390-408 (2015)
[15] Chen, K.; Zhan, H., A Green’s function method for two-dimensional reactive solute transport in a parallel fracture-matrix system, J. Contam. Hydrol., 213, 15-21 (2018)
[16] Doughty, C., Investigation of conceptual and numerical approaches for evaluating moisture, gas, chemical, and heat transport in fractured unsaturated rock, J. Contam. Hydrol., 38, 1-3, 69-106 (1999)
[17] Hu, Q.; Salve, R.; Stringfellow, WT; Wang, JS, Field tracer-transport tests in unsaturated fractured tuff, J. Contam. Hydrol., 51, 1-2, 1-12 (2001)
[18] Berkowitz, B., Characterizing flow and transport in fractured geological media: a review, Adv. Water Resour., 25, 8-12, 861-884 (2002)
[19] Hu, M.; Wang, Y.; Rutqvist, J., On continuous and discontinuous approaches for modeling groundwater flow in heterogeneous media using the numerical manifold method: model development and comparison, Adv. Water Resour., 80, 17-29 (2015)
[20] Hu, Q.; Kneafsey, TJ; Trautz, RC; Wang, JS, Tracer penetration into welded tuff matrix from flowing fractures, Vadose Zone J., 1, 1, 102-112 (2002)
[21] Wu, Y-S; Pan, L.; Pruess, K., A physically based approach for modeling multiphase fracture-matrix interaction in fractured porous media, Adv. Water Resour., 27, 9, 875-887 (2004)
[22] MacQuarrie, KT; Mayer, KU, Reactive transport modeling in fractured rock: a state-of-the-science review, Earth Sci. Rev., 72, 3-4, 189-227 (2005)
[23] Steefel, CI; Lichtner, PC, Diffusion and reaction in rock matrix bordering a hyperalkaline fluid-filled fracture, Geochim. Cosmochim. Acta, 58, 17, 3595-3612 (1994)
[24] Giambalvo, ER; Steefel, CI; Fisher, AT; Rosenberg, ND; Wheat, CG, Effect of fluid-sediment reaction on hydrothermal fluxes of major elements, eastern flank of the Juan de Fuca Ridge, Geochim. Cosmochim. Acta, 66, 10, 1739-1757 (2002)
[25] Appelo, CAJ; Wersin, P., Multicomponent diffusion modeling in clay systems with application to the diffusion of tritium, iodide, and sodium in opalinus clay, Environ. Sci. Technol., 41, 14, 5002-5007 (2007)
[26] Rolle, M.; Muniruzzaman, M.; Haberer, CM; Grathwohl, P., Coulombic effects in advection-dominated transport of electrolytes in porous media: multicomponent ionic dispersion, Geochim. Cosmochim. Acta, 120, 195-205 (2013)
[27] Tournassat, C.; Steefel, CI, Ionic transport in nano-porous clays with consideration of electrostatic effects, Rev. Mineral. Geochem., 80, 1, 287-329 (2015)
[28] Tournassat, C.; Steefel, CI, Reactive transport modeling of coupled processes in nanoporous media, Rev. Mineral. Geochem., 85, 1, 75-109 (2019)
[29] Tournassat, C.; Steefel, CI; Gimmi, T., Solving the Nernst-Planck equation in heterogeneous porous media with finite volume methods: averaging approaches at interfaces, Water Resour. Res., 56, 3, e2019WR026832 (2020)
[30] Tournassat, C., Steefel, C.I.: Modeling diffusion processes in the presence of a diffuse layer at charged mineral surfaces: a benchmark exercise. Comput. Geosci. 1-18 (2019)
[31] Wersin, P.; Gimmi, T.; Mazurek, M.; Alt-Epping, P.; Pękala, M.; Traber, D., Multicomponent diffusion in a 280 m thick argillaceous rock sequence, Appl. Geochem., 95, 110-123 (2018)
[32] Glaus, MA; Birgersson, M.; Karnland, O.; Van Loon, LR, Seeming steady-state uphill diffusion of 22Na+ in compacted montmorillonite, Environ. Sci. Technol., 47, 20, 11522-11527 (2013)
[33] Muniruzzaman, M.; Haberer, CM; Grathwohl, P.; Rolle, M., Multicomponent ionic dispersion during transport of electrolytes in heterogeneous porous media: experiments and model-based interpretation, Geochim. Cosmochim. Acta, 141, 656-669 (2014)
[34] Appelo, C.; Vinsot, A.; Mettler, S.; Wechner, S., Obtaining the porewater composition of a clay rock by modeling the in-and out-diffusion of anions and cations from an in-situ experiment, J. Contam. Hydrol., 101, 1-4, 67-76 (2008)
[35] Muniruzzaman, M.; Rolle, M., Impact of multicomponent ionic transport on pH fronts propagation in saturated porous media, Water Resour. Res., 51, 8, 6739-6755 (2015)
[36] Kirby, B.J.: Micro-and nanoscale fluid mechanics: transport in microfluidic devices. Cambridge University Press, (2010) · Zbl 1283.76002
[37] Miller, G.; Trebotich, D., An embedded boundary method for the Navier-Stokes equations on a time-dependent domain, Commun. Appl. Math. Comput. Sci., 7, 1, 1-31 (2011) · Zbl 1273.35215
[38] Hu, M.; Wang, Y.; Rutqvist, J., Development of a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method, Int. J. Numer. Anal. Methods Geomech., 39, 17, 1932-1952 (2015)
[39] Appelo, C.; Van Loon, L.; Wersin, P., Multicomponent diffusion of a suite of tracers (HTO, Cl, Br, I, Na, Sr, Cs) in a single sample of Opalinus Clay, Geochim. Cosmochim. Acta, 74, 4, 1201-1219 (2010)
[40] Tournassat, C.; Appelo, C., Modelling approaches for anion-exclusion in compacted Na-bentonite, Geochim. Cosmochim. Acta, 75, 13, 3698-3710 (2011)
[41] Tinnacher, RM; Holmboe, M.; Tournassat, C.; Bourg, IC; Davis, JA, Ion adsorption and diffusion in smectite: molecular, pore, and continuum scale views, Geochim. Cosmochim. Acta, 177, 130-149 (2016)
[42] Soler, JM; Steefel, CI; Gimmi, T.; Leupin, OX; Cloet, V., Modeling the ionic strength effect on diffusion in clay. The DR-A Experiment at Mont Terri. ACS earth and space chemistry, ACS Earth Space Chem, 3, 3, 442-451 (2019)
[43] Muniruzzaman, M.; Rolle, M., Multicomponent ionic transport modeling in physically and electrostatically heterogeneous porous media with PhreeqcRM coupling for geochemical reactions, Water Resour. Res., 55, 12, 11121-11143 (2019)
[44] Onsager, L., Reciprocal relations in irreversible processes, I. Phys. Rev., 37, 4, 405-426 (1931) · JFM 57.1168.10
[45] Onsager, L., Reciprocal relations in irreversible processes, II. Phys. Rev., 38, 12, 2265-2279 (1931) · Zbl 0004.18303
[46] Einstein, A.: Über die von der molekularkinetischen Theorie der Wärme geforderte Bewegung von in ruhenden Flüssigkeiten suspendierten Teilchen. Annalen der Physik 4 (1905) · JFM 36.0975.01
[47] Steefel, C.; Appelo, C.; Arora, B.; Jacques, D.; Kalbacher, T.; Kolditz, O.; Lagneau, V.; Lichtner, P.; Mayer, KU; Meeussen, J., Reactive transport codes for subsurface environmental simulation, Comput. Geosci., 19, 3, 445-478 (2015) · Zbl 1323.86002
[48] Davis, JA, Surface complexation modeling in aqueous geochemistry, Mineral-Water Interf Geochem., 23, 177-259 (1990)
[49] Bourg, IC; Sposito, G., Molecular dynamics simulations of the electrical double layer on smectite surfaces contacting concentrated mixed electrolyte (NaCl-CaCl2) solutions, J. Colloid Interface Sci., 360, 2, 701-715 (2011)
[50] Park, H.; Lee, J.; Kim, T., Comparison of the Nernst-Planck model and the Poisson-Boltzmann model for electroosmotic flows in microchannels, J. Colloid Interface Sci., 315, 2, 731-739 (2007)
[51] Revil, A., Leroy, P.: Constitutive equations for ionic transport in porous shales. J. Geophys. Res.: Solid Earth 109(B3) (2004)
[52] Birgersson, M.; Karnland, O., Ion equilibrium between montmorillonite interlayer space and an external solution—consequences for diffusional transport, Geochim. Cosmochim. Acta, 73, 7, 1908-1923 (2009)
[53] Gimmi, T.; Alt-Epping, P., Simulating Donnan equilibria based on the Nernst-Planck equation, Geochim. Cosmochim. Acta, 232, 1-13 (2018)
[54] Sposito, G.: The surface chemistry of natural particles. Oxford University Press on Demand, (2004)
[55] Tournassat, C., Bourg, I.C., Steefel, C.I., Bergaya, F.: Surface properties of clay minerals. In: Developments in clay science, vol. 6. pp. 5-31. Elsevier, (2015)
[56] Appelo, C., Solute transport solved with the Nernst-Planck equation for concrete pores with ‘free’ water and a double layer, Cem. Concr. Res., 101, 102-113 (2017)
[57] Steefel, CI; Lasaga, AC, A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution reactions with application to reactive flow in single phase hydrothermal systems, Am. J. Sci., 294, 5, 529-592 (1994)
[58] Balay, S., Abhyankar, S., Adams, M., Brown, J., Brune, P., Buschelman, K., Dalcin, L., Dener, A., Eijkhout, V., Gropp, W.: PETSc Users Manual. (2019)
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.