On the Selection of Random Field Evaluation Points in the p-MLQMC Method

by   Philippe Blondeel, et al.

Engineering problems are often characterized by significant uncertainty in their material parameters. A typical example coming from geotechnical engineering is the slope stability problem where the soil's cohesion is modeled as a random field. An efficient manner to account for this uncertainty is the novel sampling method called p-refined Multilevel Quasi-Monte Carlo (p-MLQMC). The p-MLQMC method uses a hierarchy of p-refined Finite Element meshes combined with a deterministic Quasi-Monte Carlo sampling rule. This combination yields a significant computational cost reduction with respect to classic Multilevel Monte Carlo. However, in previous work, not enough consideration was given how to incorporate the uncertainty, modeled as a random field, in the Finite Element model with the p-MLQMC method. In the present work we investigate how this can be adequately achieved by means of the integration point method. We therefore investigate how the evaluation points of the random field are to be selected in order to obtain a variance reduction over the levels. We consider three different approaches. These approaches will be benchmarked on a slope stability problem in terms of computational runtime. We find that the Local Nested Approach yields a speedup up to a factor five with respect to the Non-Nested approach.


page 1

page 2

page 3

page 4


h- and p-refined Multilevel Monte Carlo Methods for Uncertainty Quantification in Structural Engineering

Practical structural engineering problems are often characterized by sig...

Quantifying uncertainties on excursion sets under a Gaussian random field prior

We focus on the problem of estimating and quantifying uncertainties on t...

Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling

When solving partial differential equations with random fields as coeffi...

Multilevel Monte Carlo for uncertainty quantification in structural engineering

Practical structural engineering problems often exhibit a significant de...

Space-time multilevel quadrature methods and their application for cardiac electrophysiology

We present a novel approach which aims at high-performance uncertainty q...

Enhanced Multi-Index Monte Carlo by means of Multiple Semi-Coarsened Multigrid for Anisotropic Diffusion Problems

In many models used in engineering and science, material properties are ...

Space-time shape uncertainties in the forward and inverse problem of electrocardiography

In electrocardiography, the "classic" inverse problem consists of findin...

Please sign up or login with your details

Forgot password? Click here to reset