Multilevel Markov Chain Monte Carlo with likelihood scaling for high-resolution data assimilation
Abstract
We propose a multilevel Markov chain Monte Carlo (MCMC) method for the Bayesian inference of random field parameters in PDEs using high-resolution data. Compared to existing multilevel MCMC methods, we additionally consider level-dependent data resolution and introduce a suitable likelihood scaling to enable consistent cross-level comparisons. We theoretically show that this approach attains the same convergence rates as when using level-independent treatment of data, but at significantly reduced computational cost. Additionally, we show that assumptions of exponential covariance and log-normality of random fields, widely held in multilevel Monte Carlo literature, can be extended to a wide range of covariance structures and random fields. These results are illustrated using numerical experiments for a 2D plane stress problem, where the Young's modulus is estimated from discretisations of the displacement field.
- Publication:
-
arXiv e-prints
- Pub Date:
- January 2024
- DOI:
- 10.48550/arXiv.2401.15978
- arXiv:
- arXiv:2401.15978
- Bibcode:
- 2024arXiv240115978V
- Keywords:
-
- Mathematics - Numerical Analysis;
- 62F15;
- 35R60;
- 65C40 (Primary);
- 62M05;
- 65C05;
- 65N30 (Secondary)
- E-Print:
- 25 pages, 7 figures