Time-efficient patient-specific quantification of regional carotid artery fluid dynamics and spatial correlation with plaque burden

Numerical study on the cerebral blood flow regulation in the circle of Willis with the vascular absence and internal carotid artery stenosis

Objectives

This study explores how vascular stenosis and absence affect the regulation of cerebral blood flow in the Circle of Willis (CoW) and the hemodynamic changes downstream of the stenosis.

Materials and Methods

Forty idealized CoW models were simulated to analyze the impact of vascular absence and internal carotid artery (ICA) stenosis on hemodynamics. Inlet conditions were set using a physiological pressure waveform, and outflow boundaries were modeled using three-element Windkessel models.

Results

The absence of vessels such as RP1, LP1, RA1, or LA1 had a comparable effect on total blood flow to a 40% stenosis of the left internal carotid artery (LICA) across the entire CoW. Specifically, when LP1 and RA1 were absent with a 50% LICA stenosis, the total blood flow closely resembled that of a complete CoW with 75% LICA stenosis. In cases of proximal ICA stenosis, downstream regions showed elevated oscillatory shear index (OSI >0.2) and reduced time-averaged wall shear stress (TAWSS <1 Pa). With increasing stenosis severity, areas of high OSI shifted, and regions of low TAWSS expanded notably. At 75% stenosis, the area with TAWSS <1 Pa downstream significantly increased. Until complete occlusion, the area of low TAWSS and high OSI were maximized.

Conclusion

This study underscores how anatomical variations in the CoW, combined with ICA stenosis, impact both total cerebral blood flow and its distribution among different outlets. Moreover, it highlights the potential for increased atherosclerosis development in affected areas. Particularly notable is the finding the absence of LP1 and RA1 vessels alongside 50% LICA stenosis results in blood flow patterns similar to those seen with 75% LICA stenosis in a complete CoW, emphasizing clinical implications for the patient. Hemodynamic changes, including TAWSS and OSI, are most pronounced downstream of the stenosis especially when the stenosis rate exceeds 75%.

Quantitative evaluation using carotid ultrasonography-based high-frame-rate vector flow imaging in patients with low carotid stenosis

OBJECTIVE

To explore the role of quantitative evaluation using carotid ultrasonography (US)-based high-frame-rate vector flow (V Flow) imaging in patients with low carotid stenosis.

METHODS

This single-centre cross-sectional study consecutively recruited volunteers without carotid plaque and patients with low carotid stenosis from August 2022 to May 2023. Patients were divided into symptomatic and asymptomatic groups according to their head CT or MRI results within 8 weeks. All V Flow imaging examinations were performed using a Mindray Resona R9 US system. The wall shear stress (WSS) values, oscillatory shear index (OSI) values, and turbulence (Tur) indexes in the normal common carotid artery (CCA), normal carotid bifurcation (CB), and on the upstream and downstream surface of carotid plaque were measured. Pearson Chi-square test and Fisher exact test were used for counting data according to their type. For measurement data, independent sample t test and non-parametric rank sum test were used.

RESULTS

The results proved that patients have higher WSS values and Tur indexes of CB than volunteers, and higher WSS values were detected on the surface of the plaques in symptomatic patients. What's more, the downstream side of the plaque was more vulnerable to plaque rupture than the upstream side due to more dynamic blood flow.

CONCLUSION

Therefore, carotid US-based high-frame-rate V Flow imaging provides reliable mechanical biomarkers for assessing the haemodynamic change in patients with low stenosis. Our study may provide a new imaging tool for monitoring the progression of atherosclerosis and aiding the management of early atherosclerotic patients.

ADVANCES IN KNOWLEDGE

Our study firstly investigated the difference of V Flow parameters on the surface of carotid plaques between symptomatic and asymptomatic patients with low carotid stenosis, which is expected to provide haemodynamic information and the mechanical basis for plaque rupture.

Numerical study of hemodynamic changes in the Circle of Willis after stenosis of the internal carotid artery

BACKGROUND AND OBJECTIVES

In clinical practice a large number of patients with ischemic stroke have internal carotid artery (ICA) stenosis accompanied by Circle of Willis (CoW) stenosis. In the presence of carotid artery stenosis, CoW atherosclerosis may cause cerebral blood flow decompensation and may promote the development of ischemic stroke. The reason for the concomitant stenosis at both sites is unknown. This study investigated the hemodynamic effects of ICA stenosis on the CoW.

