Flow patterns at the stenosed carotid bifurcation: effect of concentric versus eccentric stenosis

Numerical analysis of blood flow through stenosed microvessels using a multi-phase model

Blood flow in arterioles have attracted considerable research attention due to their clinical implications. However, the fluid structure interaction between red blood cells and plasma in the blood poses formidable difficulty to the computational efforts. In this contribution, we seek to represent the red blood cells in the blood as a continuous non-Newtonian phase and construct a multi-phase model for the blood flow in microvessels. The methods are presented and validated using a channel with sudden expansion. And the resulting blood flow inside a stenosed microvessel is investigated at different inlet velocity amplitudes and hematocrits. It is show that the increase of both inlet velocity amplitude and inlet hematocrit leads to longer and thicker cell-rich layer downstream the stenosis. Besides, it is found that the maximum values of wall shear stress scales up with inlet velocity amplitudes and hematocrits. These results show the validity of the proposed computational model and provide helpful insights into blood flow behaviors inside stenosed vessels.

Time-resolved simulation of blood flow through left anterior descending coronary artery: effect of varying extent of stenosis on hemodynamics

BACKGROUND AND OBJECTIVES

Real-time blood flow variation is crucial for understanding the dynamic development of coronary atherosclerosis. The main objective of this study is to investigate the effect of varying extent of stenosis on the hemodynamic features in left anterior descending coronary artery.

METHODS

Various Computational fluid dynamics (CFD) models were constructed with patient-specific CT image data, using actual fractional flow reserve (FFR) as boundary conditions to provide a real-time quantitative description of hemodynamic properties. The hemodynamic parameters, such as the local and instantaneous wall shear stress (WSS), oscillating shear index (OSI) and relative residence time (RRT), blood flow velocity and pressure drop during various phases of cardiac cycle were provided in detail.

RESULTS

There was no evident variation in hemodynamic parameters in the cases of less than 50% stenosis while there were abrupt and dramatic changes in hemodynamics when the stenosis aggravated from 60 to 70%. Furthermore, when the stenosis was beyond 70%, there existed substantial pressure difference, WSS, and blood flow velocity in the center of the stenosis. Although OSI and RRT increased along with the aggravation of stenosis, they appeared with obvious abnormalities across all cases, even in mild stenosis.

CONCLUSION

The simulation could present a dynamic and comprehensive profile of how hemodynamic parameters vary in accordance with divergent severities of stenosis, which could serve as an effective reference for the clinicians to have a deeper insight into the pathological mechanism of coronary atherosclerosis and stenosis.

Atrial Fibrillation and Anterior Cerebral Artery Absence Reduce Cerebral Perfusion: A De Novo Hemodynamic Model

Background: Atrial fibrillation is a prevalent cardiac arrhythmia and may reduce cerebral blood perfusion augmenting the risk of dementia. We hypothesize that geometric variations in the cerebral arterial structure called the Circle of Willis (CoW) play an important role in influencing cerebral perfusion. The objective of this work was to develop a novel cardio-cerebral lumped parameter hemodynamic model to investigate the role of CoW variants on cerebral blood flow dynamics under atrial fibrillation conditions. Methods: A computational blood flow model was developed by coupling whole-body and detailed cerebral circulation descriptions, modified to represent six common variations of the CoW. Cerebral blood flow dynamics were simulated in common CoW variants, under control and imposed atrial fibrillation conditions. Risk was assessed based on the frequency of beat-wise hypoperfusion events, and sensitivity analysis was performed with respect to this model output. Results: It was found that the geometry of the CoW influenced the frequency of hypoperfusion events at different heart rates, with the variant missing a P1 segment having the highest risk. Sensitivity analysis revealed that intrinsic heart rate is most associated with the considered outcome. Conclusions: Our results suggest that CoW geometry plays an important role in influencing cerebral hemodynamics during atrial fibrillation. The presented study may assist in guiding our future clinical-imaging research.

Simulation of turbulence induced sound generation inside stenosed femoral artery models with different severities and eccentricities

BACKGROUND AND OBJECTIVES

Recent developments of low-cost, compact acoustic sensors, advanced signal processing tools and powerful computational resources allow researchers design new scoring systems for acoustic detection of arterial stenoses. In this study, numerical simulations of blood flow inside stenosed arteries are performed to understand the effect of stenosis severity and eccentricity on the turbulence induced wall pressure fluctuations and the generated sound.

METHODS

Axisymmetric and eccentric elliptic stenoses of five different severities are generated inside a 6.4 mm diameter femoral artery model. Large eddy simulations of pulsatile, non-Newtonian blood flow are performed using the open source software OpenFOAM.

RESULTS

Post-stenotic turbulence activity is found to be almost zero for 50 and 60% severities. For severities of 75% and more, turbulent kinetic energy rises significantly with increasing severity. The location of the highest turbulence activity on the vessel wall from the stenosis exit decreases with increasing severity. The maximum level of turbulent kinetic energy seen in 95% severity models is about 9 and 31 times higher than that of 87% and 75% models, respectively. Spectrum of wall pressure fluctuations show that 50 and 60% axisymmetric models are almost silent. The spectrum starts to get richer with 75% severity, and the fluctuation intensity increases with severity. Compared to the axisymmetric models, more activity is observed in the 0-150 Hz band for the 50 and 60% eccentric models. Axial extent of the acoustically active region is also longer in them. Converting wall pressure data into sound revealed that murmurs that can be considered as signs of vascular stenosis are obtained for models with 75% and higher severity.

CONCLUSIONS

Sound patterns generated from simulation results are similar to the typical sounds obtained by Doppler ultrasonography, and present distinct characters. Together with a sensor technology that can measure these sounds from within the stenosed artery, they can be processed and used for the purpose of non-invasive diagnosis. Computational fluid dynamics studies that simulate large number of cases with different stenosis severities and morphologies will play a critical role in developing the necessary sound databases, which can be used to train new diagnostic devices.

Computed Poststenotic Flow Instabilities Correlate Phenotypically With Vibrations Measured Using Laser Doppler Vibrometry: Perspectives for a Promising In Vivo Device for Early Detection of Moderate and Severe Carotid Stenosis

Early detection of asymptomatic carotid stenosis is crucial for treatment planning in the prevention of ischemic stroke. Auscultation, the current first-line screening methodology, comes with severe limitations that create urge for novel and robust techniques. Laser Doppler vibrometer (LDV) is a promising tool for inferring carotid stenosis by measuring stenosis-induced vibrations. The goal of the current study was to evaluate the feasibility of LDV for carotid stenosis detection. LDV measurements on a carotid phantom were used to validate our previously verified high-resolution computational fluid dynamics methodology, which was used to evaluate the impact of flowrate, flow split, and stenosis severity on the poststenotic intensity of flow instabilities (IFI). We evaluated sensitivity, specificity, and accuracy of using IFI for stenoses detection. Linear regression analyses showed that computationally derived pressure fluctuations correlated (R2 = 0.98) with LDV measurements of stenosis-induced vibrations. The flowrate of stenosed vessels correlated (R2 = 0.90) with the presence of poststenotic instabilities. Receiver operating characteristic analyses of power spectra revealed that the most relevant frequency bands for the detection of moderate (56-76%) and severe (86-96%) stenoses were 80-200 Hz and 0-40 Hz, respectively. Moderate stenosis was identified with sensitivity and specificity of 90%; values decreased to 70% for severe stenosis. The use of LDV as screening tool for asymptomatic stenosis can potentially provide improved accuracy of current screening methodologies for early detection. The applicability of this promising device for mass screening is currently being evaluated clinically.

