An algebraic multiscale solver with local Robin boundary value problems for flows in high-contrast media
- PDF / 3,207,562 Bytes
- 20 Pages / 547.087 x 737.008 pts Page_size
- 77 Downloads / 197 Views
An algebraic multiscale solver with local Robin boundary value problems for flows in high-contrast media Hani Akbari
· Felipe Pereira
Received: 30 September 2019 / Accepted: 1 June 2020 © Springer Nature B.V. 2020
Abstract The algebraic and semi-algebraic formulations of a new multiscale elliptic solver based on a nonoverlapping 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 twoand three-dimensional problems with coefficients exhibiting high-contrast, anisotropic, and channelized structures that are solved to reveal the good properties of the proposed schemes. Keywords Algebraic multiscale methods · Defect correction scheme · Highly heterogeneous media · Preconditioning · Robin condition
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. H. Akbari (B) Shamim Institute of HPC, Scientific Visualization and Machine Learning, Mashhad, Iran e-mail: [email protected] F. Pereira Department of Mathematical Sciences, The University of Texas at Dallas, Richardson, TX, USA e-mail: [email protected]
123
H. Akbari, F. Pereira
1 Introduction Heterogeneous and anisotropic permeability data from large reservoirs yield ill-conditioned linear systems with huge number of unknowns. The solution of these systems provides steady-state (elliptic) or transient (parabolic) pressure of oil and/or gas in petroleum reservoirs. To solve such challenging boundary value problems, several approaches have been investigated. Here we focus on multiscale methods that decompose la
Data Loading...
