We introduce and analyze a numerical strategy to approximate effective coefficients in stochastic homogenization of discrete elliptic equations. In particular, we consider the simplest case possible: An elliptic equation on the d-dimensional lattice with independent and identically distributed conductivities on the associated edges. Recent results by Otto and the author quantify the error made by approximating the homogenized coefficient by the averaged energy of a regularized corrector (with parameter T) on some box of finite size L. In this article, we replace the regularized corrector (which is the solution of a problem posed on ) by some practically computable proxy on some box of size R ≥ L, and quantify the associated additional error. In order to improve the convergence, one may also consider N independent realizations of the computable proxy, and take the empirical average of the associated approximate homogenized coefficients. A natural optimization problem consists in properly choosing T, R, L and N in order to reduce the error at given computational complexity. Our analysis is sharp and sheds some light on this question. In particular, we propose and analyze a numerical algorithm to approximate the homogenized coefficients, taking advantage of the (nearly) optimal scalings of the errors we derive. The efficiency of the approach is illustrated by a numerical study in dimension 2.
@article{M2AN_2012__46_1_1_0,
author = {Gloria, Antoine},
title = {Numerical approximation of effective coefficients in stochastic homogenization of discrete elliptic equations},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique},
volume = {46},
year = {2012},
pages = {1-38},
doi = {10.1051/m2an/2011018},
mrnumber = {2846365},
zbl = {1282.35038},
language = {en},
url = {http://dml.mathdoc.fr/item/M2AN_2012__46_1_1_0}
}
Gloria, Antoine. Numerical approximation of effective coefficients in stochastic homogenization of discrete elliptic equations. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 46 (2012) pp. 1-38. doi : 10.1051/m2an/2011018. http://gdmltest.u-ga.fr/item/M2AN_2012__46_1_1_0/
[1] , Lectures on exponential decay of solutions of second-order elliptic equations: bounds on eigenfunctions of N-body Schrödinger operators, Mathematical Notes 29. Princeton University Press, Princeton, NJ (1982). | MR 745286 | Zbl 0503.35001
[2] , and , Integral representation results for energies defined on stochastic lattices and application to nonlinear elasticity. Arch. Ration. Mech. Anal. 200 (2011) 881-943. | MR 2796134 | Zbl 1294.74056
[3] and , Approximations of effective coefficients in stochastic homogenization. Ann. Inst. H. Poincaré 40 (2004) 153-165. | Numdam | MR 2044813 | Zbl 1058.35023
[4] and , Finite volume approximation of the effective diffusion matrix: the case of independent bond disorder. Ann. Inst. H. Poincaré Probab. Statist. 39 (2003) 505-525. | Numdam | MR 1978989 | Zbl 1014.60094
[5] , Inégalité de Harnack elliptique sur les graphes. Colloq. Math. 72 (1997) 19-37. | MR 1425544 | Zbl 0871.31008
[6] , Conductivity of a two-dimensional two-phase system. Sov. Phys. JETP 32 (1971) 63-65. Russian version: Zh. Eksp. Teor. Fiz. 59 (1970) 110-5.
[7] W. E, P.B. Ming and P.W. Zhang, Analysis of the heterogeneous multiscale method for elliptic homogenization problems. J. Amer. Math. Soc. 18 (2005) 121-156. | MR 2114818 | Zbl 1060.65118
[8] , Reduction of the resonance error - Part 1: Approximation of homogenized coefficients. Math. Models Methods Appl. Sci., to appear. | MR 2826466 | Zbl 1233.35016
[9] and , An optimal variance estimate in stochastic homogenization of discrete elliptic equations. Ann. Probab. 39 (2011) 779-856. | MR 2789576 | Zbl 1215.35025
[10] and , An optimal error estimate in stochastic homogenization of discrete elliptic equations. Ann. Appl. Probab., to appear. | MR 2932541 | Zbl pre06026087
[11] and , Quantitative estimates in stochastic homogenization of linear elliptic equations. In preparation.
[12] and , A Multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys. 134 (1997) 169-189. | MR 1455261 | Zbl 0880.73065
[13] , and , Homogenization of Differential Operators and Integral Functionals. Springer-Verlag, Berlin (1994). | MR 1329546 | Zbl 0801.35001
[14] , , , and , Determination of the size of the representative volume element for random composites: statistical and numerical approach. Int. J. Sol. Struct. 40 (2003) 3647-3679. | Zbl 1038.74605
[15] , The averaging of random operators. Mat. Sb. (N.S.) 109 (1979) 188-202, 327. | MR 542557 | Zbl 0415.60059
[16] , Averaging of difference schemes. Mat. Sb. 57 (1987) 351-369. | Zbl 0639.65052
[17] , The diffusion limit for reversible jump processes on with ergodic random bond conductivities. Commun. Math. Phys. 90 (1983) 27-68. | MR 714611 | Zbl 0523.60097
[18] and , An iterative solution method for linear systems of which the coefficient matrix is a symmetric M-matrix. Math. Comp. 31 (1977) 148-162. | MR 438681 | Zbl 0349.65020
[19] and , Estimates on the variance of some homogenization problems. Preprint (1998).
[20] , Approximation of the effective conductivity of ergodic media by periodization. Probab. Theory Relat. Fields 125 (2003) 225-258. | MR 1961343 | Zbl 1040.60025
[21] and , Boundary value problems with rapidly oscillating random coefficients, in Random fields I, II (Esztergom, 1979), Colloq. Math. Soc. János Bolyai 27. North-Holland, Amsterdam (1981) 835-873. | MR 712714 | Zbl 0499.60059
[22] X. Yue and W. E, The local microscale problem in the multiscale modeling of strongly heterogeneous media: effects of boundary conditions and cell size. J. Comput. Phys. 222 (2007) 556-572. | MR 2313415 | Zbl 1158.74541
[23] , Averaging of symmetric diffusion in random medium. Sibirskii Matematicheskii Zhurnal 27 (1986) 167-180. | MR 867870 | Zbl 0614.60051