Combined effects of pulsatile flow and dynamic curvature on wall shear stress in a coronary artery bifurcation model

A three-dimensional model with simplified geometry for the branched coronary artery is presented. The bifurcation is defined by an analytical intersection of two cylindrical tubes lying on a sphere that represents an idealized heart surface. The model takes into account the repetitive variation of curvature and motion to which the vessel is subject during each cardiac cycle, and also includes the phase difference between arterial motion and blood flowrate, which may be nonzero for patients with pathologies such as aortic regurgitation. An arbitrary Lagrangian Eulerian (ALE) formulation of the unsteady, incompressible, three-dimensional Navier-Stokes equations is employed to solve for the flow field, and numerical simulations are performed using the spectral/hp element method. The results indicate that the combined effect of pulsatile inflow and dynamic geometry depends strongly on the aforementioned phase difference. Specifically, the main findings of this work show that the time-variation of flowrate ratio between the two branches is minimal (less than 5%) for the simulation with phase difference angle equal to 90 degrees, and maximal (51%) for 270 degrees. In two flow pulsatile simulation cases for fixed geometry and dynamic geometry with phase angle 270 degrees, there is a local minimum of the normalized wall shear rate amplitude in the vicinity of the bifurcation, while in other simulations a local maximum is observed.     

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.     

Inter-individual variations in wall shear stress and mechanical stress distributions at the carotid artery bifurcation of healthy humans

Fluid shear stress and mechanical wall stress may play a role in the formation of early atherosclerotic lesions, but these quantities are difficult to measure in vivo. Our objective was to quantify these parameters in normal subjects in a clinical setting, and to define regions of low wall shear stress and high mechanical stress. The right carotid bifurcations of five healthy male volunteers were investigated using a novel non-invasive technique which integrates magnetic resonance angiography, ultrasonography, tonometry and state-of-the-art computational fluid dynamics and solid mechanics models. Significant inter-subject variations in patterns as well as magnitude of wall shear stress and mechanical stress were found. In spite of individual variabilities, this study revealed that some regions of the artery wall are exposed simultaneously to low wall shear stress and high mechanical stress and that these regions correspond to areas where atherosclerotic plaque develops. The coexistence of regions of low wall shear stress and high tensile stress may be an important determinant of the formation of atheroma in human arteries.     

Wall shear stress distribution in the human carotid siphon during pulsatile flow

Wall shear stress distribution in the carotid siphon, which is a multiple curved segment of the internal carotid artery, is investigated numerically under physiological flow conditions. The computer simulation of flow through the model segment is based on the time-dependent, three-dimensional Navier-Stokes equations, solved numerically with a finite element method. The study shows the behavior of the wall shear stress-vector field and identifies the zones of high and low wall shear stress values during the cardiac cycle.     

A numerical analysis of steady flow in a three-dimensional model of the carotid artery bifurcation

A finite element approximation of steady flow in a rigid three-dimensional model of the carotid artery bifurcation is presented. A Reynolds number of 640 and a flow division ratio of about 50/50, simulating systolic flow, was used. To limit the CPU- and I/O-times needed for solving the systems of equations, a mesh-generator was developed, which gives full control over the number of elements into which the bifurcation is divided. A mini-supercomputer, based on parallel and vector processing techniques, was used to solve the system of equations. The numerical results of axial and secondary flow compare favorably with those obtained from previously performed laser-Doppler velocity measurements. Also, the influence of the Reynolds number, the flow division ratio, and the bifurcation angle on axial and secondary flow in the carotid sinus were studied in the three-dimensional model. The influence of the interventions is limited to a relatively small variation in the region with reversed axial flow, more or less pronounced C-shaped axial velocity contours, and increasing or decreasing axial velocity maxima.     

Vortex generation in pulsatile flow through arterial bifurcation models including the human carotid artery

Visualization experiments were performed to elucidate the complicated flow pattern in pulsatile flow through arterial bifurcations. Human common carotid arteries, which were made transparent, and glass-models simulating Y- and T-shaped bifurcations were used. Pulsatile flow with wave forms similar to those of arterial flow was generated with a piston pump, elastic tube, airchamber, and valves controlling the outflow resistance. Helically recirculating flow with a pattern similar to that of the horseshoe vortex produced around wall-based protuberances in circular tubes was observed in pulsatile flow through all the bifurcations used in the present study. This flow type, which we shall refer to as the horseshoe vortex, has also been demonstrated to occur at the human common carotid bifurcation in steady flow with Reynolds numbers above 100. Time-varying flows also produced the horseshoe vortex mostly during the decelerating phase. Fluid particles of dye solution approaching the bifurcation apex diverged, divided into two directions perpendicularly, and then showed helical motion representing the horseshoe vortex formation. While this helical flow was produced, the stagnation points appeared on the wall upstream of the apex. Their position was dependent upon the flow distribution ratio between the branches in the individual arteries. The region affected by the horseshoe vortex was smaller during pulsatile flow than during steady flow. Lowering the Reynolds number together with the Womersley number weakened the intensity of helical flow. A separation bubble, resulting from the divergence or wall roughness, was observed at the outer or inner wall of the branch vessels and made the flow more complicated.     

