Lotfi Azadeh, Simmons Anne, Barber Tracie
J Biomech Eng. 2016 Mar;138(3):4032502. doi: 10.1115/1.4032502.
The formation and progression of in-stent restenosis (ISR) in bifurcated vessels may vary depending on the technique used for stenting. This study evaluates the effect of a variety of mesh styles on the accuracy and reliability of computational fluid dynamics (CFD) models in predicting these regions, using an idealized stented nonbifurcated model. The wall shear stress (WSS) and the near-stent recirculating vortices are used as determinants. The meshes comprise unstructured tetrahedral and polyhedral elements. The effects of local refinement, as well as higher-order elements such as prismatic inflation layers and internal hexahedral core, have also been examined. The uncertainty associated with individual mesh style was assessed through verification of calculations using the grid convergence index (GCI) method. The results obtained show that the only condition which allows the reliable comparison of uncertainty estimation between different meshing styles is that the monotonic convergence of grid solutions is in the asymptotic range. Comparisons show the superiority of a flow-adaptive polyhedral mesh over the commonly used adaptive and nonadaptive tetrahedral meshes in terms of resolving the near-stent flow features, GCI value, and prediction of WSS. More accurate estimation of hemodynamic factors was obtained using higher-order elements, such as hexahedral or prismatic grids. Incorporating these higher-order elements, however, was shown to introduce some degrees of numerical diffusion at the transitional area between the two meshes, not necessarily translating into high GCI value. Our data also confirmed the key role of local refinement in improving the performance and accuracy of nonadaptive mesh in predicting flow parameters in models of stented artery. The results of this study can provide a guideline for modeling biofluid domain in complex bifurcated arteries stented in regards to various stenting techniques.
分叉血管内支架再狭窄(ISR)的形成和进展可能因支架置入技术的不同而有所差异。本研究使用理想化的非分叉支架模型,评估了多种网格样式对计算流体动力学(CFD)模型预测这些区域的准确性和可靠性的影响。壁面切应力(WSS)和支架附近的再循环涡流被用作决定因素。网格由非结构化四面体和多面体单元组成。还研究了局部细化以及高阶单元(如棱柱形膨胀层和内部六面体核心)的影响。通过使用网格收敛指数(GCI)方法验证计算来评估与单个网格样式相关的不确定性。所得结果表明,能够可靠比较不同网格样式之间不确定性估计的唯一条件是网格解的单调收敛处于渐近范围内。比较表明,在解析支架附近的流动特征、GCI值和WSS预测方面,流动自适应多面体网格优于常用的自适应和非自适应四面体网格。使用高阶单元(如六面体或棱柱形网格)可获得更准确的血流动力学因素估计。然而,纳入这些高阶单元在两个网格之间的过渡区域会引入一定程度的数值扩散,不一定会转化为高GCI值。我们的数据还证实了局部细化在提高非自适应网格预测支架动脉模型中流动参数的性能和准确性方面的关键作用。本研究结果可为针对各种支架置入技术的复杂分叉动脉内生物流体域建模提供指导。