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.

Effect of ventilation on systemic blood flow evaluated by echodopplercardiography

Systemic blood flow (Q) was measured by echodopplercardiography in 5 normal young adult males during apnea, eupnea and tachypnea. Measurements were made in a recumbent posture at 3-min intervals. Tachypnea was attained by doubling the respiratory frequency at eupnea at a constant tidal volume. An open glottis was maintained during apnea at the resting respiratory level. The Q values were positively correlated with the respiratory frequency, decreasing from eupnea to apnea and increasing from eupnea to tachypnea (P<0.05). These data demonstrate that echodopplercardiography, a better qualified tool for this purpose, confirms the positive and progressive effects of ventilation on systemic blood flow, as suggested by previous studies based on diverse technical approaches

Validation of transit-time flowmetry for chronic measurements of regional blood flow in resting and exercising rats

The objective of the present study was to validate the transit-time technique for long-term measurements of iliac and renal blood flow in rats. Flow measured with ultrasonic probes was confirmed ex vivo using excised arteries perfused at varying flow rates. An implanted 1-mm probe reproduced with accuracy different patterns of flow relative to pressure in freely moving rats and accurately quantitated the resting iliac flow value (on average 10.43 ± 0.99 ml/min or 2.78 ± 0.3 ml min-1 100 g body weight-1). The measurements were stable over an experimental period of one week but were affected by probe size (resting flows were underestimated by 57% with a 2-mm probe when compared with a 1-mm probe) and by anesthesia (in the same rats, iliac flow was reduced by 50-60% when compared to the conscious state). Instantaneous changes of iliac and renal flow during exercise and recovery were accurately measured by the transit-time technique. Iliac flow increased instantaneously at the beginning of mild exercise (from 12.03 ± 1.06 to 25.55 ± 3.89 ml/min at 15 s) and showed a smaller increase when exercise intensity increased further, reaching a plateau of 38.43 ± 1.92 ml/min at the 4th min of moderate exercise intensity. In contrast, exercise-induced reduction of renal flow was smaller and slower, with 18% and 25% decreases at mild and moderate exercise intensities. Our data indicate that transit-time flowmetry is a reliable method for long-term and continuous measurements of regional blood flow at rest and can be used to quantitate the dynamic flow changes that characterize exercise and recovery

Blood flow measurements in rats using four color microspheres during blockade of different vasopressor systems

The use of colored microspheres to adequately evaluate blood flow changes under different circumstances in the same rat has been validated with a maximum of three different colors due to methodological limitations. The aim of the present study was to validate the use of four different colors measuring four repeated blood flow changes in the same rat to assess the role of vasopressor systems in controlling arterial pressure (AP). Red (150,000), white (200,000), yellow (150,000), and blue (200,000) colored microspheres were infused into the left ventricle of 6 male Wistar rats 1) at rest and 2) after vasopressin (aAVP, 10 µg/kg, iv), 3) renin-angiotensin (losartan, 10 mg/kg, iv), and 4) sympathetic system blockade (hexamethonium, 20 mg/kg, iv) to determine blood flow changes. AP was recorded and processed with a data acquisition system (1-kHz sampling frequency). Blood flow changes were quantified by spectrophotometry absorption peaks for colored microsphere components in the tissues evaluated. Administration of aAVP and losartan slightly reduced the AP (-5.7 ± 0.5 and -7.8 ± 1.2 mmHg, respectively), while hexamethonium induced a 52 ± 3 mmHg fall in AP. The aAVP injection increased blood flow in lungs (78%), liver (117%) and skeletal muscle (>150%), while losartan administration enhanced blood flow in heart (126%), lungs (100%), kidneys (80%), and gastrocnemius (75%) and soleus (94%) muscles. Hexamethonium administration reduced only kidney blood flow (50%). In conclusion, four types of colored microspheres can be used to perform four repeated blood flow measurements in the same rat detecting small alterations such as changes in tissues with low blood flow.

