Sex-related differences in coronary and carotid vessel geometry, plaque composition and shear stress obtained from imaging.

Atherosclerosis manifests itself differently in men and women with respect to plaque initiation, progression and plaque composition. The observed delay in plaque progression in women is thought to be related to the hormonal status of women. Also features associated with the vulnerability of plaques to rupture seem to be less frequently present in women compared to men. Current invasive and non-invasive imaging modalities allow for visualization of plaque size, composition and high risk vulnerable plaque features. Moreover, image based modeling gives access to local shear stress and shear stress-related plaque growth. In this review, current knowledge on sex-related differences in plaque size, composition, high risk plaque features and shear stress related plaque growth in carotid and coronary arteries obtained from imaging are summarized.

The Correlation Between Wall Shear Stress and Plaque Composition in Advanced Human Carotid Atherosclerosis.

The role of wall shear stress (WSS) in atherosclerotic plaque development is evident, but the relation between WSS and plaque composition in advanced atherosclerosis, potentially resulting in plaque destabilization, is a topic of discussion. Using our previously developed image registration pipeline, we investigated the relation between two WSS metrics, time-averaged WSS (TAWSS) and the oscillatory shear index (OSI), and the local histologically determined plaque composition in a set of advanced human carotid plaques. Our dataset of 11 carotid endarterectomy samples yielded 87 histological cross-sections, which yielded 511 radial bins for analysis. Both TAWSS and OSI values were subdivided into patient-specific low, mid, and high tertiles. This cross-sectional study shows that necrotic core (NC) size and macrophage area are significantly larger in areas exposed to high TAWSS or low OSI. Local TAWSS and OSI tertile values were generally inversely related, as described in the literature, but other combinations were also found. Investigating the relation between plaque vulnerability features and different combinations of TAWSS and OSI tertile values revealed a significantly larger cap thickness in areas exposed to both low TAWSS and low OSI. In conclusion, our study confirmed previous findings, correlating high TAWSS to larger macrophage areas and necrotic core sizes. In addition, our study demonstrated new relations, correlating low OSI to larger macrophage areas, and a combination of low TAWSS and low OSI to larger cap thickness.

Plaque Composition as a Predictor of Plaque Ulceration in Carotid Artery Atherosclerosis: The Plaque At RISK Study.

BACKGROUND AND PURPOSE: Plaque ulceration is a marker of previous plaque rupture. We studied the association between atherosclerotic plaque composition at baseline and plaque ulceration at baseline and follow-up. MATERIALS AND METHODS: We included symptomatic patients with a carotid stenosis of <70% who underwent MDCTA and MR imaging at baseline (n = 180). MDCTA was repeated at 2 years (n = 73). We assessed the presence of ulceration using MDCTA. Baseline MR imaging was used to assess the vessel wall volume and the presence and volume of plaque components (intraplaque hemorrhage, lipid-rich necrotic core, and calcifications) and the fibrous cap status. Associations at baseline were evaluated with binary logistic regression and reported with an OR and its 95% CI. Simple statistical testing was performed in the follow-up analysis. RESULTS: At baseline, the prevalence of plaque ulceration was 27% (49/180). Increased wall volume (OR = 12.1; 95% CI, 3.5-42.0), higher relative lipid-rich necrotic core (OR = 1.7; 95% CI, 1.3-2.2), higher relative intraplaque hemorrhage volume (OR = 1.7; 95% CI, 1.3-2.2), and a thin-or-ruptured fibrous cap (OR = 3.4; 95% CI, 1.7-6.7) were associated with the presence of ulcerations at baseline. In 8% (6/73) of the patients, a new ulcer developed. Plaques with a new ulceration at follow-up had at baseline a larger wall volume (1.04 cm3 [IQR, 0.97-1.16 cm3] versus 0.86 cm3 [IQR, 0.73-1.00 cm3]; P = .029), a larger relative lipid-rich necrotic core volume (23% [IQR, 13-31%] versus 2% [IQR, 0-14%]; P = .002), and a larger relative intraplaque hemorrhage volume (14% [IQR, 8-24%] versus 0% [IQR, 0-5%]; P < .001). CONCLUSIONS: Large atherosclerotic plaques and plaques with intraplaque hemorrhage and lipid-rich necrotic cores were associated with plaque ulcerations at baseline and follow-up.

