A parallel solver for vortex-induced vibrations at zero mass ratio

We present a parallel strategy for the study of the vortex-induced vibrations (VIV) of a rigid, massless cylinder mounted on translational springs, in two and three dimensions. Without mass or damping, its motion is governed solely by the fluid forces and the spring restoring force, and is an immediate response to the flow. We focus on low spring stiffness, yielding high reduced velocity and maximal oscillation amplitude. A monolithic finite element method in an arbitrary Lagrangian-Eulerian framework with implicit time stepping is proposed, solving the fluid flow, cylinder dynamics, and mesh motion simultaneously. Mesh movement is modelled with a linear elastic pseudo-solid with variable Lam\'e coefficients, and the no-slip condition on the cylinder is enforced with a Lagrange multiplier coupled to the position degrees of freedom; three strategies are explored to enforce this force-position coupling in a distributed framework. To the best of our knowledge, this work is the first to propose a distributed solver for three-dimensional fluid-structure interaction at zero mass ratio. The solver is verified with manufactured solutions, and exhibits second-order accuracy in time. A numerical study of two- and three-dimensional VIV for Reynolds numbers from 100 to 300 reproduces the VIV responses obtained in previous two-dimensional studies, and predicts nonperiodic 3D oscillations starting at $Re=300$. Using the direct solver MUMPS, the proposed method scales to up to 768 compute cores.