A multiscale direct solver for the approximation of flows in high contrast porous media

Hani Akbari*, Allan Peter Engsig-Karup, Victor Ginting, Felipe Pereira

*Corresponding author for this work

Research output: Contribution to journalJournal articleResearchpeer-review

5 Downloads (Pure)


We consider a non-overlapping domain decomposition approach to approximate the solution of elliptic boundary value problems with high contrast in their coefficients. We propose a method such that initially local solutions subject to Robin boundary conditions in each primal subdomain are constructed with (locally conservative) finite element or finite volume methods. Then, a novel approach is introduced to obtain a (discontinuous) global solution in terms of linear combination of the local subdomain solutions. In the proposed algorithm the computation of local solutions for unions of subdomains are localized at nearest-neighbor subdomain boundaries, thus avoiding the solution of global interface problems. We remove discontinuities in a smoothing step that is defined on a staggered grid or dual subdomains. The resulting algorithm is naturally parallelizable and can be employed as a parallel direct solver, offering great potential for the numerical solution of large problems. In fact, subdomains can be considered small enough to fit well in GPUs and the proposed procedure can handle adaptive (in space) simulations effectively. Numerical simulations are presented and discussed. We demonstrate the effectiveness of the proposed approach with two and three dimensional high contrast and channelized coefficients, that lead to challenging approximation problems. The new procedure, although designed for parallel processing, is also of value for serial calculations.
Original languageEnglish
JournalJournal of Computational and Applied Mathematics
Pages (from-to)88-101
Number of pages14
Publication statusPublished - 2019


  • Elliptic equations
  • Combination algorithm
  • Multiscale method
  • Robin boundary condition
  • Parallel computing
  • High contrast media

Fingerprint Dive into the research topics of 'A multiscale direct solver for the approximation of flows in high contrast porous media'. Together they form a unique fingerprint.

Cite this