Development of a general method for designing microvascular networks using distribution of wall shear stress

In the present study, theoretical formulations for calculation of optimal bifurcation angle and relationship between the diameters of mother and daughter vessels using the power law model for non-Newtonian fluids are developed. The method is based on the distribution of wall shear stress in the mother and daughter vessels. Also, the effect of distribution of wall shear stress on the minimization of energy loss and flow resistance is considered. It is shown that constant wall shear stress in the mother and daughter vessels provides the minimum flow resistance and energy loss of biological flows. Moreover, the effects of different wall shear stresses in the mother and daughter branches, different lengths of daughter branches in the asymmetric bifurcations and non-Newtonian effect of biological fluid flows on the bifurcation angle and the relationship between the diameters of mother and daughter branches are considered. Using numerical simulations for non-Newtonian models such as power law and Carreau models, the effects of optimal bifurcation angle on the pressure drop and flow resistance of blood flow in the symmetric bifurcation are investigated. Numerical simulations show that optimal bifurcation angle decreases the pressure drop and flow resistance especially for bifurcations at large Reynolds number.     

Realistic non-Newtonian viscosity modelling highlights hemodynamic differences between intracranial aneurysms with and without surface blebs

Most computational fluid dynamic (CFD) simulations of aneurysm hemodynamics assume constant (Newtonian) viscosity, even though blood demonstrates shear-thinning (non-Newtonian) behavior. We sought to evaluate the effect of this simplifying assumption on hemodynamic forces within cerebral aneurysms, especially in regions of low wall shear stress, which are associated with rupture. CFD analysis was performed for both viscosity models using 3D rotational angiography volumes obtained for 26 sidewall aneurysms (12 with blebs, 12 ruptured), and parametric models incorporating blebs at different locations (inflow/outflow zone). Mean and lowest 5% values of time averaged wall shear stress (TAWSS) computed over the dome were compared using Wilcoxon rank-sum test. Newtonian modeling not only resulted in higher aneurysmal TAWSS, specifically in areas of low flow and blebs, but also showed no difference between aneurysms with or without blebs. In contrast, for non-Newtonian analysis, bleb-bearing aneurysms showed significantly lower 5% TAWSS compared to those without (p=0.005), despite no significant difference in mean dome TAWSS (p=0.32). Non-Newtonian modeling also accentuated the differences in dome TAWSS between ruptured and unruptured aneurysms (p<0.001). Parametric models further confirmed that realistic non-Newtonian viscosity resulted in lower bleb TAWSS and higher focal viscosity, especially when located in the outflow zone. The results show that adopting shear-thinning non-Newtonian blood viscosity in CFD simulations of intracranial aneurysms uncovered hemodynamic differences induced by bleb presence on aneurysmal surfaces, and significantly improved discriminant statistics used in risk stratification. These findings underline the possible implications of using a realistic model of blood viscosity in predictive computational hemodynamics.     

A study of wall shear stress in 12 aneurysms with respect to different viscosity models and flow conditions

Recent computational fluid dynamics (CFD) studies relate abnormal blood flow to rupture of cerebral aneurysms. However, it is still debated how to model blood flow with sufficient accuracy. Common assumptions made include Newtonian behaviour of blood, traction free outlet boundary conditions and inlet boundary conditions based on available literature. These assumptions are often required since the available patient specific data is usually restricted to the geometry of the aneurysm and the surrounding vasculature. However, the consequences of these assumptions have so far been inadequately addressed.     

血浆层对血管壁剪应力和剪应力梯度的影响

在细小血管中,由于血细胞明显的趋轴效应,管中的血液分为两个不同的区域,即具有血细胞的核心区和邻近管壁和血浆层。应用两相分层流模型,研究在相同的流量和管径下,当核心区中的血液分别为牛顿流体和Casson流体时,不同的血浆层厚度对细小血管壁剪应力和剪应力梯度的影响。结果表明,血浆层的存在对壁剪应力和壁剪应力梯度有较大影响,当血浆层厚度仅为血管半径的1％和3％时,壁剪应力梯度分别下降约10％和20％。     

Flow field and oscillatory shear stress in a tuning-fork-shaped model of the average human carotid bifurcation

