Schenkel Torsten, Malve Mauro, Reik Michael, Markl Michael, Jung Bernd, Oertel Herbert
Institute of Fluid Mechanics, University of Karlsruhe, Karlsruhe, Germany.
Ann Biomed Eng. 2009 Mar;37(3):503-15. doi: 10.1007/s10439-008-9627-4. Epub 2009 Jan 6.
A three-dimensional computational fluid dynamics (CFD) method has been developed to simulate the flow in a pumping left ventricle. The proposed method uses magnetic resonance imaging (MRI) technology to provide a patient specific, time dependent geometry of the ventricle to be simulated. Standard clinical imaging procedures were used in this study. A two-dimensional time-dependent orifice representation of the heart valves was used. The location and size of the valves is estimated based on additional long axis images through the valves. A semi-automatic grid generator was created to generate the calculation grid. Since the time resolution of the MR scans does not fit the requirements of the CFD calculations a third order bezier approximation scheme was developed to realize a smooth wall boundary and grid movement. The calculation was performed by a Navier-Stokes solver using the arbitrary Lagrange-Euler (ALE) formulation. Results show that during diastole, blood flow through the mitral valve forms an asymmetric jet, leading to an asymmetric development of the initial vortex ring. These flow features are in reasonable agreement with in vivo measurements but also show an extremely high sensitivity to the boundary conditions imposed at the inflow. Changes in the atrial representation severely alter the resulting flow field. These shortcomings will have to be addressed in further studies, possibly by inclusion of the real atrial geometry, and imply additional requirements for the clinical imaging processes.
已开发出一种三维计算流体动力学(CFD)方法来模拟左心室泵血时的血流情况。所提出的方法利用磁共振成像(MRI)技术来提供待模拟心室的患者特异性、随时间变化的几何形状。本研究采用了标准临床成像程序。使用了心脏瓣膜的二维随时间变化的孔口表示法。瓣膜的位置和大小是根据通过瓣膜的额外长轴图像估计的。创建了一个半自动网格生成器来生成计算网格。由于磁共振扫描的时间分辨率不符合CFD计算的要求,因此开发了一种三阶贝塞尔近似方案来实现光滑的壁边界和网格移动。计算由使用任意拉格朗日 - 欧拉(ALE)公式的纳维 - 斯托克斯求解器执行。结果表明,在舒张期,通过二尖瓣的血流形成不对称射流,导致初始涡环的不对称发展。这些流动特征与体内测量结果合理一致,但也显示出对流入处施加的边界条件具有极高的敏感性。心房表示的变化会严重改变所得的流场。这些缺点将必须在进一步的研究中加以解决,可能通过纳入真实的心房几何形状来解决,这也意味着对临床成像过程有额外的要求。