Assessment of surface roughness and blood rheology on local coronary haemodynamics: a multi-scale computational fluid dynamics study

On the nonlinear relationship between wall shear stress topology and multi-directionality in coronary atherosclerosis

BACKGROUND AND OBJECTIVE

In this paper we investigate twelve multi-directional/topological wall shear stress (WSS) derived metrics and their relationships with the formation of coronary plaques in both computational fluid dynamics (CFD) and dynamic fluid-structure interaction (FSI) frameworks. While low WSS is one of the most established biomechanical markers associated with coronary atherosclerosis progression, alone it is limited. Multi-directional and topological WSS derived metrics have been shown to be important in atherosclerosis related mechanotransduction and near-wall transport processes. However, the relationships between these twelve WSS metrics and the influence of both FSI simulations and coronary dynamics is understudied.

METHODS

We first investigate the relationships between these twelve WSS derived metrics, stenosis percentage and lesion length through a parametric, transient CFD study. Secondly, we extend the parametric study to FSI, both with and without the addition of coronary dynamics, and assess their correlations. Finally, we present the case of a patient who underwent invasive coronary angiography and optical coherence tomography imaging at two time points 18 months apart. Associations between each of the twelve WSS derived metrics in CFD, static FSI and dynamic FSI simulations were assessed against areas of positive/negative vessel remodelling, and changes in plaque morphology.

RESULTS

22-32% stenosis was the threshold beyond which adverse multi-directional/topological WSS results. Each metric produced a different relationship with changing stenoses and lesion length. Transient haemodynamics was impacted by coronary dynamics, with the topological shear variation index suppressed by up to 94%. These changes appear more critical at smaller stenosis levels, suggesting coronary dynamics could play a role in the earlier stages of atherosclerosis development. In the patient case, both dynamics and FSI vs CFD changes altered associations with measured changes in plaque morphology. An appendix of the linear fits between the various FSI- and CFD-based simulations is provided to assist in scaling CFD-based results to resemble the compliant walled characteristics of FSI more accurately.

CONCLUSIONS

These results highlight the potential for coronary dynamics to alter multi-directional/topological WSS metrics which could impact associations with changes in coronary atherosclerosis over time. These results warrant further investigation in a wider range of morphological settings and longitudinal cohort studies in the future.

Fluid-structure interaction simulation of visceral perfusion and impact of different cannulation methods on aortic dissection

Hemodynamics in aortic dissection (AD) is closely associated with the risk of aortic aneurysm, rupture, and malperfusion. Altered blood flow in patients with AD can lead to severe complications such as visceral malperfusion. In this study, we aimed to investigate the effect of cannulation flow on hemodynamics in AD using a fluid-structure interaction simulation. We developed a specific-idealized AD model that included an intimal tear in the descending thoracic aorta, a re-entry tear in the left iliac artery, and nine branches. Two different cannulation methods were tested: (1) axillary cannulation (AC) only through the brachiocephalic trunk and (2) combined axillary and femoral cannulation (AFC) through the brachiocephalic trunk and the right common iliac artery. AC was found to result in the development of a pressure difference between the true lumen and false lumen, owing to the difference in the flow rate through each lumen. This pressure difference collapsed the true lumen, disturbing blood flow to the celiac and superior mesenteric arteries. However, in AFC, the pressure levels between the two lumens were similar, and no collapse occurred. Moreover, the visceral flow was higher than that in AC. Lastly, the stiffness of the intimal flap affected the true lumen's collapse.

Combined electromechanically driven pulsating flow of nonlinear viscoelastic fluids in narrow confinements

Controlled microscale transport is at the core of many scientific and technological advancements, including medical diagnostics, separation of biomolecules, etc., and often involves complex fluids. One of the challenges in this regard is to actuate flows at small scales in an energy efficient manner, given the strong viscous forces opposing fluid motion. We try to address this issue here by probing a combined time-periodic pressure and electrokinetically driven flow of a viscoelastic fluid obeying the simplified linear Phan-Thien-Tanner model, using numerical as well as asymptotic tools, in view of the fact that oscillatory fields are less energy intensive. We establish that the interplay between oscillatory electrical and mechanical forces can lead to complex temporal mass flow rate variations with short-term bursts and peaks in the flow rate. We further demonstrate that an oscillatory pressure gradient or an electric field, in tandem with another steady actuating force can indeed change the net throughput significantly-a paradigm that is not realized in Newtonian or other simpler polymeric liquids. Our results reveal that the extent of augmentation in the flow rate strongly depends on the frequency of the imposed actuating forces along with their waveforms. We also evaluate the streaming potential resulting from an oscillatory pressure-driven flow and illustrate that akin to the volume throughput, the streaming potential also shows complex temporal variations, while its time average gets augmented in the presence of a time-periodic pressure gradient in a nonlinear viscoelastic medium.

