Non-Newtonian Endothelial Shear Stress Simulation: Does It Matter?

Patient-specific coronary endothelial shear stress (ESS) calculations using Newtonian and non-Newtonian rheological models were performed to assess whether the common assumption of Newtonian blood behavior offers similar results to a more realistic but computationally expensive non-Newtonian model. 16 coronary arteries (from 16 patients) were reconstructed from optical coherence tomographic (OCT) imaging. Pulsatile CFD simulations using Newtonian and the Quemada non-Newtonian model were performed. Endothelial shear stress (ESS) and other indices were compared. Exploratory indices including local blood viscosity (LBV) were calculated from non-Newtonian simulation data. Compared to the Newtonian results, the non-Newtonian model estimates significantly higher time-averaged ESS (1.69 (IQR 1.36)Pa versus 1.28 (1.16)Pa, p < 0.001) and ESS gradient (0.90 (1.20)Pa/mm versus 0.74 (1.03)Pa/mm, p < 0.001) throughout the cardiac cycle, under-estimating the low ESS (<1Pa) area (37.20 ± 13.57% versus 50.43 ± 14.16%, 95% CI 11.28-15.18, p < 0.001). Similar results were also found in the idealized artery simulations with non-Newtonian median ESS being higher than the Newtonian median ESS (healthy segments: 0.8238Pa versus 0.6618Pa, p < 0.001 proximal; 0.8179Pa versus 0.6610Pa, p < 0.001 distal; stenotic segments: 0.8196Pa versus 0.6611Pa, p < 0.001 proximal; 0.2546Pa versus 0.2245Pa, p < 0.001 distal) On average, the non-Newtonian model has a LBV of 1.45 times above the Newtonian model with an average peak LBV of 40-fold. Non-Newtonian blood model estimates higher quantitative ESS values than the Newtonian model. Incorporation of non-Newtonian blood behavior may improve the accuracy of ESS measurements. The non-Newtonian model also allows calculation of exploratory viscosity-based hemodynamic indices, such as local blood viscosity, which may offer additional information to detect underlying atherosclerosis.

Elevated Blood Viscosity and Microrecirculation Resulting From Coronary Stent Malapposition.

One particular complexity of coronary artery is the natural tapering of the vessel with proximal segments having larger caliber and distal tapering as the vessel get smaller. The natural tapering of a coronary artery often leads to proximal incomplete stent apposition (ISA). ISA alters coronary hemodynamics and creates pathological path to develop complications such as in-stent restenosis, and more worryingly, stent thrombosis (ST). By employing state-of-the-art computer-aided design software, generic stent hoops were virtually deployed in an idealized tapered coronary artery with decreasing malapposition distance. Pulsatile blood flow simulations were carried out using computational fluid dynamics (CFD) on these computer-aided design models. CFD results reveal unprecedented details in both spatial and temporal development of microrecirculation environments throughout the cardiac cycle (CC). Arterial tapering also introduces secondary microrecirculation. These primary and secondary microrecirculations provoke significant fluctuations in arterial wall shear stress (WSS). There has been a direct correlation with changes in WSS and the development of atherosclerosis. Further, the presence of these microrecirculations influence strongly on the local levels of blood viscosity in the vicinity of the malapposed stent struts. The observation of secondary microrecirculations and changes in blood rheology is believed to complement the wall (-based) shear stress, perhaps providing additional physical explanations for tissue accumulation near ISA detected from high resolution optical coherence tomography (OCT).

Enhanced targeted drug delivery through controlled release in a three-dimensional vascular tree.

"Controlled particle release and targeting" is a technique using particle release score map (PRSM) and transient particle release score map (TPRSM) via backtracking to determine optimal drug injection locations for achieving an enhanced target efficiency (TE). This paper investigates the possibility of targeting desired locations through an idealized but complex three-dimensional (3D) vascular tree geometry under realistic hemodynamic conditions by imposing a Poiseuille velocity profile and a Womersley velocity profile derived from cine phase contrast magnetic resonance imaging (MRI) data for steady and pulsatile simulations, respectively. The shear thinning non-Newtonian behavior of blood was accounted for by the Carreau-Yasuda model. One-way coupled Eulerian-Lagrangian particle tracking method was used to record individual drug particle trajectories. Particle size and density showed negligible influence on the particle fates. With the proposed optimal release scoring algorithm, multiple optimal release locations were determined under steady flow conditions, whereas there was one unique optimal release location under pulsatile flow conditions. The initial in silico results appear promising, showing on average 66% TE in the pulsatile simulations, warranting further studies to improve the mathematical model and experimental validation.

Numerical investigations of the haemodynamic changes associated with stent malapposition in an idealised coronary artery.

The deployment of a coronary stent near complex lesions can sometimes lead to incomplete stent apposition (ISA), an undesirable side effect of coronary stent implantation. Three-dimensional computational fluid dynamics (CFD) calculations are performed on simplified stent models (with either square or circular cross-section struts) inside an idealised coronary artery to analyse the effect of different levels of ISA to the change in haemodynamics inside the artery. The clinical significance of ISA is reported using haemodynamic metrics like wall shear stress (WSS) and wall shear stress gradient (WSSG). A coronary stent with square cross-sectional strut shows different levels of reverse flow for malapposition distance (MD) between 0mm and 0.12 mm. Chaotic blood flow is usually observed at late diastole and early systole for MD=0mm and 0.12 mm but are suppressed for MD=0.06 mm. The struts with circular cross section delay the flow chaotic process as compared to square cross-sectional struts at the same MD and also reduce the level of fluctuations found in the flow field. However, further increase in MD can lead to chaotic flow not only at late diastole and early systole, but it also leads to chaotic flow at the end of systole. In all cases, WSS increases above the threshold value (0.5 Pa) as MD increases due to the diminishing reverse flow near the artery wall. Increasing MD also results in an elevated WSSG as flow becomes more chaotic, except for square struts at MD=0.06 mm.