An MRI-based method to register patient-specific wall shear stress data to histology.

Wall shear stress (WSS), the frictional force exerted on endothelial cells by blood flow, is hypothesised to influence atherosclerotic plaque growth and composition. We developed a methodology for image registration of MR and histology images of advanced human carotid plaques and corresponding WSS data, obtained by MRI and computational fluid dynamics. The image registration method requires four types of input images, in vivo MRI, ex vivo MRI, photographs of transversally sectioned plaque tissue and histology images. These images are transformed to a shared 3D image domain by applying a combination of rigid and non-rigid registration algorithms. Transformation matrices obtained from registration of these images are used to transform subject-specific WSS data to the shared 3D image domain as well. WSS values originating from the 3D WSS map are visualised in 2D on the corresponding lumen locations in the histological sections and divided into eight radial segments. In each radial segment, the correlation between WSS values and plaque composition based on histological parameters can be assessed. The registration method was successfully applied to two carotid endarterectomy specimens. The resulting matched contours from the imaging modalities had Hausdorff distances between 0.57 and 0.70 mm, which is in the order of magnitude of the in vivo MRI resolution. We simulated the effect of a mismatch in the rigid registration of imaging modalities on WSS results by relocating the WSS data with respect to the stack of histology images. A 0.6 mm relocation altered the mean WSS values projected on radial bins on average by 0.59 Pa, compared to the output of original registration. This mismatch of one image slice did not change the correlation between WSS and plaque thickness. In conclusion, we created a method to investigate correlations between WSS and plaque composition.

Flow Patterns in Carotid Webs: A Patient-Based Computational Fluid Dynamics Study.

BACKGROUND AND PURPOSE: Carotid webs are increasingly recognized as an important cause of (recurrent) ischemic stroke in patients without other cardiovascular risk factors. Hemodynamic flow patterns induced by these lesions might be associated with thrombus formation. The aim of our study was to evaluate flow patterns of carotid webs using computational fluid dynamics. MATERIALS AND METHODS: Patients with a carotid web in the Multicenter Randomized Clinical Trial of Endovascular Treatment of Acute Ischemic Stroke in the Netherlands (MR CLEAN) were selected for hemodynamic evaluation with computational fluid dynamics models based on lumen segmentations obtained from CT angiography scans. Hemodynamic parameters, including the area of recirculation zone, time-averaged wall shear stress, transverse wall shear stress, and the oscillatory shear index, were assessed and compared with the contralateral carotid bifurcation. RESULTS: In our study, 9 patients were evaluated. Distal to the carotid webs, recirculation zones were significantly larger compared with the contralateral bifurcation (63 versus 43 mm2, P = .02). In the recirculation zones of the carotid webs and the contralateral carotid bifurcation, time-averaged wall shear stress values were comparable (both: median, 0.27 Pa; P = .30), while transverse wall shear stress and oscillatory shear index values were significantly higher in the recirculation zone of carotid webs (median, 0.25 versus 0.21 Pa; P = .02 and 0.39 versus 0.30 Pa; P = .04). At the minimal lumen area, simulations showed a significantly higher time-averaged wall shear stress in the web compared with the contralateral bifurcation (median, 0.58 versus 0.45 Pa; P = .01). CONCLUSIONS: Carotid webs are associated with increased recirculation zones and regional increased wall shear stress metrics that are associated with disturbed flow. These findings suggest that a carotid web might stimulate thrombus formation, which increases the risk of acute ischemic stroke.

High shear stress influences plaque vulnerability Part of the data presented in this paper were published in Stroke 2007;38:2379-81.