METHODS

We developed a three-dimensional/zero-dimensional (3D/0D) closed-loop geometric multiscale model of the cerebral artery to quantify the hemodynamic indicators, including time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI). Mild (<50 %), moderate (50-69 %) and severe (>69 %) ICA stenoses were established based on 3D models of cerebral arteries in two volunteers. Geometric multiscale computing models were numerically evaluated to obtain local hemodynamic changes in the CoW in order to assess the risk of stenosis in the CoW.

RESULTS

Model calculations showed that for all 3D models the A1 segment of the anterior cerebral artery (ACA) or the posterior communicating artery (PCA) within the CoW exhibited a hemodynamic environment with high OSI (>0.2) and low TAWSS (<1 Pa) when the ICA had a moderate stenosis. While in the case of mild and severe stenosis in ICA, there is no such phenomenon. The proportion of the surface area possessing high OSI and low TAWSS in the A1 segment of the ACA or in the PCA was mostly greater than 60 %, which might potentially cause the formation and development of atherosclerosis in CoW and finally lead to CoW stenosis.

CONCLUSIONS

Therefore, although moderate carotid artery stenosis may not cause ischemic stroke, it may cause hemodynamic changes in the CoW, which in turn may promote CoW stenosis and cause CoW decompensation. In clinical treatment attention should be paid not only to stenosis of the carotid arteries but also to changes in the hemodynamic environment within the CoW, in order to prevent the adverse effects of CoW stenosis.

A workflow for viewing biomedical computational fluid dynamics results and corresponding data within virtual and augmented reality environments

Researchers conducting computational fluid dynamics (CFD) modeling can spend weeks obtaining imaging data, determining boundary conditions, running simulations and post-processing files. However, results are typically viewed on a 2D display and often at one point in time thus reducing the dynamic and inherently three-dimensional data to a static image. Results from different pathologic states or cases are rarely compared in real-time, and supplementary data are seldom included. Therefore, only a fraction of CFD results are typically studied in detail, and associations between mechanical stimuli and biological response may be overlooked. Virtual and augmented reality facilitate stereoscopic viewing that may foster extraction of more information from CFD results by taking advantage of improved depth cues, as well as custom content development and interactivity, all within an immersive approach. Our objective was to develop a straightforward, semi-automated workflow for enhanced viewing of CFD results and associated data in an immersive virtual environment (IVE). The workflow supports common CFD software and has been successfully completed by novice users in about an hour, demonstrating its ease of use. Moreover, its utility is demonstrated across clinical research areas and IVE platforms spanning a range of cost and development considerations. We are optimistic that this advancement, which decreases and simplifies the steps to facilitate more widespread use of immersive CFD viewing, will foster more efficient collaboration between engineers and clinicians. Initial clinical feedback is presented, and instructional videos, manuals, templates and sample data are provided online to facilitate adoption by the community.

Computer Modeling of Carotid Endarterectomy With the Different Shape Patches and Prediction of the Atherosclerotic Plaque Formation Zones

The article describes a method for constructing geometric models of the carotid bifurcation and computer simulation of endarterectomy surgery with the patches of various configurations. The purpose of this work is to identify the areas of the greatest risk of restenosis in the constructed models and to conduct a comparative analysis of risk factors when using the patches of different widths and shapes. The method is demonstrated on a reconstructed model of a healthy vessel. Its building is based on a preoperative computed tomography study of a particular patient's affected vessel. The flow in the vessel is simulated by computational fluid dynamics using data from the patient's ultrasound Doppler velocimetry. Risk factors are assessed through the hemodynamic indices on the vessel wall associated with Wall Shear Stress. The distribution of risk zones in the healthy vessel, presumably leading to its observed lesion (plaque), is analyzed. Comparative evaluation of 10 various patches implantation results is carried out and the optimal variant is determined. The proposed method can be used to predict the hemodynamic results of surgery using patches of various sizes and shapes.

Slice-based and time-specific hemodynamic measurements discriminate carotid artery vulnerable atherosclerotic plaques

BACKGROUND AND OBJECTIVE

Hemodynamic patterns play key roles in progression of carotid vulnerable plaques. However, most of previous studies utilized maximum or averaged value of hemodynamic measurements which is not an ideal representative of hemodynamic patterns. This study aimed to investigate the association of slice-based and time-specific hemodynamic measurements with carotid vulnerable plaque using magnetic resonance (MR) vessel wall imaging and histology.

