A stable partitioned FSI algorithm for rigid bodies and incompressible flow in three dimensions

TL;DR

This paper introduces a stable, second-order accurate partitioned fluid-structure interaction algorithm for three-dimensional rigid bodies and incompressible flow, capable of handling very light or zero-mass bodies without sub-iterations.

Contribution

The paper presents the first stable, second-order accurate partitioned FSI algorithm for 3D rigid bodies with zero mass, extending previous methods with generalized added-damping tensors and surface integral quadrature.

Findings

Algorithm remains stable for bodies of any mass.

Validated with benchmark problems including sphere motion and heart valve interaction.

Demonstrated parallel performance and pressure system conditioning.

Abstract

This paper describes a novel partitioned algorithm for fluid-structure interaction (FSI) problems that couples the motion of rigid bodies and incompressible flow. This is the first partitioned algorithm that remains stable and second-order accurate, without sub-time-step iterations, for very light, and even zero-mass, bodies in three dimensions. This new added-mass partitioned (AMP) algorithm extends the previous developments in [1, 2] by generalizing the added-damping tensors to account for arbitrary three-dimensional rotations, and by employing a general quadrature for the surface integral over a rigid body to derive the discrete AMP interface condition for the fluid pressure. Stability analyses for two three-dimensional model problems show that the algorithm remains stable for bodies of any mass when applied to the relevant model problems. The resulting AMP algorithm is implemented in parallel using a moving composite grid framework to treat one or more rigid bodies in complex three-dimensional configurations. The new three-dimensional algorithm is verified and validated though several benchmark problems, including the motion of a sphere in a viscous incompressible fluid and the interaction of a bi-leaflet mechanical heart valve and a pulsating fluid. Numerical simulations confirm the predictions of the stability analysis even for complex problems, and show that the AMP algorithm remains stable, without sub-iterations, for light and even zero-mass three-dimensional rigid bodies of general shape. These benchmark problems are further used to examine the parallel performance of the algorithm and to investigate the conditioning of the linear system for the pressure including the newly derived AMP interface conditions.