Role of nitric oxide and prostaglandin in the maintenance of cortical and renal medullary blood flow

This study was undertaken in anesthetized dogs to evaluate the relative participation of prostaglandins (PGs) and nitric oxide (NO) in the maintenance of total renal blood flow (TRBF), and renal medullary blood flow (RMBF). It was hypothesized that the inhibition of NO should impair cortical and medullary circulation because of the synthesis of this compound in the endothelial cells of these two territories. In contrast, under normal conditions of perfusion pressure PG synthesis is confined to the renal medulla. Hence PG inhibition should predominantly impair the medullary circulation. The initial administration of 25 µM kg-1 min-1 NG-nitro-L-arginine methyl ester produced a significant 26% decrease in TRBF and a concomitant 34% fall in RMBF, while the subsequent inhibition of PGs with 5 mg/kg meclofenamate further reduced TRBF by 33% and RMBF by 89%. In contrast, the initial administration of meclofenamate failed to change TRBF, while decreasing RMBF by 49%. The subsequent blockade of NO decreased TRBF by 35% without further altering RMBF. These results indicate that initial PG synthesis inhibition predominantly alters the medullary circulation, whereas NO inhibition decreases both cortical and medullary flow. This latter change induced by NO renders cortical and RMBF susceptible to a further decrease by PG inhibition. However, the decrease in medullary circulation produced by NO inhibition is not further enhanced by subsequent PG inhibition.

Effects of disturbed blood flow during exercise on endothelial function: a time course analysis

This study aimed to examine the time course of endothelial function after a single handgrip exercise session combined with blood flow restriction in healthy young men. Nine participants (28±5.8 years) completed a single session of bilateral dynamic handgrip exercise (20 min with 60% of the maximum voluntary contraction). To induce blood flow restriction, a cuff was placed 2 cm below the antecubital fossa in the experimental arm. This cuff was inflated to 80 mmHg before initiation of exercise and maintained through the duration of the protocol. The experimental arm and control arm were randomly selected for all subjects. Brachial artery flow-mediated dilation (FMD) and blood flow velocity profiles were assessed using Doppler ultrasonography before initiation of the exercise, and at 15 and 60 min after its cessation. Blood flow velocity profiles were also assessed during exercise. There was a significant increase in FMD 15 min after exercise in the control arm compared with before exercise (64.09%±16.59%, P=0.001), but there was no change in the experimental arm (-12.48%±12.64%, P=0.252). FMD values at 15 min post-exercise were significantly higher for the control arm in comparison to the experimental arm (P=0.004). FMD returned to near baseline values at 60 min after exercise, with no significant difference between arms (P=0.424). A single handgrip exercise bout provoked an acute increase in FMD 15 min after exercise, returning to near baseline values at 60 min. This response was blunted by the addition of an inflated pneumatic cuff to the exercising arm.

Examining the roles of core and local temperature on forearm skin blood flow

The interaction between local and reflexive control of skin blood flow (SkBF) is unclear. This thesis isolated the roles of rectal (Tre) and local (Tloc) temperature on forearm SkBF regulation at normal and elevated body temperatures, and to investigate the interaction between local and reflexive SkBF control. While either normothermic (Tre ~37.0°C) or hyperthermic (∆Tre +1.1°C), SkBF was assessed on the dorsal aspect of each forearm in 10 participants while Tloc was manipulated in an A-B-A-B fashion between neutral (33.0°C) and hot (38.5°C). Finally, local heating to 44°C was performed to elicit maximal SkBF. Data are presented as a percentage of maximal cutaneous vascular conductance (CVC), calculated as laser-Doppler flux divided by mean arterial pressure. Tloc manipulations performed during normothermia had significantly greater effects on CVC than during hyperthermia. The decreased modification to SkBF from the Tloc changes during hyperthermia suggests that strong reflexive vasodilation attenuates local SkBF control mechanisms.

