Premium
Nonparametric density estimation for randomly perturbed elliptic problems II: Applications and adaptive modeling
Author(s) -
Estep D.,
Målqvist A.,
Tavener S.
Publication year - 2009
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.2547
Subject(s) - nonparametric statistics , partial differential equation , elliptic partial differential equation , algorithm , mathematical optimization , computer science , mathematics , probability density function , statistics , mathematical analysis
In this paper, we develop and apply an efficient adaptive algorithm for computing the propagation of uncertainty into a quantity of interest computed from numerical solutions of an elliptic partial differential equation with a randomly perturbed diffusion coefficient. The algorithm is well‐suited for problems for which limited information about the random perturbations is available and when an approximation of the probability distribution of the output is desired. We employ a nonparametric density estimation approach based on a very efficient method for computing random samples of elliptic problems described and analyzed in ( SIAM J. Sci. Comput. 2008. DOI: JCOMP‐D‐08‐00261 ). We fully develop the adaptive algorithm suggested by the analysis in that paper, discuss details of its implementation, and illustrate its behavior using a realistic data set. Finally, we extend the analysis to include a ‘modeling error’ term that accounts for the effects of the resolution of the statistical description of the random variation. We modify the adaptive algorithm to adapt the resolution of the statistical description and illustrate the behavior of the adaptive algorithm in several examples. Copyright © 2009 John Wiley & Sons, Ltd.