Shear stress of the blood at the vessel wall plays an important role in many processes in the cardiovascular system primarily focused on the regulation of vessel lumen and wall dimensions. There is ample evidence that atherosclerotic plaques are generated at low shear stress regions in the cardiovascular system, while high shear stress regions are protected. In the course of plaque progression, advanced plaques start to encroach into the lumen, and thereby start to experience high shear stress at the endothelium. Until now the consequences of high shear stress working at the endothelium of an advanced plaque are unknown. As high shear stress influences tissue regression, we hypothesised that high shear stress can destabilise the plaque by cap weakening leading to ulceration. We investigated this hypothesis in a magnetic resonance imaging (MRI) dataset of a 67-year-old woman with a plaque in the carotid artery at baseline and an ulcer at ten-month follow-up. The lumen, plaque components (lipid/necrotic core, intraplaque haemorrhage) and ulcer were reconstructed three dimensionally and the geometry at baseline was used for shear stress calculation using computational fluid dynamics. Correlation of the change in plaque composition with the shear stress at baseline showed that the ulcer was generated exclusively at the high shear stress location. In this serial MRI study we found plaque ulceration at the high shear stress location of a protruding plaque in the carotid artery. Our data suggest that high shear stress influences plaque vulnerability and therefore may become a potential parameter for predicting future events. (Neth Heart J 2008;16:280-3.).

A numerical study on the influence of vulnerable plaque composition on intravascular thermography measurements.

Intracoronary thermography is a technique that measures lumen wall temperatures for vulnerable plaque detection. In this paper the influence of vulnerable plaque composition on lumen wall temperatures was studied numerically. Concerning the vulnerable plaque heat generation, the location of the heat source and its heat production were varied. Concerning the heat transfer, the thermal properties of the lipid core and the location of the vasa vasorum were studied. The heat source location was the main determinant of the lumen wall temperature distribution. The strongest effect was noted when the heat producing macrophages were located in the shoulder region leading to focal spots of higher temperature. The maximal lumen wall temperature was mainly determined by the heat production of the macrophages and the cooling effect of blood. The insulating properties of the lipid core increased lumen wall temperatures when the heat source was located in the cap and the presence of vasa vasorum lowered the temperatures. These results show that the lumen wall temperature distribution is influenced by vulnerable plaque composition and that intracoronary thermography techniques require a high spatial resolution. To be able to couple temperature measurements to plaque vulnerability, intracoronary thermography needs to be combined with an imaging modality.

High shear stress relates to intraplaque haemorrhage in asymptomatic carotid plaques.

BACKGROUND AND AIMS: Carotid artery plaques with vulnerable plaque components are related to a higher risk of cerebrovascular accidents. It is unknown which factors drive vulnerable plaque development. Shear stress, the frictional force of blood at the vessel wall, is known to influence plaque formation. We evaluated the association between shear stress and plaque components (intraplaque haemorrhage (IPH), lipid rich necrotic core (LRNC) and/or calcifications) in relatively small carotid artery plaques in asymptomatic persons. METHODS: Participants (n = 74) from the population-based Rotterdam Study, all with carotid atherosclerosis assessed on ultrasound, underwent carotid MRI. Multiple MRI sequences were used to evaluate the presence of IPH, LRNC and/or calcifications in plaques in the carotid arteries. Images were automatically segmented for lumen and outer wall to obtain a 3D reconstruction of the carotid bifurcation. These reconstructions were used to calculate minimum, mean and maximum shear stresses by applying computational fluid dynamics with subject-specific inflow conditions. Associations between shear stress measures and plaque composition were studied using generalized estimating equations analysis, adjusting for age, sex and carotid wall thickness. RESULTS: The study group consisted of 93 atherosclerotic carotid arteries of 74 participants. In plaques with higher maximum shear stresses, IPH was more often present (OR per unit increase in maximum shear stress (log transformed) = 12.14; p = 0.001). Higher maximum shear stress was also significantly associated with the presence of calcifications (OR = 4.28; p = 0.015). CONCLUSIONS: Higher maximum shear stress is associated with intraplaque haemorrhage and calcifications.

