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.

Membrane properties of structurally modified ceramides : effects on lipid lateral distribution and sphingomyelin-interactions in artificial bilayer membranes

Ceramides comprise a class of sphingolipids that exist only in small amounts in cellular membranes, but which have been associated with important roles in cellular signaling processes. The influences that ceramides have on the physical properties of bilayer membranes reach from altered thermodynamical behavior to significant impacts on the molecular order and lateral distribution of membrane lipids. Along with the idea that the membrane physical state could influence the physiological state of a cell, the membrane properties of ceramides have gained increasing interest. Therefore, membrane phenomena related to ceramides have become a subject of intense study both in cellular as well as in artificial membranes. Artificial bilayers, the so called model membranes, are substantially simpler in terms of contents and spatio-temporal variation than actual cellular membranes, and can be used to give detailed information about the properties of individual lipid species in different environments. This thesis focuses on investigating how the different parts of the ceramide molecule, i.e., the N-linked acyl chain, the long-chain sphingoid base and the membrane-water interface region, govern the interactions and lateral distribution of these lipids in bilayer membranes. With the emphasis on ceramide/sphingomyelin(SM)-interactions, the relevance of the size of the SMhead group for the interaction was also studied. Ceramides with methylbranched N-linked acyl chains, varying length sphingoid bases, or methylated 2N (amide-nitrogen) and 3O (C3-hydroxyl) at the interface region, as well as SMs with decreased head group size, were synthesized and their bilayer properties studied by calorimetric and fluorescence spectroscopic techniques. In brief, the results showed that the packing of the ceramide acyl chains was more sensitive to methyl-branching in the mid part than in the distal end of the N-linked chain, and that disrupting the interfacial structure at the amide-nitrogen, as opposed to the C3-hydroxyl, had greater effect on the interlipid interactions of ceramides. Interestingly, it appeared that the bilayer properties of ceramides could be more sensitive to small alterations in the length of the long-chain base than what was previously reported for the N-linked acyl chain. Furthermore, the data indicated that the SM-head group does not strongly influence the interactions between SMs and ceramides. The results in this thesis illustrate the pivotal role of some essential parts of the ceramide molecules in determining their bilayer properties. The thesis provides increased understanding of the molecular aspects of ceramides that possibly affect their functions in biological membranes, and could relate to distinct effects on cell physiology.

Forecasting of wind speeds and directions with artificial neural networks

In this master’s thesis, wind speeds and directions were modeled with the aim of developing suitable models for hourly, daily, weekly and monthly forecasting. Artificial Neural Networks implemented in MATLAB software were used to perform the forecasts. Three main types of artificial neural network were built, namely: Feed forward neural networks, Jordan Elman neural networks and Cascade forward neural networks. Four sub models of each of these neural networks were also built, corresponding to the four forecast horizons, for both wind speeds and directions. A single neural network topology was used for each of the forecast horizons, regardless of the model type. All the models were then trained with real data of wind speeds and directions collected over a period of two years in the municipal region of Puumala in Finland. Only 70% of the data was used for training, validation and testing of the models, while the second last 15% of the data was presented to the trained models for verification. The model outputs were then compared to the last 15% of the original data, by measuring the mean square errors and sum square errors between them. Based on the results, the feed forward networks returned the lowest generalization errors for hourly, weekly and monthly forecasts of wind speeds; Jordan Elman networks returned the lowest errors when used for forecasting of daily wind speeds. Cascade forward networks gave the lowest errors when used for forecasting daily, weekly and monthly wind directions; Jordan Elman networks returned the lowest errors when used for hourly forecasting. The errors were relatively low during training of the models, but shot up upon simulation with new inputs. In addition, a combination of hyperbolic tangent transfer functions for both hidden and output layers returned better results compared to other combinations of transfer functions. In general, wind speeds were more predictable as compared to wind directions, opening up opportunities for further research into building better models for wind direction forecasting.

Artificial Ribonucleases: Oligonucleotides Conjugated with Metal Ion Chelates of Azacrowns

Ribonucleic acid (RNA) has many biological roles in cells: it takes part in coding, decoding, regulating and expressing of the genes as well as has the capacity to work as a catalyst in numerous biological reactions. These qualities make RNA an interesting object of various studies. Development of useful tools with which to investigate RNA is a prerequisite for more advanced research in the field. One of such tools may be the artificial ribonucleases, which are oligonucleotide conjugates that sequence-selectively cleave complementary RNA targets. This thesis is aimed at developing new efficient metal-ion-based artificial ribonucleases. On one hand, to solve the challenges related to solid-supported synthesis of metal-ion-binding conjugates of oligonucleotides, and on the other hand, to quantify their ability to cleave various oligoribonucleotide targets in a pre-designed sequence selective manner. In this study several artificial ribonucleases based on cleaving capability of metal ion chelated azacrown moiety were designed and synthesized successfully. The most efficient ribonucleases were the ones with two azacrowns close to the 3´- end of the oligonucleotide strand. Different transition metal ions were introduced into the azacrown moiety and among them, the Zn2+ ion was found to be better than Cu2+ and Ni2+ ions.

