Numerical Algorithms for Tracking Dynamic Fluid-Structure Interfaces in Embedded / Immersed Boundary Methods ∗

In any embedded/immersed boundary method, the embedded/immersed boundary needs to be tracked with respect to the non body-conforming Computational Fluid Dynamics (CFD) grid. When dealing with dynamic fluid-structure interaction problems with complex structural geometry, a robust, accurate, and efficient interface tracking algorithm becomes an important component of the computational framework. Existing methods developed in the CFD community have primarily focused on tracking closed boundaries with respect to uniform or non-uniform Cartesian CFD grids. In this work, two robust and efficient algorithms are presented for tracking embedded/immersed fluid-structure interfaces with respect to three-dimensional, arbitrary (i.e. structured or unstructured) CFD grids. First, a projection-based approach can be used when the interface is closed and fully resolved by the CFD grid, as it relies on the interface to separate the two distinct flow regimes “inside” and “outside” of it. When the surface is not closed, or is not fully resolved by the CFD grid, a slower but more general and accurate collision-based approach can be employed. Both approaches take advantage of a bounding box hierarchy to efficiently store and access the elements of the discretized interface. Moreover, for distributed flow computations, an algorithm for constructing distributed bounding box hierarchies is also presented. The performances of these algorithms are assessed in three-dimensional, dynamic fluid-structure interaction problems in the fields of aeronautics and underwater implosion. The differences between the two interface tracking algorithms in terms of capability, accuracy, and efficiency are revealed and discussed.

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

[2]  Carlos A. Felippa,et al.  Membrane triangles with corner drilling freedoms II: the ANDES element , 1992 .

[3]  Charbel Farhat,et al.  Partitioned procedures for the transient solution of coupled aeroelastic problems , 2001 .

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

[5]  Jacques Periaux,et al.  Numerical simulation and optimal shape for viscous flow by a fictitious domain method , 1995 .

[6]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

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

[8]  P. Colella,et al.  A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains , 1998 .

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

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

[11]  H. S. Udaykumar,et al.  A Sharp Interface Cartesian Grid Methodfor Simulating Flows with ComplexMoving Boundaries , 2001 .

[12]  Joachim Gudmundsson,et al.  Box-Trees and R-Trees with Near-Optimal Query Time , 2001, SCG '01.

[13]  Charbel Farhat,et al.  Aeroelastic Dynamic Analysis of a Full F-16 Configuration for Various Flight Conditions , 2003 .

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

[15]  Carlos A. Felippa,et al.  A unified formulation of small-strain corotational finite elements: I. Theory , 2005 .

[16]  E. Guendelman,et al.  Coupling water and smoke to thin deformable and rigid shells , 2005, SIGGRAPH 2005.

[17]  Ronald Fedkiw,et al.  Robust treatment of collisions, contact and friction for cloth animation , 2002, SIGGRAPH Courses.

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

[19]  HEINZ-OTTO KREISS,et al.  A Second Order Accurate Embedded Boundary Method for the Wave Equation with Dirichlet Data , 2005, SIAM J. Sci. Comput..

[20]  Ralf Deiterding,et al.  A virtual test facility for the efficient simulation of solid material response under strong shock and detonation wave loading , 2006, Engineering with Computers.

[21]  Ralf Deiterding,et al.  Large-scale fluid-structure interaction simulation of viscoplastic and fracturing thin-shells subjected to shocks and detonations , 2007 .

[22]  Rajat Mittal,et al.  A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries , 2008, J. Comput. Phys..

[23]  Ronald Fedkiw,et al.  Two-way coupling of fluids to rigid and deformable solids and shells , 2008, ACM Trans. Graph..

[24]  T. Belytschko,et al.  Robust and provably second‐order explicit–explicit and implicit–explicit staggered time‐integrators for highly non‐linear compressible fluid–structure interaction problems , 2010 .

[25]  Kevin G. Wang,et al.  Algorithms for interface treatment and load computation in embedded boundary methods for fluid and fluid–structure interaction problems , 2011 .

[26]  Ronald Fedkiw,et al.  Numerically stable fluid-structure interactions between compressible flow and solid structures , 2011, J. Comput. Phys..