Hemodynamics in the abdominal aorta: a comparison of in vitro and in vivo measurements

Multi-Objective Optimisation of a Novel Bypass Graft with a Spiral Ridge

The low long-term patency of bypass grafts is a major concern for cardiovascular treatments. Unfavourable haemodynamic conditions in the proximity of distal anastomosis are closely related to thrombus creation and lumen lesions. Modern graft designs address this unfavourable haemodynamic environment with the introduction of a helical component in the flow field, either by means of out-of-plane helicity graft geometry or a spiral ridge. While the latter has been found to lack in performance when compared to the out-of-plane helicity designs, recent findings support the idea that the existing spiral ridge grafts can be further improved in performance through optimising relevant design parameters. In the current study, robust multi-objective optimisation techniques are implemented, covering a wide range of possible designs coupled with proven and well validated computational fluid dynamics (CFD) algorithms. It is shown that the final set of suggested design parameters could significantly improve haemodynamic performance and therefore could be used to enhance the design of spiral ridge bypass grafts.

How does hemodynamics affect rupture tissue mechanics in abdominal aortic aneurysm: Focus on wall shear stress derived parameters, time-averaged wall shear stress, oscillatory shear index, endothelial cell activation potential, and relative residence time

An abdominal aortic aneurysm (AAA) is a critical health condition with a risk of rupture, where the diameter of the aorta enlarges more than 50% of its normal diameter. The incidence rate of AAA has increased worldwide. Currently, about three out of every 100,000 people have aortic diseases. The diameter and geometry of AAAs influence the hemodynamic forces exerted on the arterial wall. Therefore, a reliable assessment of hemodynamics is crucial for predicting the rupture risk. Wall shear stress (WSS) is an important metric to define the level of the frictional force on the AAA wall. Excessive levels of WSS deteriorate the remodeling mechanism of the arteries and lead to abnormal conditions. At this point, WSS-related hemodynamic parameters, such as time-averaged WSS (TAWSS), oscillatory shear index (OSI), endothelial cell activation potential (ECAP), and relative residence time (RRT) provide important information to evaluate the shear environment on the AAA wall in detail. Calculation of these parameters is not straightforward and requires a physical understanding of what they represent. In addition, computational fluid dynamics (CFD) solvers do not readily calculate these parameters when hemodynamics is simulated. This review aims to explain the WSS-derived parameters focusing on how these represent different characteristics of disturbed hemodynamics. A representative case is presented for spatial and temporal formulation that would be useful for interested researchers for practical calculations. Finally, recent hemodynamics investigations relating WSS-related parameters with AAA rupture risk assessment are presented. This review will be useful to understand the physical representation of WSS-related parameters in cardiovascular flows and how they can be calculated practically for AAA investigations.

Comparison of Prospective and Retrospective Gated 4D Flow Cardiac MR Image Acquisitions in the Carotid Bifurcation

PURPOSE

To evaluate the agreement of 4D flow cMRI-derived bulk flow features and fluid (blood) velocities in the carotid bifurcation using prospective and retrospective gating techniques.

METHODS

Prospective and retrospective ECG-gated three-dimensional (3D) cine phase-contrast cardiac MRI with three-direction velocity encoding (i.e., 4D flow cMRI) data were acquired in ten carotid bifurcations from men (n = 3) and women (n = 2) that were cardiovascular disease-free. MRI sequence parameters were held constant across all scans except temporal resolution values differed. Velocity data were extracted from the fluid domain and evaluated across the entire volume or at defined anatomic planes (common, internal, external carotid arteries). Qualitative agreement between gating techniques was performed by visualizing flow streamlines and topographical images, and statistical comparisons between gating techniques were performed across the fluid volume and defined anatomic regions.

RESULTS

Agreement in the kinematic data (e.g., bulk flow features and velocity data) were observed in the prospectively and retrospectively gated acquisitions. Voxel differences in time-averaged, peak systolic, and diastolic-averaged velocity magnitudes between gating techniques across all volunteers were 2.7%, 1.2%, and 6.4%, respectively. No significant differences in velocity magnitudes or components ([Formula: see text], [Formula: see text], [Formula: see text]) were observed. Importantly, retrospective acquisitions captured increased retrograde flow in the internal carotid artery (i.e., carotid sinus) compared to prospective acquisitions (10.4 ± 6.3% vs. 4.6 ± 5.3%; [Formula: see text] < 0.05).

CONCLUSION

Prospective and retrospective ECG-gated 4D flow cMRI acquisitions provide comparable evaluations of fluid velocities, including velocity vector components, in the carotid bifurcation. However, the increased temporal coverage of retrospective acquisitions depicts increased retrograde flow patterns (i.e., disturbed flow) not captured by the prospective gating technique.

Effects of abnormal vertebral arteries and the circle of Willis on vertebrobasilar dolichoectasia: A multi-scale simulation study

BACKGROUND

Vertebrobasilar dolichoectasia is a rare cerebrovascular disease characterized by obvious extension, dilation and tortuosity of vertebrobasilar artery, and its pathophysiological mechanism is not clear. This study focused on local hemodynamic changes in basilar arteries with typical vertebrobasilar dolichoectasia, together with unbalanced vertebral arteries and abnormal structures of the circle of Willis, through multi-scale modeling.

METHODS

Three-dimensional models of 3 types of vertebrobasilar arteries were constructed from magnetic resonance images. The first type has no vertebrobasilar dolichoectasia, the second type has vertebrobasilar dolichoectasia and balanced vertebral arteries, and the third type has vertebrobasilar dolichoectasia and unbalanced vertebral arteries. A lumped parameter model of the circle of Willis was established and coupled to these three-dimensional models.

FINDINGS

The results showed that unbalanced bilateral vertebral arteries, especially single vertebral artery deletion mutation, might associate with higher wall shear stress on anterior wall of basilar artery in patients with vertebrobasilar dolichoectasia. And unbalanced bilateral vertebral arteries would increase the blood pressure in basilar artery. Meanwhile, missing communicating arteries in the circle of Willis, especially bilateral posterior communicating arteries absences, would significantly increase blood pressure in basilar artery. The unilateral absence of posterior communicating arteries would increase differences in blood flow between the left and right posterior cerebral arteries.

INTERPRETATION

This study provided a multi-scale modeling method and some preliminary results for helping understand the role of hemodynamics in occurrence and development of vertebrobasilar dolichoectasia.

A New Method for Quantifying Abdominal Aortic Wall Shear Stress Using Phase Contrast Magnetic Resonance Imaging and the Womersley Solution

Wall shear stress (WSS) is an important mediator of cardiovascular pathologies and there is a need for its reliable evaluation as a potential prognostic indicator. The purpose of this work was to develop a method that quantifies WSS from two-dimensional (2D) phase contrast magnetic resonance (PCMR) imaging derived flow waveforms, apply this method to PCMR data acquired in the abdominal aorta of healthy volunteers, and to compare PCMR-derived WSS values to values predicted from a computational fluid dynamics (CFD) simulation. The method uses PCMR-derived flow versus time waveforms constrained by the Womersley solution for pulsatile flow in a cylindrical tube. The method was evaluated for sensitivity to input parameters, intrastudy repeatability and was compared with results from a patient-specific CFD simulation. 2D-PCMR data were acquired in the aortas of healthy men (n = 12) and women (n = 15) and time-averaged WSS (TAWSS) was compared. Agreement was observed when comparing TAWSS between CFD and the PCMR flow-based method with a correlation coefficient of 0.88 (CFD: 15.0 ± 1.9 versus MRI: 13.5 ± 2.4 dyn/cm2) though comparison of WSS values between the PCMR-based method and CFD predictions indicate that the PCMR method underestimated instantaneous WSS by 3.7 ± 7.6 dyn/cm2. We found no significant difference in TAWSS magnitude between the sexes; 8.19 ± 2.25 versus 8.07 ± 1.71 dyn/cm2, p = 0.16 for men and women, respectively.

