Optimal Local Approximation Spaces for Generalized Finite Element Methods with Application to Multiscale Problems
Abstract
The paper addresses a numerical method for solving second order elliptic partial differential equations that describe fields inside heterogeneous media. The scope is general and treats the case of rough coefficients, i.e. coefficients with values in L∞(Ω). This class of coefficients includes as examples media with micro-structure as well as media with multiple non-separated length scales. The approach taken here is based on the the generalized finite element method (GFEM) introduced in 107, and elaborated in 102, 103 and 104. The GFEM is constructed by partitioning the computational domain Ω into to a collection of preselected subsets ωi,i=1,2,..m and constructing finite dimensional approximation spaces Ψi over each subset using local information. The notion of the Kolmogorov n-width is used to identify the optimal local approximation spaces. These spaces deliver local approximations with errors that decay almost exponentially with the degrees of freedom Ni in the energy norm over ωi. The local spaces % Ψi are used within the GFEM scheme to produce a finite dimensional subspace SN of H1(Ω) which is then employed in the Galerkin method. It is shown that the error in the Galerkin approximation decays in the energy norm almost exponentially (i.e., super-algebraicly) with respect to the degrees of freedom N. When length scales "`separate" and the microstructure is sufficiently fine with respect to the length scale of the domain ωi it is shown that homogenization theory can be used to construct local approximation spaces with exponentially decreasing error in the pre-asymtotic regime.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.