Partitioned solver for strongly coupled fluid–structure interaction

• Published in: Computers & fluids 2013-01, Vol.71, p.306-319
• Main Authors: Habchi, Charbel, Russeil, Serge, Bougeard, Daniel, Harion, Jean-Luc, Lemenand, Thierry, Ghanem, Akram, Valle, Dominique Della, Peerhossaini, Hassan
• Format: Article
• Language: English
• Subjects: Aitken’s under-relaxation; Algorithms; Arbitrary Lagrangian–Eulerian; Benchmarks; Chemical and Process Engineering; Computational fluid dynamics; Elastic large deformation; Engineering Sciences; Finite volume method; Flaps; Fluid flow; Fluid-structure interaction; Fluids; Mathematical analysis; OpenFOAM; Reactive fluid environment; Solvers; St. Venant–Kirchhoff material
• ISSN: 0045-7930; 1879-0747
• DOI: 10.1016/j.compfluid.2012.11.004

► Partitioned solver using Gauss–Seidel scheme with adaptive Aitken’s relaxation. ► This solver is developed using the open-source C++ library OpenFOAM. ► Good agreement with results from the open literature is observed. ► A complex case is considered consisting on two flaps inserted in a pulsatile flow. ► For optimal material properties, a chaotic flow behavior is observed near the flaps. In this work a fluid–structure interaction solver is developed in a partitioned approach using block Gauss–Seidel implicit scheme. Finite volume method is used to discretize the fluid flow problem on a moving mesh in an arbitrary Lagrangian–Eulerian formulation and by using an adaptive time step. The pressure–velocity coupling is performed by using the PIMPLE algorithm, a combination of both SIMPLE and PISO algorithms, which permits the use of larger time steps in a moving mesh. The structural elastic deformation is analyzed in a Lagrangian formulation using the St. Venant–Kirchhoff constitutive law, for non-linear large deformations. The solid structure is discretized by the finite volume method in an iterative segregated approach. The automatic mesh motion solver is based on Laplace smoothing equation with variable mesh diffusion. The strong coupling between the different solvers and the equilibrium on the fluid–structure interface are achieved by using an iterative implicit fixed-point algorithm with dynamic Aitken’s relaxation method. The solver, which is called vorflexFoam, is developed using the open source C++ library OpenFOAM. The solver is validated on two different benchmarks largely used in the open literature. In the first one the structural deformation is induced by incompressibility. The second benchmark consists on a vortex excited elastic flap in a Von Karman vortex street. Finally, a more complex case is studied including two elastic flaps immersed in a pulsatile flow. The present solver detects accurately the interaction between the complex flow structures generated by the flaps and the effect of the flaps oscillations between each other.