Relationship between neointimal thickness and shear stress after Wallstent implantation in human coronary arteries.

BACKGROUND: In-stent restenosis by excessive intimal hyperplasia reduces the long-term clinical efficacy of coronary stents. Because shear stress (SS) is related to plaque growth in atherosclerosis, we investigated whether variations in SS distribution are related to variations in neointima formation. METHODS AND RESULTS: In 14 patients, at 6-month follow-up after coronary Wallstent implantation, 3D stent and vessel reconstruction was performed with a combined angiographic and intravascular ultrasound technique (ANGUS). The bare stent reconstruction was used to calculate in-stent SS at implantation, applying computational fluid dynamics. The flow was selected to deliver an average SS of 1.5 N/m(2). SS and neointimal thickness (Th) values were obtained with a resolution of 90 degrees in the circumferential and 2.5 mm in the longitudinal direction. For each vessel, the relationship between Th and SS was obtained by linear regression analysis. Averaging the individual slopes and intercepts of the regression lines summarized the overall relationship. Average Th was 0.44+/-0.20 mm. Th was inversely related to SS: Th=(0.59+/-0.24)-(0.08+/-0.10)xSS (mm) (P<0.05). CONCLUSIONS: These data show for the first time in vivo that the Th variations in Wallstents at 6-month follow-up are inversely related to the relative SS distribution. These findings support a hemodynamic mechanism underlying in-stent neointimal hyperplasia formation.

True 3-dimensional reconstruction of coronary arteries in patients by fusion of angiography and IVUS (ANGUS) and its quantitative validation.

BACKGROUND: True 3D reconstruction of coronary arteries in patients based on intravascular ultrasound (IVUS) may be achieved by fusing angiographic and IVUS information (ANGUS). The clinical applicability of ANGUS was tested, and its accuracy was evaluated quantitatively. METHODS AND REUSLTS: In 16 patients who were investigated 6 months after stent implantation, a sheath-based catheter was used to acquire IVUS images during an R-wave-triggered, motorized stepped pullback. First, a single set of end-diastolic biplane angiographic images documented the 3D location of the catheter at the beginning of pullback. From this set, the 3D pullback trajectory was predicted. Second, contours of the lumen or stent obtained from IVUS were fused with the 3D trajectory. Third, the angular rotation of the reconstruction was optimized by quantitative matching of the silhouettes of the 3D reconstruction with the actual biplane images. Reconstructions were obtained in 12 patients. The number of pullback steps, which determines the pullback length, closely agreed with the reconstructed path length (r=0.99). Geometric measurements in silhouette images of the 3D reconstructions showed high correlation (0.84 to 0.97) with corresponding measurements in the actual biplane angiographic images. CONCLUSIONS: With ANGUS, 3D reconstructions of coronary arteries can be successfully and accurately obtained in the majority of patients.

Shear-stress and wall-stress regulation of vascular remodeling after balloon angioplasty: effect of matrix metalloproteinase inhibition.

BACKGROUND: Constrictive vascular remodeling (VR) is the most significant component of restenosis after balloon angioplasty (PTA). Whereas in physiological conditions VR is associated with normalization of shear stress (SS) and wall stress (WS), after PTA the role of SS and WS in VR is unknown. Furthermore, whereas matrix metalloproteinase inhibition (MMPI) has been shown to modulate VR after PTA, its effect on the SS and WS control mechanisms after PTA is unknown. METHODS AND RESULTS: PTA was performed in external iliac arteries of 12 atherosclerotic Yucatan pigs, of which 6 pigs (7 vessels) received the MMPI batimastat and 6 pigs (10 vessels) served as controls. Before and after the intervention and at 6-week follow-up, intravascular ultrasound pullback was performed, allowing 3D reconstruction of the treated segment and computational fluid dynamics to calculate the media-bounded area and SS. WS was derived from the Laplace formula. Immediately after PTA, media-bounded area, WS, and SS changed by 20%, 16%, and -49%, respectively, in both groups. VR was predicted by SS and WS. In the control group, SS and WS had been normalized at follow-up with respect to the reference segment. In contrast, for the batimastat group, the SS had been normalized, but not the WS. The latter is attributed to an increase in wall area at follow-up. CONCLUSIONS: Vascular remodeling after PTA is controlled by both SS and WS. MMPI inhibited the WS control system.

