Flow pattern analysis in a highly stenotic patient-specific carotid bifurcation model using a turbulence model

The aim of this study is to investigate the blood flow pattern in carotid bifurcation with a high degree of luminal stenosis, combining in vivo magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). A newly developed two-equation transitional model was employed to evaluate wall shear stress (WSS) distribution and pressure drop across the stenosis, which are closely related to plaque vulnerability. A patient with an 80% left carotid stenosis was imaged using high resolution MRI, from which a patient-specific geometry was reconstructed and flow boundary conditions were acquired for CFD simulation. A transitional model was implemented to investigate the flow velocity and WSS distribution in the patient-specific model. The peak time-averaged WSS value of approximately 73 Pa was predicted by the transitional flow model, and the regions of high WSS occurred at the throat of the stenosis. High oscillatory shear index values up to 0.50 were present in a helical flow pattern from the outer wall of the internal carotid artery immediately after the throat. This study shows the potential suitability of a transitional turbulent flow model in capturing the flow phenomena in severely stenosed carotid arteries using patient-specific MRI data and provides the basis for further investigation of the links between haemodynamic variables and plaque vulnerability. It may be useful in the future for risk assessment of patients with carotid disease.     

Steady flow and wall compression in stenotic arteries: a three-dimensional thick-wall model with fluid-wall interactions.

Severe stenosis may cause critical flow and wall mechanical conditions related to artery fatigue, artery compression, and plaque rupture, which leads directly to heart attack and stroke. The exact mechanism involved is not well understood. In this paper a nonlinear three-dimensional thick-wall model with fluid-wall interactions is introduced to simulate blood flow in carotid arteries with stenosis and to quantify physiological conditions under which wall compression or even collapse may occur. The mechanical properties of the tube wall were selected to match a thick-wall stenosis model made of PVA hydrogel. The experimentally measured nonlinear stress-strain relationship is implemented in the computational model using an incremental linear elasticity approach. The Navier-Stokes equations are used for the fluid model. An incremental boundary iteration method is used to handle the fluid-wall interactions. Our results indicate that severe stenosis causes considerable compressive stress in the tube wall and critical flow conditions such as negative pressure, high shear stress, and flow separation which may be related to artery compression, plaque cap rupture, platelet activation, and thrombus formation. The stress distribution has a very localized pattern and both maximum tensile stress (five times higher than normal average stress) and maximum compressive stress occur inside the stenotic section. Wall deformation, flow rates, and true severities of the stenosis under different pressure conditions are calculated and compared with experimental measurements and reasonable agreement is found.     

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

Hemodynamics factors and biomechanical forces play key roles in atherogenesis, plaque development and final rupture. In this paper, we investigated the flow field and stress field for different degrees of stenoses under physiological conditions. Disease is modelled as axisymmetric cosine shape stenoses with varying diameter reductions of 30%, 50% and 70%, respectively. A simulation model which incorporates fluid-structure interaction, a turbulence model and realistic boundary conditions has been developed. The results show that wall motion is constrained at the throat by 60% for the 30% stenosis and 85% for the 50% stenosis; while for the 70% stenosis, wall motion at the throat is negligible through the whole cycle. Peak velocity at the throat varies from 1.47 m/s in the 30% stenosis to 3.2m/s in the 70% stenosis against a value of 0.78 m/s in healthy arteries. Peak wall shear stress values greater than 100 Pa were found for > or =50% stenoses, which in vivo could lead to endothelial stripping. Maximum circumferential stress was found at the shoulders of plaques. The results from this investigation suggest that severe stenoses inhibit wall motion, resulting in higher blood velocities and higher peak wall shear stress, and localization of hoop stress. These factors may contribute to further development and rupture of plaques.     

A nonlinear axisymmetric model with fluid-wall interactions for steady viscous flow in stenotic elastic tubes.

