Skip to main content
Log in

An algebraic multiscale solver with local Robin boundary value problems for flows in high-contrast media

  • Published:
Journal of Engineering Mathematics Aims and scope Submit manuscript

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11

Similar content being viewed by others

References

  1. 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

    Article  MathSciNet  MATH  Google Scholar 

  2. 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

    Article  MATH  Google Scholar 

  3. 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

    Article  MathSciNet  MATH  Google Scholar 

  4. 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

    Article  MathSciNet  MATH  Google Scholar 

  5. Wang Y, Hajibeygi H, Tchelepi HA (2016) Monotone multiscale finite volume method. Comput Geosci 20(3):509–524

    Article  MathSciNet  MATH  Google Scholar 

  6. Cortinovis D, Jenny P (2017) Zonal multiscale finite-volume framework. J Comput Phys 337:84–97

    Article  MathSciNet  MATH  Google Scholar 

  7. 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

    Article  MathSciNet  MATH  Google Scholar 

  8. Aarnes J, Hou TY (2002) Multiscale domain decomposition methods for elliptic problems with high aspect ratios. Acta Math Appl Sin 18(1):63–76

    Article  MathSciNet  MATH  Google Scholar 

  9. Efendiev Y, Galvis J, Hou TY (2013) Generalized multiscale finite element methods (GMsFEM). J Comput Phys 251:116–135

    Article  MathSciNet  MATH  Google Scholar 

  10. Madureira AL, Sarkis M (2017) Hybrid localized spectral decomposition for multiscale problems. arXiv:1706.08941v2

  11. 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

    Article  MATH  Google Scholar 

  12. Wang Y, Hajibeygi H, Tchelepi HA (2014) Algebraic multi-scale solver for flow in heterogeneous porous media. J Comput Phys 259:284–303

    Article  MathSciNet  MATH  Google Scholar 

  13. Zhou H, Tchelepi HA (2012) Two-stage algebraic multi-scale linear solver for highly heterogeneous reservoir models, SPE J., SPE-141473-PA

  14. 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

    Article  Google Scholar 

  15. 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

    Article  MathSciNet  MATH  Google Scholar 

  16. Briggs WL, Henson VE, McCormick SF (2000) A multigrid tutorial, 2nd edn. SIAM, Philadelphia

    Book  MATH  Google Scholar 

  17. 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

    Article  MathSciNet  MATH  Google Scholar 

  18. 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

  19. 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

    Article  MathSciNet  MATH  Google Scholar 

  20. Arbogast T, Pencheva G, Wheeler MF, Yotov I (2007) A multiscale mortar mixed finite element method. SIAM Multiscale Model Simul 6(1):319–346

    Article  MathSciNet  MATH  Google Scholar 

  21. 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

    Article  MathSciNet  MATH  Google Scholar 

  22. 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

    Article  MathSciNet  MATH  Google Scholar 

  23. 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

    Article  MathSciNet  Google Scholar 

  24. 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)

  25. Yousef S (2000) Iterative methods for sparse linear systems (2nd ed.)

  26. Akbari H, Quadrio N, Engsig-Karup AP (2017) Unstructured mesh generation by wavelets for multiscale methods

  27. 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

    Article  MathSciNet  MATH  Google Scholar 

  28. Christie MA, Blunt MJ (2001) Tenth SPE comparative solution project: a comparison of upscaling techniques. SPE Reserv Eval Eng 4:308–317

    Article  Google Scholar 

Download references

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

Authors

Corresponding author

Correspondence to Hani Akbari.

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10665-020-10057-4

Keywords

Navigation