Department of Biomedical Engineering, Florida International University, Miami, FL 33174.
J Biomech Eng. 2022 Nov 1;144(11). doi: 10.1115/1.4054630.
The Newtonian model has commonly been used to represent the viscosity of blood in the aorta, despite blood itself being a non-Newtonian fluid. This is justified where shear rates tend to be large. However, we hypothesized that using the Newtonian model to predict the hemodynamics on the aortic valve, particularly in those with severe calcifications, is inaccurate owing to valve leaflet geometry irregularities inducing multiple regions of low shear rates, <100 s-1, where a Newtonian model is invalid. We investigated the utility of three fluid viscosity models via quasi-static simulations: Newtonian, Carreau, and Quemada on a severely calcified aortic heart valve and compared their ability to capture important hemodynamic parameters of wall shear stress (WSS) and the oscillatory shear index (OSI). Our findings indicate that when the shear rates were large enough, >100 s-1, the use of a Newtonian model was justified. However, in spatial regions of relatively low shear rates, <100 s-1, specifically on the inner cusps of the fibrosa side of the valve, WSS calculations under a Newtonian model were found to be noticeably different when compared with their non-Newtonian, Carreau and Quemada counterparts. We hereby conclude that to facilitate more accurate computational flow simulations in severe aortic heart valve calcification, which is subjected to relatively large spatial regions of low shear (<100 s-1), a non-Newtonian model should be applied.
尽管血液本身是一种非牛顿流体,但牛顿模型通常被用于表示主动脉中的血液粘度。这在剪切率趋于较大的情况下是合理的。然而,我们假设,由于瓣叶几何形状不规则导致多个低剪切率区域(<100s-1),牛顿模型无效,因此使用牛顿模型来预测主动脉瓣上的血液动力学,特别是在严重钙化的情况下,是不准确的。我们通过准静态模拟研究了三种流体粘度模型的实用性:牛顿模型、Carreau 模型和 Quemada 模型,应用于严重钙化的主动脉心脏瓣膜,并比较了它们捕获壁面剪切应力(WSS)和振荡剪切指数(OSI)等重要血液动力学参数的能力。我们的研究结果表明,当剪切率足够大时,即>100s-1,使用牛顿模型是合理的。然而,在相对较低的剪切率空间区域,即<100s-1,特别是在瓣膜纤维状侧的内叶尖,与非牛顿模型的 Carreau 和 Quemada 相比,牛顿模型下的 WSS 计算明显不同。因此,我们得出结论,为了促进严重主动脉心脏瓣膜钙化中更准确的计算流体模拟,应应用非牛顿模型,因为这种钙化情况受到相对较大的低剪切区域(<100s-1)的影响。