Characteristics of Wall Shear Stress and Pressure of Intracranial Atherosclerosis Analyzed by a Computational Fluid Dynamics Model: A Pilot Study

Background: Although wall shear stress (WSS) and pressure play important roles in plaque vulnerability, characteristics of the two indices in intracranial atherosclerosis (ICAS) have not been fully investigated yet. This study aimed to elucidate this issue by means of establishing a non-invasive computational fluid dynamics method with time-of-flight magnetic resonance angiography (TOF-MRA) of the whole cerebral artery.

Materials and Methods: Subjects with symptomatic ICAS in the middle cerebral artery domain were enrolled, excluding those with concomitant internal carotid artery stenosis. Based on patient-specific TOF-MRA images for three-dimensional (3D) meshes and arterial blood pressure with patient-specific carotid artery ultrasonography for inlet boundary conditions, patients' three-dimensional hemodynamics were modeled by a finite element method governed by Navier-Stokes equations.

Results: Among the 55 atherosclerotic lesions analyzed by this TOF-MRA based computational fluid dynamics model, the maximum WSS (WSS_max) was most frequently detected at the apex points and the upper half of the upstream sections of the lesions, whereas the maximum pressure was most often located at the lower half of the upstream sections. As the percent stenosis increases, the relative value of WSS_max and pressure drop increased with significantly increasing steep beyond 50% stenosis. Moreover, WSS_max was found to linearly correlate with pressure drop in ICAS.

Conclusions: This study on ICAS revealed certain trends of longitudinal distribution of WSS and pressure and the influences of percent stenosis on cerebral hemodynamics, as well as the correlations between WSS and pressure drop. It represents a step forward in applying computational flow simulation techniques in studying ICAS and stroke, in a patient-specific manner.

Ultrasound vector projectile imaging for detection of altered carotid bifurcation hemodynamics during reductions in cardiac output

PURPOSE

Complex blood flow is commonly observed in the carotid bifurcation, although the factors that regulate these patterns beyond arterial geometry are unknown. The emergence of high-frame-rate ultrasound vector flow imaging allows for noninvasive, time-resolved analysis of complex hemodynamic behavior in humans, and it can potentially help researchers understand which physiological stressors can alter carotid bifurcation hemodynamics in vivo. Here, we seek to pursue the first use of vector projectile imaging (VPI), a dynamic form of vector flow imaging, to analyze the regulation of carotid bifurcation hemodynamics during experimental reductions in cardiac output induced via a physiological stressor called lower body negative pressure (LBNP).

METHODS

Seven healthy adults (age: 27 ± 4 yr, 4 men) underwent LBNP at -45 mmHg to simulate a postural hemodynamic response in a controlled environment. Using a research-grade, high-frame-rate ultrasound platform, vector flow estimation in each subject's right carotid bifurcation was performed through a multi-angle plane wave imaging (two transmission angles of 10° and -10°) formulation, and VPI cineloops were generated at a frame rate of 750 fps. Vector concentration was quantified by the resultant blood velocity vector angles within a region of interest; lower concentration indicated greater flow dispersion. Discrete concentration values during peak and late systole were compared across different segments of the carotid artery bifurcation before, and during, LBNP.

RESULTS

Vector projectile imaging revealed that external and internal carotid arteries exhibited regional hemodynamic changes during LBNP, which acted to reduce both the subject's cardiac output (Δ - 1.2 ± 0.5 L/min, -19%; P < 0.01) and peak carotid blood velocity (Δ - 6.30 ± 8.27 cm/s, -7%; P = 0.05). In these carotid artery branches, the vector concentration time trace before and during LBNP were observed to be different. The impact of LBNP on flow complexity in the two carotid artery branches showed variations between subjects.

CONCLUSIONS

Using VPI, intuitive visualization of complex hemodynamic changes can be obtained in healthy humans subjected to LBNP. This imaging tool has potential for further applications in vascular physiology to identify and quantify complex hemodynamic features in humans during different physiological stressor tests that regulate hemodynamics.

Accuracy and Precision of a Plane Wave Vector Flow Imaging Method in the Healthy Carotid Artery

The objective of the study described here was to investigate the accuracy and precision of a plane wave 2-D vector flow imaging (VFI) method in laminar and complex blood flow conditions in the healthy carotid artery. The approach was to study (i) the accuracy for complex flow by comparing the velocity field from a computational fluid dynamics (CFD) simulation to VFI estimates obtained from the scan of an anthropomorphic flow phantom and from an in vivo scan; (ii) the accuracy for laminar unidirectional flow in vivo by comparing peak systolic velocities from VFI with magnetic resonance angiography (MRA); (iii) the precision of VFI estimation in vivo at several evaluation points in the vessels. The carotid artery at the bifurcation was scanned using both fast plane wave ultrasound and MRA in 10 healthy volunteers. The MRA geometry acquired from one of the volunteers was used to fabricate an anthropomorphic flow phantom, which was also scanned using the fast plane wave sequence. The same geometry was used in a CFD simulation to calculate the velocity field. Results indicated that similar flow patterns and vortices were estimated with CFD and VFI in the phantom for the carotid bifurcation. The root-mean-square difference between CFD and VFI was within 0.12 m/s for velocity estimates in the common carotid artery and the internal branch. The root-mean-square difference was 0.17 m/s in the external branch. For the 10 volunteers, the mean difference between VFI and MRA was -0.17 m/s for peak systolic velocities of laminar flow in vivo. The precision in vivo was calculated as the mean standard deviation (SD) of estimates aligned to the heart cycle and was highest in the center of the common carotid artery (SD = 3.6% for velocity magnitudes and 4.5° for angles) and lowest in the external branch and for vortices (SD = 10.2% for velocity magnitudes and 39° for angles). The results indicate that plane wave VFI measures flow precisely and that estimates are in good agreement with a CFD simulation and MRA.

Assessment of boundary conditions for CFD simulation in human carotid artery

