SPH simulations of time-dependent Poiseuille flow at low Reynolds numbers

A working Smoothed Particle Hydrodynamics (SPH) formalism for solving the equations of motion of a viscous fluid is presented. The method is based on a standard symmetrized SPH expression for the viscous forces, which involves only first-order derivatives of the kernel through a direct evaluation of the viscous stress tensor. Therefore, the interpolation can be performed using low-order kernels of compact support without compromising the accuracy and stability of the results. In principle, the scheme is suitable for treating compressible fluids with arbitrary shear and bulk viscosities. Here, we demonstrate that when it is combined with the pressure-gradient correction proposed by Morris et al., the method is also suitable for solving the Navier-Stokes equations for incompressible flows without any further assumptions. Simulations using the method show close agreement with the analytic series solutions for plane Poiseuille and Hagen-Poiseuille flows at very low Reynolds numbers. At least for these specific tests, the results obtained are essentially independent of employing either a cubic or a quintic spline kernel.

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

[2]  L. Libersky,et al.  High strain Lagrangian hydrodynamics: a three-dimensional SPH code for dynamic material response , 1993 .

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

[4]  Anthony Peter Whitworth,et al.  A new prescription for viscosity in smoothed particle hydrodynamics. , 1996 .

[5]  H. Ruder,et al.  Smoothed Particle Hydrodynamics: Physical Viscosity and the Simulation of Accretion Disks , 1994 .

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

[7]  C. G. Hoover,et al.  Molecular dynamics, smoothed-particle applied mechanics, and irreversibility , 1994 .

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

[9]  J. Bonet,et al.  Variational and momentum preservation aspects of Smooth Particle Hydrodynamic formulations , 1999 .

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

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

[12]  H. Posch,et al.  Liquid drops and surface tension with smoothed particle applied mechanics , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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

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

[15]  J. Monaghan,et al.  SPH elastic dynamics , 2001 .

[16]  H. Posch,et al.  Viscous conducting flows with smooth-particle applied mechanics. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[17]  Larry D. Libersky,et al.  Cylindrical smoothed particle hydrodynamics , 1993 .