An incompressible smoothed particle hydrodynamics method for the motion of rigid bodies in fluids

A two-dimensional incompressible smoothed particle hydrodynamics scheme is presented for simulation of rigid bodies moving through Newtonian fluids. The scheme relies on combined usage of the rigidity constraints and the viscous penalty method to simulate rigid body motion. Different viscosity ratios and interpolation schemes are tested by simulating a rigid disc descending in quiescent medium. A viscosity ratio of 100 coupled with weighted harmonic averaging scheme has been found to provide satisfactory results. The performance of the resulting scheme is systematically tested for cases with linear motion, rotational motion and their combination. The test cases include sedimentation of a single and a pair of circular discs, sedimentation of an elliptic disc and migration and rotation of a circular disc in linear shear flow. Comparison with previous results at various Reynolds numbers indicates that the proposed method captures the motion of rigid bodies driven by flow or external body forces accurately. An ISPH scheme for the motion of rigid bodies in Newtonian fluids is presented.The scheme relies on combined use of rigidity constraints and viscous penalty.A viscosity ratio of 100 and weighted harmonic averaging was found satisfactory.The proposed scheme is easy to implement and circumvents explicit boundary conditions.The scheme is successfully tested for linear and rotational motion.

[1]  T. Inamuro,et al.  Effect of internal mass in the simulation of a moving body by the immersed boundary method , 2011 .

[2]  C. Pozrikidis.  Axisymmetric motion of a file of red blood cells through capillaries , 2005 .

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

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

[5]  A. Ladd.  Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 2. Numerical results , 1993, Journal of Fluid Mechanics.

[6]  Songdong Shao,et al.  Incompressible SPH simulation of water entry of a free‐falling object , 2009 .

[7]  P. Koumoutsakos.  MULTISCALE FLOW SIMULATIONS USING PARTICLES , 2005 .

[8]  C. W. Hirt,et al.  An Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds , 1997 .

[9]  Olivier Simonin,et al.  A Lagrangian VOF tensorial penalty method for the DNS of resolved particle-laden flows , 2014, J. Comput. Phys..

[10]  Afzal Suleman,et al.  A robust weakly compressible SPH method and its comparison with an incompressible SPH , 2012 .

[11]  Yoichiro Matsumoto,et al.  A full Eulerian finite difference approach for solving fluid-structure coupling problems , 2010, J. Comput. Phys..

[12]  Andreas Acrivos,et al.  Steady simple shear flow past a circular cylinder at moderate Reynolds numbers: a numerical solution , 1974, Journal of Fluid Mechanics.

[13]  Afzal Suleman,et al.  SPH with the multiple boundary tangent method , 2009 .

[14]  Cyrus K. Aidun,et al.  The dynamics and scaling law for particles suspended in shear flow with inertia , 2000, Journal of Fluid Mechanics.

[15]  J. Monaghan,et al.  A refined particle method for astrophysical problems , 1985 .

[16]  M. Ozbulut,et al.  A numerical investigation into the correction algorithms for SPH method in modeling violent free surface fl ows , 2013 .

[17]  Saikiran Rapaka,et al.  Flow patterns in the sedimentation of an elliptical particle , 2009, Journal of Fluid Mechanics.

[18]  Earl H. Dowell,et al.  Modeling of Fluid-Structure Interaction , 2001 .

[19]  Salvatore Marrone,et al.  An accurate SPH modeling of viscous flows around bodies at low and moderate Reynolds numbers , 2013, J. Comput. Phys..

[20]  R. Glowinski,et al.  A fictitious domain approach to the direct numerical simulation of incompressible viscous flow past moving rigid bodies: application to particulate flow , 2001 .

[21]  A. Ladd.  Numerical simulations of particulate suspensions via a discretized Boltzmann equation. Part 1. Theoretical foundation , 1993, Journal of Fluid Mechanics.

[22]  Maurizio Brocchini,et al.  A study of violent sloshing wave impacts using an improved SPH method , 2010 .

[23]  Rui Xu,et al.  Comparisons of weakly compressible and truly incompressible algorithms for the SPH mesh free particle method , 2008, J. Comput. Phys..

[24]  Mehrdad T. Manzari,et al.  A modified SPH method for simulating motion of rigid bodies in Newtonian fluid flows , 2012 .

[25]  Toshiaki Hisada,et al.  Multiphysics simulation of left ventricular filling dynamics using fluid-structure interaction finite element method. , 2004, Biophysical journal.

[26]  R. Glowinski,et al.  A new formulation of the distributed Lagrange multiplier/fictitious domain method for particulate flows , 2000 .

[27]  J. Caltagirone,et al.  Numerical modelling of solid particle motion using a new penalty method , 2005 .

