Solving elliptic boundary value problems with uncertain coefficients by the finite element method: the stochastic formulation

This work studies a linear elliptic problem with uncertainty. The introduction gives a survey of different formulations of the uncertainty and resulting numerical approximations. The major emphasis of this work is the probabilistic treatment of uncertainty, addressing the problem of solving linear e...

Full description

Saved in:
Bibliographic Details
Published in:Computer methods in applied mechanics and engineering 2005-04, Vol.194 (12-16), p.1251-1294
Main Authors: Babuška, Ivo, Tempone, Raúl, Zouraris, Georgios E.
Format: Article
Language:English
Subjects:
Adaptive methods
domain
Error control
Error estimates
Expected value
Finite elements
groundwater models cannot
k x h-version
k × h-version
Karhunen-Loeve expansion
Karhunen–Loève expansion
Monte Carlo method
p x h-version
p × h-version
partial-differential-equations
Perturbation estimates
science
sensitivity-analysis
Stochastic elliptic equation
validation
verification
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:This work studies a linear elliptic problem with uncertainty. The introduction gives a survey of different formulations of the uncertainty and resulting numerical approximations. The major emphasis of this work is the probabilistic treatment of uncertainty, addressing the problem of solving linear elliptic boundary value problems with stochastic coefficients. If the stochastic coefficients are known functions of a random vector, then the stochastic elliptic boundary value problem is turned into a parametric deterministic one with solution u(y,x), y∈Γ, x∈D, where D⊂Rd, d=1,2,3, and Γ is a high-dimensional cube. In addition, the function u is specified as the solution of a deterministic variational problem over Γ×D. A tensor product finite element method, of h-version in D and k-, or, p-version in Γ, is proposed for the approximation of u. A priori error estimates are given and an adaptive algorithm is also proposed. Due to the high dimension of Γ, the Monte Carlo finite element method is also studied here. This work compares the asymptotic complexity of the numerical methods, and shows results from numerical experiments. Comments on the uncertainty in the probabilistic characterization of the coefficients in the stochastic formulation are included.
ISSN:0045-7825
1879-2138
1879-2138
DOI:10.1016/j.cma.2004.02.026