z-logo
Premium
Numerical study in stochastic homogenization for elliptic partial differential equations: Convergence rate in the size of representative volume elements
Author(s) -
Khoromskaia Venera,
Khoromskij Boris N.,
Otto Felix
Publication year - 2020
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.2296
Subject(s) - mathematics , rate of convergence , discretization , mathematical analysis , homogenization (climate) , stiffness matrix , finite element method , biodiversity , ecology , channel (broadcasting) , electrical engineering , biology , engineering , physics , thermodynamics
Summary We describe the numerical scheme for the discretization and solution of 2D elliptic equations with strongly varying piecewise constant coefficients arising in the stochastic homogenization of multiscale composite materials. An efficient stiffness matrix generation scheme based on assembling the local Kronecker product matrices is introduced. The resulting large linear systems of equations are solved by the preconditioned conjugate gradient iteration with a convergence rate that is independent of the grid size and the variation in jumping coefficients (contrast). Using this solver, we numerically investigate the convergence of the representative volume element (RVE) method in stochastic homogenization that extracts the effective behavior of the random coefficient field. Our numerical experiments confirm the asymptotic convergence rate of systematic error and standard deviation in the size of RVE rigorously established in Gloria et al. The asymptotic behavior of covariances of the homogenized matrix in the form of a quartic tensor is also studied numerically. Our approach allows laptop computation of sufficiently large number of stochastic realizations even for large sizes of the RVE.

This content is not available in your region!

Continue researching here.

Having issues? You can contact us here