[28]  Guirong Liu,et al.  Smoothed Particle Hydrodynamics (SPH): an Overview and Recent Developments , 2010 .

[29]  S. Cummins,et al.  An SPH Projection Method , 1999 .

[30]  S. Shao,et al.  INCOMPRESSIBLE SPH METHOD FOR SIMULATING NEWTONIAN AND NON-NEWTONIAN FLOWS WITH A FREE SURFACE , 2003 .

[31]  G. Oger,et al.  Two-dimensional SPH simulations of wedge water entries , 2006, J. Comput. Phys..

[32]  Z. Feng,et al.  The immersed boundary-lattice Boltzmann method for solving fluid-particles interaction problems , 2004 .

[33]  Neelesh A. Patankar,et al.  A new mathematical formulation and fast algorithm for fully resolved simulation of self-propulsion , 2009, J. Comput. Phys..

[34]  A. Skillen,et al.  Incompressible smoothed particle hydrodynamics (SPH) with reduced temporal noise and generalised Fickian smoothing applied to body–water slam and efficient wave–body interaction , 2013 .

[35]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

[36]  Grégoire Pianet,et al.  Local penalty methods for flows interacting with moving solids at high Reynolds numbers , 2007 .

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

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

[39]  M. Yildiz,et al.  Improved Incompressible Smoothed Particle Hydrodynamics method for simulating flow around bluff bodies , 2011 .

[40]  Jeffrey F. Morris,et al.  Hydrodynamic interaction of two particles in confined linear shear flow at finite Reynolds number , 2007 .

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

[42]  M. Uhlmann.  An immersed boundary method with direct forcing for the simulation of particulate flows , 2005, 1809.08170.

[43]  Andrea Prosperetti,et al.  A Method for Particle Simulation , 2003 .

[44]  Jean-Paul Caltagirone,et al.  A numerical continuous model for the hydrodynamics of fluid particle systems , 1999 .

[45]  Zhaosheng Yu.  A DLM/FD method for fluid/flexible-body interactions , 2005 .

[46]  Rui Xu,et al.  Accuracy and stability in incompressible SPH (ISPH) based on the projection method and a new approach , 2009, J. Comput. Phys..

[47]  Xin Bian,et al.  A splitting integration scheme for the SPH simulation of concentrated particle suspensions , 2014, Comput. Phys. Commun..

[48]  Minami Yoda,et al.  The circular cylinder in simple shear at moderate Reynolds numbers: An experimental study , 2001 .

[49]  M. Yildiz,et al.  Numerical investigation of Newtonian and non-Newtonian multiphase flows using ISPH method , 2013 .

[50]  A. Colagrossi,et al.  Nonlinear water wave interaction with floating bodies in SPH , 2013 .

[51]  Daniel D. Joseph,et al.  Direct simulation of initial value problems for the motion of solid bodies in a Newtonian fluid. Part 2. Couette and Poiseuille flows , 1994, Journal of Fluid Mechanics.

[52]  A. Popel,et al.  Large deformation of red blood cell ghosts in a simple shear flow. , 1998, Physics of fluids.

[53]  J. Monaghan,et al.  SPH simulation of multi-phase flow , 1995 .

[54]  Chien-Cheng Chang,et al.  A numerical study of the motion of a neutrally buoyant cylinder in two dimensional shear flow , 2012, 1209.0805.

[55]  L. Lucy.  A numerical approach to the testing of the fission hypothesis. , 1977 .

[56]  S. Koshizuka,et al.  International Journal for Numerical Methods in Fluids Numerical Analysis of Breaking Waves Using the Moving Particle Semi-implicit Method , 2022 .

[57]  M. Ozbulut,et al.  DESCENT OF A SOLID DISK IN QUIESCENT FLUID SIMULATED USING INCOMPRESSIBLE SMOOTHED PARTICLE HYDRODYNAMICS , 2014 .

[58]  Xin Liu,et al.  An ISPH simulation of coupled structure interaction with free surface flows , 2014 .

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

[60]  Frederick Stern,et al.  Sharp interface immersed-boundary/level-set method for wave-body interactions , 2009, J. Comput. Phys..

[61]  James J. Feng,et al.  Direct simulation of initial value problems for the motion of solid bodies in a Newtonian fluid Part 1. Sedimentation , 1994, Journal of Fluid Mechanics.

[62]  Joseph P. Morris,et al.  A Study of the Stability Properties of Smooth Particle Hydrodynamics , 1996, Publications of the Astronomical Society of Australia.

[63]  Petros Koumoutsakos,et al.  An immersed boundary method for smoothed particle hydrodynamics of self-propelled swimmers , 2008, J. Comput. Phys..

[64]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.