Three-dimensional geometry of the human carotid artery

Accurate characterization of carotid artery geometry is vital to our understanding of the pathogenesis of atherosclerosis. Three-dimensional computer reconstructions based on medical imaging are now ubiquitous; however, mean carotid artery geometry has not yet been comprehensively characterized. The goal of this work was to build and study such geometry based on data from 16 male patients with severe carotid artery disease. Results of computerized tomography angiography were used to analyze the cross-sectional images implementing a semiautomated segmentation algorithm. Extracted data were used to reconstruct the mean three-dimensional geometry and to determine average values and variability of bifurcation and planarity angles, diameters and cross-sectional areas. Contrary to simplified carotid geometry typically depicted and used, our mean artery was tortuous exhibiting nonplanarity and complex curvature and torsion variations. The bifurcation angle was 36 deg?±?11 deg if measured using arterial centerlines and 15 deg?±?14 deg if measured between the walls of the carotid bifurcation branches. The average planarity angle was 11 deg?±?10 deg. Both bifurcation and planarity angles were substantially smaller than values reported in most studies. Cross sections were elliptical, with an average ratio of semimajor to semiminor axes of 1.2. The cross-sectional area increased twofold in the bulb compared to the proximal common, but then decreased 1.5-fold for the combined area of distal internal and external carotid artery. Inter-patient variability was substantial, especially in the bulb region; however, some common geometrical features were observed in most patients. Obtained quantitative data on the mean carotid artery geometry and its variability among patients with severe carotid artery disease can be used by biomedical engineers and biomechanics vascular modelers in their studies of carotid pathophysiology, and by endovascular device and materials manufacturers interested in the mean geometrical features of the artery to target the broad patient population.     

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.     

Morphological age-dependent development of the human carotid bifurcation

The unique morphology of the adult human carotid bifurcation and its sinus has been investigated extensively, but its long-term, age-dependent development has not. It is important fundamentally and clinically to understand the hemodynamics and developmental forces that play a role in remodeling of the carotid bifurcation and maturation of the sinus in association with brain maturation. This understanding can lead to better prognostication and therapy of carotid disease. We analyzed the change of sinus morphology and the angle of the carotid bifurcation in four postnatal developmental stages (Group I: 0-2 years, Group II: 3-9 years, Group III: 10-19 years, and Group IV: 20-36 years, respectively) using multiprojection digital subtraction angiograms and image post-processing techniques. The most significant findings are the substantial growth of the internal carotid artery (ICA) with age and the development of a carotid sinus at the root of the ICA during late adolescence. The bifurcation angle remains virtually unchanged from infancy to adulthood. However, the angle split between the ICA and external carotid artery (ECA) relative to the common carotid artery (CCA) undergoes significant changes. Initially, the ICA appears to emanate as a side branch. Later in life, to reduce hydraulic resistance in response to increased flow demand by the brain, the bifurcation is remodeled to a construct in which both daughter vessels are a skewed continuation of the parent artery. This study provides a new analysis method to examine the development of the human carotid bifurcation over the developmental years, despite the small and sparse database. A larger database will enable in the future a more extensive analysis such as gender or racial differences.     

Quantification of the non-planarity of the human carotid bifurcation

The carotid bifurcation has been a region of particular interest due to its predilection for clinically significant atherosclerosis. It has been shown that the vessel geometry is a major determinant of the local haemodynamic properties which are believed to be associated with the location of atherosclerotic lesions. Current knowledge of the geometry of the carotid bifurcation is insufficient and restricted to basic geometric parameters. To provide some means of quantifying the degree of complexity of the 3D shape of the bifurcation, we made an initial attempt by evaluating the non-planarity of an arterial bifurcation based upon the singular value decomposition theorem.In this paper we present our results obtained on the right carotid bifurcations of six normal subjects, each of whom was scanned twice using the 2D time-of-flight MR sequence. The acquired 2D cross sectional images were processed by using our in-house software which comprises 2D segmentation, 3D reconstruction and smoothing. The centroids of each transverse slices were determined and used as input data for the non-planarity analysis. Our results using the singular value decomposition method have demonstrated discernible differences in non-planarity among individuals. Comparisons with the planarity definition proposed by other investigators suggest that the singular value decomposition method offers more information about the linearity and planarity of the bifurcation. However, it is also realised that a single measure of non-planarity can never fully characterise a bifurcation owing to the great variety of geometries.     

Effect of common carotid artery inlet length on normal carotid bifurcation hemodynamics

