Mechanical Engineering Department, São Paulo State University (UNESP), Ilha Solteira, São Paulo 15385-000, Brazil.
Department of Mechanical Engineering, Dalhousie University, Halifax, NS B3H 4R2, Canada.
J Biomech Eng. 2021 Jul 1;143(7). doi: 10.1115/1.4050539.
When simulating blood flow in intracranial aneurysms (IAs), the Newtonian model seems to be ubiquitous. However, analyzing the results from the few studies on this subject, the doubt remains on whether it is necessary to use non-Newtonian models in computational fluid dynamics (CFD) simulations of cerebral vascular flows. The objective of this study is to investigate whether different rheology models would influence the hemodynamic parameters related to the wall shear stress (WSS) for ruptured and unruptured IA cases, especially because ruptured aneurysms normally have morphological features, such as lobular regions and blebs, that could trigger non-Newtonian phenomena in the blood flow due to low shear rates. Using CFD in an open-source framework, we simulated four ruptured and four unruptured patient-specific aneurysms to assess the influence of the blood modeling on the main hemodynamic variables associated with aneurysm formation, growth, and rupture. Results for WSS and oscillatory shear index (OSI) and their metrics were obtained using Casson and Carreau-Yasuda non-Newtonian models and were compared with those obtained using the Newtonian model. We found that all differences between non-Newtonian and the Newtonian models were consistent among all cases irrespective of their rupture status. We further found that the WSS at peak systole is overestimated by more than 50% by using the non-Newtonian models, but its metrics based on time and surface averaged values are less affected-the maximum relative difference among the cases is 7% for the Casson model. On the other hand, the surface-averaged OSI is underestimated by more than 30% by the non-Newtonian models. These results suggest that it is recommended to investigate different blood rheology models in IAs simulations when specific parameters to characterize the flow are needed, such as peak-systole WSS and OSI.
在模拟颅内动脉瘤(IA)中的血流时,牛顿模型似乎无处不在。然而,分析少数关于该主题的研究结果,仍然存在疑问,即在脑血流的计算流体动力学(CFD)模拟中是否有必要使用非牛顿模型。本研究的目的是研究不同流变学模型是否会影响与壁面切应力(WSS)相关的血流动力学参数,特别是因为破裂的动脉瘤通常具有形态特征,如小叶区域和气泡,这可能由于低剪切速率而在血流中引发非牛顿现象。我们使用开源框架中的 CFD 模拟了四个破裂和四个未破裂的患者特定动脉瘤,以评估血液建模对与动脉瘤形成、生长和破裂相关的主要血流动力学变量的影响。使用 Casson 和 Carreau-Yasuda 非牛顿模型获得了 WSS 和振荡剪切指数(OSI)及其指标的结果,并将其与牛顿模型的结果进行了比较。我们发现,无论动脉瘤是否破裂,非牛顿模型和牛顿模型之间的所有差异在所有情况下都是一致的。我们进一步发现,使用非牛顿模型会使峰值收缩期的 WSS 高估超过 50%,但其基于时间和表面平均值的指标受影响较小-在 Casson 模型中,各病例之间的最大相对差异为 7%。另一方面,非牛顿模型会使表面平均 OSI 低估超过 30%。这些结果表明,当需要特定参数来描述流动时,建议在 IA 模拟中研究不同的血液流变学模型,例如峰值收缩期 WSS 和 OSI。