On Discretization of Second-order Derivatives in Smoothed Particle Hydrodynamics

Authors: R. Fatehi, M.A. Fayazbakhsh, M.T. Manzari

Abstract:

Discretization of spatial derivatives is an important issue in meshfree methods especially when the derivative terms contain non-linear coefficients. In this paper, various methods used for discretization of second-order spatial derivatives are investigated in the context of Smoothed Particle Hydrodynamics. Three popular forms (i.e. "double summation", "second-order kernel derivation", and "difference scheme") are studied using one-dimensional unsteady heat conduction equation. To assess these schemes, transient response to a step function initial condition is considered. Due to parabolic nature of the heat equation, one can expect smooth and monotone solutions. It is shown, however in this paper, that regardless of the type of kernel function used and the size of smoothing radius, the double summation discretization form leads to non-physical oscillations which persist in the solution. Also, results show that when a second-order kernel derivative is used, a high-order kernel function shall be employed in such a way that the distance of inflection point from origin in the kernel function be less than the nearest particle distance. Otherwise, solutions may exhibit oscillations near discontinuities unlike the "difference scheme" which unconditionally produces monotone results.

Keywords: Heat conduction, Meshfree methods, Smoothed ParticleHydrodynamics (SPH), Second-order derivatives.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1083799

References:

[1] L. Lucy, "A numerical approach to the testing of fission hypothesis," Astrophysical Journal, vol. 82, pp. 1013-1020, 1977.
[2] R. Gingold and J. Monaghan, "Smoothed particle hydrodynamics: Theory and application to nonspherical stars," Astrophysical Journal, vol. 181, pp. 275-389, 1977.
[3] ÔÇöÔÇö, "Kernel estimates as a basis for general particle methods in hydrodynamics," Journal of Computational Physics, vol. 46, pp. 429- 453, 1982.
[4] J. Monaghan, "Simulating free surface flows with sph," Journal of Computational Physics, vol. 110, pp. 399-406, 1994.
[5] H. Takeda, S. Miyama, and M. Sekiya, "Numerical simulation of viscous flow by smoothed particle hydrodynamics," Progress of Theoretical Physics, vol. 92, no. 5, pp. 939-960, 1994.
[6] J. Morris, P. Fox, and Y. Zhu, "Modeling low reynolds number incompressible flows using sph," Journal of Computational Physics, vol. 136, no. 1, pp. 214-226, 1997.
[7] J. Monaghan, "Smoothed particle hydrodynamics," Reports on Progress in Physics, vol. 68, pp. 1703-1759, 2005.
[8] O. Flebbe, S. Munzel, H. Herold, H. Riffert, and H. Ruder, "Smoothed particle hydrodynamics-physical viscosity and the simulation of accretion disks," The Astrophysical Journal, vol. 431, pp. 214-226, 1994.
[9] S. Watkins, A. Bhattal, N. Francis, T. J.A., and Whitworth, "A new prescription for viscosity in smoothed particle hydrodynamics," Astronomy and Astrophysics, vol. 119, pp. 177-187, 1996.
[10] J. Jeong, M. Jhon, J. Halow, and J. van Osdol, "Smoothed particle hydrodynamics: Applications to heat conduction," Computer Physics Communications, vol. 153, no. 1, pp. 71-84, 2003.
[11] J. Bonet and T. Lok, "Variational and momentum preservation aspects of Smooth Particle Hydrodynamic formulations," Computer Methods in Applied Mechanics and Engineering, vol. 180, no. 1-2, pp. 97-115, 1999.
[12] S. Nugent and H. Posch, "Liquid drops and surface tension with smoothed particle applied mechanics," Physical Review E, vol. 62, no. 4, pp. 4968-4975, 2000.
[13] J. Bonet, S. Kulasegaram, M. Rodriguez-Paz, and M. Profit, "Variational formulation for the smooth particle hydrodynamics (SPH) simulation of fluid and solid problems," Computer Methods in Applied Mechanics and Engineering, vol. 193, no. 12-14, pp. 1245-1256, 2004.
[14] Y. Mele'an, L. Sigalotti, and A. Hasmy, "On the SPH tensile instability in forming viscous liquid drops," Computer Physics Communications, vol. 157, no. 3, pp. 191-200, 2004.
[15] H. L'opez and L. Sigalotti, "Oscillation of viscous drops with smoothed particle hydrodynamics," Physical Review E, vol. 73, no. 5, p. 51201, 2006.
[16] A. Chaniotis, D. Poulikakos, and P. Kououtsakos, "Remeshed smoothed particle hydrodynamics for the simulations of viscous and heat conducting flows," Journal of Computational Physics, vol. 182, pp. 67-90, 2002.
[17] P. Cleary, "Modelling confined multi-material heat and mass flows using sph," Applied Mathematical Modelling, vol. 22, pp. 981-993, 1998.
[18] J. Monaghan and J. Lattanzio, "A refined particle method for astrophysical problems," Astronomy and Astrophysics, vol. 149, no. 1, pp. 135-143, 1985.
[19] I. Schoenberg, "Contributions to the problem of approximation of equidistant data by analytic functions. Part A - On the problem of smoothing or graduation. A first class of analytic approximation formulas," Quarterly of Applied Mathematics, vol. 4, pp. 45-99, 1946.