Is flow in the common carotid artery fully developed?

Development and Evaluation of a Multifrequency Ultrafast Doppler Spectral Analysis (MFUDSA) Algorithm for Wall Shear Stress Measurement: A Simulation and In Vitro Study

Cardiovascular pathology is the leading cause of death and disability in the Western world, and current diagnostic testing usually evaluates the anatomy of the vessel to determine if the vessel contains blockages and plaques. However, there is a growing school of thought that other measures, such as wall shear stress, provide more useful information for earlier diagnosis and prediction of atherosclerotic related disease compared to pulsed-wave Doppler ultrasound, magnetic resonance angiography, or computed tomography angiography. A novel algorithm for quantifying wall shear stress (WSS) in atherosclerotic plaque using diagnostic ultrasound imaging, called Multifrequency ultrafast Doppler spectral analysis (MFUDSA), is presented. The development of this algorithm is presented, in addition to its optimisation using simulation studies and in-vitro experiments with flow phantoms approximating the early stages of cardiovascular disease. The presented algorithm is compared with commonly used WSS assessment methods, such as standard PW Doppler, Ultrafast Doppler, and Parabolic Doppler, as well as plane-wave Doppler. Compared to an equivalent processing architecture with one-dimensional Fourier analysis, the MFUDSA algorithm provided an increase in signal-to-noise ratio (SNR) by a factor of 4-8 and an increase in velocity resolution by a factor of 1.10-1.35. The results indicated that MFUDSA outperformed the others, with significant differences detected between the typical WSS values of moderate disease progression (p = 0.003) and severe disease progression (p = 0.001). The algorithm demonstrated an improved performance for the assessment of WSS and has potential to provide an earlier diagnosis of cardiovascular disease than current techniques allow.

Differences in Carotid Artery Geometry and Flow Caused by Body Postural Changes and Physical Exercise

Different body postures and physical exercises may lead to changes in arterial geometry and hemodynamics, which may be associated with the distribution of atherosclerosis lesions. This study was aimed at investigating potential geometric and hemodynamic changes of the carotid bifurcation in different body postures and after high-intensity interval training (HIIT) workouts. Three-dimensional vascular ultrasound (3DVUS) and Doppler ultrasound images were acquired for 21 healthy participants (aged 29 ± 6 y, 14 men and 7 women) in different body postures (sitting and three sleeping postures [supine, left lateral and right lateral]) and after physical exercises. The common carotid artery (CCA) and internal carotid artery (ICA) diameters of the left carotid artery were found to increase significantly from supine to left lateral (both p <0.05). CCA diameters (p < 0.05) and ICA/CCA diameter ratio (p < 0.01) of the left carotid artery changed significantly from supine to sitting. Significant differences in CCA peak systolic velocity (CCA PSV, p < 0.001), CCA end-diastolic velocity (CCA EDV, p < 0.001), CCA pulsatility index (CCA PI, p < 0.001) and maximum velocity-based wall shear stress at the CCA (WSS_(max) at the CCA, p < 0.001) were identified in different postures. After physical exercises, significant increases were observed in the CCA diameter (p < 0.001), CCA PSV (p < 0.001), ICA PSV (p < 0.05), WSS_(max) at the CCA (p < 0.001) and WSS_(max) at the ICA (p < 0.05), as were significantly lower values of the CCA EDV (p < 0.01) and ICA/CCA PSV ratio (p < 0.05). Side-to-side differences were also detected in different postural change scenarios and after physical exercise; more significant differences were found to occur only in the left-sided carotid artery. Significant differences were identified under postural change and after physical exercise among healthy adults, suggesting that daily activity has an effect on the carotid bifurcation. These changes may be associated with formation and development of carotid atherosclerosis. Moreover, these side differences might be severe for patients and worth further attention in clinical practice.

Carotid Doppler ultrasound for non-invasive haemodynamic monitoring: a narrative review

Objective.Accurate haemodynamic monitoring is the cornerstone in the management of critically ill patients. It guides the optimization of tissue and organ perfusion in order to prevent multiple organ failure. In the past decades, carotid Doppler ultrasound (CDU) has been explored as a non-invasive alternative for long-established invasive haemodynamic monitoring techniques. Considering the large heterogeneity in reported studies, we conducted a review of the literature to clarify the current status of CDU as a haemodynamic monitoring tool.Approach.In this article, firstly an overview is given of the equipment and workflow required to perform a CDU exam in clinical practice, the limitations and technical challenges potentially faced by the CDU sonographer, and the cerebrovascular mechanisms that may influence CDU measurement outcomes. The following chapter describes alternative techniques for non-invasive haemodynamic monitoring, detailing advantages and limitations compared to CDU. Next, a comprehensive review of the literature regarding the use of CDU for haemodynamic monitoring is presented. Furthermore, feasibility aspects, training requirements and technical developments of CDU are addressed.Main results.Based on the outcomes of these studies, we assess the applicability of CDU-derived parameters within three clinical domains (cardiac output, volume status, and fluid responsiveness), and amongst different patient groups. Finally, recommendations are provided to improve the quality and standardization of future research and clinical practice in this field.Significance.Although CDU is not yet interchangeable with invasive 'gold standard' cardiac output monitoring, the present work shows that certain CDU-derived parameters prove promising in the context of functional haemodynamic monitoring.

