Numerical analysis of pulsatile blood flow and vessel wall mechanics in different degrees of stenoses

Effect of plaque micro-watershed changes on carotid atherosclerosis

OBJECTIVE

The study aims to elucidate the mechanisms underlying plaque growth by analyzing the variations in hemodynamic parameters within the plaque region of patients' carotid arteries before and after the development of atherosclerotic lesions.

METHODS

The study enrolls 25 patients with common carotid artery stenosis and 25 with tandem carotid artery stenosis. Based on pathological analysis, three-dimensional models of the actual blood vessels before and after the lesion are constructed for two patients within a two-year period. Computational fluid dynamics is employed to conduct unsteady periodic non-Newtonian fluid numerical simulations, enabling an in-depth investigation into the changes in the micro-environment of blood flow.

RESULTS

During the systolic phase of the cardiac cycle, vortex regions are particularly prone to developing at the bifurcation point between the common carotid artery and the distal end of the internal carotid artery. In the early diastolic phase, blood reflux phenomena can be observed within the carotid artery. Towards the end of diastole, there is an expansion of vortex regions at the bifurcation point of the carotid artery. The shoulder region of initial small plaques within the blood vessel is susceptible to developing a low-speed recirculation zone, characterized by significantly reduced shear stress compared to the surrounding areas. Following vascular stenosis, the wall shear stress within the plaque domain generally increases; however, it maintains a consistent pattern of high central values and low upper shoulder values. The shear stress at the upper shoulder of the plaque of tandem carotid stenosis is below 0.4 Pa, whereas the central and lower shoulder regions exhibit shear stress exceeding 40 Pa.

CONCLUSIONS

The dynamic parameters of the blood flow micro-environment exhibit variations throughout the cardiac cycle, and temporal disparities exist in local lesions within the carotid artery. Both common and tandem carotid artery stenosis are particularly prone to developing lesions at the shoulder of initial small plaques. The micro-flow characteristics within the plaque domain undergo alterations prior to and following the onset of carotid artery disease. Furthermore, the occurrence of restenosis and rupture is associated with the location of plaque growth.

A CFD-FFT approach to hemoacoustics that enables degree of stenosis prediction from stethoscopic signals

In this paper, we identify a new (acoustic) frequency-stenosis relation whose frequencies lie within the recommended auscultation threshold of stethoscopy (< 120 Hz). We show that this relation can be used to extend the application of phonoangiography (quantifying the degree of stenosis from bruits) to widely accessible stethoscopes. The relation is successfully identified from an analysis restricted to the acoustic signature of the von Karman vortex street, which we automatically single out by means of a metric we propose that is based on an area-weighted average of the Q-criterion for the post-stenotic region. Specifically, we perform CFD simulations on internal flow geometries that represent stenotic blood vessels of different severities. We then extract their emitted acoustic signals using the Ffowcs Williams-Hawkings equation, which we subtract from a clean signal (stenosis free) at the same heart rate. Next, we transform this differential signal to the frequency domain and carefully classify its acoustic signatures per six (stenosis-)invariant flow phases of a cardiac cycle that are newly identified in this paper. We then automatically restrict our acoustic analysis to the sounds emitted by the von Karman vortex street (phase 4) by means of our Q-criterion-based metric. Our analysis of its acoustic signature reveals a strong linear relationship between the degree of stenosis and its dominant frequency, which differs considerably from the break frequency and the heart rate (known dominant frequencies in the literature). Applying our new relation to available stethoscopic data, we find that its predictions are consistent with clinical assessment. Our finding of this linear correlation is also unlike prevalent scaling laws in the literature, which feature a small exponent (i.e., low stenosis percentage sensitivity over much of the clinical range). They hence can only distinguish mild, moderate, and severe cases. Conversely, our linear law can identify variations in the degree of stenosis sensitively and accurately for the full clinical range, thus significantly improving the utility of the relevant scaling laws... Future research will investigate incorporating the vibroacoustic role of adjacent organs to expand the clinical applicability of our findings. Extending our approach to more complex 3D stenotic morphologies and including the vibroacoustic role of surrounding organs will be explored in future research to advance the clinical reach of our findings.

On the nonlinear relationship between wall shear stress topology and multi-directionality in coronary atherosclerosis

BACKGROUND AND OBJECTIVE

In this paper we investigate twelve multi-directional/topological wall shear stress (WSS) derived metrics and their relationships with the formation of coronary plaques in both computational fluid dynamics (CFD) and dynamic fluid-structure interaction (FSI) frameworks. While low WSS is one of the most established biomechanical markers associated with coronary atherosclerosis progression, alone it is limited. Multi-directional and topological WSS derived metrics have been shown to be important in atherosclerosis related mechanotransduction and near-wall transport processes. However, the relationships between these twelve WSS metrics and the influence of both FSI simulations and coronary dynamics is understudied.

METHODS

We first investigate the relationships between these twelve WSS derived metrics, stenosis percentage and lesion length through a parametric, transient CFD study. Secondly, we extend the parametric study to FSI, both with and without the addition of coronary dynamics, and assess their correlations. Finally, we present the case of a patient who underwent invasive coronary angiography and optical coherence tomography imaging at two time points 18 months apart. Associations between each of the twelve WSS derived metrics in CFD, static FSI and dynamic FSI simulations were assessed against areas of positive/negative vessel remodelling, and changes in plaque morphology.

RESULTS

22-32% stenosis was the threshold beyond which adverse multi-directional/topological WSS results. Each metric produced a different relationship with changing stenoses and lesion length. Transient haemodynamics was impacted by coronary dynamics, with the topological shear variation index suppressed by up to 94%. These changes appear more critical at smaller stenosis levels, suggesting coronary dynamics could play a role in the earlier stages of atherosclerosis development. In the patient case, both dynamics and FSI vs CFD changes altered associations with measured changes in plaque morphology. An appendix of the linear fits between the various FSI- and CFD-based simulations is provided to assist in scaling CFD-based results to resemble the compliant walled characteristics of FSI more accurately.

CONCLUSIONS

These results highlight the potential for coronary dynamics to alter multi-directional/topological WSS metrics which could impact associations with changes in coronary atherosclerosis over time. These results warrant further investigation in a wider range of morphological settings and longitudinal cohort studies in the future.

Numerical investigation of arterial stenosis location affecting hemodynamics considering microcirculation function

BACKGROUND

In recent years, arterial stenosis has become one of the serious diseases threatening people's life and health.

OBJECTIVE

The main purpose of the present study is to examine the changes of hemodynamic parameters in different stenosis locations of arteries.

METHODS

An arterial stenosis model with fluid-structure interaction and microcirculation as the outlet boundary of seepage is adopted in this paper. Considering the interaction between blood and arterial wall, a numerical simulation is carried out using the finite element method.

RESULTS

The results show that hemodynamic parameters are sensitive to the change of stenosis location. The closer to the microcirculation zone the stenosis location, the lower the blood flow velocity, pressure and the wall shear stress. In addition, the velocity trend is transformed from the gradual increase to decrease with the increasing distance away from the inlet when the stenosis location moves to the microcirculation zone.

CONCLUSION

