Abstract
An efficient and scalable numerical method for massively parallel computing of fluid-structure interaction systems has been developed for biomedical applications. The basic equation set is monolithically formulated in a full Eulerian framework, and solved on a fixed mesh. The discretization is based on a finite-difference volume-of-fluid method, corresponding to a standard numerical approach for simulating incompressible multiphase flows. Further, a novel algorithm based on an artificial compressibility method is proposed. The parameters therein are adaptively determined to guarantee the computed velocity field to be nearly incompressible. The applicability and effectiveness of the proposed method are demonstrated.