Department of Mathematical Sciences / Institutionen för matematiska vetenskaper
Permanent URI for this communityhttps://gupea-staging.ub.gu.se/handle/2077/17606
Browse
Browsing Department of Mathematical Sciences / Institutionen för matematiska vetenskaper by Subject "approximate global convergence"
Now showing 1 - 2 of 2
- Results Per Page
- Sort Options
Item A two-stage numerical procedure for an inverse scattering problem(Chalmers University of Technology and University of Gothenburg, 2015) Bondestam Malmberg, John; Department of Mathematical Sciences, Chalmers University of Technology and University of GothenburgIn this thesis we study a numerical procedure for the solution of the inverse problem of reconstructing location, shape and material properties (in particular refractive indices) of scatterers located in a known background medium. The data consist of time-resolved backscattered radar signals from a single source position. This relatively small amount of data and the ill-posed nature of the inversion are the main challenges of the problem. Mathematically, the problem is formulated as a coefficient inverse problem for a system of partial differential equations derived from Maxwell’s equations. The numerical procedure is divided into two stages. In the first stage, a good initial approximation for the unknown coefficient is computed by an approximately globally convergent algorithm. This initial approximation is refined in the second stage, where an adaptive finite element method is employed to minimize a Tikhonov functional. An important tool for the second stage is a posteriori error estimates – estimates in terms of known (computed) quantities – for the difference between the computed coefficient and the true minimizing coefficient. This thesis includes four papers. In the first two, the a posteriori error analysis required for the adaptive finite element method in the second stage is extended from the previously existing indirect error estimators to direct ones. The last two papers concern verification of the two-stage numerical procedure on experimental data. We find that location and material properties of scatterers are obtained already in the first stage, while shapes are significantly improved in the second stage.Item Efficient Adaptive Algorithms for an Electromagnetic Coefficient Inverse Problem(2017-06-08) Malmberg, John BondestamThis thesis comprises five scientific papers, all of which are focusing on the inverse problem of reconstructing a dielectric permittivity which may vary in space inside a given domain. The data for the reconstruction consist of time-domain observations of the electric field, resulting from a single incident wave, on a part of the boundary of the domain under consideration. The medium is assumed to be isotropic, non-magnetic, and non-conductive. We model the permittivity as a continuous function, and identify distinct objects by means of iso-surfaces at threshold values of the permittivity. Our reconstruction method is centred around the minimization of a Tikhonov functional, well known from the theory of ill-posed problems, where the minimization is performed in a Lagrangian framework inspired by optimal control theory for partial differential equations. Initial approximations for the regularization and minimization are obtained either by a so-called approximately globally convergent method, or by a (simpler but less rigorous) homogeneous background guess. The functions involved in the minimization are approximated with finite elements, or with a domain decomposition method with finite elements and finite differences. The computational meshes are refined adaptively with regard to the accuracy of the reconstructed permittivity, by means of an a posteriori error estimate derived in detail in the fourth paper. The method is tested with success on simulated as well as laboratory measured data.