This work proves that the stenosis location has a great influence on hemodynamics based on microcirculation function. Microcirculation is an important factor that cannot be ignored in the numerical simulation of arterial hemodynamics. The numerical results could provide the potential of clinical preconditions for disease diagnosis and treatment.

Numerical analysis of hemodynamics in pulmonary artery stenosis

BACKGROUND: Pulmonary artery stenosis is a serious threat to people’s life and health. OBJECTIVE: The hydrodynamic mechanism of pulmonary artery stenosis is investigated. METHODS: Numerical analysis of hemodynamics in pulmonary artery stenosis using computational fluid dynamics techniques is a well-established method. An idealized model of pulmonary stenosis is established, and the model is divided into main pulmonary artery, right and left pulmonary arteries, and their branches. The sections at different positions are intercepted to study the distribution trend of maximum velocity, pressure and wall shear stress. RESULTS: The numerical simulation results show that the pressure drop at both ends of the narrow is large. High velocity and wall shear stress exist in the center of stenosis, and the wall shear stress at the distal end of stenosis gradually decreases, resulting in endothelial dysfunction. CONCLUSIONS:  To some extent, this study helps clinicians make diagnosis and treatment plans in advance and improve prognosis. This method could be used in the numerical simulation of practical models.

Assessing mechanical vibration-altered wall shear stress in digital arteries

The aim of this study is to implement and validate a method for assessing acute vibration-altered Wall Shear Stress (WSS) in the proper volar digital artery of the non-exposed left forefinger when subjecting the right hand to mechanical vibration. These changes of WSS may be involved in Vibration White Finger. Hence, an experimental device was set-up to link a vibration shaker and an ultra-high frequency ultrasound scanner. The Womersley-based WSS was computed by picking up the maximum velocity from pulse Wave Doppler measurements and extracting the artery diameter from B-mode images through an in-house image processing technique. The parameters of the former method were optimised on numerical ultrasound phantoms of cylindrical and lifelike arteries. These phantoms were computed with the FIELD II and FOCUS platforms which mimicked our true ultrasound device. The Womersley-based WSS were compared to full Fluid Structure Interaction (FSI) and rigid wall models built from resonance magnetic images of a volunteer-specific forefinger artery. Our FSI model took into account the artery's surrounding tissues. The diameter computing procedure led to a bias of 4%. The Womersley-based WSS resulted in misestimating the FSI model by roughly 10% to 20%. No difference was found between the rigid wall computational model and FSI simulations. Regarding the WSS measured on a group of 20 volunteers, the group-averaged basal value was 3 Pa, while the vibration-altered WSS was reduced to 1 Pa, possibly triggering intimal hyperplasia mechanisms and leading to the arterial stenoses encountered in patients suffering from vibration-induced Raynaud's syndrome.

Hydrodynamic Approach for Revealing Venous Anastomotic Stenosis Formation Within a Dialysis Arteriovenous Graft

A conventional arteriovenous graft in patients on dialysis often leads to anastomotic stenosis, which decreases the blood flow rate and increases the risk of complications. In this study, based on hydrodynamics, the pulsatile pressure at the blood vessel graft-vein junction was investigated experimentally and numerically for revealing the causes of stenosis formation and inward remodeling. In the experiments, the pulsatile pressure and displacement at the anastomotic connection were measured at a branched collapsible tube. It was revealed that the pressure becomes negative between pressure peaks of the pulsatile flow; furthermore, tube diameter changes in accordance with the pressure pulsation. Subsequently, numerical simulations revealed that a relatively large pressure difference occurs at the anastomotic connection because of flow collision and separation as compared with the other part, and the pulsatile pressure. Therefore, it is possible that vein at an anastomotic connection may change its shape under pulsating flow. Furthermore, it was found that the pressure difference slightly increased with the anastomosis angle, but the anastomosis angle did not affect the flow rate. Clinical trials in the next step are required to reveal the causal relationship between stenosis and the pulsatile pressure, but the pulsatile flow and its pressure are likely to be one factor in stenosis and inward remodeling.

Numerical simulation of pulsatile blood flow characteristics in a multi stenosed coronary artery

BACKGROUND

Coronary artery disease is reported as one of the most common sources of death all over the world. The presence of stenosis (plaque) in the coronary arteries results in the restriction of blood supply, which leads to myocardial infarction.

OBJECTIVE

The aim of this study was to investigate the effect of multi stenosis on hemodynamics parameters in idealized coronary artery models with varying degrees of stenosis and interspace distance between the stenosis.

METHODS

A finite volume-based software package (Ansys CFX version 17.2) was employed to model the blood flow. The hemodynamic stenosis parameters of blood, such as the pressure, velocity, and wall shear stress were obtained.

RESULTS

The computed results showed that the pressure drop is maximum across the 90% area stenosis (AS). The pressure drop is increased as the distance between the proximal and distal stenosis is decreased across the proximal stenosis for the model P70_D70 during the systolic period of the cardiac cycle. A recirculation zone is formed behind the stenosis and is restricted by the occurrence of distal stenosis as the interspacing distance decreases, which could lead to further progression of stenosis in the flow-disturbed area. The wall shear stress was found to increase as the distance between the proximal and distal stenosis is increased across the distal stenosis. The maximum wall shear stress was found at 90% AS.

CONCLUSIONS

In the clinical diagnosis, an overestimation of distal stenosis severity could be possible. Furthermore, the low wall shear stress zone in between the proximal and distal stenosis may help atherosclerotic growth or merge adjacent stenosis.

Numerical Study of the Unsteady Flow in Simplified and Realistic Iliac Bifurcation Models

Cardiovascular diseases are a major cause of death and disability worldwide and they are commonly associated with the occurrence of atherosclerotic plaque deposition in the vessel walls, a process denoted as atherosclerosis. This is a chronic and progressive inflammatory disease of large-/medium-sized blood vessels that affects blood flow profiles, with the abdominal aorta and its branches being one of the locations prone to the development of this pathology, due to their curvatures and bifurcations. In this regard, the effect of flow patterns was studied and compared for both a simplified three-dimensional model of aorta bifurcation on the iliac arteries and a realistic model of iliac bifurcation, which was constructed from a computational tomography medical image. The flow patterns were analyzed in terms of velocity and wall shear stress distribution, but a special focus was given to the size and location of the recirculation zone. The simulations were performed using the Computational Fluid Dynamics software, FLUENT, taking into account the cardiac cycle profile at the infrarenal aorta. The shear stress and the velocity distribution observed for both models indicated that higher shear stress occurred along the flow divider wall (inner wall) and low shear stress occurred along the outer walls. In addition, the results demonstrated that the wall shear stress profiles were deeply affected by the transient profile of the cardiac cycle, with the deceleration phase being the most critical phase to the occurrence of backflow.

Numerical analysis of stenoses severity and aortic wall mechanics in patients with supravalvular aortic stenosis

