An incompressible multi-phase SPH method

An incompressible multi-phase SPH method is proposed. In this method, a fractional time-step method is introduced to enforce both the zero-density-variation condition and the velocity-divergence-free condition at each full time-step. To obtain sharp density and viscosity discontinuities in an incompressible multi-phase flow a new multi-phase projection formulation, in which the discretized gradient and divergence operators do not require a differentiable density or viscosity field is proposed. Numerical examples for Taylor-Green flow, capillary waves, drop deformation in shear flows and for Rayleigh-Taylor instability are presented and compared to theoretical solutions or references from literature. The results suggest good accuracy and convergence properties of the proposed method.

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

[2]  Leonardo Di G. Sigalotti,et al.  SPH simulations of time-dependent Poiseuille flow at low Reynolds numbers , 2003 .

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

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

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

[6]  J. Morris,et al.  Modeling Low Reynolds Number Incompressible Flows Using SPH , 1997 .

[7]  Leonardo Di G. Sigalotti,et al.  On the SPH tensile instability in forming viscous liquid drops , 2004 .

[8]  Nikolaus A. Adams,et al.  A multi-phase SPH method for macroscopic and mesoscopic flows , 2006, J. Comput. Phys..

[9]  J. Monaghan.  Simulating Free Surface Flows with SPH , 1994 .

[10]  C. Pozrikidis,et al.  The flow of suspensions in channels: Single files of drops , 1993 .

[11]  J. Morris.  Simulating surface tension with smoothed particle hydrodynamics , 2000 .

[12]  Jacek Pozorski,et al.  SPH computation of incompressible viscous flows , 2002 .

[13]  M E Cates,et al.  Role of inertia in two-dimensional deformation and breakdown of a droplet. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[15]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[16]  B. Jiang,et al.  SIMULATION OF TWO-FLUID FLOWS BY THE LEAST-SQUARES FINITE ELEMENT METHOD USING A CONTINUUM SURFACE TENSION MODEL , 1998 .

[17]  Geoffrey Ingram Taylor,et al.  The formation of emulsions in definable fields of flow , 1934 .

[18]  Wm. G. Hoover.  Isomorphism linking smooth particles and embedded atoms , 1998 .

[19]  Pep Español,et al.  Incompressible smoothed particle hydrodynamics , 2007, J. Comput. Phys..

[20]  J. Monaghan,et al.  Fundamental differences between SPH and grid methods , 2006, astro-ph/0610051.

[21]  Petros Koumoutsakos,et al.  Remeshed smoothed particle hydrodynamics for the simulation of viscous and heat conducting flows , 2002 .