Computational fluid dynamics (CFD) is an increasingly used method for investigation of hemodynamic parameters and their alterations under pathological conditions, which are important indicators for diagnosis of cardiovascular disease. In hemodynamic simulation models, the employment of appropriate boundary conditions (BCs) determines the computational accuracy of the CFD simulation in comparison with pressure and velocity measurements. In this study, we have first assessed the influence of inlet boundary conditions on hemodynamic CFD simulations. We selected two typical patients suspected of carotid artery disease, with mild stenosis and severe stenosis. Both patients underwent digital subtraction angiography (DSA), magnetic resonance angiography, and the invasive pressure guide wire measured pressure profile. We have performed computational experiments to (1) study the hemodynamic simulation outcomes of distributions of wall shear stress, pressure, pressure gradient and (2) determine the differences in hemodynamic performances caused by inlet BCs derived from DSA and Womersley analytical solution. Our study has found that the difference is related to the severity of the stenosis; the greater the stenosis, the more the difference ensues. Further, in our study, the two typical subjects with invasively measured pressure profile and thirty subjects with ultrasound Doppler velocimeter (UDV) measurement served as the criteria to evaluate the hemodynamic outcomes of wall shear stress, pressure, pressure gradient and velocity due to different outlet BCs based on the Windkessel model, structured-tree model, and fully developed flow model. According to the pressure profiles, the fully developed model appeared to have more fluctuations compared with the other two models. The Windkessel model had more singularities before convergence. The three outlet BCs models also showed good correlation with the UDV measurement, while the Windkessel model appeared to be slightly better ([Formula: see text]). The structured-tree model was seen to have the best performance in terms of available computational cost and accuracy. The results of our numerical simulation and the good correlation with the computed pressure and velocity with their measurements have highlighted the effectiveness of CFD simulation in patient-specific human carotid artery with suspected stenosis.

The Computational Fluid Dynamics Analyses on Hemodynamic Characteristics in Stenosed Arterial Models

Arterial stenosis plays an important role in the progressions of thrombosis and stroke. In the present study, a standard axisymmetric tube model of the stenotic artery is introduced and the degree of stenosis η is evaluated by the area ratio of the blockage to the normal vessel. A normal case (η = 0) and four stenotic cases of η = 0.25, 0.5, 0.625, and 0.75 with a constant Reynolds number of 300 are simulated by computational fluid dynamics (CFD), respectively, with the Newtonian and Carreau models for comparison. Results show that for both models, the poststenotic separation vortex length increases exponentially with the growth of stenosis degree. However, the vortex length of the Carreau model is shorter than that of the Newtonian model. The artery narrowing accelerates blood flow, which causes high blood pressure and wall shear stress (WSS). The pressure drop of the η = 0.75 case is nearly 8 times that of the normal value, while the WSS peak at the stenosis region of η = 0.75 case even reaches up to 15 times that of the normal value. The present conclusions are of generality and contribute to the understanding of the dynamic mechanisms of artery stenosis diseases.

A model of blood supply to the brain via the carotid arteries: Effects of obstructive vs. sclerotic changes

The carotid artery is one of the major supply routes of blood to the brain and a common site of vascular disease. Obstructive and sclerotic disorders within the carotid artery impact local blood flow patterns as well as overall impedance and blood supply to the brain. A lumped parameter model and an experimental in-vitro flow loop were used to study the effects of local stenosis and stiffness in the carotid artery based on a family of phantoms with different degrees of stenosis and compliance. The model also allows independent examination of the effects of downstream resistance and compliance. Mild to moderate stenosis was found to lead to minimal (∼1%) reduction in blood supply to the brain. Reduction in mean internal carotid artery (ICA) flow was statistically significant (p< 0.01) only above 70% stenosis. On the other hand, a three-fold increase in stiffness of the carotid artery, as might occur in aging, was found to lead to a modest yet statistically significant reduction (p< 0.01) in mean ICA flow. Effects of changing downstream resistance and compliance were examined. For a given pressure waveform, reduction in downstream compliance led to altered waveform shape and reduction in peak systolic flow rates where the mean flow rates were not altered. Increased downstream resistance resulted in drastic reduction in mean flow rates.

Convection-Enhanced Transport into Open Cavities : Effect of Cavity Aspect Ratio

Recirculating fluid regions occur in the human body both naturally and pathologically. Diffusion is commonly considered the predominant mechanism for mass transport into a recirculating flow region. While this may be true for steady flows, one must also consider the possibility of convective fluid exchange when the outer (free stream) flow is transient. In the case of an open cavity, convective exchange occurs via the formation of lobes at the downstream attachment point of the separating streamline. Previous studies revealed the effect of forcing amplitude and frequency on material transport rates into a square cavity (Horner in J Fluid Mech 452:199-229, 2002). This paper summarizes the effect of cavity aspect ratio on exchange rates. The transport process is characterized using both computational fluid dynamics modeling and dye-advection experiments. Lagrangian analysis of the computed flow field reveals the existence of turnstile lobe transport for this class of flows. Experiments show that material exchange rates do not vary linearly as a function of the cavity aspect ratio (A = W/H). Rather, optima are predicted for A ≈ 2 and A ≈ 2.73, with a minimum occurring at A ≈ 2.5. The minimum occurs at the point where the cavity flow structure bifurcates from a single recirculating flow cell into two corner eddies. These results have significant implications for mass transport environments where the geometry of the flow domain evolves with time, such as coronary stents and growing aneurysms. Indeed, device designers may be able to take advantage of the turnstile-lobe transport mechanism to tailor deposition rates near newly implanted medical devices.

MODELING BLOOD FLOW IN AN ECCENTRIC STENOSED ARTERY USING LARGE EDDY SIMULATION AND PARALLEL COMPUTING

Computational fluid dynamics (CFD) is an excellent computational tool to assess the hemodynamics and detailed blood-flow structure for cardiovascular applications. Modeling turbulence for cardiovascular applications can be achieved (to some extent) using available numerical models such as Reynolds average Navier–Stokes (RANS), the large eddy simulation (LES) and the direct numerical solution (DNS). In order to develop an efficient model which is as accurate as DNS and as quick as RANS, our laboratory's focus is on LES. In this study, we develop an efficient numerical model which is based on LES and structured but non-orthogonal finite volumes. Using the proposed model, the detailed flow structure and turbulent features of the blood stream in a complicated geometry is captured. The aim of this study is to model blood-flow through an eccentric stenosis accurately and quickly. The results are similar to those obtained using DNS but in a fraction of the CPU time. The computational tools implemented in this study are based on a FORTRAN based in-house code coupled with parallel computing using SHARCNET. The developed model is a significant computational tool which can be used to assess the hemodynamic properties for cardiovascular applications, e.g., prosthetic heart valves and atherosclerosis.

Vector projectile imaging: time-resolved dynamic visualization of complex flow patterns

