@article{13129, abstract = {We study the representative volume element (RVE) method, which is a method to approximately infer the effective behavior ahom of a stationary random medium. The latter is described by a coefficient field a(x) generated from a given ensemble ⟨⋅⟩ and the corresponding linear elliptic operator −∇⋅a∇. In line with the theory of homogenization, the method proceeds by computing d=3 correctors (d denoting the space dimension). To be numerically tractable, this computation has to be done on a finite domain: the so-called representative volume element, i.e., a large box with, say, periodic boundary conditions. The main message of this article is: Periodize the ensemble instead of its realizations. By this, we mean that it is better to sample from a suitably periodized ensemble than to periodically extend the restriction of a realization a(x) from the whole-space ensemble ⟨⋅⟩. We make this point by investigating the bias (or systematic error), i.e., the difference between ahom and the expected value of the RVE method, in terms of its scaling w.r.t. the lateral size L of the box. In case of periodizing a(x), we heuristically argue that this error is generically O(L−1). In case of a suitable periodization of ⟨⋅⟩ , we rigorously show that it is O(L−d). In fact, we give a characterization of the leading-order error term for both strategies and argue that even in the isotropic case it is generically non-degenerate. We carry out the rigorous analysis in the convenient setting of ensembles ⟨⋅⟩ of Gaussian type, which allow for a straightforward periodization, passing via the (integrable) covariance function. This setting has also the advantage of making the Price theorem and the Malliavin calculus available for optimal stochastic estimates of correctors. We actually need control of second-order correctors to capture the leading-order error term. This is due to inversion symmetry when applying the two-scale expansion to the Green function. As a bonus, we present a stream-lined strategy to estimate the error in a higher-order two-scale expansion of the Green function.}, author = {Clozeau, Nicolas and Josien, Marc and Otto, Felix and Xu, Qiang}, issn = {1615-3383}, journal = {Foundations of Computational Mathematics}, pages = {1305--1387}, publisher = {Springer Nature}, title = {{Bias in the representative volume element method: Periodize the ensemble instead of its realizations}}, doi = {10.1007/s10208-023-09613-y}, volume = {24}, year = {2024}, } @article{14797, abstract = {We study a random matching problem on closed compact 2 -dimensional Riemannian manifolds (with respect to the squared Riemannian distance), with samples of random points whose common law is absolutely continuous with respect to the volume measure with strictly positive and bounded density. We show that given two sequences of numbers n and m=m(n) of points, asymptotically equivalent as n goes to infinity, the optimal transport plan between the two empirical measures μn and νm is quantitatively well-approximated by (Id,exp(∇hn))#μn where hn solves a linear elliptic PDE obtained by a regularized first-order linearization of the Monge-Ampère equation. This is obtained in the case of samples of correlated random points for which a stretched exponential decay of the α -mixing coefficient holds and for a class of discrete-time Markov chains having a unique absolutely continuous invariant measure with respect to the volume measure.}, author = {Clozeau, Nicolas and Mattesini, Francesco}, issn = {1432-2064}, journal = {Probability Theory and Related Fields}, pages = {485--541}, publisher = {Springer Nature}, title = {{Annealed quantitative estimates for the quadratic 2D-discrete random matching problem}}, doi = {10.1007/s00440-023-01254-0}, volume = {190}, year = {2024}, } @article{17462, abstract = {We are interested in numerical algorithms for computing the electrical field generated by a charge distribution localized on scale l in an infinite heterogeneous correlated random medium, in a situation where the medium is only known in a box of diameter L >>l around the support of the charge. We show that the algorithm in [J. Lu, F. Otto, and L. Wang, Optimal Artificial Boundary Conditions Based on Second-Order Correctors for Three Dimensional Random Ellilptic Media, preprint, arXiv:2109.01616, 2021], suggesting optimal Dirichlet boundary conditions motivated by the multipole expansion [P. Bella, A. Giunti, and F. Otto, Comm. Partial Differential Equations, 45 (2020), pp. 561–640], still performs well in correlated media. With overwhelming probability, we obtain a convergence rate in terms of l, L, and the size of the correlations for which optimality is supported with numerical simulations. These estimates are provided for ensembles which satisfy a multiscale logarithmic Sobolev inequality, where our main tool is an extension of the semigroup estimates in [N. Clozeau, Stoch. Partial Differ. Equ. Anal. Comput., 11 (2023), pp. 1254–1378]. As part of our strategy, we construct sublinear second-order correctors in this correlated setting, which is of independent interest.}, author = {Clozeau, Nicolas and Wang, Lihan}, issn = {1540-3467}, journal = {Multiscale Modeling and Simulation}, number = {3}, pages = {973--1029}, publisher = {Society for Industrial and Applied Mathematics}, title = {{Artificial boundary conditions for random elliptic systems with correlated coefficient field}}, doi = {10.1137/23M1603819}, volume = {22}, year = {2024}, } @article{10173, abstract = {We study the large scale behavior of elliptic systems with stationary random coefficient that have only slowly decaying correlations. To this aim we analyze the so-called corrector equation, a degenerate elliptic equation posed in the probability space. In this contribution, we use a parabolic approach and optimally quantify the time decay of the semigroup. For the theoretical point of view, we prove an optimal decay estimate of the gradient and flux of the corrector when spatially averaged over a scale R larger than 1. For the numerical point of view, our results provide convenient tools for the analysis of various numerical methods.}, author = {Clozeau, Nicolas}, issn = {2194-0401}, journal = {Stochastics and Partial Differential Equations: Analysis and Computations}, pages = {1254–1378}, publisher = {Springer Nature}, title = {{Optimal decay of the parabolic semigroup in stochastic homogenization for correlated coefficient fields}}, doi = {10.1007/s40072-022-00254-w}, volume = {11}, year = {2023}, } @article{10174, abstract = {Quantitative stochastic homogenization of linear elliptic operators is by now well-understood. In this contribution we move forward to the nonlinear setting of monotone operators with p-growth. This first work is dedicated to a quantitative two-scale expansion result. Fluctuations will be addressed in companion articles. By treating the range of exponents 2≤p<∞ in dimensions d≤3, we are able to consider genuinely nonlinear elliptic equations and systems such as −∇⋅A(x)(1+|∇u|p−2)∇u=f (with A random, non-necessarily symmetric) for the first time. When going from p=2 to p>2, the main difficulty is to analyze the associated linearized operator, whose coefficients are degenerate, unbounded, and depend on the random input A via the solution of a nonlinear equation. One of our main achievements is the control of this intricate nonlinear dependence, leading to annealed Meyers' estimates for the linearized operator, which are key to the quantitative two-scale expansion result.}, author = {Clozeau, Nicolas and Gloria, Antoine}, issn = {1432-0673}, journal = {Archive for Rational Mechanics and Analysis }, number = {4}, publisher = {Springer Nature}, title = {{Quantitative nonlinear homogenization: Control of oscillations}}, doi = {10.1007/s00205-023-01895-4}, volume = {247}, year = {2023}, } @article{10175, abstract = {We study periodic homogenization by Γ-convergence of integral functionals with integrands W(x,ξ) having no polynomial growth and which are both not necessarily continuous with respect to the space variable and not necessarily convex with respect to the matrix variable. This allows to deal with homogenization of composite hyperelastic materials consisting of two or more periodic components whose the energy densities tend to infinity as the volume of matter tends to zero, i.e., W(x,ξ)=∑j∈J1Vj(x)Hj(ξ) where {Vj}j∈J is a finite family of open disjoint subsets of RN, with |∂Vj|=0 for all j∈J and ∣∣RN∖⋃j∈JVj|=0, and, for each j∈J, Hj(ξ)→∞ as detξ→0. In fact, our results apply to integrands of type W(x,ξ)=a(x)H(ξ) when H(ξ)→∞ as detξ→0 and a∈L∞(RN;[0,∞[) is 1-periodic and is either continuous almost everywhere or not continuous. When a is not continuous, we obtain a density homogenization formula which is a priori different from the classical one by Braides–Müller. Although applications to hyperelasticity are limited due to the fact that our framework is not consistent with the constraint of noninterpenetration of the matter, our results can be of technical interest to analysis of homogenization of integral functionals.}, author = {Anza Hafsa, Omar and Clozeau, Nicolas and Mandallena, Jean-Philippe}, issn = {2118-7436}, journal = {Annales mathématiques Blaise Pascal}, number = {2}, pages = {135--193}, publisher = {Université Clermont Auvergne}, title = {{Homogenization of nonconvex unbounded singular integrals}}, doi = {10.5802/ambp.367}, volume = {24}, year = {2017}, }