The oscillatory shear index (OSI) was developed based on the hypothesis that intimal hyperplasia was correlated with oscillatory shear stresses. However, the validity of the OSI was in question since the correlation between intimal thickness and the OSI at the side walls of the sinus in the Y-shaped model of the average human carotid bifurcation (Y-AHCB) was weak. The objectives of this paper are to examine whether the reason for the weak correlation lies in the deviation in geometry of Y-AHCB from real human carotid bifurcation, and whether this correlation is clearly improved in the tuning-fork-shaped model of the average human carotid bifurcation (TF-AHCB). The geometry of the TF-AHCB model was based on observation and statistical analysis of specimens from 74 cadavers. The flow fields in both models were studied and compared by using flow visualization methods under steady flow conditions and by using laser Doppler anemometer (LDA) under pulsatile flow conditions. The TF-shaped geometry leads to a more complex flow field than the Y-shaped geometry. This added complexity includes strengthened helical movements in the sinus, new flow separation zone, and directional changes in the secondary flow patterns. The results show that the OSI-values at the side walls of the sinus in the TF-shaped model were more than two times as large as those in the Y-shaped model. This study confirmed the stronger correlation between the OSI and intimal thickness in the tuning-fork geometry of human carotid bifurcation, and the TF-AHCB model is a significant improvement over the traditional Y-shaped model.     

颈动脉血管壁切应力的分析

动脉中管壁的脉动低切应力在动脉粥样硬化形成中起始动和主要的决定作用。本文比较了几种计算血管壁切应力的方法,认为采用有约束的弹性管模型计算获得的动脉壁切应力更适合于临床应用。根据检测得到的正常人和动脉硬化性脑血管病患者的颈动脉血流速度、血管管径等数据,计算两者的颈动脉壁面切应力。研究发现动脉硬化性脑血管病患者的壁面切应力比正常人显著减小。这表明,颈动脉的壁面切应力可以作为动脉硬化性脑血管疾病的早期诊断的重要参考指标。     

Hemodynamics in coronary arteries with overlapping stents

Coronary artery stenosis is commonly treated by stent placement via percutaneous intervention, at times requiring multiple stents that may overlap. Stent overlap is associated with increased risk of adverse clinical outcome. While changes in local blood flow are suspected to play a role therein, hemodynamics in arteries with overlapping stents remain poorly understood. In this study we analyzed six cases of partially overlapping stents, placed ex vivo in porcine left coronary arteries and compared them to five cases with two non-overlapping stents. The stented vessel geometries were obtained by micro-computed tomography of corrosion casts. Flow and shear stress distribution were calculated using computational fluid dynamics. We observed a significant increase in the relative area exposed to low wall shear stress (WSS<0.5 Pa) in the overlapping stent segments compared both to areas without overlap in the same samples, as well as to non-overlapping stents. We further observed that the configuration of the overlapping stent struts relative to each other influenced the size of the low WSS area: positioning of the struts in the same axial location led to larger areas of low WSS compared to alternating struts. Our results indicate that the overlap geometry is by itself sufficient to cause unfavorable flow conditions that may worsen clinical outcome. While stent overlap cannot always be avoided, improved deployment strategies or stent designs could reduce the low WSS burden.     

Can aspect ratio be used to categorize intra-aneurysmal hemodynamics?-A study of elastase induced aneurysms in rabbit

Clinical studies suggest that aneurysm aspect ratio (AR) is an important indicator of rupture likelihood. The importance of AR is hypothesized to arise from its influence on intra-aneurysmal hemodynamics. It has been conjectured that slower flow in high AR sacs leads to a cascade of biological activities that weaken the aneurysm wall (Ujiie et al.,1999). However, the connection between AR, hemodynamics and wall weakening has never been proven. Animal models of saccular aneurysms provide a venue for evaluating this conjecture. The focus of this work was to evaluate whether a commonly used elastase induced aneurysm model in rabbits is suitable for a study of this kind from a hemodynamic perspective. In particular, to assess whether hemodynamic factors in low and high AR sacs are statistically different. To achieve this objective, saccular aneurysms were created in 51 rabbits and pulsatile computational fluid dynamics (CFD) studies were performed using rabbit specific inflows. Distinct hemodynamics were found in the low AR (AR<1.8, n=25), and high AR (AR>2.2, n=18) models. A single, stable recirculation zone was present in all low AR aneurysms, whereas a second, transient recirculation zone was also found in the superior aspect of the aneurysm dome for all high AR cases. Aneurysms with AR between 1.8 and 2.2 displayed transitional flow patterns. Differences in values and distributions of hemodynamic parameters were found between low and high AR cases including time averaged wall shear stress, oscillatory shear index, relative residence time and non-dimensional inflow rate. This work lays the foundation for future studies of the dependence of growth and remodeling on AR in the rabbit model and provides a motivation for further studies of the coupling between AR and hemodynamics in human aneurysms.     