Achieving non-invasive, accurate and time-resolved imaging of vascular flow with spatiotemporal fluctuations is well acknowledged to be an ongoing challenge. In this article, we present a new ultrasound-based framework called vector projectile imaging (VPI) that can dynamically render complex flow patterns over an imaging view at millisecond time resolution. VPI is founded on three principles: (i) high-frame-rate broad-view data acquisition (based on steered plane wave firings); (ii) flow vector estimation derived from multi-angle Doppler analysis (coupled with data regularization and least-squares fitting); (iii) dynamic visualization of color-encoded vector projectiles (with flow speckles displayed as adjunct). Calibration results indicated that by using three transmit angles and three receive angles (-10°, 0°, +10° for both), VPI can consistently compute flow vectors in a multi-vessel phantom with three tubes positioned at different depths (1.5, 4, 6 cm), oriented at different angles (-10°, 0°, +10°) and of different sizes (dilated diameter: 2.2, 4.4 and 6.3 mm; steady flow rate: 2.5 mL/s). The practical merit of VPI was further illustrated through an anthropomorphic flow phantom investigation that considered both healthy and stenosed carotid bifurcation geometries. For the healthy bifurcation with 1.2-Hz carotid flow pulses, VPI was able to render multi-directional and spatiotemporally varying flow patterns (using a nominal frame rate of 416 fps or 2.4-ms time resolution). In the case of stenosed bifurcations (50% eccentric narrowing), VPI enabled dynamic visualization of flow jet and recirculation zones. These findings suggest that VPI holds promise as a new tool for complex flow analysis.

In vitro shear stress measurements using particle image velocimetry in a family of carotid artery models: effect of stenosis severity, plaque eccentricity, and ulceration

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

Simulation of a pulsatile non-Newtonian flow past a stenosed 2D artery with atherosclerosis

Atherosclerotic plaque can cause severe stenosis in the artery lumen. Blood flow through a substantially narrowed artery may have different flow characteristics and produce different forces acting on the plaque surface and artery wall. The disturbed flow and force fields in the lumen may have serious implications on vascular endothelial cells, smooth muscle cells, and circulating blood cells. In this work a simplified model is used to simulate a pulsatile non-Newtonian blood flow past a stenosed artery caused by atherosclerotic plaques of different severity. The focus is on a systematic parameter study of the effects of plaque size/geometry, flow Reynolds number, shear-rate dependent viscosity and flow pulsatility on the fluid wall shear stress and its gradient, fluid wall normal stress, and flow shear rate. The computational results obtained from this idealized model may shed light on the flow and force characteristics of more realistic blood flow through an atherosclerotic vessel.

High-frame-rate ultrasound color-encoded speckle imaging of complex flow dynamics

Realization of flow imaging at high frame rates is essential to the visualization of complex flow patterns with fast-changing spatiotemporal dynamics. In this study, we present an experimental demonstration of a novel ultrasound-based high-frame-rate flow visualization technique called color-encoded speckle imaging (CESI), which depicts flow information in a hybrid form comprising flow speckle pattern and color-encoded velocity mapping. This technique works by integrating two key principles: (i) using broad-view data acquisition schemes like plane wave compounding to obtain image data at frame rates well beyond the video display range and (ii) deriving and displaying both flow speckles and velocity estimates from the acquired broad-view image data. CESI was realized on a channel-domain ultrasound imaging research platform, and its performance was evaluated in the context of monitoring complex flow dynamics inside a carotid bifurcation flow phantom with 25% eccentric stenosis at the inlet of the internal carotid artery. Results show that, using an imaging frame rate of 2000 frames per second (based on plane wave compounding with five steering angles), CESI can effectively render flow acceleration and deceleration with visual continuity. It is also effective in depicting how stenosis-related flow disturbance events, such as flow jet formation and post-stenotic flow recirculation, evolve spatiotemporally over a pulse cycle. We anticipate that CESI can represent a rational approach to rendering flow information in ultrasound-based vascular diagnoses.

Hemodynamics analysis of patient-specific carotid bifurcation: a CFD model of downstream peripheral vascular impedance

The study of cardiovascular models was presented in this paper based on medical image reconstruction and computational fluid dynamics. Our aim is to provide a reality platform for the purpose of flow analysis and virtual intervention outcome predication for vascular diseases. By connecting two porous mediums with transient permeability at the downstream of the carotid bifurcation branches, a downstream peripheral impedance model was developed, and the effect of the downstream vascular bed impedance can be taken into consideration. After verifying its accuracy with a healthy carotid bifurcation, this model was implemented in a diseased carotid bifurcation analysis. On the basis of time-averaged wall shear stress, oscillatory shear index, and the relative residence time, fractions of abnormal luminal surface were highlighted, and the atherosclerosis was assessed from a hemodynamic point of view. The effect of the atherosclerosis on the transient flow division between the two branches because of the existence of plaque was also analysed. This work demonstrated that the proposed downstream peripheral vascular impedance model can be used for computational modelling when the outlets boundary conditions are not available, and successfully presented the potential of using medical imaging and numerical simulation to provide existing clinical prerequisites for diagnosis and therapeutic treatment.

Flow recirculation zone length and shear rate are differentially affected by stenosis severity in human coronary arteries

Flow recirculation zones and shear rate are associated with distinct pathogenic biological pathways relevant to thrombosis and atherogenesis. The interaction between stenosis severity and lesion eccentricity in determining the length of flow recirculation zones and peak shear rate in human coronary arteries in vivo is unclear. Computational fluid dynamic simulations were performed under resting and hyperemic conditions on computer-generated models and three-dimensional (3-D) reconstructions of coronary arteriograms of 25 patients. Boundary conditions for 3-D reconstructions simulations were obtained by direct measurements using a pressure-temperature sensor guidewire. In the computer-generated models, stenosis severity and lesion eccentricity were strongly associated with recirculation zone length and maximum shear rate. In the 3-D reconstructions, eccentricity increased recirculation zone length and shear rate when lesions of the same stenosis severity were compared. However, across the whole population of coronary lesions, eccentricity did not correlate with recirculation zone length or shear rate ( P = not signficant for both), whereas stenosis severity correlated strongly with both parameters ( r = 0.97, P < 0.001, and r = 0.96, P < 0.001, respectively). Nonlinear regression analyses demonstrated that the relationship between stenosis severity and peak shear was exponential, whereas the relationship between stenosis severity and recirculation zone length was sigmoidal, with an apparent threshold effect, demonstrating a steep increase in recirculation zone length between 40% and 60% diameter stenosis. Increasing stenosis severity and lesion eccentricity can both increase flow recirculation and shear rate in human coronary arteries. Flow recirculation is much more sensitive to mild changes in the severity of intermediate stenoses than is peak shear.

Transitional flow analysis in the carotid artery bifurcation by proper orthogonal decomposition and particle image velocimetry

