Computational algorithms for tracking dynamic fluid–structure interfaces in embedded boundary methods

SUMMARY A robust, accurate, and computationally efficient interface tracking algorithm is a key component of an embedded computational framework for the solution of fluid–structure interaction problems with complex and deformable geometries. To a large extent, the design of such an algorithm has focused on the case of a closed embedded interface and a Cartesian computational fluid dynamics grid. Here, two robust and efficient interface tracking computational algorithms capable of operating on structured as well as unstructured three-dimensional computational fluid dynamics grids are presented. The first one is based on a projection approach, whereas the second one is based on a collision approach. The first algorithm is faster. However, it is restricted to closed interfaces and resolved enclosed volumes. The second algorithm is therefore slower. However, it can handle open shell surfaces and underresolved enclosed volumes. Both computational algorithms exploit the bounding box hierarchy technique and its parallel distributed implementation to efficiently store and retrieve the elements of the discretized embedded interface. They are illustrated, and their respective performances are assessed and contrasted, with the solution of three-dimensional, nonlinear, dynamic fluid–structure interaction problems pertaining to aeroelastic and underwater implosion applications. Copyright © 2012 John Wiley & Sons, Ltd.

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

[2]  R. Löhner,et al.  Applied Computational Fluid Dynamics Techniques: An Introduction Based on Finite Element Methods , 2001 .

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

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

[5]  Charbel Farhat,et al.  A higher-order generalized ghost fluid method for the poor for the three-dimensional two-phase flow computation of underwater implosions , 2008, J. Comput. Phys..

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

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

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

[9]  Gregory W. Brown,et al.  Application of a three-field nonlinear fluid–structure formulation to the prediction of the aeroelastic parameters of an F-16 fighter , 2003 .

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

[11]  Rainald Lhner.  Applied Computational Fluid Dynamics Techniques , 2008 .

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

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

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

[15]  Charbel Farhat,et al.  MIXED EXPLICITlIMPLICIT TIME INTEGRATION OF COUPLED AEROELASTIC PROBLEMS: AND DISTRIBUTED SOLUTION THREE-FIELD FORMULATION, GEOMETRIC CONSERVATION , 1995 .

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

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

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

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

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

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

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

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

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

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

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

[27]  C. Farhat,et al.  Partitioned procedures for the transient solution of coupled aroelastic problems Part I: Model problem, theory and two-dimensional application , 1995 .

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