Elfverson, D;
Scheichl, R;
Weissmann, S;
Diaz De La O, FA;
(2024)
Adaptive Multilevel Subset Simulation with Selective Refinement.
SIAM/ASA Journal on Uncertainty Quantification
, 12
(3)
pp. 932-963.
10.1137/22M1515240.
Preview |
Text
Diaz De La O_22m1515240.pdf Download (3MB) | Preview |
Abstract
In this work we propose an adaptive multilevel version of subset simulation to estimate the probability of rare events for complex physical systems. Given a sequence of nested failure domains of increasing size, the rare event probability is expressed as a product of conditional probabilities. The proposed new estimator uses different model resolutions and varying numbers of samples across the hierarchy of nested failure sets. In order to dramatically reduce the computational cost, we construct the intermediate failure sets such that only a small number of expensive high-resolution model evaluations are needed, whilst the majority of samples can be taken from inexpensive low-resolution simulations. A key idea in our new estimator is the use of a posteriori error estimators combined with a selective mesh refinement strategy to guarantee the critical subset property that may be violated when changing model resolution from one failure set to the next. The efficiency gains and the statistical properties of the estimator are investigated both theoretically via shaking transformations, as well as numerically. On a model problem from subsurface flow, the new multilevel estimator achieves gains of more than a factor 60 over standard subset simulation for a practically relevant relative error of 25%.
Type: | Article |
---|---|
Title: | Adaptive Multilevel Subset Simulation with Selective Refinement |
Open access status: | An open access version is available from UCL Discovery |
DOI: | 10.1137/22M1515240 |
Publisher version: | https://doi.org/10.1137/22M1515240 |
Language: | English |
Additional information: | This version is the version of record. For information on re-use, please refer to the publisher’s terms and conditions. |
Keywords: | rare event probabilities, adaptive model hierarchies, high-dimensional problems, Markov chain Monte Carlo, shaking transformations |
UCL classification: | UCL UCL > Provost and Vice Provost Offices > UCL BEAMS UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences > Dept of Mathematics UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences > Dept of Mathematics > Clinical Operational Research Unit |
URI: | https://discovery.ucl.ac.uk/id/eprint/10194909 |
Archive Staff Only
View Item |