Adaptive Wall Shear Stress Imaging in Phantoms, Simulations and In Vivo

WSS measurement is challenging since it requires sensitive flow measurements at a distance close to the wall. The aim of this study is to develop an ultrasound imaging technique which combines vector flow imaging with an unsupervised data clustering approach that automatically detects the region close to the wall with optimally linear flow profile, to provide direct and robust WSS estimation. The proposed technique was evaluated in phantoms, mimicking normal and atherosclerotic vessels, and spatially registered Fluid Structure Interaction (FSI) simulations. A relative error of 6.7% and 19.8% was obtained for peak systolic (WSSPS) and end diastolic (WSSED) WSS in the straight phantom, while in the stenotic phantom, a good similarity was found between measured and simulated WSS distribution, with a correlation coefficient, R, of 0.89 and 0.85 for WSSPS and WSSED, respectively. Moreover, the feasibility of the technique to detect pre-clinical atherosclerosis was tested in an atherosclerotic swine model. Six swines were fed atherogenic diet, while their left carotid artery was ligated in order to disturb flow patterns. Ligated arterial segments that were exposed to low WSSPS and WSS characterized by high frequency oscillations at baseline, developed either moderately or highly stenotic plaques (p < 0.05). Finally, feasibility of the technique was demonstrated in normal and atherosclerotic human subjects. Atherosclerotic carotid arteries with low stenosis had lower WSSPS as compared to control subjects (p < 0.01), while in one subject with high stenosis, elevated WSS was found on an arterial segment, which coincided with plaque rupture site, as determined through histological examination.

Quantitative Blood Flow Measurements in the Common Carotid Artery: A Comparative Study of High-Frame-Rate Ultrasound Vector Flow Imaging, Pulsed Wave Doppler, and Phase Contrast Magnetic Resonance Imaging

V Flow is commercially developed by high-frame-rate ultrasound vector flow imaging. Compared to conventional color Doppler, V Flow is angle-independent and is capable of measuring both the magnitude and the direction of blood flow velocities. This paper aims to investigate the differences between V Flow and pulsed wave Doppler (PW) relative to phase contrast magnetic resonance imaging (PC-MRI), for the quantitative measurements of blood flow in common carotid arteries (CCA) and, consequently, to evaluate the accuracy of the new technique, V Flow. Sixty-four CCAs were measured using V Flow, PW, and PC-MRI. The maximum velocities, time-averaged mean (TAMEAN) velocities, and volume flow were measured using different imaging technologies. The mean error with standard deviation (Std), the median of absolute errors, and the r-values between V Flow and PC-MRI results for the maximum velocity, the TAMEAN velocity, and the volume flow measurements are {9.40% ± 14.91%; 11.84%; 0.84}, {21.52% ± 14.46%; 19.28%; 0.86}, and {−2.80% ± 14.01%; 10.38%; 0.7}, respectively, and are {53.44% ± 29.68%; 49.79%; 0.74}, {27.83% ± 31.60%; 23.83; 0.71}, and {21.01% ± 29.64%; 25.48%; 0.34}, respectively, for those between PW and PC-MRI. The boxplot, linear regression and Bland–Altman plots were performed for each comparison, which illustrated that the results measured via V Flow rather than via PW agreed more closely with those measured via PC-MRI.

Real-Time High Frame Rate Color Flow Mapping System

Plane wave (PW) transmission (TX) can be profitably used to improve the performance of color flow mapping (CFM) systems by increasing the autocorrelation ensemble length (EL) and/or the frame rate (FR). Although high-end scanners tend to include imaging schemes using PW TX and parallel receive beams, high frame rate (HFR) CFM has been so far experimentally implemented mostly through research platforms that transmit PWs and beamform/process the received channel data off-line. In this article, full real-time implementation of PW CFM with continuous-time clutter filtering and extended FR/EL is reported. The field-programmable gate arrays (FPGAs) and digital signal processors (DSPs) onboard the ULA-OP 256 research scanner were programmed to perform high-speed parallel beamforming and autocorrelation-based CFM processing, respectively. Different strategies were tested, in which the TX of PWs for CFM is either continuous or interleaved with the TX of packets of B-mode pulses. A fourth-order Chebyshev continuous-time high-pass filter with programmable cutoff frequency was implemented and its clutter rejection performance was positively compared with that obtained when operating on packet data. CFM FRs up to 575 were obtained. The possibility of programming the autocorrelation EL up to 64 permitted to detect flow with high sensitivity and accuracy (average relative errors down to 0.4% ± 8.4%). In vivo HFR movies are presented, showing the dynamics of flow in the common carotid artery, which highlight the presence of secondary flow components.