Arteries with high-grade stenoses may compress under physiologic conditions due to negative transmural pressure caused by high-velocity flow passing through the stenoses. To quantify the compressive conditions near the stenosis, a nonlinear axisymmetric model with fluid-wall interactions is introduced to simulate the viscous flow in a compliant stenotic tube. The nonlinear elastic properties of the tube (tube law) are measured experimentally and used in the model. The model is solved using ADINA (Automatic Dynamic Incremental Nonlinear Analysis), which is a finite element package capable of solving problems with fluid-structure interactions. Our results indicate that severe stenoses cause critical flow conditions such as negative pressure and high and low shear stresses, which may be related to artery compression, plaque cap rupture, platelet activation, and thrombus formation. The pressure filed near a stenosis has a complex pattern not seen in one-dimensional models. Negative transmural pressure as low as -24 mmHg for a 78 percent stenosis by diameter is observed at the throat of the stenosis for a downstream pressure of 30 mmHg. Maximum shear stress as a high as 1860 dyn/cm2 occurs at the throat of the stenoses, while low shear stress with reversed direction is observed right distal to the stenosis. Compressive stresses are observed inside the tube wall. The maximal principal stress and hoop stress in the 78 percent stenosis are 80 percent higher than that from the 50 percent stenosis used in our simulation. Flow rates under different pressure drop conditions are calculated and compared with experimental measurements and reasonable agreement is found for the prebuckling stage.     

Numerical modeling of pulsatile turbulent flow in stenotic vessels

Pulsatile turbulent flow in stenotic vessels has been numerically modeled using the Reynolds-averaged Navier-Stokes equation approach. The commercially available computational fluid dynamics code (CFD), FLUENT, has been used for these studies. Two different experiments were modeled involving pulsatile flow through axisymmetric stenoses. Four different turbulence models were employed to study their influence on the results. It was found that the low Reynolds number k-omega turbulence model was in much better agreement with previous experimental measurements than both the low and high Reynolds number versions of the RNG (renormalization-group theory) k-epsilon turbulence model and the standard k-epsilon model, with regard to predicting the mean flow distal to the stenosis including aspects of the vortex shedding process and the turbulent flow field. All models predicted a wall shear stress peak at the throat of the stenosis with minimum values observed distal to the stenosis where flow separation occurred.     

Numerical simulation of blood pulsatile flow in a stenosed carotid artery using different rheological models

Symmetrical 30-60% stenosis in a common carotid artery under unsteady flow condition for Newtonian and six non-Newtonian viscosity models are investigated numerically. Results show power-law model produces higher deviations, in terms of velocity and wall shear stress in comparison with other models while generalized power-law and modified-Casson models are more prone to Newtonian state. Comparing separation length of recirculation region at different critical points of cardiac cycle confirms the necessity of considering blood flow in unsteady mode. Increasing stenosis intensity causes flow patterns more disturbed downstream of the stenosis and WSS appear to develop remarkably at the stenosis throat.     

In Vitro Shear Stress Measurements Using Particle Image Velocimetry in a Family of Carotid Artery Models: Effect of Stenosis Severity,Plaque Eccentricity,and Ulceration

