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.

CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Factors predicting mortality in type 2 diabetes. With special reference to physical exercise, blood pressure, proteinuria, inflammation and P wave duration

Background: Type 2 diabetes patients have a 2-4 fold risk of cardiovascular disease (CVD) compared to the general population. In type 2 diabetes, several CVD risk factors have been identified, including obesity, hypertension, hyperglycemia, proteinuria, sedentary lifestyle and dyslipidemia. Although much of the excess CVD risk can be attributed to these risk factors, a significant proportion is still unknown. Aims: To assess in middle-aged type 2 diabetic subjects the joint relations of several conventional and non-conventional CVD risk factors with respect to cardiovascular and total mortality. Subjects and methods: This thesis is part of a large prospective, population based East-West type 2 diabetes study that was launched in 1982-1984. It includes 1,059 middle-aged (45-64 years old) participants. At baseline, a thorough clinical examination and laboratory measurements were performed and an ECG was recorded. The latest follow-up study was performed 18 years later in January 2001 (when the subjects were 63-81 years old). The study endpoints were total mortality and mortality due to CVD, coronary heart disease (CHD) and stroke. Results: Physically more active patients had significantly reduced total, CVD and CHD mortality independent of high-sensitivity C-reactive protein (hs-CRP) levels unless proteinuria was present. Among physically active patients with a hs-CRP level >3 mg/L, the prognosis of CVD mortality was similar to patients with hs-CRP levels ≤3 mg/L. The worst prognosis was among physically inactive patients with hs-CRP levels >3 mg/L. Physically active patients with proteinuria had significantly increased total and CVD mortality by multivariate analyses. After adjustment for confounding factors, patients with proteinuria and a systolic BP <130 mmHg had a significant increase in total and CVD mortality compared to those with a systolic BP between 130 and 160 mmHg. The prognosis was similar in patients with a systolic BP <130 mmHg and ≥160 mmHg. Among patients without proteinuria, a systolic BP <130 mmHg was associated with a non-significant reduction in mortality. A P wave duration ≥114 ms was associated with a 2.5-fold increase in stroke mortality among patients with prevalent CHD or claudication. This finding persisted in multivariable analyses. Among patients with no comorbidities, there was no relationship between P wave duration and stroke mortality. Conclusions: Physical activity reduces total and CVD mortality in patients with type 2 diabetes without proteinuria or with elevated levels of hs-CRP, suggesting that the anti-inflammatory effect of physical activity can counteract increased CVD morbidity and mortality associated with a high CRP level. In patients with proteinuria the protective effect was not, however, present. Among patients with proteinuria, systolic BP <130 mmHg may increase mortality due to CVD. These results demonstrate the importance of early intervention to prevent CVD and to control all-cause mortality among patients with type 2 diabetes. The presence of proteinuria should be taken into account when defining the target systolic BP level for prevention of CVD deaths. A prolongation of the duration of the P wave was associated with increased stroke mortality among high-risk patients with type 2 diabetes. P wave duration is easy to measure and merits further examination to evaluate its importance for estimation of the risk of stroke among patients with type 2 diabetes.

Modeling of high pressure pretreatment process for gold leaching

Investigation of high pressure pretreatment process for gold leaching is the objective of the present master's thesis. The gold ores and concentrates which cannot be easily treated by leaching process are called "refractory". These types of ores or concentrates often have high content of sulfur and arsenic that renders the precious metal inaccessible to the leaching agents. Since the refractory ores in gold manufacturing industry take a considerable share, the pressure oxidation method (autoclave method) is considered as one of the possible ways to overcome the related problems. Mathematical modeling is the main approach in this thesis which was used for investigation of high pressure oxidation process. For this task, available information from literature concerning this phenomenon, including chemistry, mass transfer and kinetics, reaction conditions, applied apparatus and application, was collected and studied. The modeling part includes investigation of pyrite oxidation kinetics in order to create a descriptive mathematical model. The following major steps are completed: creation of process model by using the available knowledge; estimation of unknown parameters and determination of goodness of the fit; study of the reliability of the model and its parameters.

Suodattimessa tapahtuvan kiintoaineen laskeutumisen vaikutus suodinkakun ominaisuuksiin

Laskeutus ja suodatus ovat paljon tutkittuja ja laajassa käytössä olevia mekaanisia erotusmenetelmiä. Laskeutumisen vaikutusta suodatuksen yhteydessä ei kuitenkaan ole tutkittu juurikaan. Tämän työn tavoitteena on selvittää suodattimessa ennen vakiopaine suodatusta tapahtuvan kiintoaineen laskeutumisen vaikutusta suodatuksen tuloksiin. Työn kirjallisuusosuudessa käsitellään työhön liittyvää teoriaa partikkelikokojakaumista, laskeutumisesta, suodatuksesta ja flokkulanttien käytöstä sekä tehdään yhteenvetoa aiemmasta laskeutumista ja suodatusta yhdistävästä tutkimuksesta. Koska käsitellyt aiheet ovat laajoja ja niitä on useita, joudutaan teorian esittämisessä tekemään rajauksia. Työn kokeellisessa osassa tutkitaan suodattimessa tapahtuvan laskeutuksen vaikutusta suodatukseen. Kokeet suoritetaan 20 m-% kalsiumkarbonaattilietteellä. Kokeissa tutkitaan laskeutumisen vaikutusta antamalla lietteen laskeutua suodattimessa tietyn ajan ennen suodatuksen aloitusta. Lisäksi tutkitaan flokkulanttien lisäämisen vaikutusta mahdollisiin ilmiöihin. Koetuloksista nähdään pidemmän laskeutusajan ennen suodatusta alentavan kakun keskimääräistä ominaisvastusta ja vaikutus kasvaa merkittävästi kun lietteeseen on lisätty flokkulanttia. Suodatusta edeltävä laskeutus siis helpottaa varsinaista suodatusta.