# Approximation of Quantities of Interest in Stochastic PDEs by the Random Discrete L^2 Projection on Polynomial Spaces

- Handle URI:
- http://hdl.handle.net/10754/555667
- Title:
- Approximation of Quantities of Interest in Stochastic PDEs by the Random Discrete L^2 Projection on Polynomial Spaces
- Authors:
- Abstract:
- In this work we consider the random discrete L^2 projection on polynomial spaces (hereafter RDP) for the approximation of scalar quantities of interest (QOIs) related to the solution of a partial differential equation model with random input parameters. In the RDP technique the QOI is first computed for independent samples of the random input parameters, as in a standard Monte Carlo approach, and then the QOI is approximated by a multivariate polynomial function of the input parameters using a discrete least squares approach. We consider several examples including the Darcy equations with random permeability, the linear elasticity equations with random elastic coefficient, and the Navier--Stokes equations in random geometries and with random fluid viscosity. We show that the RDP technique is well suited to QOIs that depend smoothly on a moderate number of random parameters. Our numerical tests confirm the theoretical findings in [G. Migliorati, F. Nobile, E. von Schwerin, and R. Tempone, Analysis of the Discrete $L^2$ Projection on Polynomial Spaces with Random Evaluations, MOX report 46-2011, Politecnico di Milano, Milano, Italy, submitted], which have shown that, in the case of a single uniformly distributed random parameter, the RDP technique is stable and optimally convergent if the number of sampling points is proportional to the square of the dimension of the polynomial space. Here optimality means that the weighted $L^2$ norm of the RDP error is bounded from above by the best $L^\infty$ error achievable in the given polynomial space, up to logarithmic factors. In the case of several random input parameters, the numerical evidence indicates that the condition on quadratic growth of the number of sampling points could be relaxed to a linear growth and still achieve stable and optimal convergence. This makes the RDP technique very promising for moderately high dimensional uncertainty quantification.
- KAUST Department:
- Citation:
- Approximation of Quantities of Interest in Stochastic PDEs by the Random Discrete L^2 Projection on Polynomial Spaces 2013, 35 (3):A1440 SIAM Journal on Scientific Computing
- Journal:
- Issue Date:
- 30-May-2013
- DOI:
- 10.1137/120897109
- Type:
- Article
- ISSN:
- 1064-8275; 1095-7197
- Additional Links:
- http://epubs.siam.org/doi/abs/10.1137/120897109

- Appears in Collections:
- Articles; Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division; Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division

# Full metadata record

DC Field | Value | Language |
---|---|---|

dc.contributor.author | Migliorati, G. | en |

dc.contributor.author | Nobile, F. | en |

dc.contributor.author | von Schwerin, E. | en |

dc.contributor.author | Tempone, Raul | en |

dc.date.accessioned | 2015-05-25T08:31:18Z | en |

dc.date.available | 2015-05-25T08:31:18Z | en |

dc.date.issued | 2013-05-30 | en |

dc.identifier.citation | Approximation of Quantities of Interest in Stochastic PDEs by the Random Discrete L^2 Projection on Polynomial Spaces 2013, 35 (3):A1440 SIAM Journal on Scientific Computing | en |

dc.identifier.issn | 1064-8275 | en |

dc.identifier.issn | 1095-7197 | en |

dc.identifier.doi | 10.1137/120897109 | en |

dc.identifier.uri | http://hdl.handle.net/10754/555667 | en |

dc.description.abstract | In this work we consider the random discrete L^2 projection on polynomial spaces (hereafter RDP) for the approximation of scalar quantities of interest (QOIs) related to the solution of a partial differential equation model with random input parameters. In the RDP technique the QOI is first computed for independent samples of the random input parameters, as in a standard Monte Carlo approach, and then the QOI is approximated by a multivariate polynomial function of the input parameters using a discrete least squares approach. We consider several examples including the Darcy equations with random permeability, the linear elasticity equations with random elastic coefficient, and the Navier--Stokes equations in random geometries and with random fluid viscosity. We show that the RDP technique is well suited to QOIs that depend smoothly on a moderate number of random parameters. Our numerical tests confirm the theoretical findings in [G. Migliorati, F. Nobile, E. von Schwerin, and R. Tempone, Analysis of the Discrete $L^2$ Projection on Polynomial Spaces with Random Evaluations, MOX report 46-2011, Politecnico di Milano, Milano, Italy, submitted], which have shown that, in the case of a single uniformly distributed random parameter, the RDP technique is stable and optimally convergent if the number of sampling points is proportional to the square of the dimension of the polynomial space. Here optimality means that the weighted $L^2$ norm of the RDP error is bounded from above by the best $L^\infty$ error achievable in the given polynomial space, up to logarithmic factors. In the case of several random input parameters, the numerical evidence indicates that the condition on quadratic growth of the number of sampling points could be relaxed to a linear growth and still achieve stable and optimal convergence. This makes the RDP technique very promising for moderately high dimensional uncertainty quantification. | en |

dc.relation.url | http://epubs.siam.org/doi/abs/10.1137/120897109 | en |

dc.rights | Archived with thanks to SIAM Journal on Scientific Computing | en |

dc.subject | PDE stochastic data | en |

dc.subject | discrete least squares | en |

dc.subject | polynomial approximation | en |

dc.title | en | |

dc.type | Article | en |

dc.contributor.department | Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division | en |

dc.contributor.department | Center for Uncertainty Quantification in Computational Science and Engineering (SRI-UQ) | en |

dc.identifier.journal | SIAM Journal on Scientific Computing | en |

dc.eprint.version | Publisher's Version/PDF | en |

dc.contributor.institution | MATHICSE-CSQI,École Polytechnique Fédérale de Lausanne, Lausanne CH-1015, Switzerland, and MOX-Dipartimento di Matematica “Francesco Brioschi,” Politecnico di Milano, Milano 20133, Italy | en |

dc.contributor.institution | MOX-Dipartimento di Matematica “Francesco Brioschi,” Politecnico di Milano, Milano 20133, Italy | en |

kaust.author | von Schwerin, Erik | en |

kaust.author | Tempone, Raul | en |

All Items in KAUST are protected by copyright, with all rights reserved, unless otherwise indicated.