Abstract
We propose a level set approach for elliptic inverse problems with piecewise constant coefficients. The geometry of the discontinuity of the coefficient is represented implicitly by level set functions. The inverse problem is solved using a variational augmented Lagrangian formulation with total variation regularization of the coefficient. The corresponding Euler-Lagrange equation gives the evolution equation for the level set functions and the constant values of the coefficients. We use a multiple level set representation which allows the coefficient to have multiple constant regions. Knowledge of the exact number of regionsis not required, only an upper bound is needed. Numerical experiments show that the method can recover coefficients with rather complicated geometries of discontinuities under moderate amount of noise in the observation data. The method is also robust with respect to the initial guess for the geometry of the coefficient discontinuities.
Original language | English (US) |
---|---|
Pages (from-to) | 40-66 |
Number of pages | 27 |
Journal | Journal of Computational Physics |
Volume | 193 |
Issue number | 1 |
DOIs | |
State | Published - Jan 1 2004 |
Externally published | Yes |
ASJC Scopus subject areas
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics