Multilevel Markov Chain Monte Carlo with likelihood scaling for Bayesian inversion with high-resolution observations (2401.15978v2)
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. The convergence analysis focuses on Lipschitz continuous transformations of Gaussian random fields with Mat\'ern covariance structure. 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.
- regularisation methods for finite element model updating. Mechanical Systems and Signal Processing, 12(1):47–64, Jan. 1998.
- One-year operational modal analysis of a steel bridge from high-resolution macrostrain monitoring: Influence of temperature vs. retrofitting | Elsevier Enhanced Reader. Mechanical Systems and Signal Processing, 161:107951, Apr. 2021.
- X-TMCMC: Adaptive kriging for Bayesian inverse modeling. Computer Methods in Applied Mechanics and Engineering, 289:409–428, June 2015.
- G. Baecher and J. Christian. Reliability and Statistics in Geotechnical Engineering. Wiley, 2003.
- p-Refined Multilevel Quasi-Monte Carlo for Galerkin Finite Element methods with applications in civil engineering. Algorithms, 13(5):110–139, 2020.
- J. Charrier. Strong and Weak Error Estimates for Elliptic Partial Differential Equations with Random Coefficients. SIAM Journal on Numerical Analysis, 50(1):216–246, Jan. 2012.
- Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and Its Application to Multilevel Monte Carlo Methods. SIAM Journal on Numerical Analysis, 51(1):322–352, Jan. 2013.
- Dimension-Robust MCMC in Bayesian Inverse Problems. arXiv:1803.03344 [stat], Mar. 2019. arXiv: 1803.03344.
- J. A. Christen and C. Fox. Markov chain Monte Carlo Using an Approximation. Journal of Computational and Graphical Statistics, 14(4):795–810, Dec. 2005.
- Multilevel Monte Carlo methods and applications to elliptic PDEs with random coefficients. Computing and Visualization in Science, 14(1):3–15, Jan. 2011.
- Mcmc methods for functions: Modifying old algorithms to make them faster. Statistical Science, 28(3):424–446, 2013.
- Multilevel Markov Chain Monte Carlo. SIAM Review, 61(3):509–545, Jan. 2019.
- Stabil: An educational Matlab toolbox for static and dynamic structural analysis. Computer Applications in Engineering Education, 29(5):1372–1389, 2021. _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/cae.22391.
- Finite elements for elliptic problems with stochastic coefficients. Comput. Methods Appl. Mech. Engrg., 194:205–228, 2005.
- M. Friswell and J. E. Mottershead. Finite Element Model Updating in Structural Dynamics. Springer Science & Business Media, Mar. 1995. Google-Books-ID: KyMVKqrEEZIC.
- Finite–element model updating using experimental test data: parametrization and regularization. Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 359(1778):169–186, Jan. 2001.
- Stochastic Finite Element Method: Response Statistics. In R. G. Ghanem and P. D. Spanos, editors, Stochastic Finite Elements: A Spectral Approach. Springer, New York, NY, 1991.
- M. B. Giles. Multilevel Monte Carlo Path Simulation. Operations Research, 56(3):607–617, June 2008.
- M. B. Giles. Multilevel Monte Carlo methods. Acta Numerica, 24:259–328, May 2015.
- W. Hastings. Monte-Carlo sampling methods using Markov chains and their applications. Biometrika, 57(1):97–109, 1970.
- 3D digital image correlation methods for full-field vibration measurement. Mechanical Systems and Signal Processing, 25:917–927, 2011.
- Non-planar sensing skins for structural health monitoring based on electrical resistance tomography. Computer-Aided Civil and Infrastructure Engineering, 36:1488–1507, May 2021. arXiv: 2012.04588.
- Global Weak Solutions to the Equations of Compressible Flow of Nematic Liquid Crystals in Two Dimensions. Archive for Rational Mechanics and Analysis, 214(2):403–451, Nov. 2014.
- Multilevel Sequential2 Monte Carlo for Bayesian Inverse Problems. Journal of Computational Physics, 368:154–178, Sept. 2018. arXiv: 1709.09763.
- Data Assimilation: A Mathematical Introduction, volume 62 of Texts in Applied Mathematics. Springer International Publishing, Cham, 2015.
- A multilevel Monte Carlo method for asymptotic-preserving particle schemes in the diffusive limit. Numerische Mathematik, 148(1):141–186, May 2021. arXiv: 1907.04610.
- Density filtering regularization of finite element model updating problems. Mechanical Systems and Signal Processing, 128:282–294, 2019.
- C. Schwab and R. A. Todor. Karhunen–Loève approximation of random fields by generalized fast multipole methods. Journal of Computational Physics, 217(1):100–122, Sept. 2006.
- Uncertainty Quantification in the Assessment of Progressive Damage in a 7-Story Full-Scale Building Slice. Journal of Engineering Mechanics, 139(12):1818–1830, 2013.
- On prediction error correlation in Bayesian model updating. Journal of Sound and Vibration, 332(18):4136–4152, 2013.
- Accounting for modeling errors and inherent structural variability through a hierarchical bayesian model updating approach: An overview. Sensors (Switzerland), 20(14):1–27, 2020.
- Further analysis of multilevel Monte Carlo methods for elliptic PDEs with random coefficients. Numerische Mathematik, 125(3):569–600, Nov. 2013.
- B. Titurus and M. I. Friswell. Regularization in model updating. International Journal for Numerical Methods in Engineering, 75(4):440–478, 2008. _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/nme.2257.
- A. Van Barel and S. Vandewalle. Robust Optimization of PDEs with Random Coefficients Using a Multilevel Monte Carlo Method. SIAM/ASA Journal on Uncertainty Quantification, 7(1):174–202, Nov. 2017. arXiv: 1711.02574.
- Bayesian decision analysis for the optimization of inspection and repair of spatially degrading concrete structures. Engineering Structures, 220(December 2019):111028, 2020.
- Shape features and finite element model updating from full-field strain data. International Journal of Solids and Structures, 48(11-12):1644–1657, 2011.
- Finite element model updating from full-field vibration measurement using digital image correlation. Journal of Sound and Vibration, 330(8):1599–1620, Apr. 2011.
- Consistent regularization of nonlinear model updating for damage identification. Mechanical Systems and Signal Processing, 23(6):1965–1985, Aug. 2009.