Stability for finite element discretization of some elliptic inverse parameter problems from internal data – application to elastography

by   Elie Bretin, et al.

In this article, we provide stability estimates for the finite element discretization of a class of inverse parameter problems of the form -∇·(μ S) = f in a domain Ω of ^d. Here μ is the unknown parameter to recover, the matrix valued function S and the vector valued distribution f are known. As uniqueness is not guaranteed in general for this problem, we prove a Lipschitz-type stability estimate in an hyperplane of L^2(Ω). This stability is obtained through an adaptation of the so-called discrete inf-sup constant or LBB constant to a large class of first-order differential operators. We then provide a simple and original discretization based on hexagonal finite element that satisfies the discrete stability condition and shows corresponding numerical reconstructions. The obtained algebraic inversion method is efficient as it does not require any iterative solving of the forward problem and is very general as it does not require any smoothness hypothesis for the data nor any additional information at the boundary.


page 21

page 24


Finite element analysis for a diffusion equation on a harmonically evolving domain

We study convergence of the evolving finite element semi-discretization ...

Existence, uniqueness, and approximation of a fictitious domain formulation for fluid-structure interactions

In this paper we describe a computational model for the simulation of fl...

Impact Study of Numerical Discretization Accuracy on Parameter Reconstructions and Model Parameter Distributions

Numerical models are used widely for parameter reconstructions in the fi...

Recovery by discretization corrected particle strength exchange (DC PSE) operators

A new recovery technique based on discretization corrected particle stre...

PDE-constrained optimization for electroencephalographic source reconstruction

This paper introduces a novel numerical method for the inverse problem o...

Radial complex scaling for anisotropic scalar resonance problems

The complex scaling/perfectly matched layer method is a widely spread te...

Please sign up or login with your details

Forgot password? Click here to reset