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.

Spondyloosi hevosen rintarangan etuosassa : tapausselostus

Summary: Spondylosis in the equine cranial thoracic spine : case report

Modeling the transport phenomena within arterial wall: porous media approach

The transport of macromolecules, such as low-density lipoprotein (LDL), and their accumulation in the layers of the arterial wall play a critical role in the creation and development of atherosclerosis. Atherosclerosis is a disease of large arteries e.g., the aorta, coronary, carotid, and other proximal arteries that involves a distinctive accumulation of LDL and other lipid-bearing materials in the arterial wall. Over time, plaque hardens and narrows the arteries. The flow of oxygen-rich blood to organs and other parts of the body is reduced. This can lead to serious problems, including heart attack, stroke, or even death. It has been proven that the accumulation of macromolecules in the arterial wall depends not only on the ease with which materials enter the wall, but also on the hindrance to the passage of materials out of the wall posed by underlying layers. Therefore, attention was drawn to the fact that the wall structure of large arteries is different than other vessels which are disease-resistant. Atherosclerosis tends to be localized in regions of curvature and branching in arteries where fluid shear stress (shear rate) and other fluid mechanical characteristics deviate from their normal spatial and temporal distribution patterns in straight vessels. On the other hand, the smooth muscle cells (SMCs) residing in the media layer of the arterial wall respond to mechanical stimuli, such as shear stress. Shear stress may affect SMC proliferation and migration from the media layer to intima. This occurs in atherosclerosis and intimal hyperplasia. The study of blood flow and other body fluids and of heat transport through the arterial wall is one of the advanced applications of porous media in recent years. The arterial wall may be modeled in both macroscopic (as a continuous porous medium) and microscopic scales (as a heterogeneous porous medium). In the present study, the governing equations of mass, heat and momentum transport have been solved for different species and interstitial fluid within the arterial wall by means of computational fluid dynamics (CFD). Simulation models are based on the finite element (FE) and finite volume (FV) methods. The wall structure has been modeled by assuming the wall layers as porous media with different properties. In order to study the heat transport through human tissues, the simulations have been carried out for a non-homogeneous model of porous media. The tissue is composed of blood vessels, cells, and an interstitium. The interstitium consists of interstitial fluid and extracellular fibers. Numerical simulations are performed in a two-dimensional (2D) model to realize the effect of the shape and configuration of the discrete phase on the convective and conductive features of heat transfer, e.g. the interstitium of biological tissues. On the other hand, the governing equations of momentum and mass transport have been solved in the heterogeneous porous media model of the media layer, which has a major role in the transport and accumulation of solutes across the arterial wall. The transport of Adenosine 5´-triphosphate (ATP) is simulated across the media layer as a benchmark to observe how SMCs affect on the species mass transport. In addition, the transport of interstitial fluid has been simulated while the deformation of the media layer (due to high blood pressure) and its constituents such as SMCs are also involved in the model. In this context, the effect of pressure variation on shear stress is investigated over SMCs induced by the interstitial flow both in 2D and three-dimensional (3D) geometries for the media layer. The influence of hypertension (high pressure) on the transport of lowdensity lipoprotein (LDL) through deformable arterial wall layers is also studied. This is due to the pressure-driven convective flow across the arterial wall. The intima and media layers are assumed as homogeneous porous media. The results of the present study reveal that ATP concentration over the surface of SMCs and within the bulk of the media layer is significantly dependent on the distribution of cells. Moreover, the shear stress magnitude and distribution over the SMC surface are affected by transmural pressure and the deformation of the media layer of the aorta wall. This work reflects the fact that the second or even subsequent layers of SMCs may bear shear stresses of the same order of magnitude as the first layer does if cells are arranged in an arbitrary manner. This study has brought new insights into the simulation of the arterial wall, as the previous simplifications have been ignored. The configurations of SMCs used here with elliptic cross sections of SMCs closely resemble the physiological conditions of cells. Moreover, the deformation of SMCs with high transmural pressure which follows the media layer compaction has been studied for the first time. On the other hand, results demonstrate that LDL concentration through the intima and media layers changes significantly as wall layers compress with transmural pressure. It was also noticed that the fraction of leaky junctions across the endothelial cells and the area fraction of fenestral pores over the internal elastic lamina affect the LDL distribution dramatically through the thoracic aorta wall. The simulation techniques introduced in this work can also trigger new ideas for simulating porous media involved in any biomedical, biomechanical, chemical, and environmental engineering applications.