METHODS

Thirty-two patients (mean age: 63.9±8.1 years; 25 males) with carotid atherosclerotic stenosis (≥50% stenosis) referred to carotid endarterectomy were recruited and underwent MR vessel wall imaging. Carotid plaque burden was evaluated on MR images and vulnerable plaque features including calcification, lipid-rich necrotic core, and intra-plaque hemorrhage (IPH) were identified by histology. The slice-based and time-specific hemodynamic measurements were extracted from computational fluid dynamics simulation of 3D carotid arterial model. Correlation coefficients between hemodynamic measurements and carotid plaque features were calculated and the logistic regressions with generalized estimating equation (GEE) were conducted. The value in discriminating carotid vulnerable plaque features was determined by receiver-operating-characteristic analysis.

RESULTS

Of 102 MR-histology matched slices from 32 patients, time-averaged wall shear stress (TAWSS) (r=0.263, p=0.008), oscillatory shear index (OSI) (r=-0.374, p<0.001), and peakWSS (r=0.232, p=0.019) were significantly associated with carotid IPH. The logistic regression with GEE revealed that peakWSS (OR, 1.206; 95% CI, 1.026-1.418; p, 0.023) and TAWSS (OR, 0.364, 95% CI, 0.138-0.959; p, 0.041) were significantly associated with presence of IPH after adjusting for age and BMI. In discriminating carotid IPH, the AUC of TAWSS, OSI, combined TAWSS with maximum wall thickness (MWT) and combined OSI with MWT was 0.656, 0.722, 0.761, and 0.764, respectively.

CONCLUSIONS

Slice-based and time-specific hemodynamic characteristics could effectively discriminate carotid IPH. Combination of hemodynamic measurements with carotid plaque burden might be a stronger indicator for carotid vulnerable plaque features than each measurement alone.

Model-based evaluation of local hemodynamic effects of enhanced external counterpulsation

BACKGROUND AND OBJECTIVES

The treatment benefits of enhanced external counterpulsation (EECP) heavily depends on hemodynamics. Global hemodynamics of EECP can cause blood flow redistribution in the circulatory system whereas local hemodynamic effects act on vascular endothelial cells (VECs). Local hemodynamic effects of EECP on VECs are important in the treatment of atherosclerosis, but currently cannot be not evaluated. Herein we aim to establish evaluation models of local hemodynamic effects based on the global hemodynamic indicators.

METHODS

We established 0D/3D geometric multi-scale hemodynamic models of the coronary and cerebral artery of two healthy individuals to calculate the global hemodynamic indicators and the local hemodynamic effects. Clinical EECP trials were performed to verify the accuracy of the multi-scale hemodynamic model. The global hemodynamic indicators included diastolic blood pressure/systolic blood pressure (Q = D/S), mean arterial pressure (MAP), internal carotid artery flow (ICAF) and cerebral blood flow (CBF), whereas local hemodynamic effects focused on time-averaged wall shear stress (TAWSS). The correlation between these indicators was analyzed via Pearson correlation coefficient. Significantly related indicators were selected for curve-fitting to establish evaluation models of the coronary and cerebral artery. Moreover, clinical data of a coronary heart disease patient and a cerebral ischemic stroke patient were collected to verify the effectiveness of the application of the established evaluation models to real patients.

RESULTS

For coronary artery, TAWSS was correlated to Q = D/S and ICAF (P < 0.05), whereas for cerebral artery, TAWSS was correlated to MAP and CBF (P < 0.05). The mean square error (MSE) between the evaluated values using evaluation model and the calculated values using 0D/3D model of TAWSS of the coronary and cerebral artery were 5.4% and 1.0%, respectively. The MSE of evaluation model applied to real patients was greater than that applied to healthy individuals, but within an acceptable range.

CONCLUSIONS

The presented error demonstrated validity and accuracy of the evaluation models in clinical patients. Based on the evaluation models, global hemodynamic indicators could be used to evaluate the local hemodynamic effects under the current counterpulsation mode. With TAWSS range of 4-7 Pa as the target range, EECP strategies can further be optimized.

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

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

