Thursday, January 13, 2022

UQSay #41

The forty-first UQSay seminar on UQ, DACE and related topics will take place online on Thursday afternoon, January 20, 2022.

2–3 PM — Nora Lüthen (ETH Zürich) — [slides]


Poincaré chaos expansions for derivative-enhanced surrogate modelling and sensitivity analysis

Variance-based global sensitivity analysis, and in particular Sobol' analysis, is widely adopted to determine the importance of input variables to a computational model. Sobol' indices can be computed cheaply based on spectral methods like polynomial chaos expansions (PCE). Another option is given by the recently developed Poincaré chaos expansions (PoinCE), whose orthonormal tensor-product basis is generated from the eigenfunctions of one-dimensional Poincaré differential operators.

The Poincaré differential operator is a special case of Sturm-Liouville operator and has recently been revisited for sensitivity analysis (Roustant et al. 2017). Solving the associated eigenproblem yields the Poincaré constant for a large class of one-dimensional measures with bounded support. The associated eigenfunctions form an orthonormal basis with the special (and characterizing) property that derivatives of the basis form again an orthogonal basis with respect to the same measure (Lüthen et al. 2021). The expansion of a model in terms of this basis allows the analytical computation of Sobol’ indices and derivative-based sensitivity indices (DGSM) directly from the expansion coefficients.

Furthermore, the special property of the derivatives makes PoinCE particularly well suited to account for derivative information in the computation of sensitivity indices (Roustant et al. 2020). Indeed the expansions involving either model or derivative evaluations are connected, and computations can be reused. Assuming that partial derivative evaluations of the computational model are available, we compute spectral expansions in terms of Poincaré basis functions or basis partial derivatives, respectively, by sparse regression. We show on numerical examples that the derivative-based expansions provide accurate estimates for Sobol' indices, even outperforming PCE in terms of bias and variance, and explore the performance of PoinCE as a surrogate model.

Joint work with Olivier Roustant, Fabrice Gamboa, Bertrand Iooss, Stefano Marelli and Bruno Sudret.

Ref: arXiv:2107.00394.

Organizing committee: Pierre Barbillon (MIA-Paris), Julien Bect (L2S), Nicolas Bousquet (EDF R&D), Amélie Fau (LMPS), Filippo Gatti (LMPS), Bertrand Iooss (EDF R&D), Alexandre Janon (LMO), Sidonie Lefebvre (ONERA), Didier Lucor (LISN), Emmanuel Vazquez (L2S).

Coordinator: Julien Bect (L2S).

Practical details: the seminar will be held online using Microsoft Teams.

If you want to attend this seminar (or any of the forthcoming online UQSay seminars), and if you do not already have access to the UQSay group on Teams, simply send an email and you will be invited. Please specify which email address the invitation must be sent to (this has to be the address associated with your Teams account).

You will find the link to the seminar on the "General" UQSay channel on Teams, approximately 15 minutes before the beginning.

The technical side of things: you can use Teams either directly from your web browser or using the "fat client", which is available for most platforms (Windows, Linux, Mac, Android & iOS). We strongly recommend the latter option whenever possible. Please give it a try before the seminar to anticipate potential problems.