Approaches to Multiscale Inverse Problems

Show simple item record

dc.contributor.advisor Kaipio, J en
dc.contributor.advisor Taylor, S en Nicholson, Ruanui en 2016-09-15T00:01:16Z en 2016 en
dc.identifier.uri en
dc.description.abstract Many scientific problems involve parameters such as conductivity, permeability or density which vary on multiple spatial and/or temporal scales. When such a parameter is to be estimated from noisy indirect measurements we face a challenging dichotomy: In practical situations the computational cost required to accurately take into account the small scale behaviour can be prohibitive. On the other hand, estimation of the conductivity is an ill-posed inverse problem meaning that any modelling errors, such as neglecting the small scale, which are not accounted for can render estimate useless. A typical way of reducing the computational cost is to simply neglect any behaviour at the smaller scales, enabling a coarse discretisation of the problem. However, if the effects of disregarding smaller scale characteristics along with the coarse discretisation are not taken into account the estimates attained can be misleading. In this thesis, we develop computationally feasible models to tackle the problem of estimating the multiscale distributed parameter of the Poisson equation, which is closely related to the multiscale electrical impedance tomography (EIT) problem. In the case of estimating a low dimensional representation of the conductivity (devoid of the small scale features) we show that by applying the Bayesian approximation error (BAE) approach to marginalise over small scale effects and discretisation errors we are able to substantially reduce the dimension of the problem while maintaining accurate estimates. Moreover, we show that by applying a feature extraction type modification to the BAE approach we can in some cases estimate the amplitude and correlation length of the small scale component of the conductivity. In the case of estimating the full multiscale conductivity we show that multiscale finite element methods (MsFEM) can be implemented at both the forward modelling stage and at the inversion stage, which along with a somewhat coarsened discretisation can reduce overall computational cost. By deriving a closed form expression for the Jacobian matrix which represents the rate at which the electric potential calculated using MsFEM changes with respect to the conductivity we are able to apply gradient based optimisation techniques to estimate the conductivity. The BAE approach was also implemented in this procedure to take into account any modelling discrepancy caused by the use of MsFEM and discretisation errors. Use of such a procedure leads to estimates attained by using 7 multiscale basis functions in line with those calculated using finite elements on 85 linear basis functions. en
dc.publisher ResearchSpace@Auckland en
dc.relation.ispartof PhD Thesis - University of Auckland en
dc.relation.isreferencedby UoA99264877506702091 en
dc.rights Items in ResearchSpace are protected by copyright, with all rights reserved, unless otherwise indicated. Previously published items are made available in accordance with the copyright policy of the publisher. en
dc.rights.uri en
dc.rights.uri en
dc.title Approaches to Multiscale Inverse Problems en
dc.type Thesis en Applied Mathematics en The University of Auckland en Doctoral en PhD en
dc.rights.holder Copyright: The author en
dc.rights.accessrights en
pubs.elements-id 541205 en Engineering en Engineering Science en
pubs.record-created-at-source-date 2016-09-15 en

Files in this item

Find Full text

This item appears in the following Collection(s)

Show simple item record


Search ResearchSpace