Finite Difference Approximation of Homogenization Problems for Elliptic Equations
Rafael Orive, Enrique Zuazua · Multiscale Modeling and Simulation · 2005
In this paper, the problem of the approximation by finite differences of solutions to elliptic problems with rapidly oscillating coefficients and periodic boundary conditions is considered. The mesh size is denoted by h, while $\varepsilon$ denotes the period of the rapidly oscillating coefficient. Using Bloch wave decompositions, we analyze the case where the ratio $h/\varepsilon$ is rational. We show that if $h/\varepsilon$ is kept fixed, being a rational number, even when $h,\varepsilon \to 0$, the limit of the numerical solution does not coincide with the homogenized one obtained when passing to the limit as $\varepsilon \to 0$ in the continuous problem. Explicit error estimates are given showing that, as the ratio $h/\varepsilon$ approximates an irrational number, solutions of the finite difference approximation converge to the solutions of the homogenized elliptic equation. We consider both the one-dimensional and the multidimensional case. Our analysis yields a quantitative version of previous results on numerical homogenization by Avellaneda, Hou, and Papanicolaou [RAIRO Modél. Math. Anal. Numér., 25 (1991), pp. 693--710].