Bayesian inverse problems for recovering coefficients of two scale elliptic equations

DOI10.1088/1361-6420/AAFCD6MaRDI QIDQ5382463zbMATH OpenOpenAlexWikidataFDO

Authors Viet Ha Hoang, Jia Hao Quek

Publication date 17 June 2019

Published in Inverse Problems (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1807.06258

zbMATH Keywords

Bayesian inverse problemscomplexity reductionhomogenizationMarkov chain Monte Carlomultiscale problems

Mathematics Subject Classification ID

Homogenization in context of PDEs; PDEs in media with periodic structure (35B27) Bayesian inference (62F15) Numerical analysis or methods applied to Markov chains (65C40) Numerical methods for initial value problems involving ordinary differential equations (65L05) Numerical methods for inverse problems for boundary value problems involving PDEs (65N21) Complexity and performance of numerical algorithms (65Y20)

Abstract: We consider the Bayesian inverse homogenization problem of recovering the locally periodic two scale coefficient of a two scale elliptic equation, given limited noisy information on the solution. We consider both the uniform and the Gaussian prior probability measures. We use the two scale homogenized equation whose solution contains the solution of the homogenized equation which describes the macroscopic behaviour, and the corrector which encodes the microscopic behaviour. We approximate the posterior probability by a probability measure determined by the solution of the two scale homogenized equation. We show that the Hellinger distance of these measures converges to zero when the microscale converges to zero, and establish an explicit convergence rate when the solution of the two scale homogenized equation is sufficiently regular. Sampling the posterior measure by Markov Chain Monte Carlo (MCMC) method, instead of solving the two scale equation using fine mesh for each proposal with extremely high cost, we can solve the macroscopic two scale homogenized equation. Although this equation is posed in a high dimensional tensorized domain, it can be solved with essentially optimal complexity by the sparse tensor product finite element method, which reduces the computational complexity of the MCMC sampling method substantially. We show numerically that observations on the macrosopic behaviour alone are not sufficient to infer the microstructure. We need also observations on the corrector. Solving the two scale homogenized equation, we get both the solution to the homogenized equation and the corrector. Thus our method is particularly suitable for sampling the posterior measure of two scale coefficients.

Recommendations

• A Bayesian Numerical Homogenization Method for Elliptic Multiscale Inverse Problems
• Numerical methods for multiscale inverse problems
• Numerical homogenization and model order reduction for multiscale inverse problems
• The Bayesian formulation and well-posedness of fractional elliptic inverse problems
• Fine scale uncertainty in parameter estimation for elliptic equations

Cites work

• A General Convergence Result for a Functional Related to the Theory of Homogenization
• Analytic regularity and polynomial approximation of stochastic, parametric elliptic multiscale PDEs
• Approximating infinity-dimensional stochastic Darcy's equations without uniform ellipticity
• Bayesian inverse problems for functions and applications to fluid mechanics
• Complexity analysis of accelerated MCMC methods for Bayesian inversion
• High-Dimensional Finite Elements for Elliptic Problems with Multiple Scales
• Homogenization and Two-Scale Convergence
• scientific article; zbMATH DE number 4138299 (Why is no real title available?)
• scientific article; zbMATH DE number 3948655 (Why is no real title available?)
• scientific article; zbMATH DE number 4011016 (Why is no real title available?)
• scientific article; zbMATH DE number 3596197 (Why is no real title available?)
• scientific article; zbMATH DE number 3627589 (Why is no real title available?)
• scientific article; zbMATH DE number 651015 (Why is no real title available?)
• scientific article; zbMATH DE number 1069611 (Why is no real title available?)
• scientific article; zbMATH DE number 2001584 (Why is no real title available?)
• scientific article; zbMATH DE number 195005 (Why is no real title available?)
• scientific article; zbMATH DE number 1446717 (Why is no real title available?)
• Inverse problems: a Bayesian perspective
• Multiscale convergence and reiterated homogenisation
• Multiscale modelling and inverse problems
• Numerical solution of parabolic equations in high dimensions
• Preconditioning Markov Chain Monte Carlo Simulations Using Coarse-Scale Models
• Sparse tensor discretizations of high-dimensional parametric and stochastic PDEs
• Statistical and computational inverse problems.
• Stochastic Galerkin discretization of the log-normal isotropic diffusion problem
• Strong and weak error estimates for elliptic partial differential equations with random coefficients
• Tensor product type subspace splittings and multilevel iterative methods for anisotropic problems
• The calculation of linear least squares problems
• The Geometry of Random Fields

Cited in

(5)

• Approximation of Bayesian Inverse Problems for PDEs
• Fine scale uncertainty in parameter estimation for elliptic equations
• A Bayesian Numerical Homogenization Method for Elliptic Multiscale Inverse Problems
• Kernel Methods for Bayesian Elliptic Inverse Problems on Manifolds
• Transform-based particle filtering for elliptic Bayesian inverse problems

This page was built for publication: Bayesian inverse problems for recovering coefficients of two scale elliptic equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5382463)