On applying stochastic Galerkin finite element method to electrical impedance tomography

Abstract

In this thesis, a new solution strategy based on stochastic Galerkin finite element method is introduced for the complete electrode model of electrical impedance tomography. The method allows writing an analytical approximation for the solution to the inverse problem of electrical impedance tomography in the setting of Bayesian inversion with the help of multivariate orthogonal polynomials. If the measurement setting, i.e., geometry, priors, etc., is known (well) in advance, most computations required by the introduced method can be performed and stored before the actual measurement. The formation of the approximative solution to the inverse problem, i.e., the posterior probability density, is practically free of charge once the measurements are available. Subsequently, estimates for the quantities of interest can typically be obtained by either minimizing an explicitly known polynomial or integrating a known analytical expression. In addition, some advances in the development of numerical solvers for parametric partial differential equations in the setting of generalized Polynomial Chaos and stochastic Galerkin finite element method are presented.

Väitöskirjassa esitellään uusi stokastiseen Galerkinin elementtimenetelmään perustuva ratkaisutekniikka impedanssitomografian täydelliselle elektrodimallille. Tekniikka mahdollistaa analyyttisen approksimaation kirjoittamisen impedanssitomografian käänteisongelman bayesiläiselle ratkaisulle usean muuttujan ortogonaalipolynomien avulla. Mikäli impedanssitomografian mittausasetelma tunnetaan ennakkoon, voidaan suurin osa menetelmän vaatimista laskutoimituksista suorittaa ja varastoida ennen varsinaisen mittauksen suorittamista. Mittausten jälkeen approksimaatio inversio-ongelman ratkaisulle eli niin sanotulle posterioritiheysfunktiolle voidaan muodostaa ilman raskasta laskentaa. Tämän jälkeen haluttuja ratkaisuestimaatteja voidaan tyypillisesti laskea joko minimoimalla eksplisiittisesti tunnettua polynomia tai integroimalla tunnettua analyyttistä lauseketta. Lisäksi väitöskirjassa kehitetään yleistettyyn polynomikaaokseen ja stokastiseen Galerkinin elementtimenetelmään pohjautuvia ratkaisumenetelmiä parametririippuville osittaisdifferentiaaliyhtälöille.