De Wilde David, Trachet Bram, De Meyer Guido, Segers Patrick
IBiTech-bioMMeda, Ghent University-IMinds Medical IT, Ghent, Belgium.
IBiTech-bioMMeda, Ghent University-IMinds Medical IT, Ghent, Belgium; Institute of Bioengineering, EPFL, Lausanne, Switzerland.
J Biomech. 2016 Sep 6;49(13):2741-2747. doi: 10.1016/j.jbiomech.2016.06.010. Epub 2016 Jun 15.
Low and oscillatory wall shear stresses (WSS) near aortic bifurcations have been linked to the onset of atherosclerosis. In previous work, we calculated detailed WSS patterns in the carotid bifurcation of mice using a Fluid-structure interaction (FSI) approach. We subsequently fed the animals a high-fat diet and linked the results of the FSI simulations to those of atherosclerotic plaque location on a within-subject basis. However, these simulations were based on boundary conditions measured under anesthesia, while active mice might experience different hemodynamics. Moreover, the FSI technique for mouse-specific simulations is both time- and labor-intensive, and might be replaced by simpler and easier Computational Fluid Dynamics (CFD) simulations. The goal of the current work was (i) to compare WSS patterns based on anesthesia conditions to those representing active resting and exercising conditions; and (ii) to compare WSS patterns based on FSI simulations to those based on steady-state and transient CFD simulations.
For each of the 3 computational techniques (steady state CFD, transient CFD, FSI) we performed 5 simulations: 1 for anesthesia, 2 for conscious resting conditions and 2 more for conscious active conditions. The inflow, pressure and heart rate were scaled according to representative in vivo measurements obtained from literature.
When normalized by the maximal shear stress value, shear stress patterns were similar for the 3 computational techniques. For all activity levels, steady state CFD led to an overestimation of WSS values, while FSI simulations yielded a clear increase in WSS reversal at the outer side of the sinus of the external carotid artery that was not visible in transient CFD-simulations. Furthermore, the FSI simulations in the highest locomotor activity state showed a flow recirculation zone in the external carotid artery that was not present under anesthesia. This recirculation went hand in hand with locally increased WSS reversal.
Our data show that FSI simulations are not necessary to obtain normalized WSS patterns, but indispensable to assess the oscillatory behavior of the WSS in mice. Flow recirculation and WSS reversal at the external carotid artery may occur during high locomotor activity while they are not present under anesthesia. These phenomena might thus influence plaque formation to a larger extent than what was previously assumed.
主动脉分叉处附近的低壁面剪切应力(WSS)和振荡性壁面剪切应力与动脉粥样硬化的发生有关。在之前的研究中,我们使用流固耦合(FSI)方法计算了小鼠颈动脉分叉处的详细WSS模式。随后,我们给这些动物喂食高脂饮食,并在个体水平上将FSI模拟结果与动脉粥样硬化斑块位置的结果联系起来。然而,这些模拟是基于在麻醉状态下测量的边界条件,而活跃的小鼠可能会经历不同的血流动力学。此外,用于小鼠特异性模拟的FSI技术既耗时又费力,可能会被更简单易行的计算流体动力学(CFD)模拟所取代。当前研究的目的是:(i)比较基于麻醉条件的WSS模式与代表活跃静息和运动条件的WSS模式;(ii)比较基于FSI模拟的WSS模式与基于稳态和瞬态CFD模拟的WSS模式。
对于3种计算技术(稳态CFD、瞬态CFD、FSI)中的每一种,我们都进行了5次模拟:1次用于麻醉状态,2次用于清醒静息状态,另外2次用于清醒活跃状态。根据从文献中获得的代表性体内测量数据对流入量、压力和心率进行缩放。
当用最大剪切应力值进行归一化时,3种计算技术的剪切应力模式相似。对于所有活动水平,稳态CFD导致WSS值高估,而FSI模拟在外颈动脉窦外侧产生了明显的WSS逆转增加,这在瞬态CFD模拟中不可见。此外,最高运动活动状态下的FSI模拟显示外颈动脉存在一个流动再循环区域,而在麻醉状态下不存在。这种再循环与局部增加的WSS逆转同时出现。
我们的数据表明,获得归一化的WSS模式不一定需要FSI模拟,但评估小鼠WSS的振荡行为则必不可少。外颈动脉的流动再循环和WSS逆转可能在高运动活动期间发生,而在麻醉状态下不存在。因此,这些现象可能比之前假设的对斑块形成的影响更大。