Post-stenotic core flow behavior in pulsatile flow and its effects on wall shear stress

Arteries of several species, including man, tend to adjust their diameters such that the mean wall shear stress is in the range of 10-20 dynes cm-2. Additionally, intimal thickening in the human carotid bifurcation correlates well with the reciprocal of wall shear stress as determined in model studies. The correlation indicates that wherever the local mean wall shear stress exceeds approximately 10 dynes cm-2, the artery tends to be spared from intimal thickening. However, it is not known whether mean shear stress, i.e. the time-averaged value, or the instantaneous shear stress is the appropriate correlative variable. Each of these variables suggests different mechanisms for the reaction of the artery wall to its hemodynamic environment. It is therefore important to devise means by which the effects of mean shear and pulsatile shear can be separated in the study of atherogenesis. The present investigation examines the post-stenotic flow field in Plexiglas models under pulsatile conditions approximating those in the aortas of the cynomolgus monkey, an animal often employed in atherogenesis research. Behavior of the core flow and its effects on wall shear stress are studied for stenoses of 75 and 90% area reductions using laser velocimetry. The results show that the post-stenotic field contains regions in which the mean wall shear stress is low, but the pulsatile excursions are large.(ABSTRACT TRUNCATED AT 250 WORDS)     

A numerical and experimental analysis of the flow field in a two-dimensional model of the human carotid artery bifurcation

Axial velocities were measured in an enlarged, two-dimensional, rigid model of the carotid artery bifurcation by means of a laser-Doppler anemometer, under both steady and unsteady flow conditions. Also a numerical model was developed, based on the finite element approximation of the Navier-Stokes and continuity equations. From this study it appeared that the numerically predicted velocities agree well with the experimentally obtained values. Besides, the bifurcation hardly influenced the upstream flow in the main branch (common carotid artery), high velocity gradients were observed at the divider walls of the daughter branches (internal and external carotid arteries) and large zones with reversed flow were present near the nondivider walls of these branches. For steady flow the maximal diameter of this zone at the entrance of the internal carotid artery (carotid sinus) was about 25% of the local diameter of this branch. For unsteady flow this zone was absent during the initial phase of flow acceleration and maximal at the end of flow deceleration with a maximal diameter of about 50% of the local diameter of the carotid sinus.     

颈动脉血管壁切应力的分析

动脉中管壁的脉动低切应力在动脉粥样硬化形成中起始动和主要的决定作用。本文比较了几种计算血管壁切应力的方法,认为采用有约束的弹性管模型计算获得的动脉壁切应力更适合于临床应用。根据检测得到的正常人和动脉硬化性脑血管病患者的颈动脉血流速度、血管管径等数据,计算两者的颈动脉壁面切应力。研究发现动脉硬化性脑血管病患者的壁面切应力比正常人显著减小。这表明,颈动脉的壁面切应力可以作为动脉硬化性脑血管疾病的早期诊断的重要参考指标。     

Numerical analysis of flow through a severely stenotic carotid artery bifurcation

The results of computational simulations may supplement MR and other in vivo diagnostic techniques to provide an accurate picture of the hemodynamics in particular vessels, which may help demonstrate the risks of embolism or plaque rupture posed by particular plaque deposits. In this study, a model based on an endarterectomy specimen of the plaque in a carotid bifurcation was examined. The flow conditions include steady flow at Reynolds numbers of 300, 600, and 900 as well as unsteady pulsatile flow. Both dynamic pressure and wall shear stress are very high, with shear values up to 70 N/m2, proximal to the stenosis throat in the internal carotid artery, and both vary significantly through the flow cycle. The wall shear stress gradient is also strong along the throat. Vortex shedding is observed downstream of the most severe occlusion. Two turbulence models, the Chien and Goldberg varieties of k-epsilon, are tested and evaluated for their relevance in this geometry. The Chien model better captures phenomena such as vortex shedding. The flow distal to stenosis is likely transitional, so a model that captures both laminar and turbulent behavior is needed.     

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.     

MRI and CFD studies of pulsatile flow in healthy and stenosed carotid bifurcation models