Controversy exists regarding the suitability of fully developed versus measured inlet velocity profiles for image-based computational fluid dynamics (CFD) studies of carotid bifurcation hemodynamics. Here, we attempt to resolve this by investigating the impact of the reconstructed common carotid artery (CCA) inlet length on computed metrics of "disturbed" flow. Twelve normal carotid bifurcation geometries were reconstructed from contrast-enhanced angiograms acquired as part of the Vascular Aging--The Link That Bridges Age to Atherosclerosis study (VALIDATE). The right carotid artery lumen geometry was reconstructed from its brachiocephalic origin to well above the bifurcation, and the CCA was truncated objectively at locations one, three, five, and seven diameters proximal to where it flares into the bifurcation. Relative to the simulations carried out using the full CCA, models truncated at one CCA diameter strongly overestimated the amount of disturbed flow. Substantial improvement was offered by using three CCA diameters, with only minor further improvement using five CCA diameters. With seven CCA diameters, the amounts of disturbed flow agreed unambiguously with those predicted by the corresponding full-length models. Based on these findings, we recommend that image-based CFD models of the carotid bifurcation should incorporate at least three diameters of CCA length if fully developed velocity profiles are to be imposed at the inlet. The need for imposing measured inlet velocity profiles would seem to be relevant only for those cases where the CCA is severely truncated.     

An experimental analysis of the flow field in a three-dimensional model of the human carotid artery bifurcation

Steady flow measurements were carried out in a rigid three-dimensional model of the human carotid artery bifurcation at a Reynolds number of 640 and a flow division ratio of 50/50. Both axial and secondary velocities were measured with a laser-Doppler anemometer. In the bulb opposite to the flow divider a zone with negative axial velocities was found with a maximal diameter of about 60% of the local diameter of the branch and a cross-sectional extent of about 25% of the local cross-sectional area. In the bulb the maximum axial velocity shifted towards the divider wall and at the end of the bulb an axial velocity plateau arose near the non-divider wall. Halfway through the bulb, secondary flow showed a vortex through which fluid flowed towards the divider wall near the bifurcation plane and back towards the non-divider wall near the upper walls.     

Use of factor analysis to characterize arterial geometry and predict hemodynamic risk: application to the human carotid bifurcation

The detailed geometry of atherosclerosis-prone vascular segments may influence their susceptibility by mediating local hemodynamics. An appreciation of the role of specific geometric variables is complicated by the considerable correlation among the many parameters that can be used to describe arterial shape and size. Factor analysis is a useful tool for identifying the essential features of such an inter-related data set, as well as for predicting hemodynamic risk in terms of these features and for interpreting the role of specific geometric variables. Here, factor analysis is applied to a set of 14 geometric variables obtained from magnetic resonance images of 50 human carotid bifurcations. Two factors alone were capable of predicting 12 hemodynamic metrics related to shear and near-wall residence time with adjusted squared Pearson's correlation coefficient as high as 0.54 and P-values less than 0.0001. One factor measures cross-sectional expansion at the bifurcation; the other measures the colinearity of the common and internal carotid artery axes at the bifurcation. The factors explain the apparent lack of an effect of branch angle on hemodynamic risk. The relative risk among the 50 bifurcations, based on time-average wall shear stress, could be predicted with a sensitivity and specificity as high as 0.84. The predictability of the hemodynamic metrics and relative risk is only modestly sensitive to assumptions about flow rates and flow partitions in the bifurcation.     

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.     

Pulsatile two-dimensional flow and plaque formation in a carotid artery bifurcation

Velocity and pressure fields, streamlines and wall shear stress distributions were numerically obtained for two-dimensional, steady and pulsatile flow in a carotid artery segment. Distinct regions of reverse flow near the bifurcation and wavy flow patterns in the branching channels were observed during portions of the pulse. These phenomena disappear at the end of the systolic phase of the cardiac cycle. A previously validated plaque formation model predicts that plaque sites and the local extent of atherosclerotic lesions are similar for those present on human angiograms.     

Stress-modulated collagen fiber remodeling in a human carotid bifurcation

This work concerns with the implementation of a new stress-driven remodeling model for simulating the overall structure and mechanical behavior of a human carotid bifurcation. By means of an iterative finite element based procedure collagen fiber direction and maximal principal stresses are computed. We find that the predicted fibers' architecture at the cylindrical branches and at the apex of the bifurcation correlates well with histological observations. Some insights about the mechanical response of the sinus bulb and the bifurcation apex are revealed and discussed. The results are compared with other, isotropic and orthotropic, models available in the literature.     

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.     

The geometry of unstented and stented pig common carotid artery bypass grafts

The long-term success of arterial bypass grafting with autologous saphenous veins is limited by neointimal hyperplasia (NIH), which seemingly develops preferentially at sites where hydrodynamic wall shear is low. Placement of a loose-fitting, porous stent around end-to-end, or end-to-side, autologous saphenous vein grafts on the porcine common carotid artery has been found significantly to reduce NIH, but the mechanism is unclear. In a preliminary study, we implanted autologous saphenous vein grafts bilaterally on the common carotid arteries of pigs, placing a stent around one graft and leaving the contralateral graft unstented. At sacrifice 1 month post implantation, the grafts were pressure fixed in situ and resin casts were made. Unstented graft geometry was highly irregular, with non-uniform dilatation, substantial axial lengthening, curvature, kinking, and possible long-pitch helical distortion. In contrast, stented grafts showed no major dilatation, lengthening or curvature, but there was commonly fine corrugation, occasional slight kinking or narrowing of segments, and possible long-pitch helical distortion. Axial growth of grafts against effectively tethered anastomoses could account for these changes. CFD studies are planned, using 3D MR reconstructions, on the effects of graft geometry on the flow. Abnormality of the flow could favour the development of vascular pathology, including NIH.     

