We present an efficient method for the computation of homogenized coefficients of divergence-form operators with random coefficients. The approach is based on a multiscale representation of the homogenized coefficients. We then implement the method numerically using a finite-element method with hierarchical hybrid grids, which is a semi-implicit method allowing for significant gains in memory usage and execution time. Finally, we demonstrate the efficiency of our approach on two- and three-dimensional examples, for piecewise-constant coefficients with corner discontinuities. For moderate ellipticity contrast and for a precision of a few percentage points, our method allows to compute the homogenized coefficients on a laptop computer in a few seconds, in two dimensions, or in a few minutes, in three dimensions.
Over the last decade, an intensive research effort has been devoted to developing theoretical quantitative results on stochastic homogenization. The multiscale representation of the homogenized coefficients forming the basis of the method is inspired by the renormalization'' approach to quantitative stochastic homogenization, as developed in @cite_47 @cite_38 @cite_60 @cite_15 ; see also @cite_62 for a gentle introduction to this line of research and @cite_36 for a monograph. A related approach based on the parabolic flow was put forward in @cite_37 , see also [Chapter 9] AKMbook , and will give us the most convenient statement for us to build upon here. A different approach based on concentration inequalities was put forward in @cite_42 @cite_7 @cite_27 @cite_0 @cite_55 @cite_56 , inspired by earlier insights from statistical mechanics @cite_33 @cite_10 .
@cite_38 @cite_37 @cite_62 @cite_33 @cite_7 @cite_15 @cite_60 @cite_36 @cite_55 @cite_42 @cite_56 @cite_0 @cite_27 @cite_47 @cite_10
2107101459 2010942428 2963794691 2121999212
We study quantitatively the effective large-scale behavior of discrete elliptic equations on the lattice ( Z^d ) with random coefficients. The theory of stochastic homogenization relates the random, stationary, and ergodic field of coefficients with a deterministic matrix of effective coefficients. This is done via the corrector problem, which can be viewed as a highly degenerate elliptic equation on the infinite-dimensional space of admissible coefficient fields. In this contribution we develop new quantitative methods for the corrector problem based on the assumption that ergodicity holds in the quantitative form of a Spectral Gap Estimate w.r.t. a Glauber dynamics on coefficient fields—as it is the case for independent and identically distributed coefficients. As a main result we prove an optimal decay in time of the semigroup associated with the corrector problem (i.e. of the generator of the process called “random environment as seen from the particle”). As a corollary we recover existence of stationary correctors (in dimensions (d>2 )) and prove new optimal estimates for regularized versions of the corrector (in dimensions (d 2 )). We also give a self-contained proof of a new estimate on the gradient of the parabolic, variable-coefficient Green’s function, which is a crucial analytic ingredient in our approach. As an application of these results, we prove the first (and optimal) estimates for the approximation of the homogenized coefficients by the popular periodization method in case of independent and identically distributed coefficients. 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 @math -dimensional lattice @math 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 @math ) on some box of finite size @math . In this article, we replace the regularized corrector (which is the solution of a problem posed on @math ) by some practically computable proxy on some box of size @math , and quantify the associated additional error. In order to improve the convergence, one may also consider @math 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 @math and @math 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. This work is a follow-up to our previous work "A numerical approach related to defect-type theories for some weakly random problems in homogenization" (preprint available on this archive). It extends and complements, both theoretically and experimentally, the results presented there. Under consideration is the homogenization of a model of a weakly random heterogeneous material. The material consists of a reference periodic material randomly perturbed by another periodic material, so that its homogenized behavior is close to that of the reference material. We consider laws for the random perturbations more general than in our previous work cited above. We prove the validity of an asymptotic expansion in a certain class of settings. We also extend the formal approach introduced in our former work. Our perturbative approach shares common features with a defect-type theory of solid state physics. The computational efficiency of the approach is demonstrated. This article deals with the numerical approximation of effective coefficients in stochastic homogenization of discrete linear elliptic equations. The originality of this work is the use of a well-known abstract spectral representation formula to design and analyze effective and computable approximations of the homogenized coefficients. In particular, we show that information on the edge of the spectrum of the generator of the environment viewed by the particle projected on the local drift yields bounds on the approximation error, and conversely. Combined with results by Otto and the first author in low dimension, and results by the second author in high dimension, this allows us to prove that for any dimension d ≥ 2, there exists an explicit numerical strategy to approximate homogenized coefficients which converges at the rate of the central limit theorem.