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.

Congenital upper limb defects in Finland 1993-2005

Background: Interest in limb defects has grown after the thalidomide tragedy in the 1960s. As a result, congenital malformation registries, monitoring changes in birthprevalence and defect patterns, have been established in several countries. However, there are only a few true population based studies on birth prevalence of upper limb defects. The burden of hospital care among these children, specifically in terms of the number of admissions and total time spent in hospital, is also unknown. Aims and Methods: This study is based on information gathered from the Finnish Register of Congenital malformations (FRM) and the Finnish Hospital Discharge Register (FHDR). A total of 417 children born between 1993 and 2005 with an upper limb defect were gathered from the FRM. The upper limb defects were classified using the International Federation of Societies for Surgery of the Hand -classification that enables comparison with previous and future studies. Birth and live birth prevalence, sex and side distribution, frequency of associated anomalies as well as the proportion of perinatal and infant deaths according to the different subtypes were calculated. The number of hospital admissions, days spent in hospital, number and type of surgical operations were collected from the FHDR. Special features of two subgroups, radial ray defects (RRD) and constriction band syndrome (CBS), were explored. Results: Upper limb defects were observed in 417 of 753 342 consecutive births and in 392 of 750 461 live births. Birth prevalence was 5.5 per 10 000 births and 5.2 per 10 000 live births. Multiple anomalies or a known syndrome was found in 250 cases (60%). Perinatal mortality was 139 per 1000 births and infant mortality 135 per 1000 live births (overall Finnish perinatal mortality <5 per 1000 births and infant mortality 3.7 per 1000 live births). Altogether, 138 infants had RRD and 120 (87%) of these had either a known syndrome or multiple major anomalies. The proportion of perinatal deaths in RRD group was 29% (40/138) and infant deaths 35% (43/123). Fifty-one children had CBS in upper limbs. Fifteen of these (29%) had other major anomalies associated with constriction rings. The number of hospital admissions per year of children with congenital upper limb defects was 11-fold and the time spent in hospital 13-fold as compared with the general paediatric population. Conclusions: Birth prevalence of congenital upper limb defects was 5.5 per 10 000 births and 5.2 per 10 000 live births. RRD was especially associated with other major anomalies and high mortality. Nearly one third of the children with CBS also had other major anomalies suggesting different aetiologies inside the group. The annual burden of hospital care of children with congenital upper limb defects was at least 11-fold as compared with the general paediatric population.