Computational modeling of stented coronary arteries

The application of computational fluid dynamics (CFD) and finite element analysis (FEA) has been growing rapidly in the various fields of science and technology. One of the areas of interest is in biomedical engineering. The altered hemodynamics inside the blood vessels plays a key role in the development of the arterial disease called atherosclerosis, which is the major cause of human death worldwide. Atherosclerosis is often treated with the stenting procedure to restore the normal blood flow. A stent is a tubular, flexible structure, usually made of metals, which is driven and expanded in the blocked arteries. Despite the success rate of the stenting procedure, it is often associated with the restenosis (re-narrowing of the artery) process. The presence of non-biological device in the artery causes inflammation or re-growth of atherosclerotic lesions in the treated vessels. Several factors including the design of stents, type of stent expansion, expansion pressure, morphology and composition of vessel wall influence the restenosis process. Therefore, the role of computational studies is crucial in the investigation and optimisation of the factors that influence post-stenting complications. This thesis focuses on the stent-vessel wall interactions followed by the blood flow in the post-stenting stage of stenosed human coronary artery. Hemodynamic and mechanical stresses were analysed in three separate stent-plaque-artery models. Plaque was modeled as a multi-layer (fibrous cap (FC), necrotic core (NC), and fibrosis (F)) and the arterial wall as a single layer domain. CFD/FEA simulations were performed using commercial software packages in several models mimicking the various stages and morphologies of atherosclerosis. The tissue prolapse (TP) of stented vessel wall, the distribution of von Mises stress (VMS) inside various layers of vessel wall, and the wall shear stress (WSS) along the luminal surface of the deformed vessel wall were measured and evaluated. The results revealed the role of the stenosis size, thickness of each layer of atherosclerotic wall, thickness of stent strut, pressure applied for stenosis expansion, and the flow condition in the distribution of stresses. The thicknesses of FC, and NC and the total thickness of plaque are critical in controlling the stresses inside the tissue. A small change in morphology of artery wall can significantly affect the distribution of stresses. In particular, FC is the most sensitive layer to TP and stresses, which could determine plaque’s vulnerability to rupture. The WSS is highly influenced by the deflection of artery, which in turn is dependent on the structural composition of arterial wall layers. Together with the stenosis size, their roles could play a decisive role in controlling the low values of WSS (<0.5 Pa) prone to restenosis. Moreover, the time dependent flow altered the percentage of luminal area with WSS values less than 0.5 Pa at different time instants. The non- Newtonian viscosity model of the blood properties significantly affects the prediction of WSS magnitude. The outcomes of this investigation will help to better understand the roles of the individual layers of atherosclerotic vessels and their risk to provoke restenosis at the post-stenting stage. As a consequence, the implementation of such an approach to assess the post-stented stresses will assist the engineers and clinicians in optimizing the stenting techniques to minimize the occurrence of restenosis.

Artificial reflecting structures based on metamaterials

This thesis studies metamaterial-inspired mirrors which provide the most general control over the amplitude and phase of the reflected wavefront. The goal is to explore practical possibilities in designing fully reflective electromagnetic structures with full control over reflection phase. The first part of the thesis describes a planar focusing metamirror with the focal distance less than the operating wavelength. Its practical applicability from the viewpoint of aberrations when the incident angle deviates from the normal one is verified numerically and experimentally. The results indicate that the proposed focusing metamirror can be efficiently employed in many different applications due to its advantages over other conventional mirrors. In the second part of the thesis a new theoretical concept of reflecting metasurface operation is introduced based on Huygens’ principle. This concept in contrast to known approaches takes into account all the requirements of perfect metamirror operation. The theory shows a route to improve the previously proposed metamirrors through tilting the individual inclusions of the structure at a chosen angle from normal. It is numerically tested and the results demonstrate improvements over the previous design.

Modelling the colour difference by artificial neural network

This thesis work studies the modelling of the colour difference using artificial neural network. Multilayer percepton (MLP) network is proposed to model CIEDE2000 colour difference formula. MLP is applied to classify colour points in CIE xy chromaticity diagram. In this context, the evaluation was performed using Munsell colour data and MacAdam colour discrimination ellipses. Moreover, in CIE xy chromaticity diagram just noticeable differences (JND) of MacAdam ellipses centres are computed by CIEDE2000, to compare JND of CIEDE2000 and MacAdam ellipses. CIEDE2000 changes the orientation of blue areas in CIE xy chromaticity diagram toward neutral areas, but on the whole it does not totally agree with the MacAdam ellipses. The proposed MLP for both modelling CIEDE2000 and classifying colour points showed good accuracy and achieved acceptable results.