Spatiotemporal changes of local hemodynamics and plaque components during atherosclerotic progression in rabbit

BACKGROUND AND OBJECTIVE

Recent evidence demonstrates that the atherogenic process is discontinuous. Our goal is to study changes of plaque components and local hemodynamics during atherosclerotic progression.

METHODS

The histological and immunohistochemical staining of high-fat diet rabbit aorta were evaluated at 0, 8, 10 and 12 weeks, respectively. In addition, the blood flow and LDL transport were simulated at the above four time points.

RESULTS

The plaque thickness at different characteristic regions increased at different rates. The collagen continued to increase, while the elastin, fibronectin, macrophages and smooth muscle cells increased first and then decreased. The relative surface LDL concentration decreased at 8 weeks, and then it increased first and decreased slightly. Meanwhile, the hemodynamic environment became better firstly at 8 weeks, then got slightly worse and lastly improved again.

CONCLUSIONS

The local hemodynamics and plaque components vary nonlinearly during atherosclerotic progression in rabbit aorta.

Identification of the haemodynamic environment permissive for plaque erosion

Endothelial erosion of atherosclerotic plaques is the underlying cause of approximately 30% of acute coronary syndromes (ACS). As the vascular endothelium is profoundly affected by the haemodynamic environment to which it is exposed, we employed computational fluid dynamic (CFD) analysis of the luminal geometry from 17 patients with optical coherence tomography (OCT)-defined plaque erosion, to determine the flow environment permissive for plaque erosion. Our results demonstrate that 15 of the 17 cases analysed occurred on stenotic plaques with median 31% diameter stenosis (interquartile range 28–52%), where all but one of the adherent thrombi located proximal to, or within the region of maximum stenosis. Consequently, all flow metrics related to elevated flow were significantly increased (time averaged wall shear stress, maximum wall shear stress, time averaged wall shear stress gradient) with a reduction in relative residence time, compared to a non-diseased reference segment. We also identified two cases that did not exhibit an elevation of flow, but occurred in a region exposed to elevated oscillatory flow. Our study demonstrates that the majority of OCT-defined erosions occur where the endothelium is exposed to elevated flow, a haemodynamic environment known to evoke a distinctive phenotypic response in endothelial cells.

Abnormal Flow Dynamics Result in Low Wall Shear Stress and High Oscillatory Shear Index in Abdominal Aortic Dilatation: Initial in vivo Assessment with 4D-flow MRI

PURPOSE

To characterize the non-laminar flow dynamics and resultant decreased wall shear stress (WSS) and high oscillatory shear index (OSI) of the infrarenal abdominal aortic dilatation, cardiac phase-resolved 3D phase-contrast MRI (4D-flow MRI) was performed.

METHODS

The prospective single-arm study was approved by the Institutional Review Board and included 18 subjects (median 67.5 years) with the dilated infrarenal aorta (median diameter 35 mm). 4D-flow MRI was conducted on a 1.5T MRI system. On 3D streamline images, laminar and non-laminar (i.e., vortex or helical) flow patterns were visually assessed both for the dilated aorta and for the undilated upstream aorta. Cardiac phase-resolved flow velocities, WSS and OSI, were also measured for the dilated aorta and the upstream undilated aorta.

RESULTS

Non-laminar flow represented by vortex or helical flow was more frequent and overt in the dilated aorta than in the undilated upstream aorta (P < 0.0156) with a very good interobserver agreement (weighted kappa: 0.82-1.0). The WSS was lower, and the OSI was higher on the dilated aortic wall compared with the proximal undilated segments. In mid-systole, mean spatially-averaged WSS was 0.20 ± 0.016 Pa for the dilated aorta vs. 0.68 ± 0.071 Pa for undilated upstream aorta (P < 0.0001), and OSI on the dilated aortic wall was 0.093 ± 0.010 vs. 0.041 ± 0.0089 (P = 0.013). The maximum values and the amplitudes of the WSS at the dilated aorta were inversely proportional to the ratio of dilated/undilated aortic diameter (r = -0.694, P = 0.0014).

CONCLUSION

4D-flow can characterize abnormal non-laminar flow dynamics within the dilated aorta in vivo. The wall of the infrarenal aortic dilatation is continuously and increasingly affected by atherogenic stimuli due to the flow disturbances represented by vortex or helical flow, which is reflected by lower WSS and higher OSI.

Consensus-based technical recommendations for clinical translation of renal ASL MRI

OBJECTIVES

This study aimed at developing technical recommendations for the acquisition, processing and analysis of renal ASL data in the human kidney at 1.5 T and 3 T field strengths that can promote standardization of renal perfusion measurements and facilitate the comparability of results across scanners and in multi-centre clinical studies.

METHODS

An international panel of 23 renal ASL experts followed a modified Delphi process, including on-line surveys and two in-person meetings, to formulate a series of consensus statements regarding patient preparation, hardware, acquisition protocol, analysis steps and data reporting.

RESULTS

Fifty-nine statements achieved consensus, while agreement could not be reached on two statements related to patient preparation. As a default protocol, the panel recommends pseudo-continuous (PCASL) or flow-sensitive alternating inversion recovery (FAIR) labelling with a single-slice spin-echo EPI readout with background suppression and a simple but robust quantification model.

DISCUSSION

This approach is considered robust and reproducible and can provide renal perfusion images of adequate quality and SNR for most applications. If extended kidney coverage is desirable, a 2D multislice readout is recommended. These recommendations are based on current available evidence and expert opinion. Nonetheless they are expected to be updated as more data become available, since the renal ASL literature is rapidly expanding.

Consensus-based technical recommendations for clinical translation of renal ASL MRI

Objectives

This study aimed at developing technical recommendations for the acquisition, processing and analysis of renal ASL data in the human kidney at 1.5 T and 3 T field strengths that can promote standardization of renal perfusion measurements and facilitate the comparability of results across scanners and in multi-centre clinical studies.

Methods

An international panel of 23 renal ASL experts followed a modified Delphi process, including on-line surveys and two in-person meetings, to formulate a series of consensus statements regarding patient preparation, hardware, acquisition protocol, analysis steps and data reporting.

Results

Fifty-nine statements achieved consensus, while agreement could not be reached on two statements related to patient preparation. As a default protocol, the panel recommends pseudo-continuous (PCASL) or flow-sensitive alternating inversion recovery (FAIR) labelling with a single-slice spin-echo EPI readout with background suppression and a simple but robust quantification model.

Discussion