A Pilot Study Characterizing Flow Patterns in the Thoracic Aorta of Patients With Connective Tissue Disease: Comparison to Age- and Gender-Matched Controls via Fluid Structure Interaction

Prior computational and imaging studies described changes in flow patterns for patients with Marfan syndrome, but studies are lacking for related populations. This pilot study addresses this void by characterizing wall shear stress (WSS) indices for patients with Loeys-Dietz and undifferentiated connective tissue diseases. Using aortic valve-based velocity profiles from magnetic resonance imaging as input to patient-specific fluid structure interaction (FSI) models, we determined local flow patterns throughout the aorta for four patients with various connective tissue diseases (Loeys-Dietz with the native aorta, connective tissue disease of unclear etiology with native aorta in female and male patients, and an untreated patient with Marfan syndrome, as well as twin patients with Marfan syndrome who underwent valve-sparing root replacement). FSI simulations used physiological boundary conditions and material properties to replicate available measurements. Time-averaged WSS (TAWSS) and oscillatory shear index (OSI) results are presented with localized comparison to age- and gender-matched control participants. Ascending aortic dimensions were greater in almost all patients with connective tissue diseases relative to their respective control. Differences in TAWSS and OSI were driven by local morphological differences and cardiac output. For example, the model for one twin had a more pronounced proximal descending aorta in the vicinity of the ductus ligamentum that impacted WSS indices relative to the other. We are optimistic that the results of this study can serve as a foundation for larger future studies on the connective tissue disorders presented in this article.

A Numerical Model for Simulating the Hemodynamic Effects of Enhanced External Counterpulsation on Coronary Arteries

Traditional enhanced external counterpulsation (EECP) used for the clinical treatment of patients with coronary heart disease only assesses diastolic/systolic blood pressure (Q = D/S > 1.2). However, improvement of the hemodynamic environment surrounding vascular endothelial cells of coronary arteries after long-term application of EECP is the basis of the treatment. Currently, the quantitative hemodynamic mechanism is not well understood. In this study, a standard 0D/3D geometric multi-scale model of the coronary artery was established to simulate the hemodynamic effects of different counterpulsation modes on the vascular endothelium. In this model, the neural regulation caused by counterpulsation was thoroughly considered. Two clinical trials were carried out to verify the numerical calculation model. The results demonstrated that the increase in counterpulsation pressure amplitude and pressurization duration increased coronary blood perfusion and wall shear stress (WSS) and reduced the oscillatory shear index (OSI) of the vascular wall. However, the impact of pressurization duration was the predominant factor. The results of the standard model and the two real individual models indicated that a long pressurization duration would cause more hemodynamic risk areas by resulting in excessive WSS, which could not be reflected by the change in the Q value. Therefore, long-term pressurization during each cardiac cycle therapy is not recommended for patients with coronary heart disease and clinical treatment should not just pay attention to the change in the Q value. Additional physiological indicators can be used to evaluate the effects of counterpulsation treatment.

Long-term hemodynamic mechanism of enhanced external counterpulsation in the treatment of coronary heart disease: a geometric multiscale simulation

Enhanced external counterpulsation (EECP) is a noninvasive treatment method for coronary artery atherosclerosis that acts on the vascular endothelial cells. The intracoronary hemodynamic parameters that influence long-term treatment effect are the fundamental factors for the inhibition of intimal hyperplasia, which cannot be measured in real time. In order to optimize the long-term treatment effect of coronary heart disease, it is necessary to establish a method for quantified calculation of intracoronary hemodynamic parameters during counterpulsation to research the long-term hemodynamic mechanism of EECP. A geometric multiscale model coupled by the zero-dimensional (0D) lumped parameter model and the three-dimensional (3D) model of narrow coronary artery was established for the simulation of intracoronary hemodynamic environment. The 3D model was used to calculate the hemodynamic parameters such as wall shear stress (WSS) and oscillatory shear index (OSI), while the 0D model was used to simulate the blood circulatory system. Sequential pressure was applied to calves, thighs, and buttocks module in 0D model with the consideration of vessel collapse. Hemodynamic performance was compared with clinical reports to verify the effectiveness of the method. There were significant increases of the diastolic blood pressure (DBP), coronary flow, and the area-averaged WSS during application of EECP, while OSI behind stenosis has some decrease. The waveforms of coronary flow has good similarity with the clinical measured waveforms, and the differences between calculated mean arterial pressures (MAPs) and clinical measurements were within 1%. The fundamental factor in the cure of coronary heart disease by EECP is the improvement of WSS and the decrease of OSI. Comparing with the clinical reports, the immediate hemodynamic changes demonstrate the effectiveness of model. Intracoronary hemodynamic parameters during EECP could be acquired and the method could be used to simulate the long-term treatment effect of EECP. Graphical abstract.

