Curvilinear immersed boundary method for simulating fluid structure interaction with complex 3D rigid bodies

The sharp-interface CURVIB approach of Ge and Sotiropoulos [L. Ge, F. Sotiropoulos, A Numerical Method for Solving the 3D Unsteady Incompressible Navier-Stokes Equations in Curvilinear Domains with Complex Immersed Boundaries, Journal of Computational Physics 225 (2007) 1782-1809] is extended to simulate fluid structure interaction (FSI) problems involving complex 3D rigid bodies undergoing large structural displacements. The FSI solver adopts the partitioned FSI solution approach and both loose and strong coupling strategies are implemented. The interfaces between immersed bodies and the fluid are discretized with a Lagrangian grid and tracked with an explicit front-tracking approach. An efficient ray-tracing algorithm is developed to quickly identify the relationship between the background grid and the moving bodies. Numerical experiments are carried out for two FSI problems: vortex induced vibration of elastically mounted cylinders and flow through a bileaflet mechanical heart valve at physiologic conditions. For both cases the computed results are in excellent agreement with benchmark simulations and experimental measurements. The numerical experiments suggest that both the properties of the structure (mass, geometry) and the local flow conditions can play an important role in determining the stability of the FSI algorithm. Under certain conditions unconditionally unstable iteration schemes result even when strong coupling FSI is employed. For such cases, however, combining the strong-coupling iteration with under-relaxation in conjunction with the Aitken's acceleration technique is shown to effectively resolve the stability problems. A theoretical analysis is presented to explain the findings of the numerical experiments. It is shown that the ratio of the added mass to the mass of the structure as well as the sign of the local time rate of change of the force or moment imparted on the structure by the fluid determine the stability and convergence of the FSI algorithm. The stabilizing role of under-relaxation is also clarified and an upper bound of the required for stability under-relaxation coefficient is derived.

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

[2]  Eric Haines,et al.  Point in Polygon Strategies , 1994, Graphics Gems.

[3]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[4]  P. Moin,et al.  Eddies, streams, and convergence zones in turbulent flows , 1988 .

[5]  Jung Il Choi,et al.  An immersed boundary method for complex incompressible flows , 2007, J. Comput. Phys..

[6]  Tomas Akenine-Möller,et al.  Fast, Minimum Storage Ray-Triangle Intersection , 1997, J. Graphics, GPU, & Game Tools.

[7]  Yannis Kallinderis,et al.  Unsteady Flow Structure Interaction for Incompressible Flows Using Deformable Hybrid Grids , 1998 .

[8]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[9]  F. Sotiropoulos,et al.  A hybrid Cartesian/immersed boundary method for simulating flows with 3D, geometrically complex, moving bodies , 2005 .

[10]  Wu,et al.  THE LIFT FORCE ON A VIBRATING CYLINDER IN A CURRENT , 1989 .

[11]  J. Ferziger,et al.  A ghost-cell immersed boundary method for flow in complex geometry , 2002 .

[12]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[13]  Hélène A. Simon,et al.  Vorticity dynamics of a bileaflet mechanical heart valve in an axisymmetric aorta , 2007 .

[14]  Gianni Pedrizzetti,et al.  Flow-driven opening of a valvular leaflet , 2006, Journal of Fluid Mechanics.

[15]  S. Green,et al.  Simulating the Motion of Flexible Pulp Fibres Using the Immersed Boundary Method , 1998 .

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

[17]  Jack Lemmon,et al.  A numerical simulation of mechanical heart valve closure fluid dynamics. , 2002, Journal of biomechanics.

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

[19]  H. S. Udaykumar,et al.  Two-Dimensional Dynamic Simulation of Platelet Activation During Mechanical Heart Valve Closure , 2006, Annals of Biomedical Engineering.

[20]  Robert Dillon,et al.  Modeling Biofilm Processes Using the Immersed Boundary Method , 1996 .

[21]  P. Tallec,et al.  Load and motion transfer algorithms for fluid/structure interaction problems with non-matching discrete interfaces: Momentum and energy conservation, optimal discretization and application to aeroelasticity , 1998 .

[22]  C. Williamson,et al.  Vortex-Induced Vibrations , 2004, Wind Effects on Structures.

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

[24]  R. Cheng,et al.  Three-Dimensional Fluid-Structure Interaction Simulation of Bileaflet Mechanical Heart Valve Flow Dynamics , 2004, Annals of Biomedical Engineering.

[25]  Bruce M. Irons,et al.  A version of the Aitken accelerator for computer iteration , 1969 .

[26]  R. Ho.  Algebraic Topology , 2022 .

[27]  Ekkehard Ramm,et al.  Accelerated iterative substructuring schemes for instationary fluid-structure interaction , 2001 .

[28]  H. Udaykumar,et al.  Sharp interface Cartesian grid method I: An easily implemented technique for 3D moving boundary computations , 2005 .

[29]  E. Balaras,et al.  A general reconstruction algorithm for simulating flows with complex 3D immersed boundaries on Cartesian grids , 2003 .

[30]  C. Conca,et al.  Added mass and damping in fluid-structure interaction , 1997 .

[31]  Fotis Sotiropoulos,et al.  Flow in Prosthetic Heart Valves: State-of-the-Art and Future Directions , 2005, Annals of Biomedical Engineering.

[32]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[33]  Haecheon Choi,et al.  Immersed boundary method for flow around an arbitrarily moving body , 2006, J. Comput. Phys..

[34]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .

[35]  C. Peskin,et al.  Simulation of a Flapping Flexible Filament in a Flowing Soap Film by the Immersed Boundary Method , 2002 .

[36]  Jan Vierendeels,et al.  Analysis and stabilization of fluid-structure interaction algorithm for rigid-body motion , 2005 .

[37]  P. Tallec,et al.  Fluid structure interaction with large structural displacements , 2001 .

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

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

[40]  Yannis Kallinderis,et al.  Strongly coupled flow/structure interactions with a geometrically conservative ALE scheme on general hybrid meshes , 2006, J. Comput. Phys..

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

[42]  R. Blevins,et al.  Flow-Induced Vibration , 1977 .

[43]  Z. Jane Wang,et al.  An immersed interface method for simulating the interaction of a fluid with moving boundaries , 2006, J. Comput. Phys..

[44]  W. Shyy,et al.  Computation of Solid-Liquid Phase Fronts in the Sharp Interface Limit on Fixed Grids , 1999 .

[45]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[46]  Jungwoo Kim,et al.  An immersed-boundary finite-volume method for simulations of flow in complex geometries , 2001 .

[47]  Fotis Sotiropoulos,et al.  A numerical method for solving the 3D unsteady incompressible Navier-Stokes equations in curvilinear domains with complex immersed boundaries , 2007, J. Comput. Phys..

[48]  Ajit P Yoganathan,et al.  Fluid mechanics of heart valves. , 2004, Annual review of biomedical engineering.

[49]  Turgut Sarpkaya,et al.  A critical review of the intrinsic nature of vortex-induced vibrations , 2004 .

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

[51]  Parviz Moin,et al.  The structure of two-dimensional separation , 1990, Journal of Fluid Mechanics.

[52]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[53]  Miguel Angel Fernández,et al.  A Newton method using exact jacobians for solving fluid-structure coupling , 2005 .

[54]  Boo Cheong Khoo,et al.  An immersed interface method for viscous incompressible flows involving rigid and flexible boundaries , 2006, J. Comput. Phys..