Supravalvular aortic stenosis (SVAS) is an aortic malformation characterized by a narrowing of the ascending aorta, resulting in abnormal hemodynamics and pressure drop across the stenosed region. It has been observed that the pressure drops measured from Doppler ultrasound exams often tend to be higher than those obtained from invasive cardiac catheterization. These misleadingly elevated pressure measurements may drive the decision to refer patients for surgical treatment prematurely. Considering this strong clinical association, the purpose of this work is to develop a computational modeling approach using a two-way coupled fluid-structure interaction methodology to determine an accurate prediction of trans-stenotic pressure drop and to further highlight the discrepancy between the SVAS assessment methods. Blood is modeled using Navier-Stokes equations while the aortic wall is simulated by a composite poroelastic structure to represent the three main layers of the arterial wall. The relationship between aortic wall elasticity and the blood flow conditions is examined in varying levels of stenosis, ranging from mild to severe degrees of vessel diameter narrowing. A substantial overestimation of the traditional Doppler pressure drop measurement is observed, especially for severe stenosis levels. The simulation results indicate that elasticity of the aortic wall has a relatively little effect on trans-stenotic pressure drop for the range of mild to moderate SVAS cases, but predicted to have a profound effect for severe SVAS cases. Moreover, significant sensitivity to the pressure drop across the SVAS region from stenosis severity is observed.

Effect of stenosis and dilatation on the hemodynamic parameters associated with left coronary artery

BACKGROUND AND OBJECTIVE

The main objective of the work is to examine the curvature effects of stenosis/dilatation region pertaining to left coronary artery. The hemodymamic features during the cardiac cycle is thoroughly examined.

METHODS

A numerical fluid structure interaction model incorporating multi- layered elastic artery wall, non-Newtonian blood viscosity and pulsating boundary conditions is developed. The composite arterial wall consists of a thin layer tunica intima, atheroma and a thick wall. Higher stiffness of atheroma is captured by using higher Young's modulus. The CFD and FSI models are validated with available experimental and analytical data. Computations are done with five different non-Newtonian models and arterial wall with various elasticity levels. The local and time averaged WSS, velocity contours downstream of stenosis, wall pressure and pressure drop during various phases of cardiac cycle are provided in detail.

RESULTS

The influence of non-Newtonian effects of blood viscosity is found to be significant especially at stenosis regions. The flexible wall caused wall deformation and the associated flow and pressure wave propagation affecting WSS and pressure drop compared to the rigid wall. Flow recirculation is noticed at stenosis downstream locations and its strength increases with increased severity of the stenosis. A stenosis is characterised by a sudden drop in wall pressure and a slower two stage recovery during peak velocity periods of the cardiac cycle.

CONCLUSIONS

The pressure drop, local WSS at stenosis centre, and radial velocity increase are significantly higher for stenosis cases and the effect is severe during peak diastole. The variation in hemodynamic parameters is found to be less significant for dilatation. Significantly lower WSS is noticed for the recirculation regions downstream of stenosis which can enhance the tendency for monocytes to attach to the endothelium. The radius of curvature of the stenosis is found to be the most sensitive parameter affecting the hemodynamic characteristics rather than the detailed geometry of the stenosis. The main effect of variation of artery wall stiffness is noted at recirculation regions present downstream of stenosis. The results from the study may be useful for predicting wall shear stress signatures associated with stenosis/dilatation changes and the management of specific cases.

2,3,7,8-Tetrachlorodibenzo-p-dioxin exposure disrupts development of the visceral and ocular vasculature

The aryl hydrocarbon receptor (AHR) has endogenous functions in mammalian vascular development and is necessary for mediating the toxic effects of a number of environmental contaminants. Studies in mice have demonstrated that AHR is necessary for the formation of the renal, retinal, and hepatic vasculature. In fish, exposure to the prototypic AHR agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) induces expression of the AHR biomarker cyp1a throughout the developing vasculature and produces vascular malformations in the head and heart. However, it is not known whether the vascular structures that are sensitive to loss of AHR function are also disrupted by aberrant AHR activation. Here, we report that TCDD-exposure in zebrafish disrupts development of 1) the subintestinal venous plexus (SIVP), which vascularizes the developing liver, kidney, gut, and pancreas, and 2) the superficial annular vessel (SAV), an essential component of the retinal vasculature. Furthermore, we determined that TCDD exposure increased the expression of bmp4, a key molecular mediator of SIVP morphogenesis. We hypothesize that the observed SIVP phenotypes contribute to one of the hallmarks of TCDD exposure in fish - the failure of the yolk sac to absorb. Together, our data describe novel TCDD-induced vascular phenotypes and provide molecular insight into critical factors producing the observed vascular malformations.

Blocking von Willebrand factor free thiols inhibits binding to collagen under high and pathological shear stress

BACKGROUND

Von Willebrand factor (VWF) contains a number of free thiols, the majority of which are located in its C-domains, and these have been shown to alter VWF function, However, the impact of free thiols on function following acute exposure of VWF to collagen under high and pathological shear stress has not been determined.

METHODS

VWF free thiols were blocked with N-ethylmaleimide and flow assays performed under high and pathological shear rates to determine the impact on platelet capture and collagen binding function. Atomic force microscopy (AFM) was used to probe the interaction of VWF with collagen and molecular simulations conducted to determine the effect of free thiols on the flexibility of the VWF-C4 domain.

RESULTS

Blockade of VWF free thiols reduced VWF-mediated platelet capture to collagen in a shear-dependent manner, with platelet capture virtually abolished above 5000 s-1 and in regions of stenosis in microfluidic channels. Direct visualization of VWF fibers formed under extreme pathological shear rates and analysis of collagen-bound VWF attributed the effect to altered binding of VWF to collagen. AFM measurements showed that thiol-blockade reduced the lifetime and strength of the VWF-collagen bond. Pulling simulations of the VWF-C4 domain demonstrated that with one or two reduced disulphide bonds the C4 domain has increased flexibility and the propensity to undergo free-thiol exchange.

CONCLUSIONS

We conclude that free thiols in the C-domains of VWF enhance the flexibility of the molecule and enable it to withstand high shear forces following collagen binding, demonstrating a previously unrecognized role for VWF free thiols.

Numerical modeling in arterial hemodynamics incorporating fluid-structure interaction and microcirculation

Background

The effects of arterial wall compliance on blood flow have been revealed using fluid-structure interaction in last decades. However, microcirculation is not considered in previous researches. In fact, microcirculation plays a key role in regulating blood flow. Therefore, it is very necessary to involve microcirculation in arterial hemodynamics.

Objective

The main purpose of the present study is to investigate how wall compliance affects the flow characteristics and to establish the comparisons of these flow variables with rigid wall when microcirculation is considered.

Methods

We present numerical modeling in arterial hemodynamics incorporating fluid-structure interaction and microcirculation. A novel outlet boundary condition is employed to prescribe microcirculation in an idealised model.

Results

The novel finding in this work is that wall compliance under the consideration of microcirculation leads to the increase of wall shear stress in contrast to rigid wall, contrary to the traditional result that wall compliance makes wall shear stress decrease when a constant or time dependent pressure is specified at an outlet.

Conclusions

This work provides the valuable study of hemodynamics under physiological and realistic boundary conditions and proves that wall compliance may have a positive impact on wall shear stress based on this model. This methodology in this paper could be used in real model simulations.

Coronary arteries hemodynamics: effect of arterial geometry on hemodynamic parameters causing atherosclerosis

