Numerical approximation of effective ...
Document type :
Compte-rendu et recension critique d'ouvrage
DOI :
Title :
Numerical approximation of effective coefficients in stochastic homogenization of discrete elliptic equations
Author(s) :
Journal title :
ESAIM: Mathematical Modelling and Numerical Analysis
Pages :
1-38
Publisher :
EDP Sciences
Publication date :
2012
ISSN :
0764-583X
HAL domain(s) :
Mathématiques [math]/Analyse numérique [math.NA]
English abstract : [en]
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 ...
Show more >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 $\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 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 $\Z^d$) by some practically computable proxy on some box of size $R\geq 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 arithmetic mean 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 allows us to give a clear answer to 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.Show less >
Show more >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 $\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 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 $\Z^d$) by some practically computable proxy on some box of size $R\geq 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 arithmetic mean 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 allows us to give a clear answer to 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.Show less >
Language :
Anglais
Popular science :
Non
Collections :
Source :
Files
- document
- Open access
- Access the document
- m2an11-gloria.pdf
- Open access
- Access the document
- m2an110018.pdf
- Open access
- Access the document