Pulsatile flow was studied in physiologically realistic models of a normal and a moderately stenosed (30% diameter reduction) human carotid bifurcation. Time-resolved velocity measurements were made using magnetic resonance imaging, from which wall shear stress (WSS) vectors were calculated. Velocity measurements in the inflow and outflow regions were also used as boundary conditions for a computational fluid dynamics (CFD) model. Experimental flow patterns and derived WSS vectors were compared qualitatively with the corresponding CFD predictions. In the stenosed phantom, flow in the bulb region of the "internal carotid artery" was concentrated along the outer wall, with a region of low and recirculating flow near the inner wall. In the normal phantom, the converse was found, with a low flow region near the outer wall of the bulb. Time-averaged WSS and oscillatory shear index were also markedly different for the two phantoms.     

A two-dimensional numerical analysis of unsteady flow in the carotid artery bifurcation. A comparison with three-dimensional in-vitro measurements and the influence of minor stenoses

In the present study a two-dimensional finite element model for incompressible Newtonian flow is applicated to the modelling of carotid artery flow. In earlier studies, the numerical model was validated experimentally for several flow configurations. In general the pulsatile flow is characterized by reversed flow regions at the non-divider side walls of both the internal and external carotid arteries. The unsteadiness of the flow is associated with rather complex spatial and temporal velocity distributions and leads to temporal variations of the location and length of the reversed flow regions. As a consequence, pronounced spatial and temporal variations in the wall shear stresses are found. At the non-divider side walls, wall shear stresses are relatively low and exhibits an oscillatory behaviour in space and time. At the divider side walls, wall shear stresses are relatively high and approximately follow the flow rate distribution in time. The aim of this study is not only to present two-dimensional calculations but also to compare the calculated two-dimensional velocity profiles with those from three-dimensional experiments. It is observed that in the common carotid artery and in the proximal parts of the internal and external carotid arteries, the two-dimensional numerical model provides valuable information with respect to the three-dimensional configuration. In the more distal parts of especially the internal carotid artery, deviations are found between the two-dimensional numerical and three-dimensional experimental model. These deviations can mainly be attributed to the neglect of the secondary velocity distribution in the two-dimensional model. In the two-dimensional numerical model the influence of a minor stenosis in the internal carotid artery is hardly distinguishable from a minor geometrical variation without stenosis. Full three-dimensional analyses of the influence of minor stenoses are needed to prove numerically whether in-vivo measurements of the axial velocity distribution are useful in the detection of minor stenoses.     

The frequency response of electrochemical wall shear probes in pulsatile flow

The frequency response of surface-mounted electrochemical mass transfer probes used to deduce wall shear rates has been investigated experimentally for the case of fully developed laminar pulsatile flow in a straight tube. Generally good agreement is found with the asymptotic results obtained by Lighthill's methods. The significance of the results with regard to the investigation of models of pulsatile flows of physiological interest is discussed. It is concluded that the frequency-dependent phase and amplitude corrections required to obtain accurate wall shear measurements are of such magnitudes as to render impractical the use of electrochemical probes to determine wall shear rates in these flows.     

Turbulence detection in a stenosed artery bifurcation by numerical simulation of pulsatile blood flow using the low-Reynolds number turbulence model

The pulsatile blood flow in a partially blocked artery is significantly altered as the flow regime changes through the cardiac cycle. This paper reports on the application of a low-Reynolds turbulence model for computation of physiological pulsatile flow in a healthy and stenosed carotid artery bifurcation. The human carotid artery was chosen since it has received much attention because atherosclerotic lesions are frequently observed. The Wilcox low-Re k-omega turbulence model was used for the simulation since it has proven to be more accurate in describing transition from laminar to turbulent flow. Using the FIDAP finite element code a validation showed very good agreement between experimental and numerical results for a steady laminar to turbulent flow transition as reported in a previous publication by the same authors. Since no experimental or numerical results were available in the literature for a pulsatile and turbulent flow regime, a comparison between laminar and low-Re turbulent calculations was made to further validate the turbulence model. The results of this study showed a very good agreement for velocity profiles and wall shear stress values for this imposed pulsatile laminar flow regime. To explore further the medical aspect, the calculations showed that even in a healthy or non-stenosed artery, small instabilities could be found at least for a portion of the pulse cycle and in different sections. The 40% and 55% diameter reduction stenoses did not significantly change the turbulence characteristics. Further results showed that the presence of 75% stenoses changed the flow properties from laminar to turbulent flow for a good portion of the cardiac pulse. A full 3D simulation with this low-Re-turbulence model, coupled with Doppler ultrasound, can play a significant role in assessing the degree of stenosis for cardiac patients with mild conditions.     

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.     

Wall shear in pulsatile flow

This paper deals with a mathematical attempt to determine the wall shear during normal flow of blood in the ascending and the descending thoracic aorta. A simple model is used, but the results obtained are in agreement with published experimental results for the descending thoracic aorta. It is suggested that the degree of fluctuation in the pressure gradient at a given station is the major factor in determining the level of wall shear at that point.