Coronary arteries have high curvatures, and hence, flow through them causes disturbed flow patterns, resulting in stenosis and atherosclerosis. This in turn decreases the myocardial flow perfusion, causing myocardial ischemia and infarction. Therefore, in order to understand the mechanisms of these phenomena caused by high curvatures and branching of coronary arteries, we have conducted elaborate hemodynamic analysis for both (i) idealized coronary arteries with geometrical parameters representing realistic curvatures and stenosis and (ii) patient-specific coronary arteries with stenoses. Firstly, in idealized coronary arteries with approximated realistic arterial geometry representative of their curvedness and stenosis, we have computed the hemodynamic parameters of pressure drop, wall shear stress (WSS) and wall pressure gradient (WPG), and their association with the geometrical parameters of curvedness and stenosis. Secondly, we have similarly determined the wall shear stress and wall pressure gradient distributions in four patient-specific curved stenotic right coronary arteries (RCAs), which were reconstructed from medical images of patients diagnosed with atherosclerosis and stenosis; our results show high WSS and WPG regions at the stenoses and inner wall of the arterial curves. This paper provides useful insights into the causative mechanisms of the high incidence of atherosclerosis in coronary arteries. It also provides guidelines for how simulation of blood flow in patient's coronary arteries and determination of the hemodynamic parameters of WSS and WPG can provide a medical assessment of the risk of development of atherosclerosis and plaque formation, leading to myocardial ischemia and infarction. The novelty of our paper is in our showing how in actual coronary arteries (based on their CT imaging) curvilinearity and narrowing complications affect the computed WSS and WPG, associated with risk of atherosclerosis. This is very important for cardiologists to be able to properly take care of their patients and provide remedial measures before coronary complications lead to myocardial infarctions and necessitate stenting or coronary bypass surgery. We want to go one step further and provide clinical application of our research work. For that, we are offering to cardiologists worldwide to carry out hemodynamic analysis of the medically imaged coronary arteries of their patients and compute the values of the hemodynamic parameters of WSS and WPG, so as to provide them an assessment of the risk of atherosclerosis for their patients. Graphical abstract Theme and aims: Coronary arteries have high curvatures, and hence flow through them causes disturbed flow patterns, resulting in stenosis and atherosclerosis. This in turn decreases the myocardial flow perfusion, causing myocardial ischemia and infarction. Therefore, in order to understand the mechanisms of these phenomena caused by high curvatures and branching of coronary arteries, we have conducted elaborate hemodynamic analysis for both (i) idealized coronary arteries with geometrical parameters representing curvatures and stenosis, and (ii) patient-specific coronary arteries with stenoses. Methods and results: Firstly, in idealized coronary arteries with approximated realistic arterial geometry representative of their curvedness and stenosis, we have computed the hemodynamic parameters of pressure drop, wall shear stress (WSS) and wall pressure gradient (WPG), and their association with the geometrical parameters of curvedness and stenosis. Then, we have determined the wall shear stress and wall pressure gradient distributions in four patient-specific curved stenotic right coronary arteries (RCAs), that were reconstructed from medical images of patients diagnosed with atherosclerosis and stenosis, as illustrated in Figure 1, in which the locations of the stenoses are highlighted by arrows. Figure 1: Three-dimensional CT visualization of arteries in patients with suspected coronary disease. The arteries can be seen as a combination of various curved segments with stenoses at unspecific locations highlighted by arrows. Our results show high WSS and WPG regions at the stenoses and inner wall of the arterial curves, as depicted in Figure 2. Therein, the encapsulations show (i) high WSS, and (ii) high WPG regions at the stenosis and inner wall of the arterial curves. Figure 2: WSS and WPG surface plot of realistic arteries (a), (b), (c) and (d), wherein the small squared parts are enlarged to show the detailed localized contour plots at the stenotic regions. Therein, the circular encapsulations show (i) high WSS and (ii) high WPG regions at the stenosis and inner wall of the arterial curves. Conclusion and novelty: This paper provides useful insights into the causative mechanisms of the high incidence of atherosclerosis in coronary arteries. It also provides guidelines for how simulation of blood flow in patient coronary arteries and determination of the hemodynamic parameters of WSS and WPG can provide a medical assessment of the risk of development of atherosclerosis and plaque formation, leading to myocardial ischemia and infarction. The novelty of our paper is our showing how in actual coronary arteries (based on their CT imaging), curvilinearity and narrowing complications affect the computed WSS and WPG associated with risk of atherosclerosis. This is very important for cardiologists to be able to properly take care of their patients and provide remedial measures before coronary complications lead to myocardial infarctions and necessitate stenting or coronary bypass surgery.

3D Printed Biomodels for Flow Visualization in Stenotic Vessels: An Experimental and Numerical Study

Atherosclerosis is one of the most serious and common forms of cardiovascular disease and a major cause of death and disability worldwide. It is a multifactorial and complex disease that promoted several hemodynamic studies. Although in vivo studies more accurately represent the physiological conditions, in vitro experiments more reliably control several physiological variables and most adequately validate numerical flow studies. Here, a hemodynamic study in idealized stenotic and healthy coronary arteries is presented by applying both numerical and in vitro approaches through computational fluid dynamics simulations and a high-speed video microscopy technique, respectively. By means of stereolithography 3D printing technology, biomodels with three different resolutions were used to perform experimental flow studies. The results showed that the biomodel printed with a resolution of 50 μm was able to most accurately visualize flow due to its lowest roughness values (Ra = 1.8 μm). The flow experimental results showed a qualitatively good agreement with the blood flow numerical data, providing a clear observation of recirculation regions when the diameter reduction reached 60%.

Comparison of 3D reconstruction methods based on different cardiovascular imaging: a study of multimodality reconstruction method

Coronary arterial imaging and the assessment of the severity of arterial stenoses can be achieved with several modalities classified mainly according to their invasive or noninvasive nature. These modalities can be further utilized for the 3-dimensional (3D) reconstruction of the arterial geometry. This study aims to determine the prediction performance of atherosclerotic disease progression using reconstructed arteries from three reconstruction methodologies: Quantitative Coronary Analysis (QCA), Virtual Histology Intravascular Ultrasound (VH)-IVUS-Angiography fusion method and Coronary Computed Tomography Angiography (CCTA). The accuracy of the reconstruction methods is assessed using several metrics such as Minimum lumen diameter (MLD), Reference vessel diameter (RVD), Lesion length (LL), Diameter stenosis (DS%) and the Mean wall shear stress (WSS). Five patients in a retrospective study who underwent X-ray angiography, VH-IVUS and CCTA are used for the method evaluation.

A hemodynamic model with a seepage condition and fluid-structure interactions for blood flow in arteries with symmetric stenosis

To strengthen the detailed understanding of arterial stenosis, we construct a novel hemodynamic model. Frequently used symmetric stenosis is employed in this work. Being different from a traditional model, this numerical model adopts microcirculation resistance as an outlet boundary condition, which is called a seepage condition. Meanwhile, fluid-structure interactions are used in the numerical simulation considering the interrelationship of blood and arterial wall. Our results indicate that (i) the region upstream of stenosis experiences very high pressures during cardiac cycles, (ii) pressure drops much faster as the flow moves into the stenotic region, and (iii) high flow velocities and high shear stresses occur in the post-stenosis region. This work provides evidence that there is a strong effect of the function of microcirculation on stenosis. This contributes to evaluating potential stenotic behavior in arteries and is pivotal in guiding disease treatment.