This approach is considered robust and reproducible and can provide renal perfusion images of adequate quality and SNR for most applications. If extended kidney coverage is desirable, a 2D multislice readout is recommended. These recommendations are based on current available evidence and expert opinion. Nonetheless they are expected to be updated as more data become available, since the renal ASL literature is rapidly expanding.

Electronic supplementary material

The online version of this article (10.1007/s10334-019-00800-z) contains supplementary material, which is available to authorized users.

Volumetric abdominal perfusion measurement using a pseudo-randomly sampled 3D fast-spin-echo (FSE) arterial spin labeling (ASL) sequence and compressed sensing reconstruction

PURPOSE

To improve image quality and spatial coverage for abdominal perfusion imaging by implementing an arterial spin labeling (ASL) sequence that combines variable-density 3D fast-spin-echo (FSE) with Cartesian trajectory and compressed-sensing (CS) reconstruction.

METHODS

A volumetric FSE sequence was modified to include background-suppressed pseudo-continuous ASL labeling and to support variable-density (VD) Poisson-disk sampling for acceleration. We additionally explored the benefits of center oversampling and variable outer k-space sampling. Fourteen healthy volunteers were scanned on a 3T scanner to test acceleration factors as well as the various sampling schemes described under synchronized-breathing to limit motion issues. A CS reconstruction was implemented using the BART toolbox to reconstruct perfusion-weighted ASL volumes, assessing the impact of acceleration, different reconstruction, and sampling strategies on image quality.

RESULTS

CS acceleration is feasible with ASL, and a strong renal perfusion signal could be observed even at very high acceleration rates (≈15). We have shown that ASL k-space complex subtraction was desirable before CS reconstruction. Although averaging of multiple highly accelerated images helped to reduce artifacts from physiologic fluctuations, superior image quality was achieved by interleaving of different highly undersampled pseudo-random spatial sampling patterns and using 4D-CS reconstruction. Combination of these enhancements produces high-quality ASL volumes in under 5 min.

CONCLUSIONS

High-quality isotropic ASL abdominal perfusion volumes can be obtained in healthy volunteers with a VD-FSE and CS reconstruction. This lays the groundwork for future developments toward whole abdomen free-breathing non-contrast perfusion imaging.

Estimation of wall shear stress in bypass grafts with computational fluid dynamics method

Coronary artery bypass graft (CABG) operation for coronary artery disease with different types of grafts has a large clinical application world wide. Immediately after this operation patients are usually relieved of their chest pain and have improved cardiac function. However, after a while, these bypass grafts may fail due to for example, neointimal hyperplasia or thrombosis. One of the causes for this bypass graft failure is assumed to be the blood flow with low wall shear stress. The aim of this research is to estimate the wall shear stress in a graft and thus to locate areas were wall shear stress is low. This was done with the help of a blood flow computer model. Postoperative biplane angiograms of the graft were recorded, and from these the three-dimensional geometry of the graft was reconstructed and imported into the computational fluid dynamics (CFD) program FLUENT. The stationary diastolic flow through the grafts was calculated, and the wall shear stress distribution was estimated. This procedure was carried out for one native vessel and two different types of bypass grafts. One bypass graft was a saphenous vein and the other one was a varicose saphenous vein encased in a fine, flexible metal mesh. The mesh was attached to give the graft a defined diameter. The computational results show that each graft has distinct areas of low wall shear stress. The graft with the metal mesh has an area of low wall shear stress (< 1 Pa, stationary flow), which is four times smaller than the respective areas in the other graft and in the native vessel. This is thought to be caused by the smaller and more uniform diameter of the metal mesh-reinforced graft.

Numerical simulation of haemodynamics and low-density lipoprotein transport in the rabbit aorta and their correlation with atherosclerotic plaque thickness

Two mechanisms of shear stress and mass transport have been recognized to play an important role in the development of localized atherosclerosis. However, their relationship and roles in atherogenesis are still obscure. It is necessary to investigate quantitatively the correlation among low-density lipoproteins (LDL) transport, haemodynamic parameters and plaque thickness. We simulated blood flow and LDL transport in rabbit aorta using computational fluid dynamics and evaluated plaque thickness in the aorta of a high-fat-diet rabbit. The numerical results show that regions with high luminal LDL concentration tend to have severely negative haemodynamic environments (HEs). However, for regions with moderately and slightly high luminal LDL concentration, the relationship between LDL concentration and the above haemodynamic indicators is not clear cut. Point-by-point correlation with experimental results indicates that severe atherosclerotic plaque corresponds to high LDL concentration and seriously negative HEs, less severe atherosclerotic plaque is related to either moderately high LDL concentration or moderately negative HEs, and there is almost no atherosclerotic plaque in regions with both low LDL concentration and positive HEs. In conclusion, LDL distribution is closely linked to blood flow transport, and the synergetic effects of luminal surface LDL concentration and wall shear stress-based haemodynamic indicators may determine plaque thickness.

Optimisation of a Novel Spiral-Inducing Bypass Graft Using Computational Fluid Dynamics

Graft failure is currently a major concern for medical practitioners in treating Peripheral Vascular Disease (PVD) and Coronary Artery Disease (CAD). It is now widely accepted that unfavourable haemodynamic conditions play an essential role in the formation and development of intimal hyperplasia, which is the main cause of graft failure. This paper uses Computational Fluid Dynamics (CFD) to conduct a parametric study to enhance the design and performance of a novel prosthetic graft, which utilises internal ridge(s) to induce spiral flow. This design is primarily based on the identification of the blood flow as spiral in the whole arterial system and is believed to improve the graft longevity and patency rates at distal graft anastomoses. Four different design parameters were assessed in this work and the trailing edge orientation of the ridge was identified as the most important parameter to induce physiological swirling flow, while the height of the ridge also significantly contributed to the enhanced performance of this type of graft. Building on these conclusions, an enhanced configuration of spiral graft is proposed and compared against conventional and spiral grafts to reaffirm its potential benefits.

Impact of Using Conventional Inlet/Outlet Boundary Conditions on Haemodynamic Metrics in a Subject-Specific Rabbit Aorta

Computational fluid dynamics is a tool capable of accurately measuring metrics currently used to predict the behaviour of cardiovascular diseases. This study quantifies the impact various commonly used inlet and outlet boundary conditions have on various shear rate–based haemodynamic metrics currently used for predicting the localisation of cardiovascular diseases. Simulations are conducted on an accurately represented rabbit aorta configuration and comparison has been made against available in vivo data. The boundary conditions studied include two different inlet profiles, three outlet boundary conditions, and steady-state versus pulsatile flow cases. Large variations were found in the results, particularly when using different outlet boundary conditions, and the discrepancies were evaluated both quantitatively and qualitatively. The results clearly highlight the importance of the type of boundary condition used when simulating complex cardiovascular models. By restricting the attention to the flow within the aorta and the intercostal branches, the results suggest that prescribing transient simulation and fully developed flow at the inlet are not required. Furthermore, assuming the widely accepted low wall shear stress theory of Caro, it was found that Murray’s law–based outlet boundary condition returns the most physiologically accurate results when compared to in vivo data.

Numerical Assessment of Novel Helical/Spiral Grafts with Improved Hemodynamics for Distal Graft Anastomoses

