Partitioned fluid-structure interaction simulations often use different time-integration schemes to discretize the different sub-problems. As such, the flow and structural equations can be solved with schemes that are particularly suited for each individual problem. However, using incompatible schemes, these simulations can encounter stability problems. In this research an analytical stability analysis is performed for a model of blood flow in an artery. The backward Euler scheme is used for the time discretization of the flow equations. For the structure two schemes are used: the BE scheme and the Hilber-Hughes-Taylor operator in which the numerical damping is controlled by a single parameter alpha. The influence of this parameter and some physiological parameters on the stability and the damping of the spurious modes is studied. According to this analysis, the combination of the BE and HHT scheme is stable, but the wave number, the numerical damping and the flow and structural density can affect the damping of the spurious modes considerably. To verify the analytical results, a numerical study is performed using nonlinear two-dimensional axisymmetric FSI simulations.
展开▼