Performance and Methodological Aspects in Positron Emission Tomography

Performance standards for Positron emission tomography (PET) were developed to be able to compare systems from different generations and manufacturers. This resulted in the NEMA methodology in North America and the IEC in Europe. In practices, the NEMA NU 2- 2001 is the method of choice today. These standardized methods allow assessment of the physical performance of new commercial dedicated PET/CT tomographs. The point spread in image formation is one of the factors that blur the image. The phenomenon is often called the partial volume effect. Several methods for correcting for partial volume are under research but no real agreement exists on how to solve it. The influence of the effect varies in different clinical settings and it is likely that new methods are needed to solve this problem. Most of the clinical PET work is done in the field of oncology. The whole body PET combined with a CT is the standard investigation today in oncology. Despite the progress in PET imaging technique visualization, especially quantification of small lesions is a challenge. In addition to partial volume, the movement of the object is a significant source of error. The main causes of movement are respiratory and cardiac motions. Most of the new commercial scanners are in addition to cardiac gating, also capable of respiratory gating and this technique has been used in patients with cancer of the thoracic region and patients being studied for the planning of radiation therapy. For routine cardiac applications such as assessment of viability and perfusion only cardiac gating has been used. However, the new targets such as plaque or molecular imaging of new therapies require better control of the cardiac motion also caused by respiratory motion. To overcome these problems in cardiac work, a dual gating approach has been proposed. In this study we investigated the physical performance of a new whole body PET/CT scanner with NEMA standard, compared methods for partial volume correction in PET studies of the brain and developed and tested a new robust method for dual cardiac-respiratory gated PET with phantom, animal and human data. Results from performance measurements showed the feasibility of the new scanner design in 2D and 3D whole body studies. Partial volume was corrected, but there is no best method among those tested as the correction also depends on the radiotracer and its distribution. New methods need to be developed for proper correction. The dual gating algorithm generated is shown to handle dual-gated data, preserving quantification and clearly eliminating the majority of contraction and respiration movement

Fractures in older people - incidence, predictors and consequences

In older populations, fractures are common and the consequences of fractures may be serious both for an individual and for society. However, information is scarce about the incidence, predictors and consequences of fractures in population-based unselected cohorts including both men and women and a long follow-up. The objective of this study was to analyse the incidence and predictors of fractures as well as functional decline and excess mortality due to fractures, among 482 men and 695 women aged 65 or older in the municipality of Lieto, Finland from 1991 until 2002. In analyses, Poisson’s, Cox proportional Hazards and Cumulative Logistic regression models were used for the control of several confounding variables. During the 12-year follow-up with a total of 10 040 person-years (PY), 307 (26%) persons sustained altogether 425 fractures of which 77% were sustained by women. The total incidence of fractures was 53.4 per 1000 PY (95% confidence intervals [95% CI]: 47.9 - 59.5) in women and 24.9 per 1000 PY (95% CI: 20.4 - 30.4) in men. The incidence rates of fractures at any sites and hip fractures were associated with increasing age. No significant changes in the ageadjusted incidence rates of fractures were found in either gender during the 12-year follow-up. The predictors of fractures varied by gender. In multivariate analyses, reduced handgrip strength and body mass index (BMI) lower than 30 in women and a large number of depressive symptoms in men were independent predictors of fractures. A compression fracture in one or more thoracic or upper lumbar vertebras on chest radiography at baseline was associated with subsequent fractures in both genders. Lower body fractures independently predicted both short- (0-2 years) and long-term (up to 8 years) functional decline in mobility and activities of daily living (ADL) performance during the 8-year follow-up. Upper body fractures predicted decline in ADL performance during longterm follow-up. In the 12-year follow-up, hip fractures in men (Hazard Ratio [HR] 8.1, 95% CI: 4.4-14.9) and in women (HR 3.0, 95% CI: 1.9-4.9), and fractures at the proximal humerus in men (HR 5.4, 95% CI: 1.6-17.7) were independently associated with excess mortality. In addition, leisure time inactivity in physical exercise predicted independently both functional decline and excess mortality. Fractures are common among older people posing serious individual consequences. Further studies about the effectiveness of preventing falls and fractures as well as improving care and rehabilitation after fractures are needed.

Motoneuronitautiin liittyvä hengitysvajaus: hengitystoiminnan ja energia-aineenvaihdunnan mittaaminen

