Numerical study in stochastic homogenization for elliptic PDEs: convergence rate in the size of representative volume elements
From MaRDI portal
Publication:6316374
DOI10.1002/NLA.2296arXiv1903.12227MaRDI QIDQ6316374
Venera Khoromskaia, Boris Khoromskij, Felix Otto
Publication date: 28 March 2019
Abstract: 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 CG 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 [6]. 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.
Computational methods for sparse matrices (65F50) Spectral, collocation and related methods for boundary value problems involving PDEs (65N35) Iterative numerical methods for linear systems (65F10)
This page was built for publication: Numerical study in stochastic homogenization for elliptic PDEs: convergence rate in the size of representative volume elements
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6316374)