The role of shear stress in the destabilization of vulnerable plaques and related therapeutic implications.

American Heart Association type IV plaques consist of a lipid core covered by a fibrous cap, and develop at locations of eccentric low shear stress. Vascular remodeling initially preserves the lumen diameter while maintaining the low shear stress conditions that encourage plaque growth. When these plaques eventually start to intrude into the lumen, the shear stress in the area surrounding the plaque changes substantially, increasing tensile stress at the plaque shoulders and exacerbating fissuring and thrombosis. Local biologic effects induced by high shear stress can destabilize the cap, particularly on its upstream side, and turn it into a rupture-prone, vulnerable plaque. Tensile stress is the ultimate mechanical factor that precipitates rupture and atherothrombotic complications. The shear-stress-oriented view of plaque rupture has important therapeutic implications. In this review, we discuss the varying mechanobiologic mechanisms in the areas surrounding the plaque that might explain the otherwise paradoxical observations and unexpected outcomes of experimental therapies.

The role of shear stress in the generation of rupture-prone vulnerable plaques.

Blood-flow-induced shear stress acting on the arterial wall is of paramount importance in vascular biology. Endothelial cells sense shear stress and largely control its value in a feedback-control loop by adapting the arterial dimensions to blood flow. Nevertheless, to allow for variations in arterial geometry, such as bifurcations, shear stress control is modified at certain eccentrically located sites to let it remain at near-zero levels. In the presence of risk factors for atherosclerosis, low shear stress contributes to local endothelial dysfunction and eccentric plaque build up, but normal-to-high shear stress is atheroprotective. Initially, lumen narrowing is prevented by outward vessel remodeling. Maintenance of a normal lumen and, by consequence, a normal shear stress distribution, however, prolongs local unfavorable low shear stress conditions and aggravates eccentric plaque growth. While undergoing such growth, eccentric plaques at preserved lumen locations experience increased tensile stress at their shoulders making them prone to fissuring and thrombosis. Consequent loss of the plaque-free wall by coverage with thrombus and new tissue may bring shear-stress-controlled lumen preservation to an end. This change causes shear stress to increase, which as a new condition may transform the lesion into a rupture-prone vulnerable plaque. We present a discussion of the role of shear stress, in setting the stage for the generation of rupture-prone, vulnerable plaques, and how this may be prevented.

Carotid plaque elasticity estimation using ultrasound elastography, MRI, and inverse FEA - A numerical feasibility study.

The material properties of atherosclerotic plaques govern the biomechanical environment, which is associated with rupture-risk. We investigated the feasibility of noninvasively estimating carotid plaque component material properties through simulating ultrasound (US) elastography and in vivo magnetic resonance imaging (MRI), and solving the inverse problem with finite element analysis. 2D plaque models were derived from endarterectomy specimens of nine patients. Nonlinear neo-Hookean models (tissue elasticity C1) were assigned to fibrous intima, wall (i.e., media/adventitia), and lipid-rich necrotic core. Finite element analysis was used to simulate clinical cross-sectional US strain imaging. Computer-simulated, single-slice in vivo MR images were segmented by two MR readers. We investigated multiple scenarios for plaque model elasticity, and consistently found clear separations between estimated tissue elasticity values. The intima C1 (160 kPa scenario) was estimated as 125.8 ± 19.4 kPa (reader 1) and 128.9 ± 24.8 kPa (reader 2). The lipid-rich necrotic core C1 (5 kPa) was estimated as 5.6 ± 2.0 kPa (reader 1) and 8.5 ± 4.5 kPa (reader 2). A scenario with a stiffer wall yielded similar results, while realistic US strain noise and rotating the models had little influence, thus demonstrating robustness of the procedure. The promising findings of this computer-simulation study stimulate applying the proposed methodology in a clinical setting.

