A Bayesian numerical homogenization method for elliptic multiscale inverse problems

Assyr Abdulle, Andrea Di Blasio
2020
Article

Résumé

A new strategy based on numerical homogenization and Bayesian techniques for solving multiscale inverse problems is introduced. We consider a class of elliptic problems which vary at a microscopic scale, and we aim at recovering the highly oscillatory tensor from measurements of the fine scale solution at the boundary, using a coarse model based on numerical homogenization and model order reduction. We provide a rigorous Bayesian formulation of the problem, taking into account different possibilities for the choice of the prior measure. We prove well-posedness of the effective posterior measure and, by means of G-convergence, we establish a link between the effective posterior and the fine scale model. Several numerical experiments illustrate the efficiency of the proposed scheme and confirm the theoretical findings.

Official source

https://infoscience.epfl.ch/record/255523?ln=fr

À propos de ce résultat

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Concepts associés

Chargement

Publications associées

Chargement

A new strategy based on numerical homogenization and Bayesian techniques for solvingmultiscale inverse problems is introduced. We consider a class of elliptic problems which vary ata microscopic scale, and we aim at recovering the highly oscillatory tensor from measurements ofthe fine scale solution at the boundary, using a coarse model based on numerical homogenizationand model order reduction. We provide a rigorous Bayesian formulation of the problem, takinginto account different possibilities for the choice of the prior measure. We prove well-posednessof the effective posterior measure and, by means of G-convergence, we establish a link betweenthe effective posterior and the fine scale model. Several numerical experiments illustrate theefficiency of the proposed scheme and confirm the theoretical findings.

MATHICSE2018

Chargement

Numerical homogenization and model order reduction for multiscale inverse problems

Assyr Abdulle, Andrea Di Blasio

A new numerical method based on numerical homogenization and model order reduction is introduced for the solution of multiscale inverse problems. We consider a class of elliptic problems with highly oscillatory tensors that varies on a microscopic scale. We assume that the micro structure is known and seek to recover a macroscopic scalar parametrization of the microscale tensor (e.g. volume fraction). Departing from the full fine scale model that would require mesh resolution for the forward problem down to the finest scale, we solve the inverse problem for a coarse model obtained by numerical homogenization. The input data, i.e., measurement from the Dirichlet to Neumann map, are solely based on the original fine scale model. Furthermore, reduced basis techniques are used to avoid computing effective coefficients for the forward solver at each integration point of the macroscopic mesh. Uniqueness and stability of the effective inverse problem is established based on standard assumptions for the fine scale model and a link with this latter model is established by means of G-convergence. A priori error estimates are established for our method. Numerical experiments illustrate the efficiency of the proposed scheme and confirm our theoretical finding.

2016

Chargement

A new numerical method based on numerical homogenization and model order reduction is introduced for the solution of multiscale inverse problems. We consider a class of elliptic problems with highly oscillatory tensors that varies on a microscopic scale. We assume that the micro structure is known and seek to recover a macroscopic scalar parameterization of the microscale tensor (e.g., volume fraction). Departing from the full fine-scale model, which would require mesh resolution for the forward problem down to the finest scale, we solve the inverse problem for a coarse model obtained by numerical homogenization. The input data, i.e., measurement from the Dirichlet-to-Neumann map, are solely based on the original fine-scale model. Furthermore, reduced basis techniques are used to avoid computing effective coefficients for the forward solver at each integration point of the macroscopic mesh. Uniqueness and stability of the effective inverse problem is established based on standard assumptions for the fine-scale model, and a link to this latter model is established by means of G-convergence. A priori error estimates are established for our method. Numerical experiments illustrate the efficiency of the proposed scheme and confirm our theoretical findings.

SIAM PUBLICATIONS2019

Chargement