Time consistent fluid structure interaction on collocated grids for incompressible flow

Published on Jan 1, 2016in Computer Methods in Applied Mechanics and Engineering5.763
· DOI :10.1016/J.CMA.2015.09.025
T. Gillebaart5
Estimated H-index: 5
(TU Delft: Delft University of Technology),
D.S. Blom5
Estimated H-index: 5
(TU Delft: Delft University of Technology)
+ 1 AuthorsHester Bijl28
Estimated H-index: 28
(TU Delft: Delft University of Technology)
Sources
Abstract
Abstract Consistent time integration on collocated grids for incompressible flow has been studied for static grids using the PISO method, in which the dependencies on time-step size and under-relaxation has been studied in detail. However, for moving grids a detailed description is still missing. Therefore, a step by step analysis of a time consistent fluid–structure interaction (FSI) method for incompressible flow on collocated grids is presented. The method consist of: face normal and area correction for moving grids, treatment of velocity boundary conditions for no-slip walls, time integration of structure equations and fluid force interpolation to structure. The basis of the method is the PISO method of the incompressible Navier–Stokes equations. Time consistency on static grids is shown first, after which time consistency on moving grids is described and analyzed. For moving grids consistent time integration is described in two ways: (1) constructing the face velocities from a normal and tangential component, and (2) correcting the face flux with a face normal and face area correction. For both descriptions the general formulation for the backward differencing schemes (BDF) are given and the correct behavior of the first, second and third order schemes is demonstrated by means of an academic test case (circular cavity flow). Additionally, the (force) coupling from the fluid to the structure is discussed in detail for combining a fourth order explicit Runge Kutta scheme for the structure with a BDF scheme for the fluid. Three interpolations for the fluid forces are shown, which result in either a first order or second order FSI scheme. Third order FSI is demonstrated when the third order BDF scheme is applied on both the structure and fluid equations. Also, under-relaxation for the fluid equations is considered and it is demonstrated that the order of the three BDF schemes are independent of the under-relaxation factor. Finally, the proposed method of time consistent FSI on collocated grids for incompressible flows is demonstrated by applying it to a three-dimensional flow over an elastic structure in a channel and the cylinder flap FSI benchmark case of Turek and Hron.
📖 Papers frequently viewed together
11 Citations
1990
1 Author (Laszlo Fuchs)
1 Citations
11 Citations
References26
Newest
#1V. Kazemi-Kamyab (TU Delft: Delft University of Technology)H-Index: 4
#2A.H. van Zuijlen (TU Delft: Delft University of Technology)H-Index: 14
Last. Hester Bijl (TU Delft: Delft University of Technology)H-Index: 28
view all 3 authors...
Abstract The application of a family of high order implicit Runge–Kutta time integration schemes, namely the explicit first-stage singly diagonally implicit Runge–Kutta schemes (ESDIRK), to cell-centered collocated finite volume discretization of the unsteady incompressible Navier–Stokes is considered. Although achieving computational efficiency relative to commonly used second order implicit schemes has been the motivating factor, this study focuses on temporal order analysis of the high order ...
12 CitationsSource
#1Dmitry Kolmogorov (DTU: Technical University of Denmark)H-Index: 3
#2Wen Zhong Shen (DTU: Technical University of Denmark)H-Index: 32
Last. Jens Nørkær SørensenH-Index: 66
view all 4 authors...
To increase the convergence rate of SIMPLE-like algorithms on collocated grids, a compatibility condition between mass flux interpolation methods and SIMPLE-like algorithms is presented. Results of unsteady flow computations show that the SIMPLEC algorithm, when obeying the compatibility condition, may obtain up to 35% higher convergence rate as compared to the standard SIMPLEC algorithm. Two new interpolation methods, fully compatible with the SIMPLEC algorithm, are presented and compared with ...
5 CitationsSource
#1Bradley Froehle (University of California, Berkeley)H-Index: 3
#2Per-Olof Persson (University of California, Berkeley)H-Index: 27
Abstract We present a high-order accurate scheme for coupled fluid–structure interaction problems. The fluid is discretized using a discontinuous Galerkin method on unstructured tetrahedral meshes, and the structure uses a high-order volumetric continuous Galerkin finite element method. Standard radial basis functions are used for the mesh deformation. The time integration is performed using a partitioned approach based on implicit–explicit Runge–Kutta methods. The resulting scheme fully decoupl...
36 CitationsSource
#1Philip Cardiff (UCD: University College Dublin)H-Index: 16
#2Aleksandar Karac (University of Zenica)H-Index: 12
Last. Alojz Ivankovic (UCD: University College Dublin)H-Index: 23
view all 3 authors...
Abstract This paper describes a finite volume method for orthotropic bodies with general principal material directions undergoing large strains and large rotations. The governing and constitutive relations are presented and the employed updated Lagrangian mathematical model is outlined. In order to maintain equivalence with large strain total Lagrangian methods, the constitutive stiffness tensor is updated transforming the principal material directions to the deformed configuration. Discretisati...
35 CitationsSource
In this review article, the focus is on partitioned simulation techniques for strongly coupled fluid-structure interaction problems, especially on techniques which use at least one of the solvers as a black box. First, a number of analyses are reviewed to explain why Gauss–Seidel coupling iterations converge slowly or not at all for fluid-structure interaction problems with strong coupling. This provides the theoretical basis for the fast convergence of quasi-Newton and multi-level techniques. S...
74 CitationsSource
#1Ming-Chen Hsu (UCSD: University of California, San Diego)H-Index: 39
#2Yuri Bazilevs (UCSD: University of California, San Diego)H-Index: 77
In this paper we present our aerodynamics and fluid---structure interaction (FSI) computational techniques that enable dynamic, fully coupled, 3D FSI simulation of wind turbines at full scale, and in the presence of the nacelle and tower (i.e., simulation of the "full machine"). For the interaction of wind and flexible blades we employ a nonmatching interface discretization approach, where the aerodynamics is computed using a low-order finite-element-based ALE-VMS technique, while the rotor blad...
218 CitationsSource
In this work, we present an adaptive finite element method for the numerical simulation of stationary fluid–structure interaction problems. The coupled system is given in a variational and monolithic Arbitrary Lagrangian Eulerian framework. We derive methods for goal-oriented error estimation and mesh adaptation with the dual weighted residual method. Key to applying this error estimator is the underlying canonic variational formulation of the fluid–structure interaction problem by mapping the f...
56 CitationsSource
#1Željko Tuković (University of Zagreb)H-Index: 14
#2Hrvoje Jasak (University of Zagreb)H-Index: 27
This paper describes a moving mesh interface tracking method implemented in OpenFOAM for simulating three-dimensional (3-D) incompressible and immiscible two-phase interfacial fluid flows with dominant surface tension forces. Collocated finite volume (FV) method is used for spatial discretisation of Navier–Stokes equations on moving polyhedral mesh. The mesh consists of two parts separated on interface. Fluid flow is solved on each mesh separately and coupling is accomplished in an iterative man...
110 CitationsSource
#1Joris Degroote (UGent: Ghent University)H-Index: 27
#2Robby Haelterman (Royal Military Academy)H-Index: 7
Last. Jan Vierendeels (UGent: Ghent University)H-Index: 40
view all 5 authors...
Partitioned simulations of fluid-structure interaction can be solved for the interface's position with Newton-Raphson iterations but obtaining the exact Jacobian is impossible if the solvers are ''black boxes''. It is demonstrated that only an approximate Jacobian is needed, as long as it describes the reaction to certain components of the error on the interface's position. Based on this insight, a quasi-Newton coupling algorithm with an approximation for the inverse of the Jacobian (IQN-ILS) ha...
116 CitationsSource
#1A. de Boer (TU Delft: Delft University of Technology)H-Index: 5
#2M. S. van der Schoot (TU Delft: Delft University of Technology)H-Index: 1
Last. Hester Bijl (TU Delft: Delft University of Technology)H-Index: 28
view all 3 authors...
A new mesh movement algorithm for unstructured grids is developed which is based on interpolating displacements of the boundary nodes to the whole mesh with radial basis functions (RBF's). A small system of equations, only involving the boundary nodes, has to be solved and no grid-connectivity information is needed. The method can handle large mesh deformations caused by translations, rotations and deformations, both for 2D and 3D meshes. However, the performance depends on the used RBF. The bes...
472 CitationsSource
Cited By5
Newest
An increasing number of the offshore structures in the North Sea are being operated in a life extension phase, which means that their original design life has been exceeded. As the structures age, they are deteriorating, mainly due to fatigue and corrosion. In order to ensure that the structural integrity is maintained during the life extension phase, the remaining fatigue life of the structures needs to be assessed. While simple and conservative approximations are used in fatigue design, it wil...
Source
#1Dante De Santis (Nuclear Research and Consultancy Group)H-Index: 7
#2Afaque Shams (Nuclear Research and Consultancy Group)H-Index: 13
Abstract In this work, an advanced solver for the simulation of Fluid-Structure Interaction (FSI) problems for nuclear applications is presented. The numerical framework, called as NRG-FSIFOAM, uses a Finite Volume (FV) fluid solver and a Finite Element (FE) structural solver coupled together with a partitioned approach. State-of-the-art implicit coupling algorithms are used in order to increase the accuracy and the stability of the solver for low solid-to-fluid density ratio, i.e. strongly coup...
3 CitationsSource
#1Željko Tuković (University of Zagreb)H-Index: 14
#2Martina Bukac (ND: University of Notre Dame)H-Index: 11
Last. Alojz IvankovicH-Index: 23
view all 5 authors...
This paper describes a self-contained, partitioned fluid–structure interaction solver based on a finite volume discretisation. The incompressible fluid flow is described by the Navier–Stokes equations in the arbitrary Lagrangian–Eulerian form and the solid deformation is described by the St. Venant-Kirchhoff hyperelastic model in the total Lagrangian form. Both fluid and solid are discretised in space using the second-order accurate cell-centred finite volume method, and temporal discretisation ...
3 CitationsSource
#1Željko Tuković (University of Zagreb)H-Index: 14
#2Milovan Perić (University of Duisburg-Essen)H-Index: 1
Last. Hrvoje Jasak (University of Zagreb)H-Index: 27
view all 3 authors...
A new variant of the non-iterative PISO-algorithm for the solution of implicitly discretised fluid flow equations is proposed. The governing equations are discretised in space using a second-order accurate finite volume method. A second-order accurate three-level implicit temporal discretisation scheme is used to discretise the momentum equation, where the non-linear convection term is linearised using a second-order explicit approximation of mass flux at the new time level, based on the results...
11 CitationsSource
#1Kristen RegeH-Index: 5
#2Bjørn H. HjertagerH-Index: 29
Turbulent flow around flexible structures is likely to induce structural vibrations which may eventually lead to fatigue failure. In order to assess the fatigue life of these structures, it is necessary to take the action of the flow on the structure into account, but also the influence of the vibrating structure on the fluid flow. This is achieved by performing fluid-structure interaction (FSI) simulations. In this work, we have investigated the capability of a FSI toolkit for the finite volume...
7 CitationsSource