Partitioned block-Gauss-Seidel coupling for dynamic fluid-structure interaction

This paper presents a fully coupled three-dimensional solver for the analysis of time-dependent fluid-structure interaction. A partitioned time marching algorithm is employed for the solution of the time-dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective block-Gauss-Seidel iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit dual time stepping technique with explicit multi-stage Runga-Kutta time stepping in pseudo-time and an ALE formulation for moving boundaries. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. Test cases are presented to benchmark the algorithm and to demonstrate the potential applications of this method.

[1]  Charbel Farhat,et al.  Second-order time-accurate and geometrically conservative implicit schemes for flow computations on unstructured dynamic meshes , 1999 .

[2]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[3]  J. Halleux,et al.  An arbitrary lagrangian-eulerian finite element method for transient dynamic fluid-structure interactions , 1982 .

[4]  Arthur Rizzi,et al.  Computation of inviscid incompressible flow with rotation , 1985, Journal of Fluid Mechanics.

[5]  E. Turkel,et al.  Preconditioned methods for solving the incompressible low speed compressible equations , 1987 .

[6]  V. Venkatakrishnan,et al.  IMPLICIT METHOD FOR THE COMPUTATION OF UNSTEADY FLOWS ON UNSTRUCTURED GRIDS , 1995 .

[7]  R. D. Wood,et al.  Nonlinear Continuum Mechanics for Finite Element Analysis , 1997 .

[8]  O. C. Zienkiewicz,et al.  The Finite Element Method for Solid and Structural Mechanics , 2013 .

[9]  T. Tezduyar,et al.  Fluid–structure Interaction Modeling of Aneurysmal Conditions with High and Normal Blood Pressures , 2006 .

[10]  Oubay Hassan,et al.  Parallel Multigrid Detached Eddy Simulation Algorithm for Three-Dimensional Unsteady Incompressible Flows on Unstructured Meshes , 2006 .

[11]  H. Matthies,et al.  Partitioned Strong Coupling Algorithms for Fluid-Structure-Interaction , 2003 .

[12]  J. Womersley.  Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known , 1955, The Journal of physiology.

[13]  Ning Qin,et al.  Fast dynamic grid deformation based on Delaunay graph mapping , 2006 .

[14]  Frederic Blom,et al.  A monolithical fluid-structure interaction algorithm applied to the piston problem , 1998 .

[15]  Oubay Hassan,et al.  The simulation of 3D unsteady incompressible flows with moving boundaries on unstructured meshes , 2008 .

[16]  D. Dinkler,et al.  A monolithic approach to fluid–structure interaction using space–time finite elements , 2004 .

[17]  C. Farhat,et al.  Two efficient staggered algorithms for the serial and parallel solution of three-dimensional nonlinear transient aeroelastic problems , 2000 .

[18]  Nigel P. Weatherill,et al.  Unsteady flow simulation using unstructured meshes , 2000 .

[19]  C. Peskin.  Flow patterns around heart valves: A numerical method , 1972 .

[20]  Peter Wriggers,et al.  Large strain analysis of soft biological membranes: Formulation and finite element analysis , 1996 .

[21]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[22]  Hydrostatic fluid loading in non‐linear finite element analysis , 2004 .

[23]  J. Bonet,et al.  An alternating digital tree (ADT) algorithm for 3D geometric searching and intersection problems , 1991 .

[24]  Guru P. Guruswamy,et al.  A review of numerical fluids/structures interface methods for computations using high-fidelity equations , 2002 .

[25]  Hervé Guillard,et al.  Godunov type method on non-structured meshes for three-dimensional moving boundary problems , 1994 .

[26]  Charbel Farhat,et al.  Provably second-order time-accurate loosely-coupled solution algorithms for transient nonlinear computational aeroelasticity , 2006 .

[27]  T. Tezduyar,et al.  A new strategy for finite element computations involving moving boundaries and interfaces—the deforming-spatial-domain/space-time procedure. I: The concept and the preliminary numerical tests , 1992 .

[28]  P. Thomas,et al.  Geometric Conservation Law and Its Application to Flow Computations on Moving Grids , 1979 .

[29]  Nigel P. Weatherill,et al.  A method for time accurate turbulent compressible fluid flow simulation with moving boundary components employing local remeshing , 2007 .

[30]  J. Batina.  Unsteady Euler airfoil solutions using unstructured dynamic meshes , 1989 .

[31]  D. Peric,et al.  A computational framework for fluid–structure interaction: Finite element formulation and applications , 2006 .

[32]  K. Bathe.  Finite Element Procedures , 1995 .

[33]  Wing Kam Liu,et al.  Lagrangian-Eulerian finite element formulation for incompressible viscous flows☆ , 1981 .

[34]  Charbel Farhat,et al.  A three-dimensional torsional spring analogy method for unstructured dynamic meshes , 2002 .

[35]  N. Weatherill,et al.  Efficient three‐dimensional Delaunay triangulation with automatic point creation and imposed boundary constraints , 1994 .

[36]  Charbel Farhat,et al.  Geometric conservation laws for flow problems with moving boundaries and deformable meshes, and their impact on aeroelastic computations , 1996 .

[37]  Peter Wriggers,et al.  A note on finite‐element implementation of pressure boundary loading , 1991 .

[38]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .