Pulsatile Blood Flow Simulations in Computed Tomography(CT) Scan-Based and Idealized Geometries of Human Aorta

Blood flow in human aorta is an unsteady and complex phenomenon. The complex patterns are related to the geometrical features like curvature, bends, and branching and pulsatile nature of flow from left ventricle of heart. The aim of this work was to understand the effect of aorta geometry on the flow dynamics. To achieve this, 3D realistic and idealized models of descending aorta were reconstructed from Computed Tomography (CT) images of a female patient. The geometries were reconstructed using medical image processing code. The blood flow in aorta was assumed to be laminar and incompressible and the blood was assumed to be Newtonian fluid. A time dependent pulsatile and parabolic boundary condition was deployed at inlet. Steady and unsteady blood flow simulations were performed in real and idealized geometries of descending aorta using a Finite Volume Method (FVM) code. Analysis of Wall Shear Stress (WSS) distribution, pressure distribution, and axial velocity profiles were carried out in both geometries at steady and unsteady state conditions. The results obtained in thesis work reveal that the idealization of geometry underestimates the values of WSS especially near the region with sudden change of diameter. However, the resultant pressure and velocity in idealized geometry are close to those in real geometry

CT and PET/CT Hybrid Imaging of Coronary Artery Disease

Coronary artery disease (CAD) is a chronic process that evolves over decades and may culminate in myocardial infarction (MI). While invasive coronary angiography (ICA) is still considered the gold standard of imaging CAD, non-invasive assessment of both the vascular anatomy and myocardial perfusion has become an intriguing alternative. In particular, computed tomography (CT) and positron emission tomography (PET) form an attractive combination for such studies. Increased radiation dose is, however, a concern. Our aim in the current thesis was to test novel CT and PET techniques alone and in hybrid setting in the detection and assessment of CAD in clinical patients. Along with diagnostic accuracy, methods for the reduction of the radiation dose was an important target. The study investigating the coronary arteries of patients with atrial fibrillation (AF) showed that CAD may be an important etiology of AF because a high prevalence of CAD was demonstrated within AF patients. In patients with suspected CAD, we demonstrated that a sequential, prospectively ECG-triggered CT technique was applicable to nearly 9/10 clinical patients and the radiation dose was over 60% lower than with spiral CT. To detect the functional significance of obstructive CAD, a novel software for perfusion quantification, CarimasTM, showed high reproducibility with 15O-labelled water in PET, supporting feasibility and good clinical accuracy. In a larger cohort of 107 patients with moderate 30-70% pre-test probability of CAD, hybrid PET/CT was shown to be a powerful diagnostic method in the assessment of CAD with diagnostic accuracy comparable to that of invasive angiography and fractional flow reserve (FFR) measurements. A hybrid study may be performed with a reasonable radiation dose in a vast majority of the cases, improving the performance of stand-alone PET and CT angiography, particularly when the absolute quantification of the perfusion is employed. These results can be applied into clinical practice and will be useful for daily clinical diagnosis of CAD.

Imaging of the Vulnerable Atherosclerotic Plaque. Pre-Clinical Evaluation of PET Tracers for Vascular Inflammation

Atherosclerosis is a vascular inflammatory disease causing coronary artery disease, myocardial infarct and stroke, the leading causes of death in Finland and in many other countries. The development of atherosclerotic plaques starts already in childhood and is an ongoing process throughout life. Rupture of a plaque and the following occlusion of the vessel is the main reason for myocardial infarct and stroke, but despite extensive research, the prediction of rupture remains a major clinical problem. Inflammation is considered a key factor in the vulnerability of plaques to rupture. Measuring the inflammation in plaques non-invasively is one potential approach for identification of vulnerable plaques. The aim of this study was to evaluate tracers for positron emission tomography (PET) imaging of vascular inflammation. The studies were performed with a mouse model of atherosclerosis by using ex vivo biodistribution, autoradiography and in vivo PET and computed tomography (CT). Several tracers for inflammation activity were tested and compared with the morphology of the plaques. Inflammation in the atherosclerotic plaques was evaluated as expression of active macrophages. Systematic analysis revealed that the uptake of 18F-FDG and 11C-choline, tracers for metabolic activity in inflammatory cells, was more prominent in the atherosclerotic plaques than in the surrounding healthy vessel wall. The tracer for αvβ3 integrin, 18Fgalacto- RGD, was also found to have high potential for imaging inflammation in the plaques. While 11C-PK11195, a tracer targeted to receptors in active macrophages, was shown to accumulate in active plaques, the target-to-background ratio was not found to be ideal for in vivo imaging purposes. In conclusion, tracers for the imaging of inflammation in atherosclerotic plaques can be tested in experimental pre-clinical settings to select potential imaging agents for further clinical testing. 18F-FDG, 18F-galacto-RGD and 11C-choline choline have good properties, and further studies to clarify their applicability for atherosclerosis imaging in humans are warranted.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

Design and characterization of large area position sensitive radiation detectors

Planar, large area, position sensitive silicon detectors are widely utilized in high energy physics research and in medical, computed tomography (CT). This thesis describes author's research work relating to development of such detector components. The key motivation and objective for the research work has been the development of novel, position sensitive detectors improving the performance of the instruments they are intended for. Silicon strip detectors are the key components of barrel-shaped tracking instruments which are typically the innermost structures of high energy physics experimental stations. Particle colliders such as the former LEP collider or present LHC produce particle collisions and the silicon strip detector based trackers locate the trajectories of particles emanating from such collisions. Medical CT has become a regular part of everyday medical care in all developed countries. CT scanning enables x-ray imaging of all parts of the human body with an outstanding structural resolution and contrast. Brain, chest and abdomen slice images with a resolution of 0.5 mm are possible and latest CT machines are able to image whole human heart between heart beats. The two application areas are presented shortly and the radiation detection properties of planar silicon detectors are discussed. Fabrication methods and preamplifier electronics of the planar detectors are presented. Designs of the developed, large area silicon detectors are presented and measurement results of the key operating parameters are discussed. Static and dynamic performance of the developed silicon strip detectors are shown to be very satisfactory for experimental physics applications. Results relating to the developed, novel CT detector chips are found to be very promising for further development and all key performance goals are met.

Detection of Pathologic Changes Following Traumatic Brain Injury Using Magnetic Resonance Imaging

Background: Approximately two percent of Finns have sequels after traumatic brain injury (TBI), and many TBI patients are young or middle-aged. The high rate of unemployment after TBI has major economic consequences for society, and traumatic brain injury often has remarkable personal consequences, as well. Structural imaging is often needed to support the clinical TBI diagnosis. Accurate early diagnosis is essential for successful rehabilition and, thus, may also influence the patient’s outcome. Traumatic axonal injury and cortical contusions constitute the majority of traumatic brain lesions. Several studies have shown magnetic resonance imaging (MRI) to be superior to computed tomography (CT) in the detection of these lesions. However, traumatic brain injury often leads to persistent symptoms even in cases with few or no findings in conventional MRI. Aims and methods: The aim of this prospective study was to clarify the role of conventional MRI in the imaging of traumatic brain injury, and to investigate how to improve the radiologic diagnostics of TBI by using more modern diffusion-weighted imaging (DWI) and diffusion tensor imaging (DTI) techniques. We estimated, in a longitudinal study, the visibility of the contusions and other intraparenchymal lesions in conventional MRI at one week and one year after TBI. We used DWI-based measurements to look for changes in the diffusivity of the normal-appearing brain in a case-control study. DTI-based tractography was used in a case-control study to evaluate changes in the volume, diffusivity, and anisotropy of the long association tracts in symptomatic TBI patients with no visible signs of intracranial or intraparenchymal abnormalities on routine MRI. We further studied the reproducibility of different tools to identify and measure white-matter tracts by using a DTI sequence suitable for clinical protocols. Results: Both the number and extent of visible traumatic lesions on conventional MRI diminished significantly with time. Slightly increased diffusion in the normal-appearing brain was a common finding at one week after TBI, but it was not significantly associated with the injury severity. Fractional anisotropy values, that represent the integrity of the white-matter tracts, were significantly diminished in several tracts in TBI patients compared to the control subjects. Compared to the cross-sectional ROI method, the tract-based analyses had better reproducibility to identify and measure white-matter tracts of interest by means of DTI tractography. Conclusions: As conventional MRI is still applied in clinical practice, it should be carried out soon after the injury, at least in symptomatic patients with negative CT scan. DWI-related brain diffusivity measurements may be used to improve the documenting of TBI. DTI tractography can be used to improve radiologic diagnostics in a symptomatic TBI sub-population with no findings on conventional MRI. Reproducibility of different tools to quantify fibre tracts vary considerably, which should be taken into consideration in the clinical DTI applications.

Advanced epithelial ovarian cancer – studies on preoperative [18F] FDG PET/CT and HE4 profile during primary chemotherapy

Epithelial ovarian cancer (EOC) is usually diagnosed in an advanced stage. The prognosis depends highly on the amount of the residual tumor in surgery. In patients with extensive disease, neoadjuvant chemotherapy (NACT) is used to diminish the tumor load before debulking surgery. New non-invasive methods are needed to preoperatively evaluate the disease dissemination and operability. [18F] FDG PET/CT (Positron emission tomography/computed tomography) is a promising method for cancer diagnostics and staging. The biomarker profiles during treatment can predict patient’s outcome. This prospective study included 41 EOC patients, 21 treated with primary surgery and 20 with NACT and interval surgery. The performances of preoperative contrast enhanced PET/CT (PET/ceCT) and diagnostic CT (ceCT) were compared. Perioperative visual estimation of tumor spread was studied in primary and interval surgery. The profile of the serum marker HE4 (Human epididymis 4) during primary chemotherapy was evaluated. In primary surgery, surgical findings were found to form an adequate reference standard for imaging studies. After NACT, the sensitivity for visual estimation of cancer dissemination was significantly worse. Preoperative PET/ceCT was more effective than ceCT alone in detecting extra-abdominal disease spread. The high number of supradiaphragmatic lymph node metastases detected by PET/ceCT at the time of diagnosis brings new insight in EOC spread patterns. The sensitivity of both PET/CT and ceCT remained modest in intra-abdominal areas important to operability. The HE4 profile was in concordance with the CA125 profile during primary chemotherapy. Its role in the evaluation of EOC chemotherapy response will be clarified in further studies.