Atherosclerotic disease, and the subsequent complications of thrombosis and plaque rupture, has been associated with local shear stress. In the diseased carotid artery, local variations in shear stress are induced by various geometrical features of the stenotic plaque. Greater stenosis severity, plaque eccentricity (symmetry) and plaque ulceration have been associated with increased risk of cerebrovascular events based on clinical trial studies. Using particle image velocimetry, the levels and patterns of shear stress (derived from both laminar and turbulent phases) were studied for a family of eight matched-geometry models incorporating independently varied plaque features – i.e. stenosis severity up to 70%, one of two forms of plaque eccentricity, and the presence of plaque ulceration). The level of laminar (ensemble-averaged) shear stress increased with increasing stenosis severity resulting in 2–16 Pa for free shear stress (FSS) and approximately double (4–36 Pa) for wall shear stress (WSS). Independent of stenosis severity, marked differences were found in the distribution and extent of shear stress between the concentric and eccentric plaque formations. The maximum WSS, found at the apex of the stenosis, decayed significantly steeper along the outer wall of an eccentric model compared to the concentric counterpart, with a 70% eccentric stenosis having 249% steeper decay coinciding with the large outer-wall recirculation zone. The presence of ulceration (in a 50% eccentric plaque) resulted in both elevated FSS and WSS levels that were sustained longer (∼20 ms) through the systolic phase compared to the non-ulcerated counterpart model, among other notable differences. Reynolds (turbulent) shear stress, elevated around the point of distal jet detachment, became prominent during the systolic deceleration phase and was widely distributed over the large recirculation zone in the eccentric stenoses.     

Enhanced External Counterpulsation Treatment May Intervene The Advanced Atherosclerotic Plaque Progression by Inducing The Variations of Mechanical Factors: A 3D FSI Study Based on in vivo Animal Experiment

Growing evidences suggest that long-term enhanced external counterpulsation (EECP) treatment can inhibit the initiation of atherosclerotic lesion by improving the hemodynamic environment in aortas. However, whether this kind procedure will intervene the progression of advanced atherosclerotic plaque remains elusive and causes great concern in its clinical application presently. In the current paper, a pilot study combining animal experiment and numerical simulation was conducted to investigate the acute mechanical stress variations during EECP intervention, and then to assess the possible chronic effects. An experimentally induced hypercholesterolemic porcine model was developed and the basic hemodynamic measurement was performed in vivo before and during EECP treatment. Meanwhile, A 3D fluid-structure interaction (FSI) model of blood vessel with symmetric local stenosis was developed for the numerical calculation of some important mechanical factors. The results show that EECP augmented 12.21% of the plaque wall stress (PWS), 57.72% of the time average wall shear stress (AWSS) and 43.67% of the non-dimensional wall shear stress gradient (WSSGnd) at throat site of the stenosis. We suggest that long-term EECP treatment may intervene the advanced plaque progression by inducing the significant variations of some important mechanical factors, but its proper effects will need a further research combined follow-up observation in clinic.     

Mechanical action of the blood onto atheromatous plaques: influence of the stenosis shape and morphology

The vulnerability of atheromatous plaques in the carotid artery may be related to several factors, the most important being the degree of severity of the endoluminal stenosis and the thickness of the fibrous cap. It has recently been shown that the plaque length can also affect the mechanical response significantly. However, in their study on the effect of the plaque length, the authors did not consider the variations of the plaque morphology and the shape irregularities that may exist independently of the plaque length. These aspects are developed in this paper. The mechanical interactions between the blood flow and an atheromatous plaque are studied through a numerical model considering fluid–structure interaction. The simulation is achieved using the arbitrary Lagrangian–Eulerian scheme in the COMSOL TM commercial finite element package. The stenosis severity and the plaque length are, respectively, set to 45% and 15 mm. Different shapes of the stenosis are modelled, considering irregularities made of several bumps over the plaque. The resulting flow patterns, wall shear stresses, plaque deformations and stresses in the fibrous cap reveal that the effects of the blood flow are amplified if the slope upstream stenosis is steep or if the plaque morphology is irregular with bumps. More specifically, the maximum stress in the fibrous cap is 50% larger for a steep slope than for a gentle slope. These results offer new perspectives for considering the shape of plaques in the evaluation of the vulnerability.     

Large deforming buoyant embolus passing through a stenotic common carotid artery: a computational simulation

