Convergence of quasi-optimal sparse-grid approximation of Hilbert-space-valued functions: application to random elliptic PDEs

Handle URI:
http://hdl.handle.net/10754/622185
Title:
Convergence of quasi-optimal sparse-grid approximation of Hilbert-space-valued functions: application to random elliptic PDEs
Authors:
Nobile, F.; Tamellini, L.; Tempone, Raul ( 0000-0003-1967-4446 )
Abstract:
In this work we provide a convergence analysis for the quasi-optimal version of the sparse-grids stochastic collocation method we presented in a previous work: “On the optimal polynomial approximation of stochastic PDEs by Galerkin and collocation methods” (Beck et al., Math Models Methods Appl Sci 22(09), 2012). The construction of a sparse grid is recast into a knapsack problem: a profit is assigned to each hierarchical surplus and only the most profitable ones are added to the sparse grid. The convergence rate of the sparse grid approximation error with respect to the number of points in the grid is then shown to depend on weighted summability properties of the sequence of profits. This is a very general argument that can be applied to sparse grids built with any uni-variate family of points, both nested and non-nested. As an example, we apply such quasi-optimal sparse grids to the solution of a particular elliptic PDE with stochastic diffusion coefficients, namely the “inclusions problem”: we detail the convergence estimates obtained in this case using polynomial interpolation on either nested (Clenshaw–Curtis) or non-nested (Gauss–Legendre) abscissas, verify their sharpness numerically, and compare the performance of the resulting quasi-optimal grids with a few alternative sparse-grid construction schemes recently proposed in the literature.
KAUST Department:
Applied Mathematics and Computational Science Program
Citation:
Nobile F, Tamellini L, Tempone R (2015) Convergence of quasi-optimal sparse-grid approximation of Hilbert-space-valued functions: application to random elliptic PDEs. Numerische Mathematik 134: 343–388. Available: http://dx.doi.org/10.1007/s00211-015-0773-y.
Publisher:
Springer Nature
Journal:
Numerische Mathematik
Issue Date:
30-Oct-2015
DOI:
10.1007/s00211-015-0773-y
Type:
Article
ISSN:
0029-599X; 0945-3245
Sponsors:
The authors would like to recognize the support of King Abdullah University of Science and Technology (KAUST) AEA project "Predictability and Uncertainty Quantification for Models of Porous Media" and University of Texas at Austin AEA Rnd 3 "Uncertainty quantification for predictive mobdeling of the dissolution of porous and fractured media". F. Nobile and L. Tamellini have been partially supported by the Italian grant FIRB-IDEAS (Project n. RBID08223Z) "Advanced numerical techniques for uncertainty quantification in engineering and life science problems" and by the Swiss National Science Foundation under the Project No. 140574 "Efficient numerical methods for flow and transport phenomena in heterogeneous random porous media". They also received partial support from the Center for ADvanced MOdeling Science (CADMOS). R. Tempone is a member of the KAUST SRI Center for Uncertainty Quantification in Computational Science and Engineering. We acknowledge the usage of the Matlab (R) functions patterson_rule.m by J. Burkardt (http://people.sc.fsu.edu/similar to jburkardt/m_src/patterson_rule/patterson_rule.html) for the computation of Gauss-Patterson points and lejapoints. m by M. Caliari (http://profs.sci.univr.it/similar to caliari/software/lejapoints.m) for the computation of symmetrized Leja points.
Appears in Collections:
Articles; Applied Mathematics and Computational Science Program

Full metadata record

DC FieldValue Language
dc.contributor.authorNobile, F.en
dc.contributor.authorTamellini, L.en
dc.contributor.authorTempone, Raulen
dc.date.accessioned2017-01-02T08:42:36Z-
dc.date.available2017-01-02T08:42:36Z-
dc.date.issued2015-10-30en
dc.identifier.citationNobile F, Tamellini L, Tempone R (2015) Convergence of quasi-optimal sparse-grid approximation of Hilbert-space-valued functions: application to random elliptic PDEs. Numerische Mathematik 134: 343–388. Available: http://dx.doi.org/10.1007/s00211-015-0773-y.en
dc.identifier.issn0029-599Xen
dc.identifier.issn0945-3245en
dc.identifier.doi10.1007/s00211-015-0773-yen
dc.identifier.urihttp://hdl.handle.net/10754/622185-
dc.description.abstractIn this work we provide a convergence analysis for the quasi-optimal version of the sparse-grids stochastic collocation method we presented in a previous work: “On the optimal polynomial approximation of stochastic PDEs by Galerkin and collocation methods” (Beck et al., Math Models Methods Appl Sci 22(09), 2012). The construction of a sparse grid is recast into a knapsack problem: a profit is assigned to each hierarchical surplus and only the most profitable ones are added to the sparse grid. The convergence rate of the sparse grid approximation error with respect to the number of points in the grid is then shown to depend on weighted summability properties of the sequence of profits. This is a very general argument that can be applied to sparse grids built with any uni-variate family of points, both nested and non-nested. As an example, we apply such quasi-optimal sparse grids to the solution of a particular elliptic PDE with stochastic diffusion coefficients, namely the “inclusions problem”: we detail the convergence estimates obtained in this case using polynomial interpolation on either nested (Clenshaw–Curtis) or non-nested (Gauss–Legendre) abscissas, verify their sharpness numerically, and compare the performance of the resulting quasi-optimal grids with a few alternative sparse-grid construction schemes recently proposed in the literature.en
dc.description.sponsorshipThe authors would like to recognize the support of King Abdullah University of Science and Technology (KAUST) AEA project "Predictability and Uncertainty Quantification for Models of Porous Media" and University of Texas at Austin AEA Rnd 3 "Uncertainty quantification for predictive mobdeling of the dissolution of porous and fractured media". F. Nobile and L. Tamellini have been partially supported by the Italian grant FIRB-IDEAS (Project n. RBID08223Z) "Advanced numerical techniques for uncertainty quantification in engineering and life science problems" and by the Swiss National Science Foundation under the Project No. 140574 "Efficient numerical methods for flow and transport phenomena in heterogeneous random porous media". They also received partial support from the Center for ADvanced MOdeling Science (CADMOS). R. Tempone is a member of the KAUST SRI Center for Uncertainty Quantification in Computational Science and Engineering. We acknowledge the usage of the Matlab (R) functions patterson_rule.m by J. Burkardt (http://people.sc.fsu.edu/similar to jburkardt/m_src/patterson_rule/patterson_rule.html) for the computation of Gauss-Patterson points and lejapoints. m by M. Caliari (http://profs.sci.univr.it/similar to caliari/software/lejapoints.m) for the computation of symmetrized Leja points.en
dc.publisherSpringer Natureen
dc.subject41A10en
dc.subject65C20en
dc.subject65N12en
dc.subject65N30en
dc.titleConvergence of quasi-optimal sparse-grid approximation of Hilbert-space-valued functions: application to random elliptic PDEsen
dc.typeArticleen
dc.contributor.departmentApplied Mathematics and Computational Science Programen
dc.identifier.journalNumerische Mathematiken
dc.contributor.institutionCSQI-MATHICSE, École Polytechnique Fédérale Lausanne, Station 8, Lausanne, 1015, Switzerlanden
dc.contributor.institutionDipartimento di Matematica, Università di Pavia, via Ferrata 1, Pavia, 27100, Italyen
kaust.authorTempone, Raulen
All Items in KAUST are protected by copyright, with all rights reserved, unless otherwise indicated.