Flow Complexity Estimation in Dysfunctional Arteriovenous Dialysis Fistulas using Vector Flow Imaging

Non-invasive assessment is preferred for monitoring arteriovenous dialysis fistulas (AVFs). Vector concentration assesses flow complexity, which may correlate with stenosis severity. We determined whether vector concentration could assess stenosis severity in dysfunctional AVFs. Vector concentration was estimated in four stenotic phantoms at different pulse repetition frequencies. Spectral Doppler peak velocity and vector concentration were measured in 12 patients with dysfunctional AVFs. Additionally, 5 patients underwent digital subtraction angiography (DSA). In phantoms, vector concentration exhibited an inverse relationship with stenosis severity and was less affected by aliasing in severe stenoses. In nine stenoses of 5 patients undergoing DSA, vector concentration correlated strongly with stenosis severity (first stenosis: r = -0.73, p = 0.04; other stenoses; r = -0.69, p = 0.02) and mid-stenotic diameter (first stenosis: r = 0.87, p = 0.006; other stenoses: r = 0.70, p = 0.02) as opposed to peak velocities (p > 0.05). Vector concentration is less affected by aliasing in severe stenoses and correlates with DSA in patients with dysfunctional AVF.

Wall Shear Stress Measurements Based on Ultrasound Vector Flow Imaging: Theoretical Studies and Clinical Examples

Wall shear stress (WSS) is considered as a key factor for atherosclerosis development. Previous WSS research based on pulsed wave Doppler (PWD) showed limitations in complex flows. To improve accuracy for nonlaminar flow, a commercial ultrasound vector flow imaging (UVFI)-based WSS calculation is proposed. Errors for PWD are presented theoretically when flow is not laminar. Based on this, simulations of WSS calculations between PWD and UVFI were set up for different turbulent flows. Our simulations show that UVFI has obviously better performance than PWD in WSS calculations. Wall shear stress results in different flow conditions at carotid bifurcations are described.

Effects of carotid pressure waveforms on the results of wave separation, wave intensity and reservoir pressure analysis

OBJECTIVE

Recently great attention has been paid to innovative cardiovascular biomarkers obtained from wave separation (WS), wave intensity (WI) and reservoir-wave (RW) theories. All these approaches share a requirement for pressure information. The aim of this study was to evaluate differences in WS-, WI- and RW-derived parameters obtained achieving pressure waveforms in different ways.

APPROACH

Twenty-two individuals (49  ±  17 years, 59% males) were examined. Common carotid blood flow waveforms were obtained from pulsed-wave Doppler images. Carotid pressure waveforms were achieved in four different ways: (1) with applanation tonometry, used as a reference method; (2) linear scaling from an ultrasound (US)-derived diameter curve; (3) exponential scaling from a US-derived diameter curve; and (4) linear scaling from an accelerometric-derived diameter signal. For each case, the reflection magnitude (RM) and index (RI) were obtained from the WS. The amplitude of the first positive peak (W ₁), of the second positive peak (W ₂) and of the negative peak (W _b) were calculated from the WI, while the maximum of the reservoir (maxP_r) and the excess (maxP_ex) pressure were achieved from the RW.

MAIN RESULTS

According to the intra-class coefficient values, the agreement between the standard method and all the others was excellent for the RM (linear: 0.82; exponential: 0.83; accelerometric: 0.86), RI (linear: 0.84; exponential: 0.85; accelerometric: 0.87), maxP_r (linear: 0.97; exponential: 0.96; accelerometric: 0.97) and maxP_ex (linear: 0.85; exponential: 0.87; accelerometric: 0.89), while only a fair/good level was reached for W ₁ (linear: 0.67; exponential: 0.77; accelerometric: 0.52), W ₂ (linear: 0.52; exponential: 0.69; accelerometric: 0.83) and W _b (linear: 0.60; exponential: 0.44; accelerometric: 0.50).

SIGNIFICANCE

Measuring carotid pressure waveforms with different approaches does not influence the cardiovascular parameters obtained by WS and RW; those derived by WI are affected by the carotid pressure curve employed.

Measurement of Wall Shear Stress Exerted by Flowing Blood in the Human Carotid Artery: Ultrasound Doppler Velocimetry and Echo Particle Image Velocimetry