A seepage outlet boundary condition in hemodynamics modeling

BACKGROUND

In computational fluid dynamics (CFD) models for hemodynamics applications, boundary conditions remain one of the major issues in obtaining accurate fluid flow predictions.

OBJECTIVE

As an important part of the arterial circulation, microcirculation plays important roles in many aspects, such as substance exchange, interstitial fluid generation and inverse flow. It is necessary to consider microcirculation in hemodynamics modeling. This is a methodological paper to test and validate a new type of boundary condition never applied to microcirculation before.

METHODS

In order to address this issue, we introduce microcirculation as a seepage outlet boundary condition in computational hemodynamics. Microcirculation is treated as a porous medium in this paper. Numerical comparisons of the seepage and traditional boundary conditions are made.

RESULTS

The results show that the seepage boundary condition has significant impacts on numerical simulation. Under the seepage boundary condition, the fluctuation range of the pressures progressively rises in the artery zone. The results obtained from the traditional boundary condition show that the pressure fluctuation range gradually falls. In addition, the wall shear stresses under the traditional outlet boundary condition are much higher than those under the seepage outlet boundary condition.

CONCLUSIONS

The proposed boundary condition is more suitable in hemodynamics modeling.

Structural modelling of the cardiovascular system

Computational modelling of the cardiovascular system offers much promise, but represents a truly interdisciplinary challenge, requiring knowledge of physiology, mechanics of materials, fluid dynamics and biochemistry. This paper aims to provide a summary of the recent advances in cardiovascular structural modelling, including the numerical methods, main constitutive models and modelling procedures developed to represent cardiovascular structures and pathologies across a broad range of length and timescales; serving as an accessible point of reference to newcomers to the field. The class of so-called hyperelastic materials provides the theoretical foundation for the modelling of how these materials deform under load, and so an overview of these models is provided; comparing classical to application-specific phenomenological models. The physiology is split into components and pathologies of the cardiovascular system and linked back to constitutive modelling developments, identifying current state of the art in modelling procedures from both clinical and engineering sources. Models which have originally been derived for one application and scale are shown to be used for an increasing range and for similar applications. The trend for such approaches is discussed in the context of increasing availability of high performance computing resources, where in some cases computer hardware can impact the choice of modelling approach used.

A multiscale approach for determining the morphology of endothelial cells at a coronary artery

The morphology of endothelial cells (ECs) may be an indication for determining atheroprone sites. Until now, there has been no clinical imaging technique to visualize the morphology of ECs in the arteries. The present study introduces a computational technique for determining the morphology of ECs. This technique is a multiscale simulation consisting of the artery scale and the cell scale. The artery scale is a fluid-structure interaction simulation. The input for the artery scale is the geometry of the coronary artery, that is, the dynamic curvature of the artery due to the cardiac motion, blood flow, blood pressure, heart rate, and the mechanical properties of the blood and the arterial wall, the measurements of which can be obtained for a specific patient. The results of the artery scale are wall shear stress (WSS) and cyclic strains as the mechanical stimuli of ECs. The cell scale is an inventive mass-and-spring model that is able to determine the morphological response of ECs to any combination of mechanical stimuli. The results of the multiscale simulation show the morphology of ECs at different locations of the coronary artery. The results indicate that the atheroprone sites have at least 1 of 3 factors: low time-averaged WSS, high angle of WSS, and high longitudinal strain. The most probable sites for atherosclerosis are located at the bifurcation region and lie on the myocardial side of the artery. The results also indicated that a higher dynamic curvature is a negative factor and a higher pulse pressure is a positive factor for protection against atherosclerosis.

EFFECTS OF POROSITY IN A SEEPAGE MODEL ON HEMODYNAMICS

In this paper, we numerically investigate the effects of porosity on hemodynamics using a seepage outlet boundary condition. The results with the porosities [Formula: see text] and [Formula: see text] are compared. The comparison of the results shows that there are identical pressure, velocity and wall shear stress distributions in the artery zone; however, the velocities and wall shear stresses in the microcirculation zone with [Formula: see text] are almost as twice as those with [Formula: see text]. This work suggests that the porosity changes do not provoke variations in the artery zone, but they lead to the differences of the velocities and wall shear stresses in the microcirculation zone. The results may be the explanation of syndrome X.

The Coronary Artery Running Pattern is One of the Causes of Individual Differences in the Progression of Coronary Atherosclerosis in WHHLMI Rabbits, an Animal Model for Coronary Atherosclerosis

AIMS

The relationship between the coronary artery running pattern and development of coronary lesions was re-examined herein using WHHLMI rabbits, an animal model of spontaneous coronary atherosclerosis.

METHODS

The coronary artery running pattern was analyzed using an X-ray computed tomography (CT) apparatus after perfusion. Pathological sections were prepared (Victoria blue-HE staining) at 100 µm intervals from the origin of the left circumflex artery (LCX). The severity of coronary lesions was evaluated based on cross-sectional narrowing (lesion area/inner area of the internal elastic lamina).

RESULTS

In the CT analysis, the angle of the main curvature of the LCX negatively correlated with the percentage of sections with lesions and cross-sectional narrowing. The percentage of sections with lesions was significantly higher in acute angle-type LCX than in obtuse angle-type LCX. Cross-sectional narrowing was also significantly greater in acute angle-type LCX than in obtuse angle-type LCX. The percentage of fibrous lesions was high at the proximal region of LCX, whereas that of lipid-rich lesions was high at the curvature. In 24 months age group, the percentage of sections with calcification in acute angle-type LCX was about twice that in obtuse angle-type LCX.

CONCLUSIONS

Individual differences were observed in the angle of the main curvature of the LCX, which affected the occurrence and extension of atherosclerotic lesions.

Characteristics of pulsatile flows in curved stenosed channels

Spatial and temporal variations of the hemodynamic features occur under pulsatile conditions in complex vessel geometry. Wall shear stress affected by the disturbed flow can result in endothelial cell dysfunction, which leads to atherogenesis and thrombosis. Therefore, detailed understanding of the hemodynamic characteristics in a curved stenosed channel is highly important when examining the pathological effects of hemodynamic phenomena on the progression of atherosclerosis. The present study measures the velocity fields of pulsatile flows with three different Reynolds numbers in 3D curved vessel models with stenosis using time-resolved particle image velocimetry (PIV). Three different models were cast in PDMS polymer using models made by a 3D printer with different bend angles of 0°, 10°, and 20° between the longitudinal axes at the upstream and downstream of the stenosis. To investigate the 3D flow structures, a stack of 2D velocity fields was obtained by adjusting the position of the laser sheet along the Z-direction. The structures of flow fields in the stenosed models were analyzed using the distribution of the shearing strain as well as the skewness and full width at half maximum of the velocity profile. To support experiment results, distributions of pressure and 3D vortex in the curved stenosed channels were estimated by conducting the numerical simulation. These results indicate that the curvature of the tube considerably influences the skewness of the flow, and the shear stress is intensified near the outer curvature wall due to centrifugal force. The results would be helpful in understanding the effects of geometrical factors on plaque rupture and severe cardiovascular diseases.

