In this work we focus on the numerical approximation of the solution $u$ of a linear elliptic PDE with stochastic coefficients. The problem is rewritten as a parametric PDE and the functional dependence of the solution on the parameters is approximated by multivariate polynomials. We first consider the Stochastic Galerkin method, and rely on sharp estimates for the decay of the Fourier coefficients of the spectral expansion of $u$ on an orthogonal polynomial basis to build a sequence of polynomial subspaces that features better convergence properties, in terms of error versus number of degrees of freedom, than standard choices such as Total Degree or Tensor Product subspaces. We consider then the Stochastic Collocation method, and use the previous estimates to introduce a new class of Sparse Grids, based on the idea of selecting a priori the most profitable hierarchical surpluses, that, again, features better convergence properties compared to standard Smolyak or tensor product grids. Numerical results show the effectiveness of the newly introduced polynomial spaces and sparse grids.

On the optimal polynomial approximation of Stochastic PDEs by Galerkin and Collocation methods

NOBILE, FABIO;TAMELLINI, LORENZO
2011-01-01

Abstract

In this work we focus on the numerical approximation of the solution $u$ of a linear elliptic PDE with stochastic coefficients. The problem is rewritten as a parametric PDE and the functional dependence of the solution on the parameters is approximated by multivariate polynomials. We first consider the Stochastic Galerkin method, and rely on sharp estimates for the decay of the Fourier coefficients of the spectral expansion of $u$ on an orthogonal polynomial basis to build a sequence of polynomial subspaces that features better convergence properties, in terms of error versus number of degrees of freedom, than standard choices such as Total Degree or Tensor Product subspaces. We consider then the Stochastic Collocation method, and use the previous estimates to introduce a new class of Sparse Grids, based on the idea of selecting a priori the most profitable hierarchical surpluses, that, again, features better convergence properties compared to standard Smolyak or tensor product grids. Numerical results show the effectiveness of the newly introduced polynomial spaces and sparse grids.
2011
Uncertainty Quantification; PDEs with random data; elliptic equations; multivariate polynomial approximation; best M terms approximation; Stochastic Galerkin methods; Smolyak approximation; Sparse grids; Stochastic Collocation methods.
File in questo prodotto:
File Dimensione Formato  
23-2011.pdf

Accesso riservato

: Pre-Print (o Pre-Refereeing)
Dimensione 640.36 kB
Formato Adobe PDF
640.36 kB Adobe PDF   Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11311/622567
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? ND
social impact