Simulation of forced deformable bodies interacting with two-dimensional incompressible flows: Application to fish-like swimming

We present an efficient algorithm for simulation of deformable bodies interacting with two-dimensional incompressible flows. The temporal and spatial discretizations of the Navier-Stokes equations in vorticity stream-function formulation are based on classical fourth-order Runge-Kutta and compact finite differences, respectively. Using a uniform Cartesian grid we benefit from the advantage of a new fourth-order direct solver for the Poisson equation to ensure the incompressibility constraint down to machine zero. For introducing a deformable body in fluid flow, the volume penalization method is used. A Lagrangian structured grid with prescribed motion covers the deformable body interacting with the surrounding fluid due to the hydrodynamic forces and moment calculated on the Eulerian reference grid. An efficient law for curvature control of an anguilliform fish, swimming to a prescribed goal, is proposed. Validation of the developed method shows the efficiency and expected accuracy of the algorithm for fish-like swimming and also for a variety of fluid/solid interaction problems.

[1]  P. Koumoutsakos,et al.  Optimal shapes for anguilliform swimmers at intermediate Reynolds numbers , 2013, Journal of Fluid Mechanics.

[2]  K. Namkoong,et al.  Numerical analysis of two-dimensional motion of a freely falling circular cylinder in an infinite fluid , 2008, Journal of Fluid Mechanics.

[3]  J. C. Vassilicos,et al.  A numerical strategy to combine high-order schemes, complex geometry and parallel computing for high resolution DNS of fractal generated turbulence , 2010 .

[4]  Khodor Khadra,et al.  Fictitious domain approach for numerical modelling of Navier–Stokes equations , 2000 .

[5]  Petros Koumoutsakos,et al.  Simulations of single and multiple swimmers with non-divergence free deforming geometries , 2011, J. Comput. Phys..

[6]  High-order finite-difference implementation of the immersed-boundary technique for incompressible flows☆ , 2011 .

[7]  Jae Wook Kim.  Optimised boundary compact finite difference schemes for computational aeroacoustics , 2007, J. Comput. Phys..

[8]  T. Bohr,et al.  Vortex wakes of a flapping foil , 2009, Journal of Fluid Mechanics.

[9]  Petros Koumoutsakos,et al.  Reducing the Time Complexity of the Derandomized Evolution Strategy with Covariance Matrix Adaptation (CMA-ES) , 2003, Evolutionary Computation.

[10]  Stéphane Abide,et al.  A 2D compact fourth-order projection decomposition method , 2005 .

[11]  J. C. Simo,et al.  A finite strain beam formulation. The three-dimensional dynamic problem. Part I , 1985 .

[12]  H. Brinkman.  A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles , 1949 .

[13]  P. Koumoutsakos,et al.  Simulations of optimized anguilliform swimming , 2006, Journal of Experimental Biology.

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

[15]  Nicolas Marchand,et al.  Multi-variable constrained control approach for a three-dimensional eel-like robot , 2008, 2008 IEEE/RSJ International Conference on Intelligent Robots and Systems.

[16]  Petros Koumoutsakos,et al.  C-start: optimal start of larval fish , 2012, Journal of Fluid Mechanics.

[17]  A. Roshko,et al.  Vortex formation in the wake of an oscillating cylinder , 1988 .

[18]  Kai Schneider,et al.  Simulation of confined magnetohydrodynamic flows using a pseudo-spectral method with volume penalization , 2012 .

[19]  Richard Pasquetti,et al.  A pseudo-penalization method for high Reynolds number unsteady flows , 2008 .

[20]  Chang Shu,et al.  Simulation of fish swimming and manoeuvring by an SVD-GFD method on a hybrid meshfree-Cartesian grid , 2010 .

[21]  C. Breder.  The locomotion of fishes , 1926 .

[22]  Kai Schneider,et al.  Simulation of confined magnetohydrodynamic flows with Dirichlet boundary conditions using a pseudo-spectral method with volume penalization , 2014, J. Comput. Phys..