Hemodynamic effects of enhanced external counterpulsation on cerebral arteries: a multiscale study

Background

Enhanced external counterpulsation (EECP) is an effective method for treating patients with cerebral ischemic stroke, while hemodynamics is the major contributing factor in the treatment of EECP. Different counterpulsation modes have the potential to lead to different acute and long-term hemodynamic changes, resulting in different treatment effects. However, various questions about appropriate counterpulsation modes for optimizing hemodynamic effects remain unanswered in clinical treatment.

Methods

A zero-dimensional/three-dimensional (0D/3D) geometric multiscale model of the cerebral artery was established to obtain acute hemodynamic indicators, including mean arterial pressure (MAP) and cerebral blood flow (CBF), as well as localized hemodynamic details for the cerebral artery, which includes wall shear stress (WSS) and oscillatory shear index (OSI). Counterpulsation was achieved by applying pressure on calf, thigh and buttock modules in the 0D model. Different counterpulsation modes including various pressure amplitudes and pressurization durations were applied to investigate hemodynamic responses, which impact acute and long-term treatment effects. Both vascular collapse and cerebral autoregulation were considered during counterpulsation.

Results

Variations of pressure amplitude and pressurization duration have different impacts on hemodynamic effects during EECP treatment. There were small differences in the hemodynamics when similar or different pressure amplitudes were applied to calves, thighs and buttocks. When increasing pressure amplitude was applied to the three body parts, MAP and CBF improved slightly. When pressure amplitude exceeded 200 mmHg, hemodynamic indicators almost never changed, demonstrating consistency with clinical data. However, hemodynamic indicators improved significantly with increasing pressurization duration. For pressurization durations of 0.5, 0.6 and 0.7 s, percentage increases for MAP during counterpulsation were 1.5%, 23.5% and 39.0%, for CBF were 1.2%, 23.4% and 41.6% and for time-averaged WSS were 0.2%, 43.5% and 85.0%, respectively.

Conclusions

When EECP was applied to patients with cerebral ischemic stroke, pressure amplitude applied to the three parts may remain the same. Patients may not gain much more benefit from EECP treatment by excessively increasing pressure amplitude above 200 mmHg. However, during clinical procedures, pressurization duration could be increased to 0.7 s during the cardiac circle to optimize the hemodynamics for possible superior treatment outcomes.

The Impact of Cardiac Motion on Aortic Valve Flow Used in Computational Simulations of the Thoracic Aorta

Advancements in image-based computational modeling are producing increasingly more realistic representations of vasculature and hemodynamics, but so far have not compensated for cardiac motion when imposing inflow boundary conditions. The effect of cardiac motion on aortic flow is important when assessing sequelae in this region including coarctation of the aorta (CoA) or regurgitant fraction. The objective of this investigation was to develop a method to assess and correct for the influence of cardiac motion on blood flow measurements through the aortic valve (AoV) and to determine its impact on patient-specific local hemodynamics quantified by computational fluid dynamics (CFD). A motion-compensated inflow waveform was imposed into the CFD model of a patient with repaired CoA that accounted for the distance traveled by the basal plane during the cardiac cycle. Time-averaged wall shear stress (TAWSS) and turbulent kinetic energy (TKE) values were compared with CFD results of the same patient using the original waveform. Cardiac motion resulted in underestimation of flow during systole and overestimation during diastole. Influences of inflow waveforms on TAWSS were greatest along the outer wall of the ascending aorta (AscAo) (∼30 dyn/cm2). Differences in TAWSS were more pronounced than those from the model creation or mesh dependence aspects of CFD. TKE was slightly higher for the motion-compensated waveform throughout the aortic arch. These results suggest that accounting for cardiac motion when quantifying blood flow through the AoV can lead to different conclusions for hemodynamic indices, which may be important if these results are ultimately used to predict patient outcomes.

Image-based quantification of 3D morphology for bifurcations in the left coronary artery: Application to stent design

BACKGROUND

