Huo Yunlong, Guo Xiaomei, Kassab Ghassan S
Department of Biomedical Engineering, IUPUI, Indianapolis, IN 46202, USA.
Ann Biomed Eng. 2008 May;36(5):685-99. doi: 10.1007/s10439-008-9473-4. Epub 2008 Feb 26.
There is a spatial disposition to atherosclerosis along the aorta corresponding to regions of flow disturbances. The objective of the present study is to investigate the detailed distribution of hemodynamic parameters (wall shear stress (WSS), spatial gradient of wall shear stress (WSSG), and oscillatory shear index (OSI)) in the entire length of C57BL/6 mouse aorta with all primary branches (from ascending aorta to common iliac bifurcation). The detailed geometrical parameters (e.g., diameter and length of the vessels) were obtained from casts of entire aorta and primary branches of mice. The flow velocity was measured at the inlet of ascending aorta using Doppler flowprobe in mice. The outlet pressure boundary condition was estimated based on scaling law. The continuity and Navier-Stokes equations were solved using three-dimensional finite element method (FEM). The model prediction was tested by comparing the computed flow rate with the flow rate measured just before the common iliac bifurcation, and good agreement was found. It was also found that complex flow patterns occur at bifurcations between main trunk and branches. The major branches of terminal aorta, with the highest proportion of atherosclerosis, have the lowest WSS, and the relatively atherosclerotic-prone aortic arch has much more complex WSS distribution and higher OSI value than other sites. The low WSS coincides with the high OSI, which approximately obeys a power law relationship. Furthermore, the scaling law between flow and diameter holds in the entire aorta and primary branches of mice under pulsatile blood flow conditions. This model will eventually serve to elucidate the causal relation between hemodynamic patterns and atherogenesis in KO mice.
动脉粥样硬化在主动脉上存在与血流紊乱区域相对应的空间分布。本研究的目的是调查C57BL/6小鼠主动脉全长及所有主要分支(从升主动脉到髂总动脉分叉处)的血流动力学参数(壁面切应力(WSS)、壁面切应力空间梯度(WSSG)和振荡切应力指数(OSI))的详细分布情况。详细的几何参数(如血管直径和长度)是从小鼠主动脉及主要分支的铸型中获得的。使用多普勒血流探头在小鼠升主动脉入口处测量血流速度。基于比例定律估算出口压力边界条件。采用三维有限元方法(FEM)求解连续性方程和纳维-斯托克斯方程。通过将计算得到的流量与在髂总动脉分叉前测量的流量进行比较来检验模型预测,结果发现二者吻合良好。还发现主干与分支之间的分叉处会出现复杂的血流模式。终末主动脉的主要分支动脉粥样硬化比例最高,其WSS最低,相对易发生动脉粥样硬化的主动脉弓的WSS分布比其他部位更为复杂,OSI值也更高。低WSS与高OSI同时出现,二者大致遵循幂律关系。此外,在脉动血流条件下,小鼠主动脉全长及主要分支中流量与直径之间的比例定律成立。该模型最终将有助于阐明基因敲除小鼠血流动力学模式与动脉粥样硬化发生之间的因果关系。