Systems Analysis of Thrombus Formation

The systems analysis of thrombosis seeks to quantitatively predict blood function in a given vascular wall and hemodynamic context. Relevant to both venous and arterial thrombosis, a Blood Systems Biology approach should provide metrics for rate and molecular mechanisms of clot growth, thrombotic risk, pharmacological response, and utility of new therapeutic targets. As a rapidly created multicellular aggregate with a polymerized fibrin matrix, blood clots result from hundreds of unique reactions within and around platelets propagating in space and time under hemodynamic conditions. Coronary artery thrombosis is dominated by atherosclerotic plaque rupture, complex pulsatile flows through stenotic regions producing high wall shear stresses, and plaque-derived tissue factor driving thrombin production. In contrast, venous thrombosis is dominated by stasis or depressed flows, endothelial inflammation, white blood cell-derived tissue factor, and ample red blood cell incorporation. By imaging vessels, patient-specific assessment using computational fluid dynamics provides an estimate of local hemodynamics and fractional flow reserve. High-dimensional ex vivo phenotyping of platelet and coagulation can now power multiscale computer simulations at the subcellular to cellular to whole vessel scale of heart attacks or strokes. In addition, an integrated systems biology approach can rank safety and efficacy metrics of various pharmacological interventions or clinical trial designs.

Computational estimation of the severity of coronary lesions with intravascular ultrasound images: a pilot study

The fast and accurate hemodynamic functional assessment of the coronary vasculature is of utmost importance in clinical practice due to the fact that Cardiovascular Diseases have become the leading cause of death globally. In this work we propose a novel method that combines two of the most efficient methods of hemodynamic status assessment of coronary arteries, Intravascular UtraSound and virtual Functional Assessment Index, an index that correlates well to the measured Fractional Flow Reserve. One Left Anterior Descending segment was reconstructed both in a straight manner (using only IVUS images) as well as using the actual 3D geometry of the vessel (using IvUS images combined with the respective coronary angiographic images [2]). The generated vFAI values were almost identical (Straight=0.80, 3D=0.79), presenting a relative error of 1.27%, thus proving the efficacy of the proposed method. We also calculated the Endothelial Shear Stress for the two models under rest (i.e. flow rate of 1 ml/s), observing a similar trend throughout the artery, but with a statistically important relative error of 13.49%, as expected.

Computational evaluation of smoothed particle hydrodynamics for implementing blood flow modelling through CT reconstructed arteries

Simulation of blood flow in a stenosed artery using Smoothed Particle Hydrodynamics (SPH) is a new research field, which is a particle-based method and different from the traditional continuum modelling technique such as Computational Fluid Dynamics (CFD). Both techniques harness parallel computing to process hemodynamics of cardiovascular structures. The objective of this study is to develop and test a new robust method for comparison of arterial flow velocity contours by SPH with the well-established CFD technique, and the implementation of SPH in computed tomography (CT) reconstructed arteries. The new method was developed based on three-dimensional (3D) straight and curved arterial models of millimeter range with a 25% stenosis in the middle section. In this study, we employed 1,000 to 13,000 particles to study how the number of particles influences SPH versus CFD deviation for blood-flow velocity distribution. Because further increasing the particle density has a diminishing effect on this deviation, we have determined a critical particle density of 1.45 particles/mm2 based on Reynolds number (Re = 200) at the inlet for an arterial flow simulation. Using this critical value of particle density can avoid unnecessarily big computational expenses that have no further effect on simulation accuracy. We have particularly shown that the SPH method has a big potential to be used in the virtual surgery system, such as to simulate the interaction between blood flow and the CT reconstructed vessels, especially those with stenosis or plaque when encountering vasculopathy, and for employing the simulation results output in clinical surgical procedures.

High shear stress induces atherosclerotic vulnerable plaque formation through angiogenesis

Rupture of atherosclerotic plaques causing thrombosis is the main cause of acute coronary syndrome and ischemic strokes. Inhibition of thrombosis is one of the important tasks developing biomedical materials such as intravascular stents and vascular grafts. Shear stress (SS) influences the formation and development of atherosclerosis. The current review focuses on the vulnerable plaques observed in the high shear stress (HSS) regions, which localizes at the proximal region of the plaque intruding into the lumen. The vascular outward remodelling occurs in the HSS region for vascular compensation and that angiogenesis is a critical factor for HSS which induces atherosclerotic vulnerable plaque formation. These results greatly challenge the established belief that low shear stress is important for expansive remodelling, which provides a new perspective for preventing the transition of stable plaques to high-risk atherosclerotic lesions.

Numerical analysis of the effect of turbulence transition on the hemodynamic parameters in human coronary arteries

BACKGROUND

Local hemodynamics plays an important role in atherogenesis and the progression of coronary atherosclerosis disease (CAD). The primary biological effect due to blood turbulence is the change in wall shear stress (WSS) on the endothelial cell membrane, while the local oscillatory nature of the blood flow affects the physiological changes in the coronary artery. In coronary arteries, the blood flow Reynolds number ranges from few tens to several hundreds and hence it is generally assumed to be laminar while calculating the WSS calculations. However, the pulsatile blood flow through coronary arteries under stenotic condition could result in transition from laminar to turbulent flow condition.

METHODS

In the present work, the onset of turbulent transition during pulsatile flow through coronary arteries for varying degree of stenosis (i.e., 0%, 30%, 50% and 70%) is quantitatively analyzed by calculating the turbulent parameters distal to the stenosis. Also, the effect of turbulence transition on hemodynamic parameters such as WSS and oscillatory shear index (OSI) for varying degree of stenosis is quantified. The validated transitional shear stress transport (SST) k-ω model used in the present investigation is the best suited Reynolds averaged Navier-Stokes turbulence model to capture the turbulent transition. The arterial wall is assumed to be rigid and the dynamic curvature effect due to myocardial contraction on the blood flow has been neglected.

RESULTS

Our observations shows that for stenosis 50% and above, the WSSavg, WSSmax and OSI calculated using turbulence model deviates from laminar by more than 10% and the flow disturbances seems to significantly increase only after 70% stenosis. Our model shows reliability and completely validated.

CONCLUSIONS

Blood flow through stenosed coronary arteries seems to be turbulent in nature for area stenosis above 70% and the transition to turbulent flow begins from 50% stenosis.

Validation study of a 3D-QCA coronary reconstruction method using a hybrid intravascular ultrasound and angiography reconstruction method and patient-specific Fractional Flow Reserve data

