Library Subscription: Guest
Begell Digital Portal Begell Digital Library eBooks Journals References & Proceedings Research Collections
International Journal for Uncertainty Quantification

Impact factor: 0.967

ISSN Print: 2152-5080
ISSN Online: 2152-5099

Open Access

International Journal for Uncertainty Quantification

DOI: 10.1615/Int.J.UncertaintyQuantification.2016016805
pages 429-443

A MULTIMODES MONTE CARLO FINITE ELEMENT METHOD FOR ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS WITH RANDOM COEFFICIENTS

Xiaobing Feng
Department of Mathematics, The University of Tennessee, Knoxville, Tennessee 37996, USA
Junshan Lin
Department of Mathematics and Statistics, Auburn University, Auburn, Alabama 36849, USA
Cody Lorton
Department of Mathematics and Statistics, University of West Florida, Pensacola, Floridas 32514, USA

ABSTRACT

This paper develops and analyzes an efficient numerical method for solving elliptic partial differential equations, where the diffusion coefficients are random perturbations of deterministic diffusion coefficients. The method is based upon a multimodes representation of the solution as a power series of the perturbation parameter, and the Monte Carlo technique for sampling the probability space. One key feature of the proposed method is that the governing equations for all the expanded mode functions share the same deterministic diffusion coefficient; thus an efficient direct solver by repeatedly using the LU decomposition of the discretized common deterministic diffusion operator can be employed for solving the finite element discretized linear systems. It is shown that the computational complexity of the algorithm is comparable to that of solving a few deterministic elliptic partial differential equations using the director solver. Error estimates are derived for the method, and numerical experiments are provided to test the efficiency of the algorithm and validate the theoretical results.