Arterial embolism is responsible for the death of lots of people who suffers from heart diseases. The major risk of embolism in upper limbs is that the ruptured particles are brought into the brain, thus stimulating neurological symptoms or causing the stroke. We presented a computational model using fluid-structure interactions (FSI) to investigate the physical motion of a blood clot inside the human common carotid artery. We simulated transportation of a buoyant embolus in an unsteady flow within a finite length tube having stenosis. Effects of stenosis severity and embolus size on arterial hemodynamics were investigated. To fulfill realistic nonlinear property of a blood clot, a rubber/foam model was used. The arbitrary Lagrangian-Eulerian formulation (ALE) and adaptive mesh method were used inside fluid domain to capture the large structural interfacial movements. The problem was solved by simultaneous solution of the fluid and the structure equations. Stress distribution and deformation of the clot were analyzed and hence, the regions of the embolus prone to lysis were localized. The maximum magnitude of arterial wall shear stress during embolism occurred at a short distance proximal to the throat of the stenosis. Through embolism, arterial maximum wall shear stress is more sensitive to stenosis severity than the embolus size whereas role of embolus size is more significant than the effect of stenosis severity on spatial and temporal gradients of wall shear stress downstream of the stenosis and on probability of clot lysis due to clot stresses while passing through the stenosis.     

Experimental and numerical study of pulsatile flows through stenosis: wall shear stress analysis.

Different shapes of pulsatile flows through a model of stenosis are experimentally and numerically modeled to validate both methods and to determine the wall shear stress temporal evolution downstream from the stenosis. Two-dimensional velocity measurements are performed in a 75% severity stenosis using a pulsed Doppler ultrasonic velocimeter. Finite element package is employed for the transient numerical simulations. Polynomial method, based on the experimental velocity values, is proposed to determine the wall shear stress temporal evolution. There is a good agreement between the numerical and experimental results. The wall shear stress temporal analysis shows oscillating wall shear stress values during the cycle with high wall shear stress values at the throat of about 120 dyn/cm2, and low values downstream from the stenosis of about - 2.5 dyn/cm2. The key result of the study is that the presence of the stenosis leads the artery to work in a direction which is opposite to the direction of a healthy artery.     

Direct numerical simulation of transitional flow in a stenosed carotid bifurcation

The blood flow dynamics of a stenosed, subject-specific, carotid bifurcation were numerically simulated using the spectral element method. Pulsatile inlet conditions were based on in vivo color Doppler ultrasound measurements of blood velocity. The results demonstrated the transitional or weakly turbulent state of the blood flow, which featured rapid velocity and pressure fluctuations in the post-stenotic region of the internal carotid artery (ICA) during systole and laminar flow during diastole. High-frequency vortex shedding was greatest downstream of the stenosis during the deceleration phase of systole. Velocity fluctuations had a frequency within the audible range of 100-300Hz. Instantaneous wall shear stress (WSS) within the stenosis was relatively high during systole ( approximately 25-45Pa) compared to that in a healthy carotid. In addition, high spatial gradients of WSS were present due to flow separation on the inner wall. Oscillatory flow reversal and low pressure were observed distal to the stenosis in the ICA. This study predicts the complex flow field, the turbulence levels and the distribution of the biomechanical stresses present in vivo within a stenosed carotid artery.     

Effect of geometrical assumptions on numerical modeling of coronary blood flow under normal and disease conditions

Shear stress plays a pivotal role in pathogenesis of coronary heart disease. The spatial and temporal variation in hemodynamics of blood flow, especially shear stress, is dominated by the vessel geometry. The goal of the present study was to investigate the effect of 2D and 3D geometries on the numerical modeling of coronary blood flow and shear stress distribution. We developed physiologically realistic 2D and 3D models (with similar geometries) of the human left coronary artery under normal and stenosis conditions (30%, 60%, and 80%) using PROE (WF 3). Transient blood flows in these models were solved using laminar and turbulent (k-ω) models using a computational fluid dynamics solver, FLUENT (v6.3.26). As the stenosis severity increased, both models predicted a similar pattern of increased shear stress at the stenosis throat, and in recirculation zones formed downstream of the stenosis. The 2D model estimated a peak shear stress value of 0.91, 2.58, 5.21, and 10.09 Pa at the throat location under normal, 30%, 60%, and 80% stenosis severity. The peak shear stress values at the same location estimated by the 3D model were 1.41, 2.56, 3.15, and 13.31 Pa, respectively. The 2D model underestimated the shear stress distribution inside the recirculation zone compared with that of 3D model. The shear stress estimation between the models diverged as the stenosis severity increased. Hence, the 2D model could be sufficient for analyzing coronary blood flow under normal conditions, but under disease conditions (especially 80% stenosis) the 3D model was more suitable.     

