A fictitious domain/distributed Lagrange multiplier based fluid–structure interaction scheme with hierarchical B-Spline grids

• Published in: Computer methods in applied mechanics and engineering 2016-04, Vol.301, p.1-27
• Main Authors: Kadapa, C., Dettmer, W.G., Perić, D.
• Format: Article
• Language: English
• Subjects: Computational fluid dynamics; Cylinders; Distributed Lagrange multipliers; Fictitious domain methods; Flow-induced vibrations; Fluid flow; Fluid-structure interaction; Hierarchical B-Splines; Immersed boundary methods; Lagrange multipliers; Mathematical analysis; Mathematical models; Navier-Stokes equations
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2015.12.023

We present a numerical scheme for fluid–structure interaction based on hierarchical B-Spline grids and fictitious domain/distributed Lagrange multipliers. The incompressible Navier–Stokes equations are solved over a Cartesian grid discretised with B-Splines. The fluid grid near the immersed solids is refined locally using hierarchical B-Splines. The immersed solid is modelled as geometrically-exact beam discretised with standard linear Lagrange shape functions. The kinematic constraint at the fluid–solid interface is enforced with distributed Lagrange multipliers. The unconditionally-stable and second-order accurate generalised-α method is used for integration in time for both the fluid and solid domains. A fully-implicit and fully-coupled solution scheme is developed by using the Newton–Raphson method to solve the non-linear system of equations obtained with Galerkin weak formulation. First, the spatial and temporal convergence of the proposed scheme is assessed by studying steady and unsteady flow past a fixed cylinder. Then, the scheme is applied to several benchmark problems to demonstrate the efficiency and robustness of the proposed scheme. The results obtained with the present scheme are compared with the reference values. •We present a new numerical scheme for fluid-structure interaction based on B-Spline cartesian grids. The fluid grid near the immersed solids can be refined locally using hierarchical B-Splines.•The solid is assumed to be infinitely thin strcture and is modelled using geometrically-exact beam formulation.•A second-order and unconditionally stable generalised-α scheme is used for time integration of both the fluid and solid domains. This allows us use of larger time increments without losing accuracy.•The proposed scheme has been applied to several fluid-flexible body interaction examples involving moderate and large deformations of the structure as well as simulating dynamically complex bi-stable phenomenon observed in filaments and swimming organisms.