Incompressible smoothed particle hydrodynamics

We present a smoothed particle hydrodynamic model for incompressible fluids. As opposed to solving a pressure Poisson equation in order to get a divergence-free velocity field, here incompressibility is achieved by requiring as a kinematic constraint that the volume of the fluid particles is constant. We use Lagrangian multipliers to enforce this restriction. These Lagrange multipliers play the role of non-thermodynamic pressures whose actual values are fixed through the kinematic restriction. We use the SHAKE methodology familiar in constrained molecular dynamics as an efficient method for finding the non-thermodynamic pressure satisfying the constraints. The model is tested for several flow configurations.

[1]  W. Press,et al.  Numerical Recipes - Example Book (FORTRAN) - Second Edition , 1992 .

[2]  William H. Press,et al.  Numerical Recipes in Fortran 77 , 1992 .

[3]  Gregory I. Sivashinsky,et al.  Weak turbulence in periodic flows , 1985 .

[4]  Harald A. Posch,et al.  Simulation of two-dimensional Kolmogorov flow with smooth particle applied mechanics , 1997 .

[5]  L. Libersky,et al.  Smoothed Particle Hydrodynamics: Some recent improvements and applications , 1996 .

[6]  Pep Español,et al.  Smoothed dissipative particle dynamics. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[7]  Fabrice Colin,et al.  Computing a null divergence velocity field using smoothed particle hydrodynamics , 2006, J. Comput. Phys..

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

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

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

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

[12]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[13]  Roland W. Lewis,et al.  A variational formulation based contact algorithm for rigid boundaries in two-dimensional SPH applications , 2004 .

[14]  J. Monaghan.  SPH without a Tensile Instability , 2000 .

[15]  Nikolaus A. Adams,et al.  Angular-momentum conservative smoothed particle dynamics for incompressible viscous flows , 2006 .

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

[17]  S. Miyama,et al.  Numerical Simulation of Viscous Flow by Smoothed Particle Hydrodynamics , 1994 .

[18]  J. W. Humberston.  Classical mechanics , 1980, Nature.

[19]  R. Tanner,et al.  SPH simulations of transient viscoelastic flows at low Reynolds number , 2005 .

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

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

[22]  Jean-Pierre Hansen,et al.  Brownian dynamics with constraints , 2001 .

[23]  Dominique Laurence,et al.  Incompressible separated flows simulations with the smoothed particle hydrodynamics gridless method , 2005 .