Stochastic spectral methods are numerical techniques for approximating solutions to partial differential equations with random parameters. In this work, we present and examine the multi-element probabilistic collocation method (ME-PCM), which is a generalized form of the probabilistic collocation me
Multi-element probabilistic collocation method in high dimensions
β Scribed by Jasmine Foo; George Em Karniadakis
- Publisher
- Elsevier Science
- Year
- 2010
- Tongue
- English
- Weight
- 1021 KB
- Volume
- 229
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
β¦ Synopsis
We combine multi-element polynomial chaos with analysis of variance (ANOVA) functional decomposition to enhance the convergence rate of polynomial chaos in high dimensions and in problems with low stochastic regularity. Specifically, we employ the multi-element probabilistic collocation method MEPCM [1] and so we refer to the new method as MEPCM-A. We investigate the dependence of the convergence of MEPCM-A on two decomposition parameters, the polynomial order l and the effective dimension m, with m ( N, and N the nominal dimension. Numerical tests for multi-dimensional integration and for stochastic elliptic problems suggest that m P l for monotonic convergence of the method. We also employ MEPCM-A to obtain error bars for the piezometric head at the Hanford nuclear waste site under stochastic hydraulic conductivity conditions. Finally, we compare the cost of MEPCM-A against Monte Carlo in several hundred dimensions, and we find MEPCM-A to be more efficient for up to 600 dimensions for a specific multi-dimensional integration problem involving a discontinuous function.
π SIMILAR VOLUMES
This paper outlines the use of non-conforming (discontinuous) elements in the collocation boundary element method for solving two-dimensional potential and Poisson type problems. The roots of an orthogonal polynomial (shifted Jacobi polynomial) are used as the collocation points. This results in inc
This paper presents a multi-dimensional particle tracking technique for applying the Lagrangian-Eulerian finite element method to solve transport equations in transient-state simulations. In the Lagrangian-Eulerian approach, the advection term is handled in the Lagrangian step so that the associated
The presence of singularities in the integral operators of the boundary element methods requires that the density functions must satisfy certain continuity requirements if the displacements and stresses are to be bounded. Quite often the continuity conditions, particularly on the derivatives of the
A research code has been written to solve an elliptic system of coupled non-linear partial differential equations of conservation form on a rectangularly shaped three-dimensional domain. The code uses the method of collocation of Gauss points with tricubic Hermite piecewise continuous polynomial bas