In the present work, numerical simulations were conducted for a typical end-to-side distal graft anastomosis to assess the effects of inducing secondary flow, which is believed to remove unfavourable flow environment. Simulations were carried out for four models, generated based on two main features of 'out-of-plane helicity' and 'spiral ridge' in the grafts as well as their combination. Following a qualitative comparison against in vitro data, various mean flow and hemodynamic parameters were compared and the results showed that helicity is significantly more effective in inducing swirling flow in comparison to a spiral ridge, while their combination could be even more effective. In addition, the induced swirling flow was generally found to be increasing the wall shear stress and reducing the flow stagnation and particle residence time within the anastomotic region and the host artery, which may be beneficial to the graft longevity and patency rates. Finally, a parametric study on the spiral ridge geometrical features was conducted, which showed that the ridge height and the number of spiral ridges have significant effects on inducing swirling flow, and revealed the potential of improving the efficiency of such designs.

Role of Pulse Pressure and Geometry of Primary Entry Tear in Acute Type B Dissection Propagation

The hemodynamic and geometric factors leading to propagation of acute Type B dissections are poorly understood. The objective is to elucidate whether geometric and hemodynamic parameters increase the predilection for aortic dissection propagation. A pulse duplicator set-up was used on porcine aorta with a single entry tear. Mean pressures of 100 and 180 mmHg were used, with pulse pressures ranging from 40 to 200 mmHg. The propagation for varying geometric conditions (%circumference of the entry tear: 15–65%, axial length: 0.5–3.2 cm) were tested for two flap thicknesses (1/3rd and 2/3rd of the thickness of vessel wall, respectively). To assess the effect of pulse and mean pressure on flap dynamics, the %true lumen (TL) cross-sectional area of the entry tear were compared. The % circumference for propagation of thin flap (47 ± 1%) was not significantly different (p = 0.14) from thick flap (44 ± 2%). On the contrary, the axial length of propagation for thin flap (2.57 ± 0.15 cm) was significantly different (p < 0.05) from the thick flap (1.56 ± 0.10 cm). TL compression was observed during systolic phase. For a fixed geometry of entry tear (%circumference = 39 ± 2%; axial length = 1.43 ± 0.13 cm), mean pressure did not have significant (p = 0.84) effect on flap movement. Increase in pulse pressure resulted in a significant change (p = 0.02) in %TL area (52 ± 4%). The energy acting on the false lumen immediately before propagation was calculated as 75 ± 9 J/m² and was fairly uniform across different specimens. Pulse pressure had a significant effect on the flap movement in contrast to mean pressure. Hence, mitigation of pulse pressure and restriction of flap movement may be beneficial in patients with type B acute dissections.

Abdominal Aortic Hemodynamics in Intermittent Claudication Patients at Rest and during Dynamic Pedaling Exercise

BACKGROUND

Lower-extremity exercise has been shown to eliminate adverse hemodynamics conditions, such as low and oscillating blood flow and wall shear stress, in the abdominal aortas of healthy young and older adults.

METHODS

We use cine phase-contrast magnetic resonance imaging and a custom MRI-compatible exercise cycle to quantify hemodynamic changes because of pedaling exercise in patients diagnosed with intermittent claudication.

RESULTS

With only an average heart increase of 35 ± 18% and exercise workload of 36 ± 16 watts, the patients experienced approximately 3- and 6-fold increases in blood flow, and 4- and 16-fold increases in wall shear stress at the supraceliac and infrarenal aortic locations, respectively. Also, all oscillations in flow and shear stress at rest were eliminated with exercise.

CONCLUSIONS

Claudication patients experience 3- to 4-fold lower oscillations in flow and shear stress at rest as compared with healthy age-matched controls, likely because of reduced distal arterial compliance as a result of distal atherosclerosis. The magnitude of flow and shear oscillatory indices may be good indicators of distal arterial compliance and health, and may provide predictive power for the efficacy of focal interventions.

A comparative review of the hemodynamics and pathogenesis of cerebral and abdominal aortic aneurysms: lessons to learn from each other

OBJECTIVE

Cerebral aneurysms (CAs) and abdominal aortic aneurysms (AAAs) are degenerative vascular pathologies that manifest as abnormal dilations of the arterial wall. They arise with different morphologies in different types of blood vessels under different hemodynamic conditions. Although treated as different pathologies, we examine common pathways in their hemodynamic pathogenesis in order to elucidate mechanisms of formation.

MATERIALS AND METHODS

A systematic review of the literature was performed. Current concepts on pathogenesis and hemodynamics were collected and compared.

RESULTS

CAs arise as saccular dilations on the cerebral arteries of the circle of Willis under high blood flow, high wall shear stress (WSS), and high wall shear stress gradient (WSSG) conditions. AAAs arise as fusiform dilations on the infrarenal aorta under low blood flow, low, oscillating WSS, and high WSSG conditions. While at opposite ends of the WSS spectrum, they share high WSSG, a critical factor in arterial remodeling. This alone may not be enough to initiate aneurysm formation, but may ignite a cascade of downstream events that leads to aneurysm development. Despite differences in morphology and the structure, CAs and AAAs share many histopathological and biomechanical characteristics. Endothelial cell damage, loss of elastin, and smooth muscle cell loss are universal findings in CAs and AAAs. Increased matrix metalloproteinases and other proteinases, reactive oxygen species, and inflammation also contribute to the pathogenesis of both aneurysms.

CONCLUSION

Our review revealed similar pathways in seemingly different pathologies. We also highlight the need for cross-disciplinary studies to aid in finding similarities between pathologies.

Pro-atherogenic shear stress and HIV proteins synergistically upregulate cathepsin K in endothelial cells

Major advances in highly active antiretroviral therapies (HAART) have extended the lives of people living with HIV, but there still remains an increased risk of death by cardiovascular diseases (CVD). HIV proteins have been shown to contribute to cardiovascular dysfunction with effects on the different cell types that comprise the arterial wall. In particular, HIV-1 transactivating factor (Tat) has been shown to bind to endothelial cells inducing a range of responses that contribute to vascular dysfunction. It is well established that hemodynamics also play an important role in endothelial cell mediated atherosclerotic development. When exposed to low or oscillatory shear stress, such as that found at branches and bifurcations, endothelial cells contribute to proteolytic vascular remodeling by upregulating cathepsins, potent elastases and collagenases that contribute to altered biomechanics and plaque formation. Mechanisms to understand the influence of Tat on shear stress mediated vascular remodeling have not been fully elucidated. Using an in vivo HIV-Tg mouse model and an in vitro cone and plate shear stress bioreactor to actuate physiologically relevant pro-atherogenic or atheroprotective shear stress on human aortic endothelial cells, we have shown synergism between HIV proteins and pro-atherogenic shear stress to increase endothelial cell expression of the powerful protease cathepsin K, and may implicate this protease in accelerated CVD in people living with HIV.

To establish whether McKenzie lumbar flexion and extension mobility exercises performed in lying affect central as well as systemic hemodynamics: a crossover experimental study

OBJECTIVE

Examine systemic and central hemodynamic responses following McKenzie lumbar flexion and extension mobility exercises performed in lying (FIL and EIL).

DESIGN

Crossover experimental study.

SETTING

Clinical laboratory.

PARTICIPANTS

Healthy male volunteers (n=25) (mean(SD) age: 28(3)years; range 21 to 34).

INTERVENTIONS

