OpenMEEG: opensource software for quasistatic bioelectromagnetics

Background: Interpreting and controlling bioelectromagnetic phenomena require realistic physiological models and accurate numerical solvers. A semi-realistic model often used in practise is the piecewise constant conductivity model, for which only the interfaces have to be meshed. This simplified model makes it possible to use Boundary Element Methods. Unfortunately, most Boundary Element solutions are confronted with accuracy issues when the conductivity ratio between neighboring tissues is high, as for instance the scalp/skull conductivity ratio in electro-encephalography. To overcome this difficulty, we proposed a new method called the symmetric BEM, which is implemented in the OpenMEEG software. The aim of this paper is to presen

Models and methods for computational electromagnetic dosimetry

The interaction between electromagnetic fields and the human body is a very complicated issue. In most cases it is not possible to measure accurately the electrical response of the human body to external sources. Because of this computational methods are used as an aid when determining the safety levels for human exposure to electromagnetic fields. In this thesis the field distribution caused by various sources is determined in different parts of the human body using detailed and anatomically correct human body models and computational methods. Both the distribution of electric currents induced by low-frequency magnetic fields and the absorption of radio-frequency fields are studied. The accuracy and reliability of the models and methods used is verified by comparing the acquired results to known closed-form solutions and calibration measurements. The obtained results can be utilised in the reliability analysis of computational methods used in electromagnetic dosimetry. Furthermore, some of the results are needed, for example, in the safety guidelines of medical personnel working close to magnetic resonance imaging scanners

Fast and Efficient Formulations for Electroencephalography-Based Neuroimaging Strategies

L'abstract è presente nell'allegato / the abstract is in the attachmen

Doctor of Philosophy

dissertationInverse Electrocardiography (ECG) aims to noninvasively estimate the electrophysiological activity of the heart from the voltages measured at the body surface, with promising clinical applications in diagnosis and therapy. The main challenge of this emerging technique lies in its mathematical foundation: an inverse source problem governed by partial differential equations (PDEs) which is severely ill-conditioned. Essential to the success of inverse ECG are computational methods that reliably achieve accurate inverse solutions while harnessing the ever-growing complexity and realism of the bioelectric simulation. This dissertation focuses on the formulation, optimization, and solution of the inverse ECG problem based on finite element methods, consisting of two research thrusts. The first thrust explores the optimal finite element discretization specifically oriented towards the inverse ECG problem. In contrast, most existing discretization strategies are designed for forward problems and may become inappropriate for the corresponding inverse problems. Based on a Fourier analysis of how discretization relates to ill-conditioning, this work proposes refinement strategies that optimize approximation accuracy o f the inverse ECG problem while mitigating its ill-conditioning. To fulfill these strategies, two refinement techniques are developed: one uses hybrid-shaped finite elements whereas the other adapts high-order finite elements. The second research thrust involves a new methodology for inverse ECG solutions called PDE-constrained optimization, an optimization framework that flexibly allows convex objectives and various physically-based constraints. This work features three contributions: (1) fulfilling optimization in the continuous space, (2) formulating rigorous finite element solutions, and (3) fulfilling subsequent numerical optimization by a primal-dual interiorpoint method tailored to the given optimization problem's specific algebraic structure. The efficacy o f this new method is shown by its application to localization o f cardiac ischemic disease, in which the method, under realistic settings, achieves promising solutions to a previously intractable inverse ECG problem involving the bidomain heart model. In summary, this dissertation advances the computational research of inverse ECG, making it evolve toward an image-based, patient-specific modality for biomedical research