The Immersed Structural Potential Method for haemodynamic applications
The Immersed Structural Potential Method for haemodynamic applications
复制标题
DOI：
10.1016/j.jcp.2010.08.005
复制
发表时间：
2010-11-01
影响因子：
4.1
通讯作者：
Hassan, O.
中科院分区：
文献类型：
--
作者：
Gil, A. J.;Carreno, A. Arranz;Hassan, O.
在本文中，一种新的流体-结构相互作用浸没计算方法，基于原始的浸没边界方法（IBM）[1]，概述了模拟心血管现象，特别是心脏瓣膜相关问题的最终目标。这种浸入技术的主要特点是表示任何可变形或刚体浸入不可压缩的粘性流场作为Navier-Stokes方程中的动量强迫源。浸入式配方中的一些缺点仍然需要进一步的研究和改进，包括由插值/扩散过程引起的过度数值扩散，需要包括更准确地描述心血管组织的性质的现实粘弹性复合材料本构模型，以及需要更有效地捕获在流体-结构界面处开发的应力。遵循与原始IBM相同的原理，本文导出了一种更复杂的公式，即“浸没结构势法（ISPM）”。介绍的方法提出了一种替代方法来计算等效的流体-结构相互作用力在流体网格，占一个复杂的粘弹性纤维增强本构模型，以更好地描述力学的心血管组织和利用一种新的时间积分方法的变形梯度张量的计算，确保符合不可压缩性约束。一系列的数值例子，以证明这种新方法的鲁棒性和适用性。(C)2010年爱思唯尔公司All rights reserved.
In this paper, a new fluid-structure interaction immersed computational methodology, based upon the original Immersed Boundary Method (IBM) [1] is outlined with the final aim of modelling cardiovascular phenomena, specifically, heart valve related problems. The principal characteristic of such immersed techniques is the representation of any deformable or rigid body immersed within an incompressible viscous flow field as a momentum forcing source in the Navier-Stokes equations. A number of shortcomings within the immersed formulation still require further investigation and improvement, including the excessive numerical diffusion caused by the interpolation/spreading process, the need to include realistic viscoelastic composite constitutive models describing more accurately the nature of cardiovascular tissues and also the need to capture more effectively stresses developed at the fluid-structure interface. By following the same philosophy as the original IBM, a more sophisticated formulation is derived in this paper, the "Immersed Structural Potential Method (ISPM)". The method introduced presents an alternative approach to compute the equivalent fluid-structure interaction forces at the fluid mesh, accounts for a sophisticated viscoelastic fibre-reinforced constitutive model to better describe the mechanics of cardiovascular tissues and utilises a novel time-integration methodology for the computation of the deformation gradient tensor which ensures compliance with the incompressibility constraint. A series of numerical examples will be presented in order to demonstrate the robustness and applicability of this new methodology. (C) 2010 Elsevier Inc. All rights reserved.