Temperature distribution in atherosclerotic coronary arteries: influence of plaque geometry and flow (a numerical study).

Intravascular coronary thermography is a method that may detect vulnerable, atherosclerotic plaques and is currently evaluated in a clinical setting. Active macrophages or enzymatic heat releasing processes in vulnerable plaques may act as heat sources. To better understand the parameters of influence on thermographic measurements, numerical simulations have been performed on a model of a coronary artery segment containing a heat source. Heat source parameters and flow were varied to study their influence on temperatures at the lumen wall. Maximal temperature differences at the lumen wall increased when the source volume increased and they differ with the source geometry. The simulations showed that blood flow acts as a coolant to the lumen wall. Blood flow decreased maximal temperatures depending on the source geometry, source volume and the maximal flow velocity. Influence of flow was highest for circumferentially extended sources, up to a factor 3.7, and lowest for longitudinally extended sources, down to a factor 1.9. When cap thickness increased, maximal temperatures decreased and the influence of flow increased. This study shows that correct interpretation of intravascular thermographic measurements requires data on the flow and on the morphologic characteristics of the atherosclerotic plaque.

Coronary stent implantation changes 3-D vessel geometry and 3-D shear stress distribution.

Mechanisms of in-stent restenosis are not fully understood. Shear stress is known to play a role in plaque and thrombus formation and is sensitive to changes in regional vessel geometry. Hence, we evaluated the regional changes in 3-D geometry and shear stress induced by stent placement in coronary arteries of pigs.Methods. 3-D reconstruction was performed, applying a combined angiographic and IVUS technique (ANGUS), from seven Wallstents (diameter 3.5 (n=3) and 5mm (n=4)), which were implanted in seven coronary arteries of five pigs. This 3-D geometry was used to calculate locally the curvature, while the shear stress distribution was obtained by computational fluid dynamics. Local changes in shear stress were obtained at the entrance and exit of the stent for baseline (0. 65+/-0.22 ml/s) and hyperemic flow (2.60+/-0.86 ml/s) conditions. Results. After stent implantation, the curvature increased by 121% at the entrance and by 100% at the exit of the stent, resulting in local changes in shear stress. In general, at the entrance of the stent local maxima in shear stress were generated, while at the exit both local maxima and minima in shear stress were observed (p<0.05). Additionally, the shear stress at the entrance and exit of the stent were correlated with the local curvature (r: 0.30-0.84).Conclusion. Stent implantation changes 3-D vessel geometry in such a way that regions with decreased and increased shear stress occur close to the stent edges. These changes might be related to the asymmetric patterns of in-stent restenosis.

On the IVUS plaque volume error in coronary arteries when neglecting curvature.

Plaque volume determined by common linear 3-D IVUS analysis systems will show under- or overestimation in curved vessel segments because these systems approximate the true 3-D transducer pull-back trajectory by a straight line. We developed a mathematical model that showed that the error is primarily dependent on the curvature of the pull-back trajectory and not on vessel tortuosity. Furthermore, we measured this error in vivo in the coronary arteries of 15 patients, comparing the plaque volume using a true 3-D reconstruction method with that of the linear approach. The in vivo plaque volume error ranged from 2.3% to -1.2% for 15 coronary segments with lengths ranging from 38.8 to 89.1 mm (62.2 +/- 13 mm). The volume error introduced by linear 3-D IVUS analysis systems is dependent on the curvature of the pull-back trajectory. The error measured in vivo was small and inversely related to segment length.

Effect of catheter placement on 3-D velocity profiles in curved tubes resembling the human coronary system.