Effects of Cerebral Blood Flow and PETCO2 on Cognitive Function During Passive Heat Stress

This thesis tested whether cognitive performance during passive heat stress may be affected by changes in cerebrovascular variables as opposed to strictly thermally-induced changes. A pharmacological reduction in cerebral blood flow (CBF) using indomethacin along with a hypocapnia-induced CBF reduction during passive heat stress (Tre ~1.5°C above baseline) were used to investigate any cerebrovascular-mediated changes in cognitive performance. Repeated measures analysis of variance indicated that One-Touch Stockings of Cambridge (OTS) performance was not affected by a significant reduction in CBF during passive heat stress. More specifically, OTS accuracy measures did not change as a result of either a reduction in CBF or increasing passive heat stress. However, it was found that OTS response time indices improved with increasing passive heat stress independent of CBF changes. In conclusion, a significant reduction in CBF does not cause additional changes in performance of an executive functioning task during severe passive heat stress.

The Influence of Cerebral Blood Flow and Carbon Dioxide on Neuromuscular Responses During Environmental Stress

Although reductions in cerebral blood flow (CBF) may be implicated in the development of central fatigue during environmental stress, the contribution from hypocapnia-induced reductions in CBF versus reductions in CBF per se has yet to be isolated. The current research program examined the influence of CBF, with and without consequent hypocapnia, on neuromuscular responses during hypoxia and passive heat stress. To this end, neuromuscular responses, as indicated by motor evoked potentials (MEP), maximal M-wave (Mmax) and cortical voluntary activation (cVA) of the flexor carpi radialis muscle during isometric wrist flexion, was assessed in three separate projects: 1) hypocapnia, independent of concomitant reductions in CBF; 2) altered CBF during severe hypoxia and; 3) thermal hyperpnea-mediated reductions in CBF, independent of hypocapnia. All projects employed a custom-built dynamic end-tidal forcing system to control end-tidal PCO2 (PETCO2), independent of the prevailing environmental conditions, and cyclooxygenase inhibition using indomethacin (Indomethacin, 1.2 mg·Kg-1) to selectively reduce CBF (estimated using transcranial Doppler ultrasound) without changes in PETCO2. A primary finding of the present research program is that the excitability of the corticospinal tract is inherently sensitive to changes in PaCO2, as demonstrated by a 12% increase in MEP amplitude in response to moderate hypocapnia. Conversely, CBF mediated reductions in cerebral O2 delivery appear to decrease corticospinal excitability, as indicated by a 51-64% and 4% decrease in MEP amplitude in response to hypoxia and passive heat stress, respectively. The collective evidence from this research program suggests that impaired voluntary activation is associated with reductions in CBF; however, it must be noted that changes in cVA were not linearly correlated with changes in CBF. Therefore, other factors independent of CBF, such as increased perception of effort, distress or discomfort, may have contributed to the reductions in cVA. Despite the functional association between reductions in CBF and hypocapnia, both variables have distinct and independent influence on the neuromuscular system. Therefore, future studies should control or acknowledge the separate mechanistic influence of these two factors.

Effect of hormone replacement therapy on retinal and optic nerve head blood flow and topography in postmenopausal women, and retinal tissue perfusion in ovariectomized rats

Thèse numérisée par la Division de la gestion de documents et des archives de l'Université de Montréal.

Evidence for a reversal with age in the pattern of near-wall blood flow around aortic branches

The changes that occur with age in the distribution of atherosclerotic lesions around arterial branch points challenge accepted theories relating disease to haemodynamic stresses. We investigated whether flow near branch points changes with age in a way that can account for the different lesion distributions. Flow around 20 branches from immature and mature aortas was investigated by examining the length:width ratio and orientation of endothelial nuclei; these properties depend on the magnitude and direction of near-wall flows, respectively. There were significant changes in the pattern of nuclear shape with age, consistent with a reversal in the pattern of shear around branches. In control regions away from branches, there were no such changes. The role of haemodynamic stresses in atherogenesis may require re-evaluation in the light of these results. (C) 2003 Elsevier Ireland Ltd. All rights reserved.