Vascular endothelial cells lining the arteries are sensitive to wall shear stress (WSS) exerted by flowing blood. An important component of the pathophysiology of vascular diseases, WSS is commonly estimated by centerline ultrasound Doppler velocimetry (UDV). However, the accuracy of this method is uncertain. We have previously validated the use of a novel, ultrasound-based, particle image velocimetry technique (echo PIV) to compute 2-D velocity vector fields, which can easily be converted into WSS data. We compared WSS data derived from UDV and echo PIV in the common carotid artery of 27 healthy participants. Compared with echo PIV, time-averaged WSS was lower using UDV (28 ± 35%). Echo PIV revealed that this was due to considerable spatiotemporal variation in the flow velocity profile, contrary to the assumption that flow is steady and the velocity profile is parabolic throughout the cardiac cycle. The largest WSS underestimation by UDV was found during peak systole (118 ± 16%) and the smallest during mid-diastole (4.3± 46%). The UDV method underestimated WSS for the accelerating and decelerating systolic measurements (68 ± 30% and 24 ± 51%), whereas WSS was overestimated for end-diastolic measurements (-44 ± 55%). Our data indicate that UDV estimates of WSS provided limited and largely inaccurate information about WSS and that the complex spatiotemporal flow patterns do not fit well with traditional assumptions about blood flow in arteries. Echo PIV-derived WSS provides detailed information about this important but poorly understood stimulus that influences vascular endothelial pathophysiology.

Effects of inflow velocity profile on two-dimensional hemodynamic analysis by ordinary and ultrasonic-measurement-integrated simulations

Two-dimensional ultrasonic-measurement-integrated (2D-UMI) simulation correctly reproduces hemodynamics even with an inexact inflow velocity distribution. This study aimed to investigate which is superior, a two-dimensional ordinary (2D-O) simulation with an accurate inflow velocity distribution or a 2D-UMI simulation with an inaccurate one. 2D-O and 2D-UMI simulations were performed for blood flow in a carotid artery with four upstream velocity boundary conditions: a velocity profile with backprojected measured Doppler velocities (condition A), and velocity profiles with a measured Doppler velocity distribution, a parabolic one, and a uniform one, magnitude being obtained by inflow velocity estimation (conditions B, C, and D, respectively). The error of Doppler velocity against the measurement data was sensitive to the inflow velocity distribution in the 2D-O simulation, but not in the 2D-UMI simulation with the inflow velocity estimation. Among the results in conditions B, C, and D, the error in the worst 2D-UMI simulation with condition D was 31 % of that in the best 2D-O simulation with condition B, implying the superiority of the 2D-UMI simulation with an inaccurate inflow velocity distribution over the 2D-O simulation with an exact one. Condition A resulted in a larger error than the other conditions in both the 2D-O and 2D-UMI simulations.

Association of common carotid artery Doppler-determined dicrotic notch velocity with the left ventricular ejection fraction

OBJECTIVES

The appearance of the dicrotic notch on blood pressure tracings is associated with impaired cardiac function. Common carotid artery waveforms have similar fiduciary markers, yet they have not been related to cardiac function. We studied associations of common carotid artery dicrotic notch velocities with the left ventricular ejection fraction (LVEF) determined by echocardiography.

METHODS

We conducted a retrospective study of 37 patients who had cardiac echocardiography and carotid Doppler evaluations within 1 day of each other. The LVEF was determined by the biplane modified Simpson rule. Doppler parameters were measured from tracings of the left common carotid artery 4 cm from the flow divider. Linear regression and stepwise multivariable linear regression models were used to evaluate any association between the LVEF and the following variables: age, sex, peak systolic velocity (PSV), end-diastolic velocity (EDV), dicrotic notch velocity, rise time (EDV to PSV), resistive index, and cardiac cycle length.

RESULTS

The dicrotic notch velocity was the only variable associated with the LVEF (P = .028) in a bivariate analyses. A backward selection stepwise multivariable equation predicting the LVEF had the dicrotic notch (P = .001) and resistive index (P = .01) as significant predictors, whereas the cardiac cycle length (P = .08) and PSV (P = .08) were borderline not significant. Model goodness of fit was R(2) = 0.37 (P = .004).

CONCLUSIONS

Dicrotic notch velocities measured from common carotid artery Doppler waveforms are associated with the LVEF and might offer some clinical value in selected cases.

A framework for the co-registration of hemodynamic forces and atherosclerotic plaque components