Based on alternating assignment of either FIL or EIL to participants, three sets of the first exercise (10, 15, 20 repetitions) were performed with 5-minute rest between sets; after 15-minute rest, the protocol was repeated for the other exercise.

MAIN OUTCOME MEASURES

Systemic hemodynamic parameters included heart rate (HR), and systolic and diastolic blood pressures (SBP, DBP). Central hemodynamic parameters included abdominal aortic diameter (AD), peak systolic velocity (PSV/AD), end diastolic velocity (EDV/AD) and resistive index (RI). Measures recorded after each exercise set.

RESULTS

FIL RPP at baseline was 9.1 (1.4), after 20 repetitions 18.3 (2.5), mean difference 8.9 (95% confidence interval (CI) 7.9 to 9.8) compared to EIL at baseline 9.1 (1.5), after 20 repetitions 13.0 (3.1), mean difference 4.1 (95% CI 3.3 to 5.0). FIL RI at baseline was 0.78 (0.03), after 20 repetitions 0.87 (0.03), mean difference 0.08 (95% CI 0.06 to 0.10) compared to EIL at baseline 0.78 (0.03), after 20 repetitions 0.83 (0.03), mean difference 0.05 (95% CI 0.04 to 0.07).

CONCLUSIONS

Although 10 repetitions of FIL and EIL may be regarded as safe, our findings support screening patients with lifestyle risk factors, and cautioning about adhering to recommended repetition number given associated increased work of the heart. The extent of AD mechanical perturbation remains unclear.

Involvement of the renin-angiotensin system in abdominal and thoracic aortic aneurysms

Aortic aneurysms are relatively common maladies that may lead to the devastating consequence of aortic rupture. AAAs (abdominal aortic aneurysms) and TAAs (thoracic aortic aneurysms) are two common forms of aneurysmal diseases in humans that appear to have distinct pathologies and mechanisms. Despite this divergence, there are numerous and consistent demonstrations that overactivation of the RAS (renin-angiotensin system) promotes both AAAs and TAAs in animal models. For example, in mice, both AAAs and TAAs are formed during infusion of AngII (angiotensin II), the major bioactive peptide in the RAS. There are many proposed mechanisms by which the RAS initiates and perpetuates aortic aneurysms, including effects of AngII on a diverse array of cell types and mediators. These experimental findings are complemented in humans by genetic association studies and retrospective analyses of clinical data that generally support a role of the RAS in both AAAs and TAAs. Given the lack of a validated pharmacological therapy for any form of aortic aneurysm, there is a pressing need to determine whether the consistent findings on the role of the RAS in animal models are translatable to humans afflicted with these diseases. The present review compiles the recent literature that has shown the RAS as a critical component in the pathogenesis of aortic aneurysms.

Correlation between local hemodynamics and lesion distribution in a novel aortic regurgitation murine model of atherosclerosis

Following surgical induction of aortic valve regurgitation (AR), extensive atherosclerotic plaque development along the descending thoracic and abdominal aorta of Ldlr⁻/⁻ mice has been reported, with distinct spatial distributions suggestive of a strong local hemodynamic influence. The objective of this study was to test, using image-based computational fluid dynamics (CFD), whether this is indeed the case. The lumen geometry was reconstructed from micro-CT scanning of a control Ldlr⁻/⁻ mouse, and CFD simulations were carried out for both AR and control flow conditions derived from Doppler ultrasound measurements and literature data. Maps of time-averaged wall shear stress magnitude (TAWSS), oscillatory shear index (OSI) and relative residence time (RRT) were compared against the spatial distributions of plaque stained with oil red O, previously acquired in a group of AR and control mice. Maps of OSI and RRT were found to be consistent with plaque distributions in the AR mice and the absence of plaque in the control mice. TAWSS was uniformly lower under control vs. AR flow conditions, suggesting that levels (> 100 dyn/cm²) exceeded those required to alone induce a pro-atherogenic response. Simulations of a straightened CFD model confirmed the importance of anatomical curvature for explaining the spatial distribution of lesions in the AR mice. In summary, oscillatory and retrograde flow induced in the AR mice, without concomitant low shear, may exacerbate or accelerate lesion formation, but the distinct anatomical curvature of the mouse aorta is responsible for the spatial distribution of lesions.

Role of hemodynamic shear stress in cardiovascular disease

Atherosclerosis is the main cause of morbidity and mortality in the Western world. Inflammation and blood flow alterations are new markers emerging as possible determinants for the development of atherosclerotic lesions. In particular, blood flow exerts a shear stress on vessel walls that alters cell physiology. Shear stress arises from the friction between two virtual layers of a fluid and is induced by the difference in motion and viscosity between these layers. Regions of the arterial tree with uniform geometry are exposed to a unidirectional and constant flow, which determines a physiologic shear stress, while arches and bifurcations are exposed to an oscillatory and disturbed flow, which determines a low shear stress. Atherosclerotic lesions develop mainly in areas of low shear stress, while those exposed to a physiologic shear stress are protected. The presence of areas of the arterial tree with different wall shear stress may explain, in part, the different localization of atherosclerotic lesions in both coronary and extracoronary arteries. The measurement of this parameter may help in identifying atherosclerotic plaques at higher risk as well as in evaluating the efficacy of different pharmacological interventions. Moreover, an altered shear stress is associated with the occurrence of both aortic and intracranial aneurysms, possibly leading to their growth and rupture. Finally, the evaluation of shear stress may be useful for predicting the risk of developing restenosis after coronary and peripheral angioplasty and for devising a coronary stent with a strut design less thrombogenic and more conducive to endothelization.

Patient-based abdominal aortic aneurysm rupture risk prediction with fluid structure interaction modeling

Elective repair of abdominal aortic aneurysm (AAA) is warranted when the risk of rupture exceeds that of surgery, and is mostly based on the AAA size as a crude rupture predictor. A methodology based on biomechanical considerations for a reliable patient-specific prediction of AAA risk of rupture is presented. Fluid-structure interaction (FSI) simulations conducted in models reconstructed from CT scans of patients who had contained ruptured AAA (rAAA) predicted the rupture location based on mapping of the stresses developing within the aneurysmal wall, additionally showing that a smaller rAAA presented a higher rupture risk. By providing refined means to estimate the risk of rupture, the methodology may have a major impact on diagnostics and treatment of AAA patients.

Aortic regurgitation dramatically alters the distribution of atherosclerotic lesions and enhances atherogenesis in mice

OBJECTIVE

Hemodynamics plays a critical role in atherogenesis, but the association between flow pattern and preferential localization of lesion is not fully understood. We developed a mouse model of aortic valve regurgitation (AR) to change the aortic flow pattern and observed the effects on plaque formation.

METHODS AND RESULTS

High-frequency Doppler ultrasound imaging of 10 untreated C57BL/6J mice and 6 sham-treated low-density lipoprotein receptor-deficient (Ldlr(-/-)) mice revealed consistent antegrade blood flow throughout the aorta and oscillatory flow only along the lesser curvature of the aortic arch. Catheter-induced AR in 7 Ldlr(-/-) mice produced various degrees of diastolic retrograde flow throughout the aorta. After the mice were fed a cholesterol-enriched diet for 6 weeks, the burden of atherosclerotic lesions was increased 6-fold, with the naturally plaque-resistant descending aorta becoming susceptible. The AR severity correlated positively with the lesion burden in the descending thoracic and abdominal aorta but negatively with the lesions in the ascending aorta and aortic arch.

