A partitioned coupling approach for dynamic fluid–structure interaction with applications to biological membranes

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 to the time-dependent coupled discretized problem, thus enabling the use of highly developed, robust and well-tested solvers for each field. Coupling of the fields is achieved through a conservative transfer of information at the fluid-structure interface. An implicit coupling is achieved when the solutions to the fluid and structure subproblems are cycled at each time step until convergence is reached. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit dual time-stepping technique with an explicit multistage Runga-Kutta time stepping in pseudo-time and arbitrary Lagrangian-Eulerian formulation for the 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 optimization method. Various test cases are presented for benchmarking and to demonstrate the potential applications of this method.

[1]  F P T Baaijens,et al.  A three-dimensional computational analysis of fluid-structure interaction in the aortic valve. , 2003, Journal of biomechanics.

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

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

[4]  Peter Hansbo,et al.  The characteristic streamline diffusion method for the time-dependent incompressible Navier-Stokes equations , 1992 .

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

[6]  J. Z. Zhu,et al.  The finite element method , 1977 .

[7]  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 .

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

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

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

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

[12]  Tayfun E. Tezduyar,et al.  Fluid-structure interactions of a cross parachute: Numerical simulation , 2001 .

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

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

[15]  E. Ramm,et al.  Artificial added mass instabilities in sequential staggered coupling of nonlinear structures and incompressible viscous flows , 2007 .

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

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

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

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

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

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

[22]  Fabio Nobile,et al.  Added-mass effect in the design of partitioned algorithms for fluid-structure problems , 2005 .

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

[24]  Thomas J. R. Hughes,et al.  Computational investigations in vascular disease , 1996 .

[25]  T. Hughes,et al.  Space-time finite element methods for elastodynamics: formulations and error estimates , 1988 .

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

[27]  David Farrell,et al.  Immersed finite element method and its applications to biological systems. , 2006, Computer methods in applied mechanics and engineering.

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

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

[30]  H. Saunders.  Book Reviews : NUMERICAL METHODS IN FINITE ELEMENT ANALYSIS K.-J. Bathe and E.L. Wilson Prentice-Hall, Inc, Englewood Cliffs, NJ , 1978 .

[31]  C. Farhat,et al.  Mixed explicit/implicit time integration of coupled aeroelastic problems: Three‐field formulation, geometric conservation and distributed solution , 1995 .

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

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

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

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

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

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

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

[39]  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.

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

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

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

[43]  Charbel Farhat,et al.  The discrete geometric conservation law and the nonlinear stability of ALE schemes for the solution of flow problems on moving grids , 2001 .