Evidence for a reversal with age in the pattern of near-wall blood flow around aortic branches

The changes that occur with age in the distribution of atherosclerotic lesions around arterial branch points challenge accepted theories relating disease to haemodynamic stresses. We investigated whether flow near branch points changes with age in a way that can account for the different lesion distributions. Flow around 20 branches from immature and mature aortas was investigated by examining the length:width ratio and orientation of endothelial nuclei; these properties depend on the magnitude and direction of near-wall flows, respectively. There were significant changes in the pattern of nuclear shape with age, consistent with a reversal in the pattern of shear around branches. In control regions away from branches, there were no such changes. The role of haemodynamic stresses in atherogenesis may require re-evaluation in the light of these results.

A time-invariant visco-elastic windkessel model relating blood flow and blood volume

The difference between the rate of change of cerebral blood volume (CBV) and cerebral blood flow (CBF) following stimulation is thought to be due to circumferential stress relaxation in veins (Mandeville, J.B., Marota, J.J.A., Ayata, C., Zaharchuk, G., Moskowitz, M.A., Rosen, B.R., Weisskoff, R.M., 1999. Evidence of a cerebrovascular postarteriole windkessel with delayed compliance. J. Cereb. Blood Flow Metab. 19, 679–689). In this paper we explore the visco-elastic properties of blood vessels, and present a dynamic model relating changes in CBF to changes in CBV. We refer to this model as the visco-elastic windkessel (VW) model. A novel feature of this model is that the parameter characterising the pressure–volume relationship of blood vessels is treated as a state variable dependent on the rate of change of CBV, producing hysteresis in the pressure–volume space during vessel dilation and contraction. The VW model is nonlinear time-invariant, and is able to predict the observed differences between the time series of CBV and that of CBF measurements following changes in neural activity. Like the windkessel model derived by Mandeville, J.B., Marota, J.J.A., Ayata, C., Zaharchuk, G., Moskowitz, M.A., Rosen, B.R., Weisskoff, R.M., 1999. Evidence of a cerebrovascular postarteriole windkessel with delayed compliance. J. Cereb. Blood Flow Metab. 19, 679–689, the VW model is primarily a model of haemodynamic changes in the venous compartment. The VW model is demonstrated to have the following characteristics typical of visco-elastic materials: (1) hysteresis, (2) creep, and (3) stress relaxation, hence it provides a unified model of the visco-elastic properties of the vasculature. The model will not only contribute to the interpretation of the Blood Oxygen Level Dependent (BOLD) signals from functional Magnetic Resonance Imaging (fMRI) experiments, but also find applications in the study and modelling of the brain vasculature and the haemodynamics of circulatory and cardiovascular systems.

Model estimation of cerebral hemodynamics between blood flow and volume changes: A data-based modeling approach

It is well known that there is a dynamic relationship between cerebral blood flow (CBF) and cerebral blood volume (CBV). With increasing applications of functional MRI, where the blood oxygen-level-dependent signals are recorded, the understanding and accurate modeling of the hemodynamic relationship between CBF and CBV becomes increasingly important. This study presents an empirical and data-based modeling framework for model identification from CBF and CBV experimental data. It is shown that the relationship between the changes in CBF and CBV can be described using a parsimonious autoregressive with exogenous input model structure. It is observed that neither the ordinary least-squares (LS) method nor the classical total least-squares (TLS) method can produce accurate estimates from the original noisy CBF and CBV data. A regularized total least-squares (RTLS) method is thus introduced and extended to solve such an error-in-the-variables problem. Quantitative results show that the RTLS method works very well on the noisy CBF and CBV data. Finally, a combination of RTLS with a filtering method can lead to a parsimonious but very effective model that can characterize the relationship between the changes in CBF and CBV.