CONCLUSIONS

This AR model is valuable for elucidating the relationship between hemodynamics and predisposition of the artery wall to atherosclerosis, because of the significant alterations in local flow patterns and the conversion of large regions in the descending aorta from lesion resistant to lesion prone.

In vivo assessment of blood flow patterns in abdominal aorta of mice with MRI: implications for AAA localization

Abdominal aortic aneurysms (AAA) localize in the infrarenal aorta in humans, while they are found in the suprarenal aorta in mouse models. It has been shown previously that humans experience a reversal of flow during early diastole in the infrarenal aorta during each cardiac cycle. This flow reversal causes oscillatory wall shear stress (OWSS) to be present in the infrarenal aorta of humans. OWSS has been linked to a variety of proatherogenic and proinflammatory factors. The presence of reverse flow in the mouse aorta is unknown. In this study we investigated blood flow in mice, using phase-contrast magnetic resonance (PCMR) imaging. We measured blood flow in the suprarenal and infrarenal abdominal aorta of 18 wild-type C57BL/6J mice and 15 apolipoprotein E (apoE)-/- mice. Although OWSS was not directly evaluated, results indicate that, unlike humans, there is no reversal of flow in the infrarenal aorta of wild-type or apoE-/- mice. Distensibility of the mouse aortic wall in both the suprarenal and infrarenal segments is higher than reported values for the human aorta. We conclude that normal mice do not experience the reverse flow in the infrarenal aorta that is observed in humans.

A numerical study on the flow of blood and the transport of LDL in the human aorta: the physiological significance of the helical flow in the aortic arch

It has been proposed that a mass transfer phenomenon called concentration polarization of low-density lipoproteins (LDLs) may occur in the arterial system and is likely involved in the localization of atherogenesis. To test the hypothesis that concentration polarization of LDL may be suppressed by the helical flow pattern in the human aorta, hence sparing the ascending aorta from atherosclerosis, the effects of aortic torsion, branching, curvature, and taper on blood flow and LDL transport in the lumen were simulated numerically under steady-state flow conditions using four aorta models constructed based on in vivo MRI slices. The results showed that it was the aortic torsion that induced the helical flow in the aortic arch, stabilizing the flow of blood in the aorta, and compensated the adverse effects of the aortic curvature on blood flow and LDL transport. The helical flow reduced the luminal surface LDL concentration in the aortic arch and probably played a role in suppressing severe polarization of LDL at the entrances of the three branches on the arch, hence, protecting them from atherogenesis. The taper of the aorta was another important feature of the aorta that further stabilized the flow of blood and delayed the attenuation of the helical flow, making it move beyond the arch and into the beginning part of the descending aorta. The results therefore may account for why the ascending aorta and the arch are relatively free of atherosclerosis.

CFD and PTV steady flow investigation in an anatomically accurate abdominal aortic aneurysm

There is considerable interest in computational and experimental flow investigations within abdominal aortic aneurysms (AAAs). This task stipulates advanced grid generation techniques and cross-validation because of the anatomical complexity. The purpose of this study is to examine the feasibility of velocity measurements by particle tracking velocimetry (PTV) in realistic AAA models. Computed tomography and rapid prototyping were combined to digitize and construct a silicone replica of a patient-specific AAA. Three-dimensional velocity measurements were acquired using PTV under steady averaged resting boundary conditions. Computational fluid dynamics (CFD) simulations were subsequently carried out with identical boundary conditions. The computational grid was created by splitting the luminal volume into manifold and nonmanifold subsections. They were filled with tetrahedral and hexahedral elements, respectively. Grid independency was tested on three successively refined meshes. Velocity differences of about 1% in all three directions existed mainly within the AAA sack. Pressure revealed similar variations, with the sparser mesh predicting larger values. PTV velocity measurements were taken along the abdominal aorta and showed good agreement with the numerical data. The results within the aneurysm neck and sack showed average velocity variations of about 5% of the mean inlet velocity. The corresponding average differences increased for all velocity components downstream the iliac bifurcation to as much as 15%. The two domains differed slightly due to flow-induced forces acting on the silicone model. Velocity quantification through narrow branches was problematic due to decreased signal to noise ratio at the larger local velocities. Computational wall pressure and shear fields are also presented. The agreement between CFD simulations and the PTV experimental data was confirmed by three-dimensional velocity comparisons at several locations within the investigated AAA anatomy indicating the feasibility of this approach.

Patient-specific modeling of cardiovascular mechanics

Advances in numerical methods and three-dimensional imaging techniques have enabled the quantification of cardiovascular mechanics in subject-specific anatomic and physiologic models. Patient-specific models are being used to guide cell culture and animal experiments and test hypotheses related to the role of biomechanical factors in vascular diseases. Furthermore, biomechanical models based on noninvasive medical imaging could provide invaluable data on the in vivo service environment where cardiovascular devices are employed and on the effect of the devices on physiologic function. Finally, patient-specific modeling has enabled an entirely new application of cardiovascular mechanics, namely predicting outcomes of alternate therapeutic interventions for individual patients. We review methods to create anatomic and physiologic models, obtain properties, assign boundary conditions, and solve the equations governing blood flow and vessel wall dynamics. Applications of patient-specific models of cardiovascular mechanics are presented, followed by a discussion of the challenges and opportunities that lie ahead.

Numerical implementation of viscoelastic blood flow in a simplified arterial geometry

The influence of the non-Newtonian stress-strain relation of blood on the oscillatory shear index (OSI) and mean wall shear stress (WSS) are described. The unsteady linear 1D momentum equation is solved for a viscoelastic fluid, with six elements of the Maxwell type and one dashpot element connected in parallel. A novel numerical algorithm is described which uses the upwind finite difference method to solve the equation of momentum. Results obtained by using a finite difference approach show significantly higher values of OSI when blood is assumed to be a viscoelastic fluid compared with those of simplified Newtonian fluid model. The calculation of OSI in human normal conditions for the Newtonian fluid differs in 12% (if alpha=0.02) from the results obtained from using the viscoelastic model.

Mesoscopic simulations of systolic flow in the human abdominal aorta

The complex nature of blood flow in the human arterial system is still gaining more attention, as it has become clear that cardiovascular diseases localize in regions of complex geometry and complex flow fields. In this article, we demonstrate that the lattice Boltzmann method can serve as a mesoscopic computational hemodynamic solver. We argue that it may have benefits over the traditional Navier-Stokes techniques. The accuracy of the method is tested by studying time-dependent systolic flow in a 3D straight rigid tube at typical hemodynamic Reynolds and Womersley numbers as an unsteady flow benchmark. Simulation results of steady and unsteady flow in a model of the human aortic bifurcation reconstructed from magnetic resonance angiography, are presented as a typical hemodynamic application.

Fluid-structure interaction in abdominal aortic aneurysms: effects of asymmetry and wall thickness

BACKGROUND

Abdominal aortic aneurysm (AAA) is a prevalent disease which is of significant concern because of the morbidity associated with the continuing expansion of the abdominal aorta and its ultimate rupture. The transient interaction between blood flow and the wall contributes to wall stress which, if it exceeds the failure strength of the dilated arterial wall, will lead to aneurysm rupture. Utilizing a computational approach, the biomechanical environment of virtual AAAs can be evaluated to study the affects of asymmetry and wall thickness on this stress, two parameters that contribute to increased risk of aneurysm rupture.