The estimation of the severity of coronary lesions is of utmost importance in today's clinical practice, since Cardiovascular diseases often have fatal consequences. The most efficient method to estimate the severity of a lesion is the calculation of the Fractional Flow Reserve. The necessary use of a pressure wire, however, makes this method invasive and strenuous for the patient. In this work, we present a novel 3-Dimensional Quantitative Coronary Analysis coronary reconstruction method and a framework for the computation of the virtual Functional Assessment Index (vFAI). In a dataset of 5 coronary arterial segments, we use the aforementioned method to reconstruct them in 3D, and compare them to the respective 3D models reconstructed from our already validated hybrid IVUS-angiography reconstruction method [2]. The obtained results indicate a high correlation between the two methods in terms of the calculated FFR values, presenting a difference of 3.19% in the worst case scenario. Furthermore, when compared to the actual FFR values that derive from a pressure wire, the differences were statistically insignificant.

FSI SIMULATION OF A HEALTHY CORONARY BIFURCATION FOR STUDYING THE MECHANICAL STIMULI OF ENDOTHELIAL CELLS UNDER DIFFERENT PHYSIOLOGICAL CONDITIONS

Atherosclerosis is a world-spread and well-known disease. This disease strongly relates to the endothelial cells (ECs) function. Normally, the endothelial cells align in the flow direction in the atheroprotected sites; however, in the case of atheroprone sites these cells orient randomly. The mechanical stimuli such as wall shear stress and strains could determine the morphology and function of the endothelial cells. In the present study, we numerically simulated the left main coronary artery (LCA) and its branches to left anterior descending (LAD) and left circumflex coronary (LCX) artery using fluid–structure interaction (FSI) modeling. The results were presented as longitudinal and circumferential strains of ECs as well as wall shear stress. Wide ranges of heart rate, cardiac motion, systolic and diastolic pressures were considered and their effects on mechanical stimuli were described in detail. The results showed that these factors could greatly influence the risk of atherosclerosis and the location of atherosclerotic lesions.

Computational assessment of the fractional flow reserve from intravascular ultrasound and coronary angiography data: a pilot study

Cardiovascular disease is one of the primary causes of morbidity and mortality around the globe. Thus, the diagnosis of critical lesions in coronary arteries is of utmost importance in clinical practice. One useful and efficient method to assess the functional severity of one or multiple lesions in a coronary artery is the calculation of the fractional flow reserve (FFR). In the current work, we present a method which allows the calculation of the FFR value computationally, without the use of a pressure wire and the induction of hyperemia, using intravascular ultrasound (IVUS) and biplane angiography images for three-dimensional (3D) coronary artery reconstruction and measurements of the volumetric flow rate derived from angiographic sequences. The simulated FFR values were compared to the invasively measured FFR values in 7 cases, presenting high correlation (r=0.85) and good agreement (mean difference=0.002). FFR assessment without employing a pressure wire and the induction of hyperemia is feasible using 3D reconstructed coronary artery models from angiographic and IVUS data coupled with computational fluid dynamics.

A comparative study of three speckle reducing methods for intima-media thickness ultrasound images

BACKGROUND

Ultrasonic evaluation of intima-media thickness (IMT) is an early marker of assessing the development of atherosclerosis and determining cardiovascular risk. To attain the best possible diagnosis, it is essential that medical images be clear, sharp and without noise and artifacts.

OBJECTIVES

Comparison of speckle reducing anisotropic diffusion (SRAD), discrete (DTD) and continuum topological derivative (CTD) on B-mode ultrasound images of common carotid and brachial arteries throughout the cardiac cycle.

PATIENTS AND METHODS

In a cross-sectional design, an examination was performed on forty-two human subjects with a mean age of 44 ± 6 years from April 2013 to June 2013. This study was approved by the ethics committees of Kashan University of Medical Sciences and Beheshti Hospital. An ultrasonic examination of common carotid and brachial arteries of forty-two human subjects was performed. The program was designed in MATLAB software to extract consecutive B-mode images and apply region of interest (ROI) on the IMT of the common carotid and brachial arteries. Then, three different noise reduction filters with the Canny edge detection were used in ROI separately. Finally, the program measured the image quality metrics.

RESULTS

According to values of eleven different image quality metrics (mentioned in the main text), there was a significant difference between CTD, DTD and SRAD filters with the Canny edge detection status in the common carotid and brachial arteries throughout the cardiac cycle (all P values < 0.001). For example, peak signal to noise ratios (PSNR) using CTD, DTD and SRAD filters were 95.43 ± 0.64, 88.86 ± 0.82 and 73.02 ± 0.20 in common carotid and 96.39 ± 1.25, 92.58 ± 0.11 and 88.27 ± 0.63 in brachial arteries, respectively (both P values < 0.001).

CONCLUSIONS

By measuring image quality metrics, this study showed that DTD and CTD filters with the Canny edge detection respectively, are better than SRAD filter with the Canny detection for speckle suppression and details preservation in both arteries in the ultrasound images.

A Computational Model for Biomechanical Effects of Arterial Compliance Mismatch

Background. Compliance mismatch is a negative factor and it needs to be considered in arterial bypass grafting. Objective. A computational model was employed to investigate the effects of arterial compliance mismatch on blood flow, wall stress, and deformation. Methods. The unsteady blood flow was assumed to be laminar, Newtonian, viscous, and incompressible. The vessel wall was assumed to be linear elastic, isotropic, and incompressible. The fluid-wall interaction scheme was constructed using the finite element method. Results. The results show that there are identical wall shear stress waveforms, wall stress, and strain waveforms at different locations. The comparison of the results demonstrates that wall shear stresses and wall strains are higher while wall stresses are lower at the more compliant section. The differences promote the probability of intimal thickening at some locations. Conclusions. The model is effective and gives satisfactory results. It could be extended to all kinds of arteries with complicated geometrical and material factors.

Patient-specific simulation of coronary artery pressure measurements: an in vivo three-dimensional validation study in humans

Pressure measurements using finite element computations without the need of a wire could be valuable in clinical practice. Our aim was to compare the computed distal coronary pressure values with the measured values using a pressure wire, while testing the effect of different boundary conditions for the simulation. Eight coronary arteries (lumen and outer vessel wall) from six patients were reconstructed in three-dimensional (3D) space using intravascular ultrasound and biplane angiographic images. Pressure values at the distal and proximal end of the vessel and flow velocity values at the distal end were acquired with the use of a combo pressure-flow wire. The 3D lumen and wall models were discretized into finite elements; fluid structure interaction (FSI) and rigid wall simulations were performed for one cardiac cycle both with pulsatile and steady flow in separate simulations. The results showed a high correlation between the measured and the computed coronary pressure values (coefficient of determination [r²] ranging between 0.8902 and 0.9961), while the less demanding simulations using steady flow and rigid walls resulted in very small relative error. Our study demonstrates that computational assessment of coronary pressure is feasible and seems to be accurate compared to the wire-based measurements.

Assessing the hemodynamic influence between multiple lesions in a realistic right coronary artery segment: A computational study

