Redaelli A, Boschetti F, Inzoli F
Dipartimento di Bioingegneria, Politecnico di Milano, Italy.
Comput Biol Med. 1997 May;27(3):233-47. doi: 10.1016/s0010-4825(97)00006-1.
In this paper we present a new method for the assignment of pulsatile velocity profiles as input boundary conditions in finite element models of arteries. The method is based on the implementation of the analytical solution for developed pulsatile flow in a rigid straight tube. The analytical solution provides the fluid dynamics of the region upstream from the fluid domain to be investigated by means of the finite element approach. In standard fluid dynamics finite element applications, the inlet developed velocity profiles are achieved assuming velocity boundary conditions to be easily implementable-such as flat or parabolic velocity profiles-applied to a straight tube of appropriate length. The tube is attached to the inflow section of the original fluid domain so that the flow can develop fully. The comparison between the analytical solution and the traditional numerical approach indicates that the analytical solution has some advantages over the numerical one. Moreover, the results suggest that subroutine employment allows a consistent reduction in solving time especially for complex fluid dynamic model, and significantly decreases the storage and memory requirements for computations.
在本文中,我们提出了一种新方法,用于在动脉有限元模型中指定脉动速度剖面作为输入边界条件。该方法基于刚性直管中充分发展的脉动流解析解的实现。该解析解通过有限元方法提供了待研究流体域上游区域的流体动力学。在标准流体动力学有限元应用中,假设速度边界条件易于实现(例如平坦或抛物线速度剖面),并应用于适当长度的直管,从而获得入口处充分发展的速度剖面。该直管连接到原始流体域的流入部分,以便流动能够充分发展。解析解与传统数值方法的比较表明,解析解比数值方法具有一些优势。此外,结果表明,使用子程序可以显著减少求解时间,特别是对于复杂的流体动力学模型,并且大大降低了计算所需的存储和内存要求。