Ensemble Kalman filter for multiscale inverse problems

Abstract

We present a novel algorithm based on the ensemble Kalman filter to solve inverse problems involving multiscale elliptic partial differential equations. Our method is based on numerical homogenization and finite element discretization and allows to recover a highly oscillatory tensor from measurements of the multiscale solution in a computationally inexpensive manner. The properties of the approximate solution are analysed with respect to the multiscale and discretization parameters, and a convergence result is shown to hold. A reinterpretation of the solution from a Bayesian perspective is provided, and convergence of the approximate conditional posterior distribution is proved with respect to the Wasserstein distance. A numerical experiment validates our methodology, with a particular emphasis on modelling error and computational cost.

References in corpus (2)

• The Ensemble Kalman Filter for Inverse Problems
• Iterative Updating of Model Error for Bayesian Inversion

Cited by in corpus (4)

• A probabilistic finite element method based on random meshes: Error estimators and Bayesian inverse problems
• Physics-informed neural networks for learning the homogenized coefficients of multiscale elliptic equations
• Fourier Series-Based Approximation of Time-Varying Parameters in Ordinary Differential Equations
• On the performance of multi-fidelity and reduced-dimensional neural emulators for inference of physiological boundary conditions