Efficiency and Accuracy of Micro‐Macro Models for Mineral Dissolution. Issue 8 (24th August 2020)

• Record Type: Journal Article
• Authors: Gärttner, Stephan; Frolkovič, Peter; Knabner, Peter; Ray, Nadja
• Is Part Of: Water resources research. Volume 56:Issue 8(2020)
• Page Start: n/a
• Page End: n/a
• Publication Date: 2020-08-24
• Subjects: Hydrology -- Periodicals; 333.91
• DOI: 10.1029/2020WR027585
• Languages: English
• ISSNs: 0043-1397

Abstract

Abstract: Micro‐macro models for dissolution processes are derived from detailed pore‐scale models applying upscaling techniques. They consist of flow and transport equations at the scale of the porous medium (macroscale). Both include averaged time‐ and space‐dependent coefficient functions (permeability, porosity, reactive surface, and effective diffusion). These are in turn explicitly computed from the time‐ and space‐dependent geometry of unit cells and by means of auxiliary cell problems defined therein (microscale). The explicit geometric structure is characterized by a level set. For its evolution, information from the transport equations solutions is taken into account (micro‐macro scales). A numerical scheme is introduced, which is capable of evaluating such complex settings. For the level‐set equation a second‐order scheme is applied, which enables us to accurately determine the dynamic reactive surface. Local mesh refinement methods are applied to evaluate Stokes type cell problems using P 2 / P 1 elements and a Uzawa type linear solver. Applications of our permeability solver to scenarios involving static and evolving geometries are presented. Furthermore, macroscopic flow and transport equations are solved applying mixed finite elements. Finally, adaptive strategies to overcome the computational burden are discussed. We apply our approach to the dissolution of an array of dolomite grains in the micro‐macro context and validate our numerical scheme. Key Points: We investigate reactive flow and transport problems for mineral dissolution A fully coupled nonlinear micro‐macro model is considered including evolving microstructures by means of level‐set methods Taking advantage of parallel computing and a tailored adaptivity scheme, computation times are comparable to purely macroscopic simulations