On moment estimation from polynomial chaos expansion models

Polynomial Chaos Expansions (PCEs) offer an efficient alternative to assess the statistical properties of a model output taking into account the statistical properties of several uncertain model inputs, particularly, under the restriction of probing the forward model as little as possible. The use of PCE knows a steady increase in recent literature with applications spanning from mere system analysis to robust control design and optimization. The principle idea is to model the output as a series expansion, with the expansion basis depending only on the stochastic variables. The basis is then chosen so that it is better suited to support the uncertainty propagation. Once the expansion coefficients have been identified, accessing the statistical moments benefits the linearity of the expansion and the properties of the moment operator. As a result, the first two statistical moments of the model output can be computed easily using well known analytical expressions. For high-order moments, analytical expressions also exist but are inefficient to evaluate. In this letter we present three strategies to efficiently calculate high-order moments from the expansion model and provide an empirical study of the associated computation time supporting any potential users in making an informed choice