Accuracy of vascular tortuosity measures using computational modelling

Severe coronary tortuosity has previously been linked to low shear stresses at the luminal surface, yet this relationship is not fully understood. Several previous studies considered different tortuosity metrics when exploring its impact of on the wall shear stress (WSS), which has likely contributed to the ambiguous findings in the literature. Here, we aim to analyze different tortuosity metrics to determine a benchmark for the highest correlating metric with low time-averaged WSS (TAWSS). Using Computed Tomography Coronary Angiogram (CTCA) data from 127 patients without coronary artery disease, we applied all previously used tortuosity metrics to the left main coronary artery bifurcation, and to its left anterior descending and left circumflex branches, before modelling their TAWSS using computational fluid dynamics (CFD). The tortuosity measures included tortuosity index, average absolute-curvature, root-mean-squared (RMS) curvature, and average squared-derivative-curvature. Each tortuosity measure was then correlated with the percentage of vessel area that showed a < 0.4 Pa TAWSS, a threshold associated with altered endothelial cell cytoarchitecture and potentially higher disease risk. Our results showed a stronger correlation between curvature-based versus non-curvature-based tortuosity measures and low TAWSS, with the average-absolute-curvature showing the highest coefficient of determination across all left main branches (p < 0.001), followed by the average-squared-derivative-curvature (p = 0.001), and RMS-curvature (p = 0.002). The tortuosity index, the most widely used measure in literature, showed no significant correlation to low TAWSS (p = 0.86). We thus recommend the use of average-absolute-curvature as a tortuosity measure for future studies.

Numerical modelling of blood rheology and platelet activation through a stenosed left coronary artery bifurcation

Coronary bifurcations are prone to atherosclerotic plaque growth, experiencing regions of reduced wall shear stress (WSS) and increased platelet adhesion. This study compares effects across different rheological approaches on hemodynamics, combined with a shear stress exposure history model of platelets within a stenosed porcine bifurcation. Simulations used both single/multiphase blood models to determine which approach best predicts phenomena associated with atherosclerosis and atherothrombosis. A novel Lagrangian platelet tracking model was used to evaluate residence time and shear history of platelets indicating likely regions of thrombus formation. Results show a decrease in area of regions with pathologically low time-averaged WSS with the use of multiphase models, particularly in a stenotic bifurcation. Significant non-Newtonian effects were observed due to low-shear and varying hematocrit levels found on the outer walls of the bifurcation and distal to the stenosis. Platelet residence time increased 11% in the stenosed artery, with exposure times to low-shear sufficient for red blood cell aggregation (>1.5 s). increasing the risk of thrombosis. This shows stenotic artery hemodynamics are inherently non-Newtonian and multiphase, with variations in hematocrit (0.163-0.617) and elevated vorticity distal to stenosis (+15%) impairing the function of the endothelium via reduced time-averaged WSS regions, rheological properties and platelet activation/adhesion.

Computational fluid dynamics of the right atrium: Assessment of modelling criteria for the evaluation of dialysis catheters

Central venous catheters are widely used in haemodialysis therapy, having to respect design requirements for appropriate performance. These are placed within the right atrium (RA); however, there is no prior computational study assessing different catheter designs while mimicking their native environment. Here, a computational fluid dynamics model of the RA, based on realistic geometry and transient physiological boundary conditions, was developed and validated. Symmetric, split and step catheter designs were virtually placed in the RA and their performance was evaluated by: assessing their interaction with the RA haemodynamic environment through prediction of flow vorticity and wall shear stress (WSS) magnitudes (1); and quantifying recirculation and tip shear stress (2). Haemodynamic predictions from our RA model showed good agreement with the literature. Catheter placement in the RA increased average vorticity, which could indicate alterations of normal blood flow, and altered WSS magnitudes and distribution, which could indicate changes in tissue mechanical properties. All designs had recirculation and elevated shear stress values, which can induce platelet activation and subsequently thrombosis. The symmetric design, however, had the lowest associated values (best performance), while step design catheters working in reverse mode were associated with worsened performance. Different tip placements also impacted on catheter performance. Our findings suggest that using a realistically anatomical RA model to study catheter performance and interaction with the haemodynamic environment is crucial, and that care needs to be given to correct tip placement within the RA for improved recirculation percentages and diminished shear stress values.