Blood flow instabilities in the carotid artery bifurcation have been highly correlated to clot formation and mobilization resulting in ischemic stroke. In this work, PIV-measured flow velocities in normal and stenosed carotid artery bifurcation models were analyzed by means of proper orthogonal decomposition (POD). Through POD analysis, transition to more complex flow was visualized and quantified for increasing stenosis severity. While no evidence of transitional flow was seen in the normal model, the 50%-stenosed model started to show characteristics of transitional flow, which became highly evident in the 70% model, with greatest manifestation during the systolic phase of the cardiac cycle. By means of a model comparison, we demonstrate two quantitative measures of the flow complexity through the power-law decay slope of the energy spectrum and the global entropy. The more complex flow in the 70%-stenosed model showed a flatter slope of energy decay (-0.91 compared to -1.34 for 50% stenosis) and higher entropy values (0.26 compared to 0.17). Finally, the minimum temporal resolution required for POD analysis of carotid artery flow was found to be 100 Hz when determined through a more typical energy-mode convergence test, as compared to 400 Hz based on global entropy values.

Comparison of carotid plaque ulcer detection using contrast-enhanced and time-of-flight MRA techniques

BACKGROUND AND PURPOSE

Ulceration in carotid plaque is a risk indicator for ischemic stroke. Our aim was to compare plaque ulcer detection by standard TOF and CE-MRA techniques and to identify factors that influence its detection.

MATERIALS AND METHODS

Carotid MR imaging scans were acquired on 2066 participants in the ARIC study. We studied the 600 thickest plaques. TOF-MRA, CE-MRA, and black-blood MR images were analyzed together to define ulcer presence (plaque surface niche ≥2 mm in depth). Sixty ulcerated arteries were detected. These arteries were randomly assigned, along with 40 nonulcerated plaques from the remaining 540, for evaluation of ulcer presence by 2 neuroradiologists. Associations between ulcer detection and ulcer characteristics, including orientation, location, and size, were determined and explored by CFD modeling.

RESULTS

One CE-MRA and 3 TOF-MRAs were noninterpretable and excluded. Of 71 ulcers in 56 arteries, readers detected an average of 39 (55%) on both TOF-MRA and CE-MRA, 26.5 (37.5%) only on CE-MRA, and 1 (1.5%) only on TOF-MRA, missing 4.5 (6%) ulcers by both methods. Ulcer detection by TOF-MRA was associated with its orientation (distally pointing versus perpendicular: OR = 5.57 [95% CI, 1.08-28.65]; proximally pointing versus perpendicular: OR = 0.21 [95% CI, 0.14-0.29]); location relative to point of maximum stenosis (distal versus isolevel: OR = 5.17 [95% CI, 2.10-12.70]); and neck-to-depth ratio (OR = 1.96 [95% CI, 1.11-3.45]) after controlling for stenosis and ulcer volume.

CONCLUSIONS

CE-MRA detects more ulcers than TOF-MRA in carotid plaques. Missed ulcers on TOF-MRA are influenced by ulcer orientation, location relative to point of maximum stenosis, and neck-to-depth ratio.

Targeted particle tracking in computational models of human carotid bifurcations

A significant and largely unsolved problem of computational fluid dynamics (CFD) simulation of flow in anatomically relevant geometries is that very few calculated pathlines pass through regions of complex flow. This in turn limits the ability of CFD-based simulations of imaging techniques (such as MRI) to correctly predict in vivo performance. In this work, I present two methods designed to overcome this filling problem, firstly, by releasing additional particles from areas of the flow inlet that lead directly to the complex flow region ("preferential seeding") and, secondly, by tracking particles both "downstream" and "upstream" from seed points within the complex flow region itself. I use the human carotid bifurcation as an example of complex blood flow that is of great clinical interest. Both idealized and healthy volunteer geometries are investigated. With uniform seeding in the inlet plane (in the common carotid artery (CCA)) of an idealized bifurcation geometry, approximately half the particles passed through the internal carotid artery (ICA) and half through the external carotid artery. However, of those particles entering the ICA, only 16% passed directly through the carotid bulb region. Preferential seeding from selected regions of the CCA was able to increase this figure to 47%. In the second method, seeding of particles within the carotid bulb region itself led to a very high proportion (97%) of pathlines running from CCA to ICA. Seeding of particles in the bulb plane of three healthy volunteer carotid bifurcation geometries led to much better filling of the bulb regions than by particles seeded at the inlet alone. In all cases, visualization of the origin and behavior of recirculating particles led to useful insights into the complex flow patterns. Both seeding methods produced significant improvements in filling the carotid bulb region with particle tracks compared with uniform seeding at the inlet and led to an improved understanding of the complex flow patterns. The methods described may be combined and are generally applicable to CFD studies of fluid and gas flow and are, therefore, of relevance in hemodynamics, respiratory mechanics, and medical imaging science.

Fluid Structure Interaction With Contact Surface Methodology for Evaluation of Endovascular Carotid Implants for Drug-Resistant Hypertension Treatment

Drug-resistant hypertensive patients may be treated by mechanical stimulation of stretch-sensitive baroreceptors located in the sinus of carotid arteries. To evaluate the efficacy of endovascular devices to stretch the carotid sinus such that the induced strain might trigger baroreceptors to increase action potential firing rate and thereby reduce systemic blood pressure, numerical simulations were conducted of devices deployed in subject-specific carotid models. Two models were chosen—a typical physiologic carotid and a diminutive atypical physiologic model representing a clinically worst case scenario—to evaluate the effects of device deployment in normal and extreme cases, respectively. Based on the anatomical dimensions of the carotids, two different device sizes were chosen out of five total device sizes available. A fluid structure interaction (FSI) simulation methodology with contact surface between the device and the arterial wall was implemented for resolving the stresses and strains induced by device deployment. Results indicate that device deployment in the carotid sinus of the physiologic model induces an increase of 2.5% and 7.5% in circumferential and longitudinal wall stretch, respectively, and a maximum of 54% increase in von Mises arterial stress at the sinus wall baroreceptor region. The second device, deployed in the diminutive carotid model, induces an increase of 6% in both circumferential and longitudinal stretch and a 50% maximum increase in von Mises stress at the sinus wall baroreceptor region. Device deployment has a minimal effect on blood-flow patterns, indicating that it does not adversely affect carotid bifurcation hemodynamics in the physiologic model. In the smaller carotid model, deployment of the device lowers wall shear stress at sinus by 16% while accelerating flow entering the external carotid artery branch. Our FSI simulations of carotid arteries with deployed device show that the device induces localized increase in wall stretch at the sinus, suggesting that this will activate baroreceptors and subsequently may control hypertension in drug-resistant hypertensive patients, with no consequential deleterious effects on the carotid sinus hemodynamics.

Computational fluid dynamic characterization of carotid bifurcation stenosis in patient-based geometries