Numerical investigation of the non-Newtonian blood flow in a bifurcation model with a non-planar branch

The non-Newtonian fluid flow in a bifurcation model with a non-planar daughter branch is investigated by using finite element method to solve the three-dimensional Navier–Stokes equations coupled with a non-Newtonian constitutive model, in which the shear thinning behavior of the blood fluid is incorporated by the Carreau–Yasuda model. The objective of this study is to investigate the influence of the non-Newtonian property of fluid as well as of curvature and out-of-plane geometry in the non-planar daughter vessel on wall shear stress (WSS) and flow phenomena. In the non-planar daughter vessel, the flows are typified by the skewing of the velocity profile towards the outer wall, creating a relatively low WSS at the inner wall. In the downstream of the bifurcation, the velocity profiles are shifted towards the flow divider. The low WSS is found at the inner walls of the curvature and the lateral walls of the bifurcation. Secondary flow patterns that swirl fluid from the inner wall of curvature to the outer wall in the middle of the vessel are also well documented for the curved and bifurcating vessels. The numerical results for the non-Newtonian fluid and the Newtonian fluid with original Reynolds number and the corresponding rescaled Reynolds number are presented. Significant difference between the non-Newtonian flow and the Newtonian flow is revealed; however, reasonable agreement between the non-Newtonian flow and the rescaled Newtonian flow is found. Results of this study support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.     

What is needed to make low-density lipoprotein transport in human aorta computational models suitable to explore links to atherosclerosis? Impact of initial and inflow boundary conditions

Personalized computational hemodynamics (CH) is a promising tool to clarify/predict the link between low density lipoproteins (LDL) transport in aorta, disturbed shear and atherogenesis. However, CH uses simplifying assumptions that represent sources of uncertainty. In particular, modelling blood-side to wall LDL transfer is challenged by the cumbersomeness of protocols needed to obtain reliable LDL concentration profile estimations. This paucity of data is limiting the establishment of rigorous CH protocols able to balance the trade-offs among the variety of in vivo data to be acquired, and the accuracy required by biological/clinical applications.In this study, we analyze the impact of LDL concentration initialization (initial conditions, ICs) and inflow boundary conditions (BCs) on CH models of LDL blood-to-wall transfer in aorta. Technically, in an image-based model of human aorta, two different inflow BCs are generated imposing subject-specific inflow 3D PC-MRI measured or idealized (flat) velocity profiles. For each simulated BC, four different ICs for LDL concentration are applied, imposing as IC the LDL distribution resulting from steady-state simulations with average conditions, or constant LDL concentration values.Based on CH results, we conclude that: (1) the imposition of realistic 3D velocity profiles as inflow BC reduces the uncertainty affecting the representation of LDL transfer; (2) different LDL concentration ICs lead to markedly different patterns of LDL transfer.Given that it is not possible to verify in vivo the proper LDL concentration initialization to be applied, we suggest to carefully set and unambiguously declare the imposed BCs and LDL concentration IC when modelling LDL transfer in aorta, in order to obtain reproducible and ultimately comparable results among different laboratories.     

Peak systolic or maximum intra-aneurysmal hemodynamic condition? Implications on normalized flow variables

Background: CFD has been used to assess intra-aneurysmal hemodynamics. Nevertheless, the lack of patient-specific flow information has triggered the possibility of implementing a wide variety of physiological flow conditions. Due to these uncertainties in the patient flow conditions, the normalization of the intra-aneurysmal hemodynamics is generally conducted.