Improved strategies for stent-based treatment of coronary artery disease at bifurcations require a greater understanding of artery morphology.

OBJECTIVE

We developed a workflow to quantify morphology in the left main coronary (LMCA), left anterior descending (LAD), and left circumflex (LCX) artery bifurcations.

METHODS

Computational models of each bifurcation were created for 55 patients using computed tomography images in 3D segmentation software. Metrics including cross-sectional area, length, eccentricity, taper, curvature, planarity, branching law parameters, and bifurcation angles were assessed using open-sources software and custom applications. Geometric characterization was performed by comparison of means, correlation, and linear discriminant analysis (LDA).

RESULTS

Differences between metrics suggest dedicated or multistent approaches should be tailored for each bifurcation. For example, the side branch of the LCX (i.e., obtuse marginal; OM) was longer than that of the LMCA (i.e., LCXprox) and LAD (i.e., first diagonal; D1). Bifurcation metrics for some locations (e.g., LMCA Finet ratio) provide results and confidence intervals agreeing with prior findings, while revised metric values are presented for others (e.g., LAD and LCX). LDA revealed several metrics that differentiate between artery locations (e.g., LMCA vs. D1, LMCA vs. OM, LADprox vs. D1, and LCXprox vs. D1).

CONCLUSIONS

These results provide a foundation for elucidating common parameters from healthy coronary arteries and could be leveraged in the future for treating diseased arteries. Collectively the current results may ultimately be used for design iterations that improve outcomes following implantation of future dedicated bifurcation stents. © 2015 Wiley Periodicals, Inc.

MRI-based biomechanical parameters for carotid artery plaque vulnerability assessment

Carotid atherosclerotic plaques are a major cause of ischaemic stroke. The biomechanical environment to which the arterial wall and plaque is subjected to plays an important role in the initiation, progression and rupture of carotid plaques. MRI is frequently used to characterize the morphology of a carotid plaque, but new developments in MRI enable more functional assessment of carotid plaques. In this review, MRI based biomechanical parameters are evaluated on their current status, clinical applicability, and future developments. Blood flow related biomechanical parameters, including endothelial wall shear stress and oscillatory shear index, have been shown to be related to plaque formation. Deriving these parameters directly from MRI flow measurements is feasible and has great potential for future carotid plaque development prediction. Blood pressure induced stresses in a plaque may exceed the tissue strength, potentially leading to plaque rupture. Multi-contrast MRI based stress calculations in combination with tissue strength assessment based on MRI inflammation imaging may provide a plaque stress-strength balance that can be used to assess the plaque rupture risk potential. Direct plaque strain analysis based on dynamic MRI is already able to identify local plaque displacement during the cardiac cycle. However, clinical evidence linking MRI strain to plaque vulnerability is still lacking. MRI based biomechanical parameters may lead to improved assessment of carotid plaque development and rupture risk. However, better MRI systems and faster sequences are required to improve the spatial and temporal resolution, as well as increase the image contrast and signal-to-noise ratio.

Immersive visualization for enhanced computational fluid dynamics analysis

Modern biomedical computer simulations produce spatiotemporal results that are often viewed at a single point in time on standard 2D displays. An immersive visualization environment (IVE) with 3D stereoscopic capability can mitigate some shortcomings of 2D displays via improved depth cues and active movement to further appreciate the spatial localization of imaging data with temporal computational fluid dynamics (CFD) results. We present a semi-automatic workflow for the import, processing, rendering, and stereoscopic visualization of high resolution, patient-specific imaging data, and CFD results in an IVE. Versatility of the workflow is highlighted with current clinical sequelae known to be influenced by adverse hemodynamics to illustrate potential clinical utility.

Comparison of morphological and rheological conditions between conventional and eversion carotid endarterectomy using computational fluid dynamics – a pilot study

Purpose:

To compare postoperative morphological and rheological conditions after eversion carotid endarterectomy versus conventional carotid endarterectomy using computational fluid dynamics.

Basic methods:

Hemodynamic metrics (velocity, wall shear stress, time-averaged wall shear stress and temporal gradient wall shear stress) in the carotid arteries were simulated in one patient after conventional carotid endarterectomy and one patient after eversion carotid endarterectomy by computational fluid dynamics analysis based on patient specific data.

Principal findings:

Systolic peak of the eversion carotid endarterectomy model showed a gradually decreased pressure along the stream path, the conventional carotid endarterectomy model revealed high pressure (about 180 Pa) at the carotid bulb. Regions of low wall shear stress in the conventional carotid endarterectomy model were much larger than that in the eversion carotid endarterectomy model and with lower time-averaged wall shear stress values (conventional carotid endarterectomy: 0.03–5.46 Pa vs. eversion carotid endarterectomy: 0.12–5.22 Pa).

Conclusions:

Computational fluid dynamics after conventional carotid endarterectomy and eversion carotid endarterectomy disclosed differences in hemodynamic patterns. Larger studies are necessary to assess whether these differences are consistent and might explain different rates of restenosis in both techniques.

Computational fluid dynamics in coronary artery disease

Computational fluid dynamics (CFD) is a widely used method in mechanical engineering to solve complex problems by analysing fluid flow, heat transfer, and associated phenomena by using computer simulations. In recent years, CFD has been increasingly used in biomedical research of coronary artery disease because of its high performance hardware and software. CFD techniques have been applied to study cardiovascular haemodynamics through simulation tools to predict the behaviour of circulatory blood flow in the human body. CFD simulation based on 3D luminal reconstructions can be used to analyse the local flow fields and flow profiling due to changes of coronary artery geometry, thus, identifying risk factors for development and progression of coronary artery disease. This review aims to provide an overview of the CFD applications in coronary artery disease, including biomechanics of atherosclerotic plaques, plaque progression and rupture; regional haemodynamics relative to plaque location and composition. A critical appraisal is given to a more recently developed application, fractional flow reserve based on CFD computation with regard to its diagnostic accuracy in the detection of haemodynamically significant coronary artery disease.

Wall shear stress distribution at the carotid bifurcation: influence of eversion carotid endarterectomy

OBJECTIVES

To test the feasibility of four-dimensional (4D) flow MRI to quantify the systolic wall shear stress (WSSsystole) and oscillatory shear index (OSI) in high-grade internal carotid artery (ICA) stenosis before and after endarterectomy (CEA).

METHODS

Twenty patients with ≥60 % ICA stenosis were prospectively and consequently included. Four-dimensional flow MRI was used to measure individual time-resolved 3D blood flow velocities. Segmental WSSsystole and OSI were derived at eight wall segments in analysis planes positioned along the ICA, common (CCA) and external carotid artery (ECA).

RESULTS

Regional WSSsystole of all patients decreased after CEA (P < 0.05). Changes were most prominent at the ICA bulb but remained unchanged in the CCA and ECA. OSI was significantly lower after CEA in the lateral vessel walls (P < 0.05). For analysis planes at the stenosis in- and outlet, a reduction of mean WSSsystole by 32 % and 52 % (P < 0.001) and OSI distal to the stenosis (40 %, P = 0.01) was found after CEA.

CONCLUSIONS

Our findings show the potential of in vivo 4D flow MRI to quantify haemodynamic changes in wall shear stress even in patients with complex flow conditions.

The effect of statins on high-risk atherosclerotic plaque associated with low endothelial shear stress

PURPOSE OF REVIEW

Low endothelial shear stress (ESS) plays an important role in the progression and severity of atherosclerotic lesions. As 3-hydroxy-3-methylglutaryl coenzyme A reductase inhibitors (statins) appear to stabilize plaque, it would be valuable to understand how statins affect the nature of lesions in the proatherogenic and proinflammatory environment of low ESS and the effect of statins on that atherosclerotic process. The purpose of this review is to summarize the relationship among low ESS, high-risk plaque and statins.

RECENT FINDINGS

Low ESS is a critically important determinant of plaque development and progression to high-risk plaques with large necrotic lipid core, intensive inflammation and thin fibrous cap. In addition to the proatherogenic phenotypic switching in areas of low ESS, local LDL cholesterol concentrations are also increased in areas of low ESS, which exacerbates the local atherogenic process. In experimental models, statins appear to reduce the inflammation in lesions associated with low ESS and reduce the atherosclerotic phenotype even in these high-risk prone vascular areas.

SUMMARY

The relationship between low ESS and statins has not been fully investigated, but the available data underscore the vasculoprotective effect of statins. Understanding the mechanisms whereby statins reduce the atherogenic and inflammatory phenotype resulting from a low ESS environment would provide new insights to design strategies to prevent regional formation of high-risk, inflamed plaques likely to rupture and cause an adverse clinical event.

