Subdivision-stabilised immersed b-spline finite elements for moving boundary flows

An immersed finite element method is presented to compute flows with complex moving boundaries on a fixed Cartesian grid. The viscous, incompressible fluid flow equations are discretized with b-spline basis functions. The two-scale relation for b-splines is used to implement an elegant and efficient technique to satisfy the LBB condition. On non-grid-aligned fluid domains and at moving boundaries, the boundary conditions are enforced with a consistent penalty method as originally proposed by Nitsche. In addition, a special extrapolation technique is employed to prevent the loss of numerical stability in presence of arbitrarily small cut-cells. The versatility and accuracy of the proposed approach is demonstrated by means of convergence studies and comparisons with previous experimental and computational investigations.

A subdivision-based implementation of the hierarchical b-spline finite element method

A novel technique is presented to facilitate the implementation of hierarchical b-splines and their interfacing with conventional finite element implementations. The discrete interpretation of the two-scale relation, as common in subdivision schemes, is used to establish algebraic relations between the basis functions and their coefficients on different levels of the hierarchical b-spline basis. The subdivision projection technique introduced allows us first to compute all element matrices and vectors using a fixed number of same-level basis functions. Their subsequent multiplication with subdivision matrices projects them, during the assembly stage, to the correct levels of the hierarchical b-spline basis. The proposed technique is applied to convergence studies of linear and geometrically nonlinear problems in one, two and three space dimensions. © 2012 Elsevier B.V.

A fixed-grid b-spline finite element technique for fluid-structure interaction

We present a fixed-grid finite element technique for fluid-structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b-spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision-stabilisation technique is used to ensure inf-sup stability. The beam equations are discretised with b-splines and the shell equations with subdivision basis functions, both leading to a rotation-free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet-Robin partitioning scheme, and the fluid equations are solved with a pressure-correction method. Auxiliary techniques employed for improving numerical robustness include the level-set based implicit representation of the structure interface on the fluid grid, a cut-cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. © 2013 John Wiley & Sons, Ltd.