Local hemodynamic forces, such as wall shear stress (WSS), are thought to trigger cellular and molecular mechanisms that determine atherosclerotic plaque vulnerability to rupture. Magnetic resonance imaging has emerged as a powerful tool to characterize human carotid atherosclerotic plaque composition and morphology, and to identify plaque features shown to be key determinants of plaque vulnerability. Image-based computational fluid dynamics has allowed researchers to obtain time-resolved WSS information of atherosclerotic carotid arteries. A deeper understanding of the mechanisms of initiation and progression of atherosclerosis can be obtained through the comparison of WSS and plaque composition and morphology. To date, however, advance in knowledge has been limited greatly due to the lack of a reliable infrastructure to perform such analysis. The aim of this study is to establish a framework that will allow for the co-registration and analysis of the three-dimensional distribution of WSS and plaque components and morphology. The use of this framework will lead to future studies targeted to determining the role of WSS in atherosclerotic plaque progression and vulnerability.

A mathematical evaluation of hemodynamic parameters after carotid eversion and conventional patch angioplasty

Carotid endarterectomy has a long history in stroke prevention, yet controversy remains concerning optimal techniques. Two methods frequently used are endarterectomy with patch angioplasty (CEAP) and eversion endarterectomy (CEE). The objective of this study was to compare hemodynamics-related stress and strain distributions between arteries repaired using CEAP and CEE. Mathematical models were based on in vivo three-dimensional arterial geometry, pulsatile velocity profiles, and intraluminal pressure inputs obtained from 16 patients with carotid artery disease. These data were combined with experimentally derived nonlinear, anisotropic carotid artery mechanical properties to create fluid-structure interaction models of CEAP and CEE. These models were then used to calculate hemodynamic parameters thought to promote recurrent disease and restenosis. Combining calculations of stress and strain into a composite risk index, called the integral abnormality factor, allowed for an overall comparison between CEAP and CEE. CEE demonstrated lower mechanical stresses in the arterial wall, whereas CEAP straightened the artery and caused high stress and strain concentrations at the suture-artery interface. CEAP produced a larger continuous region of oscillatory, low-shear, vortical flow in the carotid bulb. There was a more than two-fold difference in the integral abnormality factor, favoring CEE. In conclusion, in a realistically simulated carotid artery, fluid-structure interaction modeling demonstrated CEE to produce less mechanical wall stress and improved flow patterns compared with CEAP. Clinical validation with larger numbers of individual patients will ultimately be required to support modeling approaches to help predict arterial disease progression and comparative effectiveness of reconstruction methods and devices.

Effect of velocity profile skewing on blood velocity and volume flow waveforms derived from maximum Doppler spectral velocity

Given evidence that fully developed axisymmetric flow may be the exception rather than the rule, even in nominally straight arteries, maximum velocity (V(max)) can lie outside the Doppler sample volume (SV). The link between V(max) and derived quantities, such as volume flow (Q), may therefore be more complex than commonly thought. We performed idealized virtual Doppler ultrasound on data from image-based computational fluid dynamics (CFD) models of the normal human carotid artery and investigated how velocity profile skewing and choice of sample volume affected V(max) waveforms and derived Q variables, considering common assumptions about velocity profile shape (i.e., Poiseuille or Womersley). Severe velocity profile skewing caused substantial errors in V(max) waveforms when using a small, centered SV, although peak V(max) was reliably detected; errors with a long SV covering the vessel diameter were orientation dependent but lower overall. Cycle-averaged Q calculated from V(max) was typically within ±15%, although substantial skewing and use of a small SV caused 10%-25% underestimation. Peak Q derived from Womersley's theory was generally accurate to within ±10%. V(max) pulsatility and resistance indexes differed from Q-based values, although the Q-based resistance index could be predicted reliably. Skewing introduced significant error into V(max)-derived Q waveforms, particularly during mid-to-late systole. Our findings suggest that errors in the V(max) and Q waveforms related to velocity profile skewing and use of a small SV, or orientation-dependent errors for a long SV, could limit their use in wave analysis or for constructing characteristic or patient-specific flow boundary conditions for model studies.

Common carotid and brachial artery hemodynamic alterations in periodontal disease

AIM

In patients affected by periodontal disease, hypertension and systemic inflammation might cause an arterial hemodynamic derangement; this, in turn, can act as a mediator of the atherogenic process often seen in these patients. This study aimed at a comprehensive hemodynamic evaluation in periodontal patients.

METHODS

Fourty-eight subjects participating to a cardiovascular prevention programme were enrolled. Periodontitis, classical risk factors for atherosclerosis, and shear and tensile forces in both carotid and brachial arteries were evaluated. Calculated periodontal indexes were plaque, gingival and pocket deep (PD) indexes. Simple and multiple regression analyses were performed. Afterwards, 30 of them with normal PD index were compared with 30 carefully - matched patients with periodontitis.

RESULTS

Brachial and carotid parietal tension were significantly associated with periodontal indexes, especially PD-Sum, in both simple (r = 0.42, p < 0.001 for carotid artery and r = 0.36, p < 0.02 for brachial artery) and multiple regression analyses. Shear stress gave similar results. In case-control analysis, shear stress was lower by 15% and 30%, respectively, in carotid and brachial artery in patients with high PD; common carotid parietal tension was higher. Arterial stiffness resulted not associated with periodontitis.

CONCLUSIONS

Periodontal disease is associated to a complex atherosclerotic prone hemodynamic derangement, particularly in large elastic arteries.

Errors in the estimation of wall shear stress by maximum Doppler velocity

OBJECTIVE

Wall shear stress (WSS) is an important parameter with links to vascular (dys)function. Difficult to measure directly, WSS is often inferred from maximum spectral Doppler velocity (Vmax) by assuming fully-developed flow, which is valid only if the vessel is long and straight. Motivated by evidence that even slight/local curvatures in the nominally straight common carotid artery (CCA) prevent flow from fully developing, we investigated the effects of velocity profile skewing on Vmax-derived WSS.

METHODS

Velocity profiles, representing different degrees of skewing, were extracted from the CCA of image-based computational fluid dynamics (CFD) simulations carried out as part of the VALIDATE study. Maximum velocities were calculated from idealised sample volumes and used to estimate WSS via fully-developed (Poiseuille or Womersley) velocity profiles, for comparison with the actual (i.e. CFD-derived) WSS.

RESULTS

For cycle-averaged WSS, mild velocity profile skewing caused ±25% errors by assuming Poiseuille or Womersley profiles, while severe skewing caused a median error of 30% (maximum 55%). Peak systolic WSS was underestimated by ~50% irrespective of skewing with Poiseuille; using a Womersley profile removed this bias, but ±30% errors remained. Errors were greatest in late systole, when skewing was most pronounced. Skewing also introduced large circumferential WSS variations: ±60%, and up to ±100%, of the circumferentially averaged value.

CONCLUSION

Vmax-derived WSS may be prone to substantial variable errors related to velocity profile skewing, and cannot detect possibly large circumferential WSS variations. Caution should be exercised when making assumptions about velocity profile shape to calculate WSS, even in vessels usually considered long and straight.

Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.

Human common carotid wall shear stress as a function of age and gender: a 12-year follow-up study

Hemodynamic forces play a role in the development of atherosclerosis. Their variations with age have been assessed in cross-sectional, but not longitudinal, studies. The aim of the present study was to investigate in both sexes the age-dependent change in wall shear stress and arterial stiffness in subjects studied twice 12 years apart. Forty-eight subjects (15 women and 33 men) were studied twice 12 years apart. Subjects underwent blood viscosity measurement and echo-Doppler of carotid arteries, from which the intima-media thickness (IMT) was measured and the wall shear stress and Peterson's elastic modulus were calculated. Blood viscosity increased in both sexes, more markedly in women (+13.2%) than men (+7.2%). Common carotid diameter increased in both sexes, but in men (+7.4%) more than in women (+5.5%). Peak and mean velocity decreased at follow-up by 10.7% and 18.9% in women and by 14.2% and 18.1% in men. Peak and mean shear stress significantly decreased in men by 13.0% and 17.5%, respectively, while in women only the mean shear stress was reduced (-11.8%). The IMT of the common carotid artery increased by 29% in women and 20% in men. Arterial stiffness significantly increased (+74.5% in women and +28.0% in men). The percent change in mean shear stress was significantly and inversely associated with the percent change in IMT. The data of this study show that, in a middle-aged population observed for almost 12 years, the mean shear stress decreases significantly in both sexes, while peak shear stress decreases significantly only in men. The change in mean shear stress is inversely associated with changes in IMT. Arterial stiffness, on the other hand, increases with aging.

Comparison of non-invasive methods for measurement of local pulse wave velocity using FSI-simulations and in vivo data

In the search for better predictors of cardiovascular events, pulse wave velocity (PWV) has gained considerable interest. We compared three single-location methods to locally estimate PWV based on simultaneous measurements of pressure (P), velocity (U) or arterial diameter (D): the PU, ln(D)U and QA-method. First, the performance of these methods was analyzed using 3D fluid-structure interaction simulations (FSI) in a tube and patient-specific carotid artery. We demonstrated that the outcome was dependent on whether the methods were tested in the ideal conditions of a 3D axisymmetrical and reflection-free tube or in the more realistic setup of a carotid artery. The three single-location PWV methods performed similarly in the tube (4.29 m/s for PU, 4.44 m/s for ln(D)U and 4.38 m/s for QA) while the carotid data showed that the PU-method dramatically overestimates PWV (9.16 m/s), and the ln(D)U and QA-method underestimate (3.86 and 3.84 m/s, respectively). The erroneously high wavespeeds from the PU-method were attributed to wave reflections, which was confirmed by measurements in 37 healthy adults. This in vivo study showed similar discrepancies between the 3 single-location techniques as present in the carotid simulations, with the difference between the PU- and ln(D)U-method related to the magnitude of wave reflection.