Hemodynamic forces play a role in determining endothelial cell (EC) phenotype and influence vascular remodeling. We present a lesion-based computational fluid dynamic (CFD) pilot analysis to understand the complex spatial and temporal hemodynamic changes that prevail in patients with high-grade carotid artery stenosis (CS). High-resolution three-dimensional (3D) rotational angiography datasets were acquired in eight patients, and used to generate computational meshes. CFD analysis was carried out implementing realistic shear-dependent viscosity for blood. The mean wall shear stress (WSS) within the stenosis region was 107 ± 73 dyn/cm(2) rapidly followed by direction reversal and lower oscillating values in the recirculation zone at a mean of 19 ± 14 dyn/cm(2). WSS vectors exhibited complex dynamic directional and amplitude oscillations not seen in healthy segments, along with time-dependent convergence and divergence strips during the cardiac cycle. The spatial gradient of WSS revealed an elevated average magnitude at the throat of the stenosis of 1425 ± 1012 dyn/cm(3). In conclusion, patient-based CFD analysis of CS predicts a complex hemodynamic environment with large spatial WSS variations that occur very rapidly over short distances. Our results improve estimates of the flow changes and forces at the vessel wall in CS and the link between hemodynamic changes and stenosis pathophysiology.

A computational simulation of the effect of hemodilution on oxygen transport in middle cerebral artery vasospasm

Cerebral vasospasm after aneurysmal subarachnoid hemorrhage is a potentially severe sequel. The induction of hypertension, hypervolemia, and hemodilution is advocated for vasospasm, but it is unclear whether hemodilution confers any benefit. A finite element model of oxygen transport in the proximal middle cerebral artery (MCA) was used to evaluate the complex relationship among hematocrit, viscosity, oxygen content, and blood flow in the setting of vasospasm. A single-phase non-Newtonian finite element model based on three-dimensional incompressible Navier-Stokes equations was constructed of the M1 segment. The model was solved at vessel stenoses ranging from 0% to 90% and hematocrit from 0.2 to 0.6. A small area of poststenotic recirculation was seen with mild (30%) stenosis. Poststenotic eddy formation was noted with more severe (60% to 90%) stenosis. Volumetric flow was inversely related to hematocrit at mild stenosis (0% to 30%). With near-complete stenosis (90%), a paradoxical increase in flow was seen with increasing hematocrit. Oxygen transport across the segment was related to hematocrit at all levels of stenosis with increasing oxygen transport despite a reduction in blood flow, suggesting that with clinically significant vasospasm in the MCA, hemodilution does not improve oxygen transport, but to the contrary, that ischemia may be worsened.

Analysis of peripheral artery velocity tracing in a porcine model

BACKGROUND

The aim of the study was to trace the peripheral artery velocity with ultrasound in pigs and provide inference on diagnosis of the type, location and severity of vascular diseases.

MATERIALS AND METHODS

Limb tightening, adrenaline administration and arterial wall pinching were performed independently in six pigs, and then the evolution of the external iliac artery or femoral artery velocity tracing were monitored.

RESULTS