Coronary artery disease is the primary cause of morbidity and mortality worldwide. Therefore, detailed assessment of lesions in the coronary vasculature is critical in current clinical practice. Fractional flow reserve (FFR) has been proven as an efficient method for assessing the hemodynamic severity of a coronary stenosis. However, functional assessment of a coronary segment with multiple stenoses (≥ 2) remains complex for guiding the strategy of percutaneous coronary intervention due to the hemodynamic interplay between adjacent stenoses. In this work, we created four 3-dimensional (3D) arterial models that derive from a healthy patient-specific right coronary artery segment. The initial healthy model was reconstructed using fusion of intravascular ultrasound (IVUS) and biplane angiographic patient data. The healthy 3D model presented a measured FFR value of 0.96 (pressure-wire) and a simulated FFR value of 0.98. We then created diseased models with two artificial sequential stenoses of 90% lumen area reduction or with the proximal and distal stenosis separately. We calculated the FFR value for each case: 0.65 for the case with the two stenoses, 0.73 for the case with the distal stenosis and 0.90 for the case with the proximal stenosis. This leads to the conclusion that although both stenoses had the same degree of lumen area stenosis, there was a large difference in hemodynamic severity, thereby indicating that angiographic lumen assessment by itself is often not adequate for accurate assessment of coronary lesions.

Effect of residual stress on peak cap stress in arteries

Vulnerable plaques are a subset of atherosclerotic plaques that are prone to rupture when high stresses occur in the cap. The roles of residual stress, plaque morphology, and cap stiffness on the cap stress are not completely understood. Here, arteries are modeled within the framework of nonlinear elasticity as incompressible cylindrical structures that are residually stressed through differential growth. These structures are assumed to have a nonlinear, anisotropic, hyperelastic response to stresses in the media and adventitia layers and an isotropic response in the intima and necrotic layers. The effect of differential growth on the peak stress is explored in a simple, concentric geometry and it is shown that axial differential growth decreases the peak stress in the inner layer. Furthermore, morphological risk factors are explored. The peak stress in residually stressed cylinders is not greatly affected by changing the thickness of the intima. The thickness of the necrotic layer is shown to be the most important morphological feature that affects the peak stress in a residually stressed vessel.

Arterial pressure and flow wave analysis using time-domain 1-D hemodynamics

We reviewed existing methods for analyzing, in the time domain, physical mechanisms underlying the patterns of blood pressure and flow waveforms in the arterial system. These are wave intensity analysis and separations into several types of waveforms: (i) forward- and backward-traveling, (ii) peripheral and conduit, or (iii) reservoir and excess. We assessed the physical information provided by each method and showed how to combine existing methods in order to quantify contributions to numerically generated waveforms from previous cardiac cycles and from specific regions and properties of the numerical domain: the aortic root, arterial bifurcations and tapered vessels, peripheral reflection sites, and the Windkessel function of the aorta. We illustrated our results with numerical examples involving generalized arterial stiffening in a distributed one-dimensional model or localized changes in the model parameters due to a femoral stenosis, carotid stent or abdominal aortic aneurysm.

WALL SHEAR STRESSES IN A FLUID–STRUCTURE INTERACTION MODEL OF PULSE WAVE PROPAGATION

In our prior paper, a fluid–structure interaction model of pulse wave propagation, called the elastic tube model, has been developed. The focus of this paper is wall shear stress (WSS) in this model and the effects of different parameters, including rigid walls, wall thickness, and internal radius. The unsteady flow was assumed to be laminar, Newtonian and incompressible, and the vessel wall to be linear-elastic isotropic, and incompressible. A fluid–structure interaction scheme is constructed using a finite element method. The results demonstrate the elastic tube plays an important role in WSS distributions of wave propagation. It is shown that there is a time delay between the WSS waveforms at different locations in the elastic tube model while the time delay cannot be observed clearly in the rigid tube model. Compared with the elastic tube model, the increase of the wall thickness makes disturbed WSS distributions, however WSS values are increased greatly due to the decrease of the internal radius. The results indicate that the effects of different parameters on WSS distributions are significant. The proposed model gives valid results.

Effect of postural changes on normal and stenosed common carotid artery using FSI

Gravity associated with postural changes has a strong bearing on haemodynamics of blood flow in arteries. Its effect on stenosed cases has not been widely investigated. In the present study, variation observed in blood flow during postural changes is investigated for different conditions like standing, sleeping and head-down position. A fluid structure interaction study is carried out for idealized normal and 75% eccentric and concentric stenosed common carotid normal artery. The results clearly indicate the effects of altered gravity on flow conditions. It was found to be very significant during head-down position and demonstrated very high arterial blood pressure in stenosed common carotid when compared with normal carotid.

Quantification of Arterial Wall Inhomogeneity Size, Distribution, and Modulus Contrast Using FSI Numerical Pulse Wave Propagation

Changes in aortic wall material properties, such as stiffness, have been shown to accompany onset and progression of various cardiovascular pathologies. Pulse Wave velocity (PWV) and propagation along the aortic wall have been shown to depend on the wall stiffness (i.e. stiffer the wall, higher the PWV), and can potentially enhance the noninvasive diagnostic techniques. Conventional clinical methods involve a global examination of the pulse traveling between femoral and carotid arteries, to provide an average PWV estimate. Such methods may not prove effective in detecting wall focal changes as entailed by a range of cardiovascular diseases. A two-way-coupled fluid-structure interaction (FSI) simulation study of pulse wave propagation along inhomogeneous aortas with focal stiffening and softening has previously proved the model reliable. In this study, simulations are performed on inhomogeneous aortic walls with hard inclusions of different numbers, size and modulus in order to further characterize the effects of focal hardening on pulse wave propagation. Spatio-temporal maps of the wall displacement were used to analyze the regional pulse wave propagations and velocities. The findings showed that the quantitative markers -such as PWVs and r² s on the pre-inclusion forward, reflected and post-inclusion waves, and the width of the standing wave- as well as qualitative markers -such as diffracted reflection zone versus single reflection wave- allow the successful and reliable distinction between the changes in inclusion numbers, size and modulus. Future studies are needed to incorporate the wall softening and physiologically-relevant wall inhomogeneities such as those seen in calcifications or aneurysms.

Blood flow vectoring control in aortic arch using full and partial clamps

BACKGROUND

Early diagnosis and treatment of aneurysm plays an important role in reducing the mortality risk of rupture. The aneurysm is a complex phenomenon and caused by different reasons, such as arteriosclerosis and heredity. In addition, pressure and Wall Shear Stress are two known factors influencing the establishment of an aneurysm. The aim of this study is to investigate the effect of using a full or partial clamp to control the blood flow streamlines and hence the location of stress concentration in a clean configuration of aorta. The main question is how to control the stresses distribution in order to reduce the possibility of aneurysm growth with less negative effects on the other sides.

METHODS AND RESULTS

A simple form of aortic arch with three branches is considered to simulate the effect of changing blood flow streamlines directions. A parameter study has been performed on the main characteristics of clamp, i.e. size, location, and the percentage of coverage. The Shear Stress Transport model is employed to simulate steady-state Newtonian blood flow when the Reynolds number is about 6500. Simulations are conducted using the commercial CFD solver ANSYS Fluent. The obtained results show that the location of clamp is more effective than the size. It is also found that increasing the depth of clamp has a negative impact on mean velocity field and hence on stress concentration.

CONCLUSION

The present results demonstrate that the Blood Flow Vectoring Control (BFVC) can change the main form of flow streamlines and consequently the distributions of pressure and Wall Shear Stress. A partial clamp leads to better results.