Computational fluid dynamics in cardiovascular disease

Computational fluid dynamics (CFD) is a mechanical engineering field for analyzing fluid flow, heat transfer, and associated phenomena, using computer-based simulation. CFD is a widely adopted methodology for solving complex problems in many modern engineering fields. The merit of CFD is developing new and improved devices and system designs, and optimization is conducted on existing equipment through computational simulations, resulting in enhanced efficiency and lower operating costs. However, in the biomedical field, CFD is still emerging. The main reason why CFD in the biomedical field has lagged behind is the tremendous complexity of human body fluid behavior. Recently, CFD biomedical research is more accessible, because high performance hardware and software are easily available with advances in computer science. All CFD processes contain three main components to provide useful information, such as pre-processing, solving mathematical equations, and post-processing. Initial accurate geometric modeling and boundary conditions are essential to achieve adequate results. Medical imaging, such as ultrasound imaging, computed tomography, and magnetic resonance imaging can be used for modeling, and Doppler ultrasound, pressure wire, and non-invasive pressure measurements are used for flow velocity and pressure as a boundary condition. Many simulations and clinical results have been used to study congenital heart disease, heart failure, ventricle function, aortic disease, and carotid and intra-cranial cerebrovascular diseases. With decreasing hardware costs and rapid computing times, researchers and medical scientists may increasingly use this reliable CFD tool to deliver accurate results. A realistic, multidisciplinary approach is essential to accomplish these tasks. Indefinite collaborations between mechanical engineers and clinical and medical scientists are essential. CFD may be an important methodology to understand the pathophysiology of the development and progression of disease and for establishing and creating treatment modalities in the cardiovascular field.

Carotid plaque regression following 6-month statin therapy assessed by 3T cardiovascular magnetic resonance: comparison with ultrasound intima media thickness

BACKGROUND

Cardiovascular magnetic resonance (CMR) allows volumetric carotid plaque measurement that has advantage over 2-dimensional ultrasound (US) intima-media thickness (IMT) in evaluating treatment response. We tested the hypothesis that 6-month statin treatment in patients with carotid plaque will lead to plaque regression when measured by 3 Tesla CMR but not by IMT.

METHODS

Twenty-six subjects (67 ± 2 years, 7 females) with known carotid plaque (> 1.1 mm) and coronary or cerebrovascular atherosclerotic disease underwent 3T CMR (T1, T2, proton density and time of flight sequences) and US at baseline and following 6 months of statin therapy (6 had initiation, 7 had increase and 13 had maintenance of statin dosing). CMR plaque volume (PV) was measured in the region 12 mm below and up to 12 mm above carotid flow divider using software. Mean posterior IMT in the same region was measured. Baseline and 6-month CMR PV and US IMT were compared. Change in lipid rich/necrotic core (LR/NC) and calcification plaque components from CMR were related to change in PV.

RESULTS

Low-density lipoprotein cholesterol decreased (86 ± 6 to 74 ± 4 mg/dL, p = 0.046). CMR PV decreased 5.8 ± 2% (1036 ± 59 to 976 ± 65 mm3, p = 0.018). Mean IMT was unchanged (1.12 ± 0.06 vs. 1.14 ± 0.06 mm, p = NS). Patients with initiation or increase of statins had -8.8 ± 2.8% PV change (p = 0.001) while patients with maintenance of statin dosing had -2.7 ± 3% change in PV (p = NS). There was circumferential heterogeneity in CMR plaque thickness with greatest thickness in the posterior carotid artery, in the region opposite the flow divider. Similarly there was circumferential regional difference in change of plaque thickness with significant plaque regression in the anterior carotid region in region of the flow divider. Change in LR/NC (R = 0.62, p = 0.006) and calcification (R = 0.45, p = 0.03) correlated with PV change.

CONCLUSIONS

Six month statin therapy in patients with carotid plaque led to reduced plaque volume by 3T CMR, but ultrasound posterior IMT did not show any change. The heterogeneous spatial distribution of plaque and regional differences in magnitude of plaque regression may explain the difference in findings and support volumetric measurement of plaque. 3T CMR has potential advantage over ultrasound IMT to assess treatment response in individuals and may allow reduced sample size, duration and cost of clinical trials of plaque regression.