With the increase of the extents of hindlimb tightening, peak systolic velocity (PSV) of ipsilateral external iliac artery turned from 36.33±1.77 cm/s to 59.72±2.67 cm/s, minimum post-principal wave velocity (MPV from 13.68±1.11 cm/s to -7.48±0.82 cm/s, peak diastolic velocity (PDV) from 19.31±0.86 cm/s to 8.98±0.45 cm/s, and, end diastolic velocity (EDV) from 13.2±0.45 cm/s to 0. With the increase of the dose of the epinephrine injection, PSV increased from 36.33±1.77 cm/s to 43.97±2.15 cm/s but then decreased to 35.43±3.01 cm/s, and MPV negatively increased to -23.53±0.82 cm/s after decreasing from 13.68±1.11 cm/s to 0. PDV and EDV gradually decreased to zero. With the increase of the stenosis severity in the abdominal aortic wall pinching, PSV was reduced and had a linearly negative correlation with the stenosis severity (R=0.983, R2=0.967). MPV gradually increased, and its direction reversed when the stenosis severity increased, then diminished when the blood flow was occluded by more than 2/3.

CONCLUSIONS

The formation of peripheral artery velocity is the result of concurrent effects of cardiac ejection, vascular resistance, effective circulating blood volume and elastic recoil. Vascular resistance exerts pronounced effects on the diastolic waveform, and the occurrence of backward wave indicates that the downstream circulation resistance significantly increases.

3-D flow characterization and shear stress in a stenosed carotid artery bifurcation model using stereoscopic PIV technique

The carotid artery bifurcation is a common site of atherosclerosis which is a major leading cause of ischemic stroke. The impact of stenosis in the atherosclerotic carotid artery is to disturb the flow pattern and produce regions with high shear rate, turbulence, and recirculation, which are key hemodynamic factors associated with plaque rupture, clot formation, and embolism. In order to characterize the disturbed flow in the stenosed carotid artery, stereoscopic PIV measurements were performed in a transparent model with 50% stenosis under pulsatile flow conditions. Simulated ECG gating of the flowrate waveform provides external triggering required for volumetric reconstruction of the complex flow patterns. Based on the three-component velocity data in the lumen region, volumetric shear-stress patterns were derived.

Computational simulations and experimental studies of 3D phase-contrast imaging of fluid flow in carotid bifurcation geometries

PURPOSE

To evaluate the use of computational fluid dynamic (CFD)-based magnetic resonance imaging (MRI) simulations to predict the image appearance and velocity measurement of fluid flow in human carotid bifurcation geometries, and to compare the results with images from experimental MRI studies.

MATERIALS AND METHODS

Simulated particle paths were calculated from available CFD datasets of normal and moderately stenosed carotid bifurcation geometries. An MRI simulator based on the spin isochromat method was used to generate images corresponding to a 3D phase-contrast sequence with velocity encoding in three orthogonal directions. The resulting images were compared qualitatively with experimental MRI scans of the corresponding physical models.

RESULTS

The simulations predicted the main features observed in experimental studies, such as the low image intensity in regions of complex flow and the position and bright appearance of the jet in the stenosed bifurcation. Simulated velocity images also agreed well with experimental results. The effects of sequence parameters such as repetition time (TR) and echo time (TE) were readily demonstrated by the simulations.

CONCLUSION

CFD-based MRI simulations can be used to predict the appearance of MRI images of regions of physiological flow, and may be useful in the development of improved pulse sequences for flow measurement.

Introduction to the biomechanics of carotid plaque pathogenesis and rupture: review of the clinical evidence

The management of patients with asymptomatic carotid disease is currently under debate and new methods are warranted for better risk stratification. The role of the biomechanical properties of the atherosclerotic arterial wall together with the effect of different stress types in plaque destabilisation has only been recently investigated. PubMed and Scopus databases were reviewed. There is preliminary clinical evidence demonstrating that the analysis of the combined effect of the various types of biomechanical stress acting on the carotid plaque may help us to identify the vulnerable plaque. At present, MRI and two-dimensional ultrasound are combined with fluid-structure interaction techniques to produce maps of the stress variation within the carotid wall, with increased cost and complexity. Stress wall analysis can be a useful tool for carotid plaque evaluation; however, further research and a multidisciplinary approach are deemed as necessary for further development in this direction.

Flow patterns in carotid bifurcation models using pulsed Doppler ultrasound: effect of concentric vs. eccentric stenosis on turbulence and recirculation

Hemodynamics play a significant role in stroke risk, where thrombus formation may be accelerated in regions of slow or recirculating flow, high shear and increased turbulence. An in vitro investigation was performed with pulsed Doppler ultrasound (DUS) using the complete spectral data to investigate the three-dimensional (3-D) distribution of advanced parameters that may have potential for making a more specific in vivo diagnosis of carotid disease and stroke risk. The effect of stenosis symmetry and the potential of DUS spectral parameters for visualizing regions of recirculation or turbulence were explored. DUS was used to map pulsatile flow in four model geometries representing two different plaque symmetries (eccentricity) and two stenosis severities (mild, severe). Qualitative comparisons were made with flow patterns visualized using digital particle imaging. Color-encoded maps of DUS spectral parameters (mean velocity, spectral-broadening index and turbulence intensity) clearly distinguished regions of slow or recirculating flow and disturbed or turbulent flow. Distinctly different flow patterns resulted from stenoses of equal severity but different eccentricity. Noticeable differences were seen in both the size and location of recirculation zones and in the paths of high-velocity jets. Highly elevated levels of turbulence intensity were seen distal to severe stenosis. Results demonstrated the importance of plaque shape, which is typically not considered in standard diagnosis, in addition to stenosis severity.

In vivo Doppler ultrasound quantification of turbulence intensity using a high-pass frequency filter method

The objective of this investigation was to implement a high-pass frequency filter method to analyze Doppler ultrasound velocity waveforms and quantify turbulence intensity (TI) in vivo. Doppler velocity data were analyzed using two techniques, based on either ensemble averaging or high-pass frequency domain filtering of the periodic waveforms. The accuracy and precision of TI measurements were determined with controlled in vitro experiments, using a pulsatile-flow model of a stenosed carotid bifurcation. The high-pass filter technique was also applied in vivo to determine whether this technique could successfully distinguish between pertinent hemodynamic sites within the carotid artery bifurcation. Twenty-five seconds of Doppler audio data were acquired at three sites (common carotid artery [CCA], internal carotid artery [ICA] stenosis and distal ICA) within 10 human carotid arteries, and repeated three times. Doppler velocity data were analyzed using a ninth-order high-pass Butterworth filter with a 12-Hz inflection point. TI measured within the CCA and distal ICA was found to be significantly different (p < 0.0001) for moderate to nearly occluded carotid artery classifications. Also, TI measured within the distal ICA increased with stenosis severity, with the ability to distinguish between each stenosis class (p < 0.05). This investigation demonstrated the ability to precisely quantify TI using a conventional Doppler ultrasound machine in human subjects, without interfering with normal clinical protocols.

Analyzing transient turbulence in a stenosed carotid artery by proper orthogonal decomposition

High-resolution three-dimensional simulations (involving 100 million degrees of freedom) were employed to study transient turbulent flow in a carotid arterial bifurcation with a stenosed internal carotid artery (ICA). The geometrical model was reconstructed from MRI images, and in vivo velocity measurements were incorporated in the simulations to provide inlet and outlet boundary conditions. Due to the high degree of the ICA occlusion and the variable flow rate, a transitional and intermittent flow between laminar and turbulent states was established. Time- and space-window proper orthogonal decomposition (POD) was applied to quantify the different flow regimes in the occluded artery. A simplified version of the POD analysis that utilizes 2D slices only--more appropriate in the clinical setting--was also investigated.

Ultrasound simulation of complex flow velocity fields based on computational fluid dynamics

In this work, a simulation environment for the development of flow-related ultrasound algorithms is presented. Ultrasound simulations of realistic Doppler signals require accurate modeling of blood flow. Instead of using analytically described flow behavior, complex blood movement can be derived from velocity fields obtained with computational fluid dynamics (CFD). By further modeling blood as a collection of point scatterers, resulting RF-signals can be efficiently retrieved using an existing ultrasound simulation model. The main aim of this paper is to elaborate on creating CFDbased phantoms for ultrasound simulations. The coupling of a computed flow field with an ultrasound model offers flexible control of flow and ultrasound imaging parameters, beneficial for improving and developing imaging algorithms. The proposed method was validated in a straight tube with a stationary parabolic velocity profile and further demonstrated by an eccentrically stenosis carotid bifurcation. The estimated flow velocities are in good agreement with the CFD reference, both for color flow imaging and pulsed-wave doppler simulations. The presented method can also be extended to include wall mechanics simulations in future work.

In vitro Doppler ultrasound investigation of turbulence intensity in pulsatile flow with simulated cardiac variability

An in vitro investigation of turbulence intensity (TI) associated with a severe carotid stenosis in the presence of physiological cardiac variability is described. The objective of this investigation was to determine if fluctuations due to turbulence could be quantified with conventional Doppler ultrasound (DUS) in the presence of normal physiological cycle-to-cycle cardiac variability. An anthropomorphic model of a 70% stenosed carotid bifurcation was used in combination with a programmable flow pump to generate pulsatile flow with a mean flow rate of 6 mL/s. Utilizing the pump, we studied normal, nonrepetitive cycle-to-cycle cardiac variability (+/-3.9%) in flow, as well as waveform shapes with standard deviations equal to 0, 2 and 3 times the normal variation. Eighty cardiac cycles of Doppler data were acquired at two regions within the model, representing either laminar or turbulent flow; each measurement was repeated six times. Turbulence intensity values were found to be 11 times higher (p < 0.001), on average, in the turbulent region than in the laminar region, with a mean difference of 24 cm/s. Twenty cardiac cycles were required for confidence in TI values. In conclusion, these results indicate that it is possible to quantify TI in vitro, even in the presence of normal and exaggerated cycle-to-cycle cardiac variability.

Doppler ultrasound compatible plastic material for use in rigid flow models

A technique for the rapid but accurate fabrication of multiple flow phantoms with variations in vascular geometry would be desirable in the investigation of carotid atherosclerosis. This study demonstrates the feasibility and efficacy of implementing numerically controlled direct-machining of vascular geometries into Doppler ultrasound (DUS)-compatible plastic for the easy fabrication of DUS flow phantoms. Candidate plastics were tested for longitudinal speed of sound (SoS) and acoustic attenuation at the diagnostic frequency of 5 MHz. Teflon was found to have the most appropriate SoS (1376 +/- 40 m s(-1) compared with 1540 m s(-1) in soft tissue) and thus was selected to construct a carotid bifurcation flow model with moderate eccentric stenosis. The vessel geometry was machined directly into Teflon using a numerically controlled milling technique. Geometric accuracy of the phantom lumen was verified using nondestructive micro-computed tomography. Although Teflon displayed a higher attenuation coefficient than other tested materials, Doppler data acquired in the Teflon flow model indicated that sufficient signal power was delivered throughout the depth of the vessel and provided comparable velocity profiles to that obtained in the tissue-mimicking phantom. Our results indicate that Teflon provides the best combination of machinability and DUS compatibility, making it an appropriate choice for the fabrication of rigid DUS flow models using a direct-machining method.

Eccentric stenosis of the carotid artery associated with ipsilateral cerebrovascular events

BACKGROUND AND PURPOSE

Eccentric stenosis of the coronary artery is associated with plaque disruption and acute coronary syndrome. The purpose of the present study was to determine whether eccentric stenosis of the carotid artery contributes to cerebrovascular events.

MATERIALS AND METHODS

Of 6859 patients with vascular diseases who underwent duplex carotid ultrasonography, we studied 512 internal carotid arteries in 441 patients who had a maximum area stenosis at or more than 70%, which corresponds with approximately 50% or more by the NASCET method. The maximal (A) and minimal wall thicknesses (B) were measured on cross-sectional sonography images, and an eccentricity index was calculated using the following formula: (A - B)/A. Arteries in the lowest quartile of the eccentricity index (<0.69) were defined as having a concentric stenosis, whereas the others were defined as having eccentric stenosis. The underlying clinical characteristics and plaque morphologies, as well as the occurrence of ipsilateral ischemic stroke or transient ischemic attack in the preceding year, were compared between patients with eccentric and concentric stenosis.

RESULTS

Patient characteristics and plaque morphology were similar between the 2 groups. Cerebrovascular events occurred more frequently ipsilaterally to the artery with eccentric stenosis (13.5%) than to the artery with concentric stenosis (5.5%; P = .013); the difference was more evident when cerebrovascular events of presumed carotid arterial origin were assessed (P = .005). After adjusting for risk factors and plaque morphology, eccentric stenosis was independently related to the presence of recent cerebrovascular events (odds ratio = 2.76; 95% confidence interval = 1.19-6.40).

CONCLUSIONS

In patients with an area carotid stenosis of 70% or more, eccentric plaque was associated with a significantly increased incidence of ipsilateral cerebrovascular events compared with patients with concentric stenosis.

PREDICTION OF COMPLEX FLOW PATTERNS IN INTRACRANIAL ATHEROSCLEROTIC DISEASE USING COMPUTATIONAL FLUID DYNAMICS

OBJECTIVE

Although carotid and vertebral intracranial atherosclerotic disease (ICAD) can lead to both hemodynamic insufficiency and thromboembolism, its fluid dynamic properties remain undefined because of its intricate features and complex three-dimensional geometry. We used computational fluid dynamic (CFD) analysis to model the hemodynamics of symptomatic ICAD lesions.

METHODS

Nine ICAD lesions (six carotid, two vertebral, one middle cerebral) underwent high-resolution catheter-based digital rotational angiography. The reconstructed three-dimensional volumes of the target lesions were segmented and used to generate hybrid computational meshes. Dynamic pulsatile CFD analysis was performed using a non-Newtonian shear-dependent model of blood's viscosity.

RESULTS

CFD results revealed complex flow patterns within ICAD lesions with midstenotic shear rates of greater than 19,000/s, sufficiently high to induce high-shear platelet activation. Vorticity and helicity within the stenoses were followed by sudden deceleration with formation of vortex cores. Pressure gradients were significant mostly at greater than 75% stenosis with a mean time-averaged drop of 27.2 ±17.8 mmHg. Unlike the smoothly-varying helicity imparted by the three-dimensional anatomy of the intracranial circulation, poststenotic regions of ICAD lesions showed significant and rapidly fluctuating helicity and vorticity patterns, which may contribute to the propagation of platelets activated by the high shear region within the stenosis throat. Stent angioplasty restored the hemodynamic profile of ICAD lesions to within contralateral controls.

CONCLUSION

Patient-based symptomatic ICAD lesions studied using CFD analysis appear to harbor a hemodynamically pathological environment that favors the activation, aggregation and distal embolization of platelets and is reversed by endovascular stent angioplasty.

Turbulence modeling in three-dimensional stenosed arterial bifurcations

Under normal healthy conditions, blood flow in the carotid artery bifurcation is laminar. However, in the presence of a stenosis, the flow can become turbulent at the higher Reynolds numbers during systole. There is growing consensus that the transitional k-omega model is the best suited Reynolds averaged turbulence model for such flows. Further confirmation of this opinion is presented here by a comparison with the RNG k-epsilon model for the flow through a straight, nonbifurcating tube. Unlike similar validation studies elsewhere, no assumptions are made about the inlet profile since the full length of the experimental tube is simulated. Additionally, variations in the inflow turbulence quantities are shown to have no noticeable affect on downstream turbulence intensity, turbulent viscosity, or velocity in the k-epsilon model, whereas the velocity profiles in the transitional k-omega model show some differences due to large variations in the downstream turbulence quantities. Following this validation study, the transitional k-omega model is applied in a three-dimensional parametrically defined computer model of the carotid artery bifurcation in which the sinus bulb is manipulated to produce mild, moderate, and severe stenosis. The parametric geometry definition facilitates a powerful means for investigating the effect of local shape variation while keeping the global shape fixed. While turbulence levels are generally low in all cases considered, the mild stenosis model produces higher levels of turbulent viscosity and this is linked to relatively high values of turbulent kinetic energy and low values of the specific dissipation rate. The severe stenosis model displays stronger recirculation in the flow field with higher values of vorticity, helicity, and negative wall shear stress. The mild and moderate stenosis configurations produce similar lower levels of vorticity and helicity.

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

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

Real-time numerical simulation of Doppler ultrasound in the presence of nonaxial flow

Numerical simulations of Doppler ultrasound (DUS) relying on computational fluid dynamics (CFD) models of nonaxial flow have traditionally employed detailed (but computationally intensive) models of the DUS physics, or have sacrificed much of the physics in the interest of computational or conceptual simplicity. In this paper, we present a compromise between these extremes, with the objective of simulating the essential characteristics of DUS spectrograms in a real-time manner. Specifically, a precomputed pulsatile CFD velocity field is interrogated at some number, N, of discrete points distributed spatially within a sample volume of prescribed geometry and power distribution and temporally within a prescribed sampling window. Intrinsic spectral broadening is accounted for by convolving each of the point velocities with a semiempirical broadening function. Real-time performance is facilitated through the use of an efficient algorithm for interpolating the unstructured CFD data. A spherical sample volume with Gaussian power distribution, N = 1000 sampling points, and quadratic broadening function are shown to be adequate for simulating, at frame rates of 86 Hz on a 1.5 GHz desktop workstation, realistic-looking spectrograms at representative locations within a stenosed carotid bifurcation model. Via qualitative comparisons with matched in vitro data, these simulated spectrograms are shown to mimic the distinctive spectral envelopes, broadening and power characteristics associated with common carotid, stenotic jet and poststenotic recirculating flows. We conclude that the complex interaction between Doppler ultrasound and complicated clinically relevant blood flow dynamics can be simulated in real time via this relatively straightforward semiempirical approach.