Wall shear stress and near-wall flows in the stenosed femoral artery

Stenotic artery hemodynamics are often characertised by metrics including oscillatory shear index (OSI) and residence time (RT). This analysis was conducted to clarify the link between the near-wall flow behaviour and these resultant flow metrics. A computational simulation was conducted of a stenosed femoral artery, with an idealised representative geometry and a physiologically realistic inlet profile. The overall flow behaviour was characterised through consideration of the axial flow, which was non-dimensionalised against mean flow velocity. The OSI and RT metrics, which are a useful indicator of likely atherosclerotic sites, were explained through a discussion of the WSS values at different time points, the velocity behaviour and velocity profiles, with a particular focus on the near-wall behaviour which influences wall shear stress and the transient evolution of the wall shear stress. While, the stenosis throat experiences high values of wall shear stress, the smooth flow through this contracted region results in low variation in wall shear stress vectors and limited opportunity for any particle stasis. However, regions were noted distal and proximal (though to a lesser extent), where the change in recirculation zones over the cycle created highly elevated regions of both OSI and RT.     

Long time evolution of atherosclerotic plaques

The evolution of atherosclerosis in general, and the influence of wall shear stress on the growth of atherosclerotic plaques in particular, is an intricate phenomenon which is still only partly understood. We therefore propose a qualitative mathematical model which consists of a number of ordinary differential equations for the concentrations of the most relevant constituents of the atherosclerotic plaque. These equations were studied both for the case that the wall shear stress is a parameter (model A), and for the case in which the plaque evolution is coupled to the blood flow (model B) which results in a time dependent wall shear stress. We find that both models exhibit a class of marginally stable equilibria, all reflecting states in which the plaque only grows for a short period of time after a perturbation. The uncoupled model A, however, shows bi-stability between this class of equilibria and another equilibrium state in which the plaque experiences unlimited growth in time, if the LDL cholesterol intake exceeds a threshold value. In model B the bi-stability vanishes, but we find that there is still a critical value of the LDL cholesterol intake beyond which the lumen radius drastically decreases. We show that this decrease is quite sensitive to the value of the wall shear stress.     

Laser Doppler anemometer measurements of pulsatile flow in a model carotid bifurcation

Hemodynamics at the human carotid bifurcation is important to the understanding of atherosclerotic plaque initiation and progression as well as to the diagnosis of clinically important disease. Laser Doppler anemometry was performed in a large scale model of an average human carotid. Pulsatile waveforms and physiologic flow divisions were incorporated. Disturbance levels and shear stresses were computed from ensemble averages of the velocity waveform measurements. Flow in the common carotid was laminar and symmetric. Flow patterns in the sinus, however, were complex and varied considerably during the cycle. Strong helical patterns and outer wall flow separation waxed and waned during each systole. The changing flow patterns resulted in an oscillatory shear stress at the outer wall ranging from -13 to 9 dyn cm-2 during systole with a time-averaged mean of only -0.5 dyn cm-2. This contrasts markedly with an inner wall shear stress range of 17-50, (mean 26) dyn cm-2. The region of transient separation was confined to the carotid sinus outer wall with no reverse velocities detected in the distal internal carotid. Notable disturbance velocities were also time-dependent, occurring only during the deceleration phase of systole and the beginning of diastole. The present pulsatile flow studies have aided in identifying hemodynamic conditions which correlate with early intimal thickening and predict the physiologic level of flow disturbances in the bulb of undiseased internal carotid arteries.     