On the shape of the common carotid artery with implications for blood velocity profiles

Clinical and engineering studies typically assume that the common carotid artery (CCA) is straight enough to assume fully developed flow, yet recent studies have demonstrated the presence of skewed velocity profiles. Toward elucidating the influence of mild vascular curvatures on blood flow patterns and atherosclerosis, this study aimed to characterize the three-dimensional shape of the human CCA. The left and right carotid arteries of 28 participants (63 ± 12 years) in the VALIDATE (Vascular Aging--The Link that Bridges Age to Atherosclerosis) study were digitally segmented from 3D contrast-enhanced magnetic resonance angiograms, from the aortic arch to the carotid bifurcation. Each CCA was divided into nominal cervical and thoracic segments, for which curvatures were estimated by least-squares fitting of the respective centerlines to planar arcs. The cervical CCA had a mean radius of curvature of 127 mm, corresponding to a mean lumen:curvature radius ratio of 1:50. The thoracic CCA was significantly more curved at 1:16, with the plane of curvature tilted by a mean angle of 25° and rotated close to 90° with respect to that of the cervical CCA. The left CCA was significantly longer and slightly more curved than the right CCA, and there was a weak but significant increase in CCA curvature with age. Computational fluid dynamic simulations carried out for idealized CCA geometries derived from these and other measured geometric parameters demonstrated that mild cervical curvature is sufficient to prevent flow from fully-developing to axisymmetry, independent of the degree of thoracic curvature. These findings reinforce the idea that fully developed flow may be the exception rather than the rule for the CCA, and perhaps other nominally long and straight vessels.

Estimation of blood velocity, volumetric flow and wall shear rate using Doppler ultrasound

Commercial ultrasound systems can make a number of measurements related to haemodynamics which are relevant to clinical practice and to clinical research. These include maximum velocity, volumetric flow and wall shear rate. Using appropriate protocols, measurements can be made averaged over the cardiac cycle, or as a function of time through the cardiac cycle. Maximum velocity underpins most of these measurements. Maximum velocity is overestimated as a result of geometric spectral broadening, by typically up to 30%, but by much larger amounts as the angle approaches 90°. Though not used in clinical practice, a simple correction technique using a string phantom can substantially reduce these errors. For volumetric flow and wall shear rate, methods such as specialist multi-gate ultrasound systems, magnetic resonance imaging and image guided modelling are available. Before resorting to these more complex methods users might consider that, with care and attention to procedure, high quality information may be obtained using commercial ultrasound systems. Manufacturers could make more use of the colour flow image for quantification of velocity, and adopt vector Doppler techniques.

In vivo three-dimensional blood velocity profile shapes in the human common, internal, and external carotid arteries

OBJECTIVE

True understanding of carotid bifurcation pathophysiology requires a detailed knowledge of the hemodynamic conditions within the arteries. Data on carotid artery hemodynamics are usually based on simplified, computer-based, or in vitro experimental models, most of which assume that the velocity profiles are axially symmetric away from the carotid bulb. Modeling accuracy and, more importantly, our understanding of the pathophysiology of carotid bifurcation disease could be considerably improved by more precise knowledge of the in vivo flow properties within the human carotid artery. The purpose of this work was to determine the three-dimensional pulsatile velocity profiles of human carotid arteries.

METHODS

Flow velocities were measured over the cardiac cycle using duplex ultrasonography, before and after endarterectomy, in the surgically exposed common (CCA), internal (ICA), and external (ECA) carotid arteries (n = 16) proximal and distal to the stenosis/endarterectomy zone. These measurements were linked to a standardized grid across the flow lumina of the CCA, ICA, and ECA. The individual velocities were then used to build mean three-dimensional pulsatile velocity profiles for each of the carotid artery branches.

RESULTS

Pulsatile velocity profiles in all arteries were asymmetric about the arterial centerline. Posterior velocities were higher than anterior velocities in all arteries. In the CCA and ECA, velocities were higher laterally, while in the ICA, velocities were higher medially. Pre- and postendarterectomy velocity profiles were significantly different. After endarterectomy, velocity values increased in the common and internal and decreased in the external carotid artery.

CONCLUSIONS

The in vivo hemodynamics of the human carotid artery are different from those used in most current computer-based and in vitro models. The new information on three-dimensional blood velocity profiles can be used to design models that more closely replicate the actual hemodynamic conditions within the carotid bifurcation. Such models can be used to further improve our understanding of the pathophysiologic processes leading to stroke and for the rational design of medical and interventional therapies.

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.

Carotid bifurcation hemodynamics in older adults: effect of measured versus assumed flow waveform

