A parallel monolithic algorithm for the numerical simulation of large-scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large-scale fluid structure interaction problems. The governing incompressible Navier-Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian-Eulerian formulation-based side-centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant-Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large-scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one-level restricted additive Schwarz preconditioner with a block-incomplete factorization within each partitioned sub-domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd. A novel FSI algorithm is proposed for the large-scale simulation of fluid-structure interaction problems in a fully coupled form. A special attention is given to satisfy both the local and global discrete geometric conservation law (DGCL) in order to conserve the total fluid volume/mass in machine precision. Large-scale numerical results are presented for several classical FSI benchmark problems.