METHODS

Ten virtual aneurysm models were created with five different asymmetry parameters ranging from beta = 0.2 to 1.0 and either a uniform or variable wall thickness to study the flow and wall dynamics by means of fully coupled fluid-structure interaction (FSI) analyses. The AAA wall was designed to have a (i) uniform 1.5 mm thickness or (ii) variable thickness ranging from 0.5-1.5 mm extruded normally from the boundary surface of the lumen. These models were meshed with linear hexahedral elements, imported into a commercial finite element code and analyzed under transient flow conditions. The method proposed was then compared with traditional computational solid stress techniques on the basis of peak wall stress predictions and cost of computational effort.

RESULTS

The results provide quantitative predictions of flow patterns and wall mechanics as well as the effects of aneurysm asymmetry and wall thickness heterogeneity on the estimation of peak wall stress. These parameters affect the magnitude and distribution of Von Mises stresses; varying wall thickness increases the maximum Von Mises stress by 4 times its uniform thickness counterpart. A pre-peak systole retrograde flow was observed in the AAA sac for all models, which is due to the elastic energy stored in the compliant arterial wall and the expansion force of the artery during systole.

CONCLUSION

Both wall thickness and geometry asymmetry affect the stress exhibited by a virtual AAA. Our results suggest that an asymmetric AAA with regional variations in wall thickness would be exposed to higher mechanical stresses and an increased risk of rupture than a more fusiform AAA with uniform wall thickness. Therefore, it is important to accurately reproduce vessel geometry and wall thickness in computational predictions of AAA biomechanics.

Simultaneous measurement of blood and myocardial velocity in the rat heart by phase contrast MRI using sparseq-space sampling

PURPOSE

To measure cardiac blood flow patterns and ventricular wall velocities through the cardiac cycle in anesthetized Wistar Kyoto (WKY) rats.

MATERIALS AND METHODS

A gradient-echo cine pulse sequence incorporating pulsed field gradients (PFGs) provided phase contrast (PC) motion encoding. We achieved a range of velocity sensitivity that was sufficient to measure simultaneously the large flow velocities within the cardiac chambers and aortic outflow tract (up to 70 cm s(-1) during systole), and the comparatively small velocities of the cardiac wall (0-3 cm s(-1)). A scheme of sparsely sampling q-space combined with a probability-based method of velocity calculation permitted such measurements along three orthogonal axes, and yielded velocity vector maps in all four chambers of the heart and the aorta, in both longitudinal and transverse sections, for up to 12 time-points in the cardiac cycle.

RESULTS

Left ventricular systole was associated with a symmetrical laminar flow pattern along the cardiac axis, with no appearance of turbulence. In contrast, blood showed a swirling motion within the right ventricle (RV) in the region of the pulmonary outflow tract. During left ventricular diastole a plume of blood entered the left ventricle (LV) from the left atrium. The ventricular flow patterns could also be correlated with measurements of left ventricular wall motion. The greatest velocities of the ventricular walls occurred in the transverse cardiac plane and were maximal during diastolic refilling. The cardiac wall motion in the longitudinal axis demonstrated a caudal-apical movement that may also contribute to diastolic refilling.

CONCLUSION

The successful measurements of blood and myocardial velocity during normal myocardial function may be extended to quantify pathological cardiac changes in animal models of human cardiac disease.

The relationship between wall shear stress distributions and intimal thickening in the human abdominal aorta

PURPOSE

The goal of this work was to determine wall shear stress (WSS) patterns in the human abdominal aorta and to compare these patterns to measurements of intimal thickness (IT) from autopsy samples.

METHODS

The WSS was experimentally measured using the laser photochromic dye tracer technique in an anatomically faithful in vitro model based on CT scans of the abdominal aorta in a healthy 35-year-old subject. IT was quantified as a function of circumferential and axial position using light microscopy in ten human autopsy specimens.

RESULTS

The histomorphometric analysis suggests that IT increases with age and that the distribution of intimal thickening changes with age. The lowest WSS in the flow model was found on the posterior wall inferior to the inferior mesenteric artery, and coincided with the region of most prominent IT in the autopsy samples. Local geometrical features in the flow model, such as the expansion at the inferior mesenteric artery (common in younger individuals), strongly influenced WSS patterns. The WSS was found to correlate negatively with IT (r2 = 0.3099; P = 0.0047).

CONCLUSION

Low WSS in the abdominal aorta is co-localized with IT and may be related to atherogenesis. Also, rates of IT in the abdominal aorta are possibly influenced by age-related geometrical changes.

Cyclic flow characteristics in an idealized asymmetric abdominal aortic aneurysm model

The flow characteristics and the corresponding hydrodynamic stability in an idealized asymmetric abdominal aortic aneurysm (AAA) model have been investigated using a laser Doppler anemometer. A rectified sine waveform was used to simulate aortic flow conditions (Re(delta) = 806 and alpha = 12.2). The flow around the distal neck of the AAA model undergoes transition and becomes turbulent for a fraction of time shortly after the commencement of the deceleration phases at every flow cycle while the rest of the flow inside the model stayed laminar throughout the cycle. As a result of non-symmetric vortical structure development inside the model, the distribution of turbulent shear stresses was found to be highly uneven along the radial direction of the model; this is in contrast to results found by the present authors in the symmetrical AAA model. The maximum turbulent shear stress found at the straight side of the distal neck are four times more than the maximum turbulent shear stress measured at the most dilated side of the distal neck. One of the interesting biological implications of the results is that the outward dilation of the arterial wall may be a physiological response to avoid the high turbulent shear load from the momentarily turbulent blood flow.

Flow changes in the aorta associated with the deployment of a AAA stent graft

The purpose of this study was to investigate the hemodynamic implications of a proximal shift in the aortic bifurcation that results from abdominal aortic aneurysm (AAA) stent graft deployment. A flow model was constructed in which an anatomically accurate model of the aorta was subjected to physiologic pulsatile flow. The model included the celiac, superior mesenteric, left and right renal arteries. The aortic bifurcation, leading to the right and left iliac arteries was included, as well as the lumbar curvature. Flow simulations were performed under resting and mild exercise conditions with and without a Cordis AAA stent graft deployed. Flow patterns were visualized with dye injection and recorded onto video. The flow rates through the iliac and renal arteries were continuously monitored using ultrasonic flowmeters. Flow visualization revealed that flow disturbances at the level of the renal arteries were slightly increased with the deployment of the stent graft. The orientation of the endolegs within the aorta had no perceptible effect on these disturbances. Under mild exercise conditions, very little flow disturbance was observed. In conclusion, there are slight changes in flow disturbance near the renal arteries due to stent graft deployment, but these changes would not be expected to have significant clinical implications.

The influence of out-of-plane geometry on pulsatile flow within a distal end-to-side anastomosis

