A convergence study of a new partitioned fluid-structure interaction algorithm based on fictitious mass and damping

We develop, analyze and validate a new method for simulating fluid-structure interactions (FSIs), which is based on fictitious mass and fictitious damping in the structure equation. We employ a partitioned method for the fluid and structure motions in conjunction with sub-iteration and Aitken relaxation. In particular, the use of such fictitious parameters requires sub-iterations in order to reduce the induced error in addition to the local temporal truncation error. To this end, proper levels of tolerance for terminating the sub-iteration procedure have been obtained in order to recover the formal order of temporal accuracy. For the coupled FSI problem, these fictitious terms have a significant effect, leading to better convergence rate and hence substantially smaller number of sub-iterations. Through analysis we identify the proper range of these parameters, which we then verify by corresponding numerical tests. We implement the method in the context of spectral element discretization, which is more sensitive than low-order methods to numerical instabilities arising in the explicit FSI coupling. However, the method we present here is simple and general and hence applicable to FSI based on any other discretization. We demonstrate the effectiveness of the method in applications involving 2D vortex-induced vibrations (VIV) and in 3D flexible arteries with structural density close to blood density. We also present 3D results for a patient-specific aneurysmal flow under pulsatile flow conditions examining, in particular, the sensitivity of the results on different values of the fictitious parameters. (C) 2011 Elsevier Inc. All rights reserved.