Numerical approximation of effective coefficients in stochastic homogenization of discrete elliptic equations

Antoine Gloria

ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique (2012)

  • Volume: 46, Issue: 1, page 1-38
  • ISSN: 0764-583X

Abstract

top
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 d 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 parameterT) on some box of finite size L. In this article, we replace the regularized corrector (which is the solution of a problem posed on d ) by some practically computable proxy on some box of sizeR≥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.

How to cite

top

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 46.1 (2012): 1-38. <http://eudml.org/doc/273104>.

@article{Gloria2012,
abstract = {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 $\mathbb \{Z\}^d$ 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 parameterT) on some box of finite size L. In this article, we replace the regularized corrector (which is the solution of a problem posed on $\mathbb \{Z\}^d$) by some practically computable proxy on some box of sizeR≥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.},
author = {Gloria, Antoine},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique},
keywords = {stochastic homogenization; effective coefficients; difference operator; numerical method; error estimate},
language = {eng},
number = {1},
pages = {1-38},
publisher = {EDP-Sciences},
title = {Numerical approximation of effective coefficients in stochastic homogenization of discrete elliptic equations},
url = {http://eudml.org/doc/273104},
volume = {46},
year = {2012},
}

TY - JOUR
AU - Gloria, Antoine
TI - Numerical approximation of effective coefficients in stochastic homogenization of discrete elliptic equations
JO - ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique
PY - 2012
PB - EDP-Sciences
VL - 46
IS - 1
SP - 1
EP - 38
AB - 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 $\mathbb {Z}^d$ 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 parameterT) on some box of finite size L. In this article, we replace the regularized corrector (which is the solution of a problem posed on $\mathbb {Z}^d$) by some practically computable proxy on some box of sizeR≥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.
LA - eng
KW - stochastic homogenization; effective coefficients; difference operator; numerical method; error estimate
UR - http://eudml.org/doc/273104
ER -

References

top
  1. [1] S. Agmon, 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). Zbl0503.35001MR745286
  2. [2] R. Alicandro, M. Cicalese and A. Gloria, Integral representation results for energies defined on stochastic lattices and application to nonlinear elasticity. Arch. Ration. Mech. Anal.200 (2011) 881–943. Zbl1294.74056MR2796134
  3. [3] A. Bourgeat and A. Piatnitski, Approximations of effective coefficients in stochastic homogenization. Ann. Inst. H. Poincaré40 (2004) 153–165. Zbl1058.35023MR2044813
  4. [4] P. Caputo and D. Ioffe, Finite volume approximation of the effective diffusion matrix: the case of independent bond disorder. Ann. Inst. H. Poincaré Probab. Statist.39 (2003) 505–525. Zbl1014.60094MR1978989
  5. [5] T. Delmotte, Inégalité de Harnack elliptique sur les graphes. Colloq. Math.72 (1997) 19–37. Zbl0871.31008MR1425544
  6. [6] A. Dykhne, 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. [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. Zbl1060.65118MR2114818
  8. [8] A. Gloria, Reduction of the resonance error – Part 1: Approximation of homogenized coefficients. Math. Models Methods Appl. Sci., to appear. Zbl1233.35016MR2826466
  9. [9] A. Gloria and F. Otto, An optimal variance estimate in stochastic homogenization of discrete elliptic equations. Ann. Probab.39 (2011) 779–856. Zbl1215.35025MR2789576
  10. [10] A. Gloria and F. Otto, An optimal error estimate in stochastic homogenization of discrete elliptic equations. Ann. Appl. Probab., to appear. Zbl06026087MR2932541
  11. [11] A. Gloria and F. Otto, Quantitative estimates in stochastic homogenization of linear elliptic equations. In preparation. Zbl06540728
  12. [12] T.Y. Hou and X.H. Wu, A Multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys.134 (1997) 169–189. Zbl0880.73065MR1455261
  13. [13] V.V. Jikov, S.M. Kozlov and O.A. Oleinik, Homogenization of Differential Operators and Integral Functionals. Springer-Verlag, Berlin (1994). Zbl0801.35001MR1329546
  14. [14] T. Kanit, S. Forest, I. Galliet, V. Mounoury and D. Jeulin, Determination of the size of the representative volume element for random composites: statistical and numerical approach. Int. J. Sol. Struct.40 (2003) 3647–3679. Zbl1038.74605
  15. [15] S.M. Kozlov, The averaging of random operators. Mat. Sb. (N.S.) 109 (1979) 188–202, 327. Zbl0415.60059MR542557
  16. [16] S.M. Kozlov, Averaging of difference schemes. Mat. Sb.57 (1987) 351–369. Zbl0639.65052
  17. [17] R. Künnemann, The diffusion limit for reversible jump processes on d with ergodic random bond conductivities. Commun. Math. Phys.90 (1983) 27–68. Zbl0523.60097MR714611
  18. [18] J.A. Meijerink and H.A. van der Vorst, An iterative solution method for linear systems of which the coefficient matrix is a symmetric M-matrix. Math. Comp.31 (1977) 148–162. Zbl0349.65020MR438681
  19. [19] A. Naddaf and T. Spencer, Estimates on the variance of some homogenization problems. Preprint (1998). Zbl0871.35010
  20. [20] H. Owhadi, Approximation of the effective conductivity of ergodic media by periodization. Probab. Theory Relat. Fields125 (2003) 225–258. Zbl1040.60025MR1961343
  21. [21] G.C. Papanicolaou and S.R.S. Varadhan, 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. Zbl0499.60059MR712714
  22. [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. Zbl1158.74541MR2313415
  23. [23] V.V. Yurinskii, Averaging of symmetric diffusion in random medium. Sibirskii Matematicheskii Zhurnal27 (1986) 167–180. Zbl0614.60051MR867870

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.