Finite element approximations of Green function \(G_{x_{0}}^{\epsilon}\) based on the method of multiscale asymptotic expansions (Q2490201)

The authors propose two kinds of effective finite element algorithms to obtain numerical approximations of Green function \(G_{x_0}^{\varepsilon}\) of the following elliptic boundary-value problem with oscillating coefficients: \[ \frac{\partial}{\partial x_i}\left(a^{ij}(\frac{x}{\varepsilon}) \frac{\partial G_{x_0}^{\varepsilon}}{\partial x_j} \right) = \delta(x - x_0), \quad \text{in}\quad \Omega \subseteq \mathbb R^2; \quad G_{x_0}^{\varepsilon} = 0, \quad \text{on}\quad \partial\Omega, \] where \(\delta(x_0 - x_0) = +\infty\); \(\delta(x - x_0) = 0, x \neq x_0\); \(\int_{\Omega} \delta(x - x_0)\,dx = 1,\) using the method of multiscale asymptotic expansions. This work is based on the article by \textit{W.-M. He} and \textit{J.-Z. Cui} [IMA J. Appl. Math. 70, 241--269 (2005; Zbl 1085.35024)]. Pointwise error estimates are presented, and their computational cost (memory and CPU time) are analyzed. Two numerical experiments are carried out to validate the theoretical results.