No Arabic abstract
Invasive intracranial electroencephalography (iEEG) or electrocorticography (ECoG) measures electrical potential directly on the surface of the brain, and, combined with numerical modeling, can be used to inform treatment planning for epilepsy surgery. Accurate solution of the iEEG or ECoG forward problem, which is a crucial prerequisite for solving the inverse problem in epilepsy seizure onset localization, requires accurate representation of the patients brain geometry and tissue electrical conductivity after implantation of electrodes. However, implantation of subdural grid electrodes causes the brain to deform, which invalidates preoperatively acquired image data. Moreover, postoperative MRI is incompatible with implanted electrodes and CT has insufficient range of soft tissue contrast, which precludes both MRI and CT from being used to obtain the deformed postoperative geometry. In this paper, we present a biomechanics-based image warping procedure using preoperative MRI for tissue classification and postoperative CT for locating implanted electrodes to perform non-rigid registration of the preoperative image data to the postoperative configuration. We solve the iEEG forward problem on the predicted postoperative geometry using the finite element method (FEM) which accounts for patient-specific inhomogeneity and anisotropy of tissue conductivity. Results for the simulation of a current source in the brain show large differences in electrical potential predicted by the models based on the original images and the deformed images corresponding to the brain geometry deformed by placement of invasive electrodes. Computation of the leadfield matrix also showed significant differences between the different models. The results suggest that significant improvements in source localization accuracy may be realized by the application of the proposed modeling methodology.
Knowledge of appropriate local fiber architecture is necessary to simulate patient-specific electromechanics in the human heart. However, it is not yet possible to reliably measure in-vivo fiber directions, especially in human atria. Thus, we present a method which defines the fiber architecture in arbitrarily shaped atria using image registration and reorientation methods based on atlas atria with fibers predefined from detailed histological observations. Thereby, it is possible to generate detailed fiber families in every new patient-specific geometry in an automated, time-efficient process. We demonstrate the good performance of the image registration and fiber definition on ten differently shaped human atria. Additionally, we show that characteristics of the electrophysiological activation pattern which appear in the atlas atria also appear in the patients atria. We arrive at analogous conclusions for coupled electro-mechano-hemodynamical computations.
Augmented Reality is used in Image Guided surgery (AR IG) to fuse surgical landmarks from preoperative images into a video overlay. Physical simulation is essential to maintaining accurate position of the landmarks as surgery progresses and ensuring patient safety by avoiding accidental damage to vessels etc. In liver procedures, AR IG simulation accuracy is hampered by an inability to model stiffness variations unique to the patients disease. We introduce a novel method to account for patient specific stiffness variation based on Magnetic Resonance Elastography (MRE) data. To the best of our knowledge we are the first to demonstrate the use of in-vivo biomechanical data for AR IG landmark placement. In this early work, a comparative evaluation of our MRE data driven simulation and the traditional method shows clinically significant differences in accuracy during landmark placement and motivates further animal model trials.
In its permanent quest of mechanobiological homeostasis, our vascula-ture significantly adapts across multiple length and time scales in various physiological and pathological conditions. Computational modeling of vascular growth and remodeling (G&R) has significantly improved our insights of the mechanobio-logical processes of diseases such as hypertension or aneurysms. However, patient-specific computational modeling of ascending thoracic aortic aneurysm (ATAA) evolution, based on finite-element models (FEM), remains a challenging scientific problem with rare contributions, despite the major significance of this topic of research. Challenges are related to complex boundary conditions and geometries combined with layer-specific G&R responses. To address these challenges, in the current paper, we employed the constrained mixture model (CMM) to model the arterial wall as a mixture of different constituents such as elastin, collagen fiber families and smooth muscle cells (SMCs). Implemented in Abaqus as a UMAT, this first patient-specific CMM-based FEM of G&R in human ATAA was first validated for canonical problems such as single-layer thick-wall cylindrical and bi-layer thick-wall toric arterial geometries. Then it was used to predict ATAA evolution for a patient-specific aortic geometry, showing that the typical shape of an ATAA can be simply produced by elastin proteolysis localized in regions of deranged hemodymanics. The results indicate a transfer of stress to the adventitia by elastin loss and continuous adaptation of the stress distribution due to change of ATAA shape. Moreover, stress redistribution leads to collagen deposition where the maximum elastin mass is lost, which in turn leads to stiffening of the arterial wall. As future work, the predictions of this G&R framework will be validated on datasets of patient-specific ATAA geometries followed up over a significant number of years.
An arteriovenous fistula, created by artificially connecting segments of a patients vasculature, is the preferred way to gain access to the bloodstream for kidney dialysis. The increasing power and availability of supercomputing infrastructure means that it is becoming more realistic to use simulations to help identify the best type and location of a fistula for a specific patient. We describe a 3D fistula model that uses the lattice Boltzmann method to simultaneously resolve blood flow in patient-specific arteries and veins. The simulations conducted here, comprising vasculatures of the whole forearm, demonstrate qualified validation against clinical data. Ongoing research to further encompass complex biophysics on realistic time scales will permit the use of human-scale physiological models for basic and clinical medicine.
A patient-specific fluid-structure interaction (FSI) model of a phase-contrast magnetic resonance angiography (PC-MRA) imaged arteriovenous fistula is presented. The numerical model is developed and simulated using a commercial multiphysics simulation package where a semi-implicit FSI coupling scheme combines a finite volume method blood flow model and a finite element method vessel wall model. A pulsatile mass-flow boundary condition is prescribed at the artery inlet of the model, and a three-element Windkessel model at the artery and vein outlets. The FSI model is freely available for analysis and extension. This work shows the effectiveness of combining a number of stabilisation techniques to simultaneously overcome the added-mass effect and optimise the efficiency of the overall model. The PC-MRA data, fluid model, and FSI model results show almost identical flow features in the fistula; this applies in particular to a flow recirculation region in the vein that could potentially lead to fistula failure.