Dinh DÅ©ng

Dinh DÅ©ng

We establish a sparsity in terms of $\ell_p$-summability and weighted $\ell_2$-summability for the coefficients of the Laguerre generalized piecewise-polynomial chaos expansion of solutions to parametric elliptic PDEs with log-Laplace random inputs. From the sparsity, we derive convergence rates for semi-discrete approximations with respect to parametric variables. These rates are valid for sparse-grid, piecewise-polynomial interpolations and the generated quadratures, and to related extended least-squares approximations and generated quadratures.

Dinh Dũng, Van Kien Nguyen, Viet Ha Hoang

We propose a novel method for establishing the sparsity of the coefficients of the Laguerre generalized polynomial chaos expansion of solutions to parametric elliptic PDEs with log-gamma inputs on $\mathbb{R}_+^\infty$. The established sparsity is quantified by $\ell_p$-summability and weighted $\ell_2$-summability of the coefficients. Building on these sparsity results, we derive convergence rates for semi-discrete approximations in the parametric variables. These rates apply to sparse-grid polynomial interpolations, extended least-squares approximations and the associated semi-discrete quadrature rules. Moreover, a counterpart of our method for parametric elliptic PDEs with log-normal inputs yields a significant improvement in the sufficient condition for $\ell_p$-summability when the component functions in the log-normal representation of the parametric diffusion coefficients have global support, compared with results obtained in prior works.

Felix Bartel, Dinh DÅ©ng

We prove convergence rates of linear sampling recovery of functions in abstract Bochner spaces satisfying weighted summability of their generalized polynomial chaos expansion coefficients. The underlying algorithm is a function-valued extension of the least squares method widely used and thoroughly studied in scalar-valued function recovery. We apply our theory to two core problems in Computational Uncertainty Quantification. First, we address non-intrusive approximations of solutions to parametric elliptic or parabolic PDEs with log-normal or affine inputs using a finite set of particular solvers. Second, we consider approximating infinite-dimensional holomorphic functions that arise as solutions to more general parametric PDEs with Gaussian random field inputs. Our framework allows a unified treatment of the log-normal and affine input models and yields substantial improvements in the state of the art of these problems. More specifically, we obtain convergence rates that improve known results by a polynomial factor in the log-normal case and by a logarithmic factor in the affine case.

Dinh Dũng

We establish convergence rates for a fully discrete, multi-level, linear collocation method solving parametric elliptic PDEs on bounded polygonal domains with log-normal inputs. The method uses a finite set of function evaluations in the spatial-parametric domain. Compared with the best-known fully discrete collocation rates, these rates are significantly improved and, up to logarithmic factors, match the rates of best n-term approximations. The results follow from applying general multi-level linear sampling recovery theory in abstract Bochner spaces -- via extended least-squares -- to infinite-dimensional holomorphic functions. The abstract multi-level recovery in Bochner spaces guarantees yield the improved rates when specialized to the parametric PDE setting.

Dinh DÅ©ng, Helmut Harbrecht, Van Kien Nguyen et al.

In this article, we consider the solution to elliptic diffusion problems on a class of random domains obtained by log-Gaussian random homothety of the unit disk respectively an annulus. We model the problem under consideration and verify the existence and uniqueness of the random solution by path-wise pullback to the nominal unit disk respectively annulus. We prove the analytic regularity of the solution with respect to the random input parameter. We consider the numerical approximation of the random diffusion problem by means of continuous, piecewise linear Lagrangian Galerkin Finite Elements with numerical quadrature in the nominal domain, and by sparse grid interpolation and quadrature of Gauss-Hermite Smolyak and Quasi-Monte Carlo type in the parameter domain. The theoretical findings are complemented by numerical results.