Novel measurement techniques based on intravenous ultrasound (IVUS) technology ('IVUS-Flowmetry') require the location of a catheter inside the coronary bed. The present study quantifies disturbances in the 3-D velocity profile induced by catheter placement inside a tube, applying computational fluid dynamics. Two curved, circular meshes (radius K = 0.025 m and K = 0.035 m) with and without a catheter inside the lumen were applied. The catheter was located at the inner curve, the outer curve and at the top position. Boundary conditions were: no slip on the wall, zero stress at the outlet, uniform inflow with entrance velocities of 0.1, 0.2 and 0.4 m/s. Curvature-associated centrifugal forces shifted the maximal velocity to the outer curve and introduced two symmetrical vortices. Additional catheter placement redistributed the 3-D axial velocity field away from the catheter, which was accompanied by the appearance of multiple low-strength vortices. In addition, peak axial velocity increased, peak secondary velocities decreased, axial pressure drop increased and shear stress increased. Flow calculations simulated to resemble IVUS-based flowmetry changed by only 1% after considering secondary velocity. In conclusion, placement of a catheter inside a curved tube resembling the human coronary system changes the velocity field and reduces secondary patterns. The present study supports the usefulness of catheter-based flowmetry during resting flow conditions. During hyperemic flow conditions, flow measurements might be accompanied by large axial pressure drops because the catheter, itself, might act as a significant stenosis.

Geometry-based pressure drop prediction in mildly diseased human coronary arteries.

Pressure drop (△p) estimations in human coronary arteries have several important applications, including determination of appropriate boundary conditions for CFD and estimation of fractional flow reserve (FFR). In this study a △p prediction was made based on geometrical features derived from patient-specific imaging data. Twenty-two mildly diseased human coronary arteries were imaged with computed tomography and intravascular ultrasound. Each artery was modelled in three consecutive steps: from straight to tapered, to stenosed, to curved model. CFD was performed to compute the additional △p in each model under steady flow for a wide range of Reynolds numbers. The correlations between the added geometrical complexity and additional △p were used to compute a predicted △p. This predicted △p based on geometry was compared to CFD results. The mean △p calculated with CFD was 855±666Pa. Tapering and curvature added significantly to the total △p, accounting for 31.4±19.0% and 18.0±10.9% respectively at Re=250. Using tapering angle, maximum area stenosis and angularity of the centerline, we were able to generate a good estimate for the predicted △p with a low mean but high standard deviation: average error of 41.1±287.8Pa at Re=250. Furthermore, the predicted △p was used to accurately estimate FFR (r=0.93). The effect of the geometric features was determined and the pressure drop in mildly diseased human coronary arteries was predicted quickly based solely on geometry. This pressure drop estimation could serve as a boundary condition in CFD to model the impact of distal epicardial vessels.

The influence of axial image resolution on atherosclerotic plaque stress computations.

Biomechanical models are used extensively to study risk factors, such as peak stresses, for vulnerable atherosclerotic plaque rupture. Typically, 3D patient-specific arterial models are reconstructed by interpolating between cross sectional contour data which have a certain axial sampling, or image, resolution. The influence of the axial sampling resolution on computed stresses, as well as the comparison of 3D with 2D simulations, is quantified in this study. A set of histological data of four atherosclerotic human coronary arteries was used which were reconstructed in 3D with a high sampling (HS) and low sampling (LS) axial resolution, and 4 slices were treated separately for 2D simulations. Stresses were calculated using finite element analysis (FEA). High stresses were found in thin cap regions and regions of thin vessel walls, low stresses were found inside the necrotic cores and media and adventitia layers. Axial sampling resolution was found to have a minor effect on general stress distributions, peak plaque/cap stress locations and the relationship between peak cap stress and minimum cap thickness. Axial sampling resolution did have a profound influence on the error in computed magnitude of peak plaque/cap stresses (±15.5% for HS vs. LS geometries and ±24.0% for HS vs. 2D geometries for cap stresses). The findings of this study show that axial under sampling does not influence the qualitative stress distribution significantly but that high axially sampled 3D models are needed when accurate computation of peak stress magnitudes is required.