Abstract
The algebraic and semi-algebraic formulations of a new multiscale elliptic solver based on a non-overlapping domain decomposition method are presented. By algebraic we mean that all information needed to implement the proposed procedure can be extracted from the underlying fine grid finite volume linear system. In addition to the entries of this linear system, the proposed semi-algebraic procedure will also use the coefficients of the elliptic equation at subdomain boundaries. Initially we construct multiscale basis functions (or local solutions) subject to Robin boundary conditions. Although the implementation of Robin conditions in the formulation of the multiscale method uses the coefficients of the elliptic equation of interest, we modify it such that only entries of the finite volume linear system appear in the calculation of local solutions. A linear combination of local solutions gives a (discontinuous at subdomain boundaries) solution that is refined in a smoothing step. To this end we propose an algebraic scheme to remove discontinuities that needs a staggered (or dual) coarse grid. By iterating with the defect correction scheme, we construct an effective two-stage preconditioner that combines both local solutions and the smoothing step; the final result is obtained in a small number of iterations. Our focus in this work is the presentation of the algebraic and semi-algebraic formulations of the multiscale elliptic solver, which are carefully explained, and in verifying their accuracy and efficiency in the solution of challenging problems. We consider two- and three-dimensional problems with coefficients exhibiting high-contrast, anisotropic, and channelized structures that are solved to reveal the good properties of the proposed schemes.
Similar content being viewed by others
References
Kippe V, Aarnes J, Lie K (2008) A comparison of multiscale methods for elliptic problems in porous media flow. Comput Geosci 12(3):377–398
Jenny P, Lee SH, Tchelepi HA (2003) Multi-scale finite-volume method for elliptic problems in subsurface flow simulation. J Comput Phys 187(1):47–67
Jenny P, Lee SH, Tchelepi HA (2005) Adaptive multiscale finite-volume method for multiphase flow and transport in porous media. Multiscale Model Simul 3(1):50–64
Hajibeygi H, Karvounis D, Jenny P (2011) A hierarchical fracture model for the iterative multiscale finite volume method. J Comput Phys 230(24):8729–8743
Wang Y, Hajibeygi H, Tchelepi HA (2016) Monotone multiscale finite volume method. Comput Geosci 20(3):509–524
Cortinovis D, Jenny P (2017) Zonal multiscale finite-volume framework. J Comput Phys 337:84–97
Hou T, Wu XH (1997) A multiscale finite element method for elliptic problems in composite materials and porous media. J Comput Phys 134:169–189
Aarnes J, Hou TY (2002) Multiscale domain decomposition methods for elliptic problems with high aspect ratios. Acta Math Appl Sin 18(1):63–76
Efendiev Y, Galvis J, Hou TY (2013) Generalized multiscale finite element methods (GMsFEM). J Comput Phys 251:116–135
Madureira AL, Sarkis M (2017) Hybrid localized spectral decomposition for multiscale problems. arXiv:1706.08941v2
Hughes TJ, Feijoo R, Mazzei L, Quincy JB (1998) The variational multi-scale method: a paradigm for computational mechanics. Comput Methods Appl Mech Eng 166(1–2):3–24
Wang Y, Hajibeygi H, Tchelepi HA (2014) Algebraic multi-scale solver for flow in heterogeneous porous media. J Comput Phys 259:284–303
Zhou H, Tchelepi HA (2012) Two-stage algebraic multi-scale linear solver for highly heterogeneous reservoir models, SPE J., SPE-141473-PA
Francisco A, Ginting V, Pereira F, Rigelo J (2014) Design and implementation of a multiscale mixed method based on a non overlapping domain decomposition procedure. Math Comput Simul 99:125–138
Akbari H, Engsig-Karup AP, Ginting V, Pereira F (2019) A multiscale direct solver for the approximation of flows in high contrast porous media. J Comput Appl Math 359:88–101
Briggs WL, Henson VE, McCormick SF (2000) A multigrid tutorial, 2nd edn. SIAM, Philadelphia
Guiraldello RT, Ausas RF, Sousa FS, Pereira F, Buscaglia GC (2018) The multiscale Robin coupled method for flows in porous media. J Comput Phys 355:1–21
Guiraldello RT, Ausas RF, Sousa FS, Pereira F, Buscaglia GC (2018) Interface spaces for the multiscale Robin coupled method in reservoir simulation. Math Comput Simul 164(C):103–119
Douglas J, Paes-Leme PJ, Roberts JE, Wang JP (1993) A parallel iterative procedure applicable to the approximate solution of second order partial differential equations by mixed finite element methods. Numer Math 65(1):95–108
Arbogast T, Pencheva G, Wheeler MF, Yotov I (2007) A multiscale mortar mixed finite element method. SIAM Multiscale Model Simul 6(1):319–346
Harder C, Paredes D, Valentin F (2013) A family of multiscale hybrid-mixed finite element methods for the Darcy equation with rough coefficients. J Comput Phys 245:107–130
Amini S, Toutounian F, Gachpazan M (2018) The block CMRH method for solving non-symmetric linear systems with multiple right-hand sides. J Comput Appl Math 337:166–174
Amini S, Toutounian F (2018) Weighted and flexible versions of block CMRH method for solving non-symmetric linear systems with multiple right-hand sides. Comput Math Appl 76:2011–2021
Abreu E, Ferraz P, Pereira F, Santo A, Santos L, Sousa FS (2020) A recursive formulation of multiscale mixed methods and its parallel implementation (in press)
Yousef S (2000) Iterative methods for sparse linear systems (2nd ed.)
Akbari H, Quadrio N, Engsig-Karup AP (2017) Unstructured mesh generation by wavelets for multiscale methods
Guiraldello RT, Ausas RF, Sousa FS, Pereira F, Buscaglia GC (2020) Velocity postprocessing schemes for multiscale mixed methods applied to contaminant transport in subsurface flows. Comput Geosci 24:1141–1161
Christie MA, Blunt MJ (2001) Tenth SPE comparative solution project: a comparison of upscaling techniques. SPE Reserv Eval Eng 4:308–317
Acknowledgements
FP was funded in part by NSF-DMS 1514808, a Science Without Borders/CNPq-Brazil grant and UT Dallas. He also wishes to thank the hospitality of ICMC/USP in Brazil during the summer of 2019 and funding from a CAPES/Print grant where a revision of this word was carried out.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
H. Akbari: This work was written, in part, while H. A. was a postdoctoral research associate at the Department of Infrastructure Engineering, The University of Melbourne, Parkville, Australia.
Rights and permissions
About this article
Cite this article
Akbari, H., Pereira, F. An algebraic multiscale solver with local Robin boundary value problems for flows in high-contrast media. J Eng Math 123, 109–128 (2020). https://doi.org/10.1007/s10665-020-10057-4
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10665-020-10057-4