Tavoitteet: Tämän tutkimussarjan tavoitteena oli tutkia hengitystoiminnan sekä energia-aineen¬vaihdunnan muutoksia motoneuronitautia (amyotrofinen lateraaliskleroosi, ALS) sairastavilla potilailla. Erityisenä mielenkiinnon kohteena olivat kotihoitoon soveltuvan hengityslaitteen vai¬kutus elinajan ennusteeseen sekä hengitysvajauksen etenemistä kuvaavien keuhkotoimintakokei¬den arviointi ALS-potilailla, epäsuoran kalorimetrian mittaustarkkuus ja perusaineenvaihdunnan (PAV) suuruus kajoavaa hengityslaitetta käyttävillä ALS-potilailla. Aineisto ja menetelmät: Kajoamattoman hengityslaitteen käytön ja iän vaikutusta elinajan en¬nusteeseen arvioitiin 84:llä ja hengitystoiminnan muutoksia 42 ALS-potilaalla. Epäsuoran kalo¬rimetrian mittaustarkkuutta kajoamatonta hengityslaitetta käytettäessä arvioitiin hereillä olevilla 12 vapaaehtoisella mieshenkilöllä. PAV:n suuruutta arvioitiin viidellä kajoavaa hengityslaitetta käyttävällä ALS-potilaalla. Osatöistä kaksi ensimmäistä olivat luonteeltaan havainnoivia (retros¬pektiivisiä) ja kaksi viimeistä seurantatutkimuksia (prospektiivisia). Tulokset: Alle 65-vuotiailla ALS-potilailla ei havaittu eroa elinajan ennusteessa kajoamaton¬ta hengityslaitetta käyttävien ja käyttämättömien potilaiden välillä. Sen sijaan yli 65-vuotiail¬la ALS-potilailla elinajan ennuste piteni merkittävästi kajoamatonta hengityslaitetta käyttävillä potilailla (elinaika diagnoosin jälkeen 22 vs. 8 kk, Hazard Ratio = 0.25, 95 % luottamusväli 0.11 – 0.55, p <0.001). ALS-potilailla, joilla kajoamaton hengityslaite katsottiin tarpeelliseksi kuuden kuukauden kuluessa diagnoosihetkestä, hengitystiheys osoittautui diagnoosihetkellä mer-kittävästi kiihtyneeksi (21/min) ja rintakehän liike merkittävästi alentuneeksi (2.9 cm) verrattuna ALS-potilaisiin, joille kajoamaton hengityslaite katsottiin tarpeelliseksi myöhemmin (16/min ja 4.0 cm). Kajoamattoman hengityslaitehoidon aikana keskimääräinen mitattu PAV vapaaehtoisilla miehillä oli 1858 kcal/vrk kun PAV ilman hengityslaitetta oli 1852 kcal/vrk, p = 0.8. Kajoavaa hengityslaitehoitoa käyttävien viiden ALS-potilaan keskimääräinen PAV vastaavalla mittausase¬telmalla mitattaessa oli 1130 kcal/vrk, kun vastaava PAV laskettuna viidellä eri laskentakaavalla oli 1700 kcal/vrk, p < 0.001. Johtopäätökset: Yli 65-vuotiailla ALS-potilailla, jotka eivät sopeutuneet kajoamattomaan hen¬gityslaitehoitoon, oli nelinkertainen riski menehtyä aiemmin kuin kajoamattomaan hengityslai¬tehoitoon sopeutuneilla ALS-potilailla. Hengitystiheys osoittautui merkittävästi kiihtyneeksi ja rintakehän liike alentuneeksi ALS-potilailla, joille kajoamaton hengityslaitehoito katsottiin ai¬heelliseksi kuuden kuukauden kuluessa diagnoosihetkestä. Kajoamattoman hengityslaitehoidon aikana mitattu PAV ei poikennut mitatusta PAV:sta itsenäisen hengityksen aikana. Näin ollen epäsuoraa kalorimetriamenetelmää voidaan käyttää luotettavasti PAV:n määrittämiseen käytet¬täessä samanaikaisesti kotihoitoon soveltuvaa hengityslaitehoitoa. Elämää ylläpitävää kajoavaa hengityslaitehoitoa käyttävien ALS-potilaiden PAV oli merkittävästi hidastunut laskennallisella menetelmällä arvioituun PAV verrattuna.