Recent work has illuminated differences in carotid artery blood flow rate dynamics of older versus young adults. To what degree flow waveform shape, and indeed the use of measured versus assumed flow rates, affects the simulated hemodynamics of older adult carotid bifurcations has not been elucidated. Image-based computational fluid dynamics models of N=9 normal, older adult carotid bifurcations were reconstructed from magnetic resonance angiography. Subject-specific hemodynamics were computed by imposing each individual's inlet and outlet flow rates measured by cine phase-contrast magnetic resonance imaging or by imposing characteristic young and older adult flow waveform shapes adjusted to cycle-averaged flow rates measured or allometrically scaled to the inlet and outlet areas. Despite appreciable differences in the measured versus assumed flow conditions, the locations and extents of low wall shear stress and elevated relative residence time were broadly consistent; however, the extent of elevated oscillatory shear index was substantially underestimated, more by the use of assumed cycle-averaged flow rates than the assumed flow waveform shape. For studies of individual vessels, use of a characteristic flow waveform shape is likely sufficient, with some benefit offered by scaling to measured cycle-averaged flow rates. For larger-scale studies of many vessels, ranking of cases according to presumed hemodynamic or geometric risk is robust to the assumed flow conditions.

Outflow Conditions for Image-Based Hemodynamic Models of the Carotid Bifurcation: Implications for Indicators of Abnormal Flow

Computational fluid dynamics (CFD) models have become very effective tools for predicting the flow field within the carotid bifurcation, and for understanding the relationship between local hemodynamics, and the initiation and progression of vascular wall pathologies. As prescribing proper boundary conditions can affect the solutions of the equations governing blood flow, in this study, we investigated the influence to assumptions regarding the outflow boundary conditions in an image-based CFD model of human carotid bifurcation. Four simulations were conducted with identical geometry, inlet flow rate, and fluid parameters. In the first case, a physiological time-varying flow rate partition at branches along the cardiac cycle was obtained by coupling the 3D model of the carotid bifurcation at outlets with a lumped-parameter model of the downstream vascular network. Results from the coupled model were compared with those obtained by imposing three fixed flow rate divisions (50/50, 60/40, and 70/30) between the two branches of the isolated 3D model of the carotid bifurcation. Three hemodynamic wall parameters were considered as indicators of vascular wall dysfunction. Our findings underscore that the overall effect of the assumptions done in order to simulate blood flow within the carotid bifurcation is mainly in the hot-spot modulation of the hemodynamic descriptors of atherosusceptible areas, rather than in their distribution. In particular, the more physiological, time-varying flow rate division deriving from the coupled simulation has the effect of damping wall shear stress (WSS) oscillations (differences among the coupled and the three fixed flow partition models are up to 37.3% for the oscillating shear index). In conclusion, we recommend to adopt more realistic constraints, for example, by coupling models at different scales, as in this study, when the objective is the outcome prediction of alternate therapeutic interventions for individual patients, or to test hypotheses related to the role of local fluid dynamics and other biomechanical factors in vascular diseases.

Endothelial shear stress: a critical determinant of arterial remodeling and arterial stiffness in humans--a carotid 3.0-T MRI study

BACKGROUND

Low endothelial shear stress (ESS) elicits endothelial dysfunction. However, the relationship between ESS and arterial remodeling and arterial stiffness is unknown in humans. We developed a 3.0-T MRI protocol to evaluate the contribution of ESS to arterial remodeling and stiffness.

METHODS AND RESULTS

Fifteen young (aged 26 ± 3 years) and 15 older (aged 57 ± 3 years) healthy volunteers as well as 15 patients with cardiovascular disease (aged 63 ± 10 years) were enrolled. Phase-contrast MRI of the common carotid arteries was used to derive ESS data from the spatial velocity gradients close to the arterial wall. ESS measurements were performed on 3 occasions and showed excellent reproducibility (intraclass correlation coefficient, 0.79). Multiple linear regression analysis accounting for age and blood pressure revealed that ESS was an independent predictor of the following response variables: carotid wall thickness (regression coefficient [b], -0.19 mm(2) per N/m(2); P=0.02), lumen area (b, -15.5 mm(2) per N/m(2); P<0.001), and vessel size (b, -24.0 mm(2) per N/m(2); P<0.001). Segments of the artery wall exposed to lower ESS were significantly thicker than segments exposed to higher ESS within the same artery (P=0.009). Furthermore, ESS was associated with arterial compliance, accounting for age, blood pressure, and wall thickness (b, -0.003 mm(2)/mm Hg per N/m(2); P=0.04).

CONCLUSIONS

Our carotid MRI data show that ESS is an important determinant of arterial remodeling and arterial stiffness in humans. The data warrant further studies to evaluate use of carotid ESS as a noninvasive tool to improve the understanding of individual cardiovascular disease risk and to assess novel drug therapies in cardiovascular disease prevention.