Computation of higher-order moments of generalized polynomial chaos expansions
Autor: | Béatrice Faverjon, Éric Savin |
---|---|
Rok vydání: | 2017 |
Předmět: |
Numerical Analysis
Polynomial Polynomial chaos Applied Mathematics 010102 general mathematics Gauss MathematicsofComputing_NUMERICALANALYSIS General Engineering 01 natural sciences Askey scheme 010101 applied mathematics ComputingMethodologies_SYMBOLICANDALGEBRAICMANIPULATION Orthogonal polynomials Applied mathematics 0101 mathematics Uncertainty quantification Polynomial expansion Mathematics Taylor expansions for the moments of functions of random variables |
Zdroj: | International Journal for Numerical Methods in Engineering. 111:1192-1200 |
ISSN: | 0029-5981 |
DOI: | 10.1002/nme.5505 |
Popis: | Because of the complexity of fluid flow solvers, non-intrusive uncertainty quantification techniques have been developed in aerodynamic simulations in order to compute the quantities of interest required in an optimization process, for example. The objective function is commonly expressed in terms of moments of these quantities, such as the mean, standard deviation, or even higher-order moments. Polynomial surrogate models based on polynomial chaos expansions have often been implemented in this respect. The original approach of uncertainty quantification using polynomial chaos is however intrusive. It is based on a Galerkin-type formulation of the model equations to derive the governing equations for the polynomial expansion coefficients. Third-order, indeed fourth-order moments of the polynomials are needed in this analysis. Besides, both intrusive and non-intrusive approaches call for their computation provided that higher-order moments of the quantities of interest need be post-processed. In most applications they are evaluated by Gauss quadratures, and eventually stored for use throughout the computations. In this paper analytical formulas are rather considered for the moments of the continuous polynomials of the Askey scheme, so that they can be evaluated by quadrature-free procedures instead. Matlab codes have been developed for this purpose and tested by comparisons with Gauss quadratures. |
Databáze: | OpenAIRE |
Externí odkaz: |