On the influence of vessel planarity on local hemodynamics at the human carotid bifurcation

The human carotid artery bifurcation is a complex, three-dimensional structure exhibiting non-planarity and both in- and out-of-plane curvature. The aim of this study was to determine the relative importance of vessel planarity, a potential geometric risk factor for atherogenesis, in determining the local hemodynamics. A combination of computational fluid dynamics and magnetic resonance imaging was used to reconstruct the subject-specific hemodynamics for three subjects. Planar models were then constructed by translating the centroids of the lumen contours onto a plane defined by the centroids of the vessel branches near the bifurcation apex. A novel "patching" technique was used to convert the continuous arterial surfaces into contiguous but discrete patches according to an objective scheme, making it possible to compare the original and planar models without the need for registration and warping. Results suggest that the planarity of the vessel has a relatively minor effect on the spatial distribution of mean and oscillatory wall shear stress. Out-of-plane curvature was, however, found to have a marked influence on the extent and magnitude of these hemodynamic variables. We conclude that vessel curvature - whether in- or out-of-plane - rather than planarity may deserve further scrutiny as a potential geometric risk for atherogenesis.     

Three-dimensional numerical analysis of pulsatile flow and wall shear stress in the carotid artery bifurcation.

To analyse the pulsatile flow field and the mechanical stresses in a three-dimensional carotid artery bifurcation model, computer simulation is applied. The approximation of the Navier-Stokes equations uses a pressure correction finite element method. Numerical results are presented for axial and secondary flow velocity and wall shear stresses with special emphasis on the fluid dynamics in the carotid sinus. This region is of major interest because it is affected preferentially by lesions. Detailed local flow studies as carried out here should lead to a further insight into the mechanisms of atherogenesis. The flow conditions used in the study were chosen according to Ku et al. (Arteriosclerosis 5, 293-302, 1985). The results of this numerical analysis agree in the essential features with their experimental results.     

Pulsatile non-Newtonian flow characteristics in a three-dimensional human carotid bifurcation model.

Numerical analysis of flow phenomena and wall shear stresses in the human carotid artery bifurcation has been carried out using a three-dimensional geometrical model. The primary aim of this study is the detailed discussion of non-Newtonian flow velocity and wall shear stress during the pulse cycle. A comparison of non-Newtonian and Newtonian results is also presented. The applied non-Newtonian behavior of blood is based on measured dynamic viscosity. In the foreground of discussion are the flow characteristics in the carotid sinus. The investigation shows complex flow patterns especially in the carotid sinus where flow separation occurs at the outer wall throughout the systolic deceleration phase. The changing sign of the velocity near the outer sinus wall results in oscillating shear stress during the pulse cycle. At the outer wall of the sinus at maximum diameter level the shear stress ranges from -1.92 N/m2 to 1.22 N/m2 with a time-averaged value of 0.04 N/m2. At the inner wall of the sinus at maximum diameter level the shear stress range is from 1.16 N/m2 to 4.18 N/m2 with a mean of 1.97 N/m2. The comparison of non-Newtonian and Newtonian results indicates unchanged flow phenomena and rather minor differences in the basic flow characteristics.     

Hemodynamics in the carotid artery bifurcation: a comparison between numerical simulations and in vitro MRI measurements

The presence of atherosclerotic plaques has been shown to be closely related to the vessel geometry. Studies on postmortem human arteries and on the experimental animal show positive correlation between the presence of plaque thickness and low shear stress, departure of unidirectional flow and regions of flow separation and recirculation. Numerical simulations of arterial blood flow and direct blood flow velocity measurements by magnetic resonance imaging (MRI) are two approaches for the assessment of arterial blood flow patterns. In order to verify that both approaches give equivalent results magnetic resonance velocity data measured in a compliant anatomical carotid bifurcation model were compared to the results of numerical simulations performed for a corresponding computational vessel model. Cross sectional axial velocity profiles were calculated and measured for the midsinus and endsinus internal carotid artery. At both locations a skewed velocity profile with slow velocities at the outer vessel wall, medium velocities at the side walls and high velocities at the flow divider (inner) wall were observed. Qualitative comparison of the axial velocity patterns revealed no significant differences between simulations and in vitro measurements. Even quantitative differences such as for axial peak flow velocities were less than 10%. Secondary flow patterns revealed some minor differences concerning the form of the vortices but maximum circumferential velocities were in the same range for both methods.     

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.