UCL Discovery
UCL home » Library Services » Electronic resources » UCL Discovery

Adaptive Multilevel Subset Simulation with Selective Refinement

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. Green open access

[thumbnail of Diaz De La O_22m1515240.pdf]
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
Downloads since deposit
18Downloads
Download activity - last month
Download activity - last 12 months
Downloads by country - last 12 months

Archive Staff Only

View Item View Item