[23]  Angelo Iollo,et al.  Modeling and simulation of fish-like swimming , 2010, J. Comput. Phys..

[24]  Ayman Belkhiri,et al.  Modélisation dynamique de la locomotion compliante : Application au vol battant bio-inspiré de l'insecte , 2013 .

[25]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

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

[27]  G. Carey,et al.  High‐order compact scheme for the steady stream‐function vorticity equations , 1995 .

[28]  Dmitry Kolomenskiy,et al.  A Fourier spectral method for the Navier-Stokes equations with volume penalization for moving solid obstacles , 2009, J. Comput. Phys..

[29]  Jung Hee Seo,et al.  A high-order immersed boundary method for acoustic wave scattering and low-Mach number flow-induced sound in complex geometries , 2011, J. Comput. Phys..

[30]  G. Taylor.  Stability of a Viscous Liquid Contained between Two Rotating Cylinders , 1923 .

[31]  R. Hirsh,et al.  Higher order accurate difference solutions of fluid mechanics problems by a compact differencing technique , 1975 .

[32]  Paolo Orlandi,et al.  Fluid Flow Phenomena: A Numerical Toolkit , 1999 .

[33]  Frédéric Boyer,et al.  Macro-continuous computed torque algorithm for a three-dimensional eel-like robot , 2006, IEEE Transactions on Robotics.

[34]  Diego Rossinelli,et al.  GPU accelerated simulations of bluff body flows using vortex particle methods , 2010, J. Comput. Phys..

[35]  Philippe Angot,et al.  A penalization method to take into account obstacles in incompressible viscous flows , 1999, Numerische Mathematik.

[36]  Patrick Bontoux,et al.  Optimisation of Hermitian methods for Navier-Stokes equations in the vorticity and stream-function formulation , 1980 .

[37]  D. Pritchard.  Sinking inside the box , 2013, Journal of Fluid Mechanics.

[38]  Babak Hejazialhosseini,et al.  Reinforcement Learning and Wavelet Adapted Vortex Methods for Simulations of Self-propelled Swimmers , 2014, SIAM J. Sci. Comput..

[39]  Frederick G. Shuman,et al.  NUMERICAL METHODS IN WEATHER PREDICTION: II. SMOOTHING AND FILTERING* , 1957 .

[40]  Gilles Carbou,et al.  Boundary layer for a penalization method for viscous incompressible flow , 2003, Advances in Differential Equations.

[41]  Kai Schneider,et al.  An efficient algorithm for simulation of forced deformable bodies interacting with incompressible flows; Application to fish swimming , 2014 .

[42]  Mathieu Coquerelle,et al.  ARTICLE IN PRESS Available online at www.sciencedirect.com Journal of Computational Physics xxx (2008) xxx–xxx , 2022 .

[43]  Bendiks Jan Boersma,et al.  A 6th order staggered compact finite difference method for the incompressible Navier-Stokes and scalar transport equations , 2011, J. Comput. Phys..

[44]  William H. Press,et al.  Numerical recipes , 1990 .

[45]  Williams,et al.  Self-propelled anguilliform swimming: simultaneous solution of the two-dimensional navier-stokes equations and Newton's laws of motion , 1998, The Journal of experimental biology.

[46]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[47]  H. Fasel,et al.  A high-order immersed interface method for simulating unsteady incompressible flows on irregular domains , 2005 .

[48]  Kai Schneider,et al.  Two-dimensional simulation of the fluttering instability using a pseudospectral method with volume penalization , 2013 .

[49]  F. Sotiropoulos,et al.  Immersed boundary methods for simulating fluid-structure interaction , 2014 .

[50]  Kai Schneider,et al.  Approximation of the Laplace and Stokes operators with Dirichlet boundary conditions through volume penalization: a spectral viewpoint , 2012, Numerische Mathematik.