Multiscale finite element methods for high-contrast problems using local spectral basis functions

Yalchin R. Efendiev, Juan Galvis, Xiao-Hui Wu

Research output: Contribution to journalArticlepeer-review

174 Scopus citations

Abstract

In this paper we study multiscale finite element methods (MsFEMs) using spectral multiscale basis functions that are designed for high-contrast problems. Multiscale basis functions are constructed using eigenvectors of a carefully selected local spectral problem. This local spectral problem strongly depends on the choice of initial partition of unity functions. The resulting space enriches the initial multiscale space using eigenvectors of local spectral problem. The eigenvectors corresponding to small, asymptotically vanishing, eigenvalues detect important features of the solutions that are not captured by initial multiscale basis functions. Multiscale basis functions are constructed such that they span these eigenfunctions that correspond to small, asymptotically vanishing, eigenvalues. We present a convergence study that shows that the convergence rate (in energy norm) is proportional to (H/Λ*)1/2, where Λ* is proportional to the minimum of the eigenvalues that the corresponding eigenvectors are not included in the coarse space. Thus, we would like to reach to a larger eigenvalue with a smaller coarse space. This is accomplished with a careful choice of initial multiscale basis functions and the setup of the eigenvalue problems. Numerical results are presented to back-up our theoretical results and to show higher accuracy of MsFEMs with spectral multiscale basis functions. We also present a hierarchical construction of the eigenvectors that provides CPU savings. © 2010.
Original languageEnglish (US)
Pages (from-to)937-955
Number of pages19
JournalJournal of Computational Physics
Volume230
Issue number4
DOIs
StatePublished - Feb 2011
Externally publishedYes

Fingerprint Dive into the research topics of 'Multiscale finite element methods for high-contrast problems using local spectral basis functions'. Together they form a unique fingerprint.

Cite this