Exploiting locality in sparse polynomial approximation of parametric elliptic PDEs and application to parameterized domains

Institute for Mathematics, Astrophysics and Particle Physics, Radboud University, Nijmegen, The Netherlands

^* Corresponding author: w.vanharten@science.ru.nl

Received: 11 August 2023
Accepted: 15 June 2024

Abstract

This work studies how the choice of the representation for parametric, spatially distributed inputs to elliptic partial differential equations (PDEs) affects the efficiency of a polynomial surrogate, based on Taylor expansion, for the parameter-to-solution map. In particular, we show potential advantages of representations using functions with localized supports. As model problem, we consider the steady-state diffusion equation, where the diffusion coefficient and right-hand side depend smoothly but potentially in a highly nonlinear way on a parameter y ∈ [−1, 1]^N. Following previous work for affine parameter dependence and for the lognormal case, we use pointwise instead of norm-wise bounds to prove ℓ^p-summability of the Taylor coefficients of the solution. As application, we consider surrogates for solutions to elliptic PDEs on parametric domains. Using a mapping to a nominal configuration, this case fits in the general framework, and higher convergence rates can be attained when modeling the parametric boundary via spatially localized functions. The theoretical results are supported by numerical experiments for the parametric domain problem, illustrating the efficiency of the proposed approach and providing further insight on numerical aspects. Although the methods and ideas are carried out for the steady-state diffusion equation, they extend easily to other elliptic and parabolic PDEs.