A Monte Carlo approach to computing stiffness matrices arising in polynomial chaos approximations
Numerical Analysis
2017-09-12 v2
Abstract
We use a Monte Carlo method to assemble finite element matrices for polynomial Chaos approximations of elliptic equations with random coefficients. In this approach, all required expectations are approximated by a Monte Carlo method. The resulting methodology requires dealing with sparse block-diagonal matrices instead of block-full matrices. This leads to the solution of a coupled system of elliptic equations where the coupling is given by a Kronecker product matrix involving polynomial evaluation matrices. This generalizes the Classical Monte Carlo approximation and Collocation method for approximating functionals of solutions of these equations.
Cite
@article{arxiv.1704.06339,
title = {A Monte Carlo approach to computing stiffness matrices arising in polynomial chaos approximations},
author = {Juan Galvis and O. Andres Cuervo},
journal= {arXiv preprint arXiv:1704.06339},
year = {2017}
}