We present an experimental and computational investigation of time-varying flow in an idealized fully occluded 45 degrees distal end-to-side anastomosis. Two geometric configurations are assessed, one where the centerlines of host and bypass vessels lie within a plane, and one where the bypass vessel is deformed out of the plane of symmetry, respectively, termed planar and non-planar. Flow experiments were conducted by magnetic resonance imaging in rigid wall models and computations were performed using a high order spectral/hp algorithm. Results indicate a significant change in the spatial distribution of wall shear stress and a reduction of the time-averaged peak wall shear stress magnitude by 10% in the non-planar model as compared to the planar configuration. In the planar geometry the stagnation point follows a straight-line path along the host artery bed with a path length of 0.8 diameters. By contrast in the non-planar case the stagnation point oscillates about a center that is located off the symmetry plane intersection with the host artery bed wall, and follows a parabolic path with a 0.7 diameter longitudinal and 0.5 diameter transverse excursion. A definition of the oscillatory shear index (OSI) is introduced that varies between 0 and 0.5 and that accounts for a continuous range of wall shear stress vector angles. In both models, regions of elevated oscillatory shear were spatially associated with regions of separated or oscillating stagnation point flow. The mean oscillatory shear magnitude (considering sites where OSI>0.1) in the non-planar geometry was reduced by 22% as compared to the planar configuration. These changes in the dynamic behavior of the stagnation point and the oscillatory shear distribution introduced by out-of-plane graft curvature may influence the localization of vessel wall sites exposed to physiologically unfavorable flow conditions.

Fluid mechanics of vascular systems, diseases, and thrombosis

The cardiovascular system is an internal flow loop with multiple branches circulating a complex liquid. The hallmarks of blood flow in arteries are pulsatility and branches, which cause wall stresses to be cyclical and nonuniform. Normal arterial flow is laminar, with secondary flows generated at curves and branches. Arteries can adapt to and modify hemodynamic conditions, and unusual hemodynamic conditions may cause an abnormal biological response. Velocity profile skewing can create pockets in which the wall shear stress is low and oscillates in direction. Atherosclerosis tends to localize to these sites and creates a narrowing of the artery lumen--a stenosis. Plaque rupture or endothelial injury can stimulate thrombosis, which can block blood flow to heart or brain tissues, causing a heart attack or stroke. This small lumen and elevated shear rate in a stenosis create conditions that accelerate platelet accumulation and occlusion. The relationship between thrombosis and fluid mechanics is complex, especially in the post-stenotic flow field. New convection models have been developed to predict clinical from platelet thrombosis in diseased arteries. Future hemodynamic studies should address the complex mechanics of flow-induced, large-scale wall motion and convection of semisolid particles and cells in flowing blood.

The influence of out-of-plane geometry on the flow within a distal end-to-side anastomosis

This paper describes a computational and experimental investigation of flow in a proto-type model geometry of a fully occluded 45 deg distal end-to-side anastomosis. Previous investigations have considered a similar configuration where the centerlines of the bypass and host vessels lie within a plane, thereby producing a plane of symmetry within the flow. We have extended these investigations by deforming the bypass vessel out of the plane of symmetry, thereby breaking the symmetry of the flow and producing a nonplanar geometry. Experimental data were obtained using magnetic resonance imaging of flow within perspex models and computational data were obtained from simulations using a high-order spectral/hp element method. We found that the nonplanar three-dimensional flow notably alters the distribution of wall shear stress at the bed of the anastomosis, reducing the peak wall shear stress peak by approximately 10 percent when compared with the planar model. Furthermore, an increase in the absolute flux of velocity into the occluded region, proximal to the anastomosis, of 80 percent was observed in the nonplanar geometry when compared with the planar geometry.

Development of velocity profiles and retrograde flow in the porcine abdominal aorta under different haemodynamic conditions

Low and/or oscillating wall shear stresses are related to the development of atherosclerosis and this oscillation is influenced by changes in basic haemodynamics (exercise). The objective of this study was to provide in vivo data on the development of velocity profiles and oscillating blood velocities in the abdominal aorta under varying haemodynamic conditions. Six anaesthetized, 90-kg pigs were used in the study. Abdominal aortic velocity profiles across the anterior-posterior diameter were acquired at different axial positions using 10 MHz pulsed Doppler ultrasound. Measurements were obtained under normal conditions and during cardiac pacing up to 170 beats/min. Velocity profiles were obtained during heart rates ranging between 58 and 169 beats/min, and during flow rates ranging between 0.57 and 2.89 l/min. Main outcome measures included minimum velocities, frequency index, shape of velocity profiles (velocity distribution index), Reynolds' numbers, and Womersley's frequency parameter. Velocity profiles were blunted, with lowest velocities at the distal posterior vessel wall. Multiple regression analysis showed the development of velocity profiles to be inversely correlated with the pulsatility index, Womersley's frequency parameter and the mean Reynolds' number (r = 0.89, p < 0.0005). Minimum velocities were negatively correlated with the PI, Womersley's frequency parameter and positively with the mean Reynolds' number (r = 0.94, p < 10(-8)). Retrograde velocities (and hence oscillating wall shear stresses) were present at mean Reynolds' number < 1000. The oscillation of blood velocities at the wall in the porcine abdominal aorta was highly dependent on general haemodynamics (i.e. flow, heart rate and vessel diameter as expressed in the Reynolds' numbers and Womersley's frequency parameters). The velocity profiles in the abdominal aorta were found to be far from parabolic. These findings have important implications for the understanding and future modelling of the complex haemodynamics in the abdominal aorta and their relation to the development of atherosclerotic disease.

Effect of exercise on hemodynamic conditions in the abdominal aorta

PURPOSE

The beneficial effect of exercise in the retardation of the progression of cardiovascular disease is hypothesized to be caused, at least in part, by the elimination of adverse hemodynamic conditions, including flow recirculation and low wall shear stress. In vitro and in vivo investigations have provided qualitative and limited quantitative information on flow patterns in the abdominal aorta and on the effect of exercise on the elimination of adverse hemodynamic conditions. We used computational fluid mechanics methods to examine the effects of simulated exercise on hemodynamic conditions in an idealized model of the human abdominal aorta.

METHODS

A three-dimensional computer model of a healthy human abdominal aorta was created to simulate pulsatile aortic blood flow under conditions of rest and graded exercise. Flow velocity patterns and wall shear stress were computed in the lesion-prone infrarenal aorta, and the effects of exercise were determined.

RESULTS

A recirculation zone was observed to form along the posterior wall of the aorta immediately distal to the renal vessels under resting conditions. Low time-averaged wall shear stress was present in this location, along the posterior wall opposite the superior mesenteric artery and along the anterior wall between the superior and inferior mesenteric arteries. Shear stress temporal oscillations, as measured with an oscillatory shear index, were elevated in these regions. Under simulated light exercise conditions, a region of low wall shear stress and high oscillatory shear index remained along the posterior wall immediately distal to the renal arteries. Under simulated moderate exercise conditions, all the regions of low wall shear stress and high oscillatory shear index were eliminated.

CONCLUSION

This numeric investigation provided detailed quantitative data on the effect of exercise on hemodynamic conditions in the abdominal aorta. Our results indicated that moderate levels of lower limb exercise are necessary to eliminate the flow reversal and regions of low wall shear stress in the abdominal aorta that exist under resting conditions. The lack of flow reversal and increased wall shear stress during exercise suggest a mechanism by which exercise may promote arterial health, namely with the elimination of adverse hemodynamic conditions.