An Adaptive Sparse Grid Algorithm for Elliptic PDEs with Lognormal Diffusion Coefficient

Fabio Nobile, Lorenzo Tamellini, Francesco Tesei, Raul Tempone

Research output: Chapter in Book/Report/Conference proceedingConference contribution

29 Scopus citations


In this work we build on the classical adaptive sparse grid algorithm (T. Gerstner and M. Griebel, Dimension-adaptive tensor-product quadrature), obtaining an enhanced version capable of using non-nested collocation points, and supporting quadrature and interpolation on unbounded sets. We also consider several profit indicators that are suitable to drive the adaptation process. We then use such algorithm to solve an important test case in Uncertainty Quantification problem, namely the Darcy equation with lognormal permeability random field, and compare the results with those obtained with the quasi-optimal sparse grids based on profit estimates, which we have proposed in our previous works (cf. e.g. Convergence of quasi-optimal sparse grids approximation of Hilbert-valued functions: application to random elliptic PDEs). To treat the case of rough permeability fields, in which a sparse grid approach may not be suitable, we propose to use the adaptive sparse grid quadrature as a control variate in a Monte Carlo simulation. Numerical results show that the adaptive sparse grids have performances similar to those of the quasi-optimal sparse grids and are very effective in the case of smooth permeability fields. Moreover, their use as control variate in a Monte Carlo simulation allows to tackle efficiently also problems with rough coefficients, significantly improving the performances of a standard Monte Carlo scheme.
Original languageEnglish (US)
Title of host publicationLecture Notes in Computational Science and Engineering
PublisherSpringer Nature
Number of pages30
ISBN (Print)9783319282602
StatePublished - Mar 17 2016


Dive into the research topics of 'An Adaptive Sparse Grid Algorithm for Elliptic PDEs with Lognormal Diffusion Coefficient'. Together they form a unique fingerprint.

Cite this