Measured wall shear stress distribution pattern upstream and downstream of a unilateral stenosis by an electrochemical method

An electrochemical surface shear stress measurement was applied to a model of unilateral arterial stenosis. The unilateral stenosis model was made up of a removable stenosis plug, in an electrochemical shear stress measurement test section with 100 cathodes. Three dimensional wall shear stress distribution was measured under steady flow field. At a relatively low Reynolds number, Re = 270, there was a characteristic high and low wall shear distribution pattern downstream of the unilateral stenosis. There were also remarkable high shear stress areas on the opposite wall up- and downstream, and both side walls of the stenosis upstream. It was clearly shown that detailed three dimensional structure of the flow field must be studied in order to correlate it to pathological findings.     

A negative correlation between human carotid atherosclerotic plaque progression and plaque wall stress: in vivo MRI-based 2D/3D FSI models

It is well accepted that atherosclerosis initiation and progression correlate positively with low and oscillating flow wall shear stresses (FSS). However, this mechanism cannot explain why advanced plaques continue to grow under elevated FSS conditions. In vivo magnetic resonance imaging (MRI)-based 2D/3D multi-component models with fluid-structure interactions (FSI, 3D only) for human carotid atherosclerotic plaques were introduced to quantify correlations between plaque progression measured by wall thickness increase (WTI) and plaque wall (structure) stress (PWS) conditions. A histologically validated multi-contrast MRI protocol was used to acquire multi-year in vivo MRI images. Our results using 2D models (200-700 data points/patient) indicated that 18 out of 21 patients studied showed significant negative correlation between WTI and PWS at time 2 (T2). The 95% confidence interval for the Pearson correlation coefficient is (-0.443,-0.246), p<0.0001. Our 3D FSI model supported the 2D correlation results and further indicated that combining both plaque structure stress and flow shear stress gave better approximation results (PWS, T2: R(2)=0.279; FSS, T1: R(2)=0.276; combining both: R(2)=0.637). These pilot studies suggest that both lower PWS and lower FSS may contribute to continued plaque progression and should be taken into consideration in future investigations of diseases related to atherosclerosis.     

To what extent does a minimal atherosclerotic plaque alter the arterial wall shear stress distribution? A model study by an electrochemical method

An electrochemical surface shear stress measurement was applied to a model of very thin unilateral arterial stenosis (height of 1/8 of the model pipe diameter with very smooth surface). Three dimensional wall shear stress distribution was measured under steady flow field from a relatively low Reynolds number, Re = 270, to a high Reynolds number, Re = 1200. There was a characteristic high and low wall shear distribution pattern around the stenosis. There were also remarkable high shear stress areas on the opposite wall and both side walls of the stenosis. It was clearly shown that three dimensional structure of the flow field, hence, the wall shear stress distribution, is affected by a minimal change on the arterial wall.     

Hemodynamics analysis of a stenosed carotid bifurcation and its plaque-mitigating design

Considering transient two-dimensional laminar flow in a diseased carotid artery segment with realistic inlet and outflow conditions, detailed velocity profiles, pressure fields, wall shear stress distributions and coupled, localized plaque formations have been simulated. The type of outflow boundary condition influences to a certain degree the extent of plaque build-up, which in turn reduces "disturbed flow" phenomena such as flow separations, recirculation zones, and wavy flow patterns in the artery branches during portions of the pulse. Based on computer experiments varying key geometric factors, a plaque-mitigating design of a carotid artery bifurcation has been proposed. Elimination of the carotid bulb, a smaller bifurcation angle, lower area ratios, and smooth wall curvatures generated a design with favorable hemodynamics parameters, leading to reduced plaque build-up by factors of 10 and 2 in the internal carotid and in the external carotid, respectively.