Error estimation in smoothed particle hydrodynamics and a new scheme for second derivatives

Several schemes for discretization of first and second derivatives are available in Smoothed Particle Hydrodynamics (SPH). Here, four schemes for approximation of the first derivative and three schemes for the second derivative are examined using a theoretical analysis based on Taylor series expansion both for regular and irregular particle distributions. Estimation of terms in the truncation errors shows that only the renormalized (the first-order consistent) scheme has acceptable convergence properties to approximate the first derivative. None of the second derivative schemes has the first-order consistency. Therefore, they converge only when the particle spacing decreases much faster than the smoothing length of the kernel function. In addition, using a modified renormalization tensor, a new SPH scheme is presented for approximating second derivatives that has the property of first-order consistency. To assess the computational performance of the proposed scheme, it is compared with the best available schemes when applied to a 2D heat equation. The numerical results show at least one order of magnitude improvement in accuracy when the new scheme is used. In addition, the new scheme has higher-order convergence rate on regular particle arrangements even for the case of only four particles in the neighborhood of each particle.

[1]  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.

[2]  Jean-Paul Vila,et al.  SPH Renormalized Hybrid Methods for Conservation Laws: Applications to Free Surface Flows , 2005 .

[3]  A. Tielens.  Surface chemistry of deuterated molecules , 1983 .

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

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

[6]  Jean-Paul Vila,et al.  ON PARTICLE WEIGHTED METHODS AND SMOOTH PARTICLE HYDRODYNAMICS , 1999 .

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

[8]  M. Lastiwka,et al.  Truncation error in mesh‐free particle methods , 2006 .

[9]  P. Raviart,et al.  A particle method for first-order symmetric systems , 1987 .

[10]  M. Lastiwka,et al.  Robustness and accuracy of SPH formulations for viscous flow , 2009 .

[11]  Wei Ge,et al.  Is standard symmetric formulation always better for smoothed particle hydrodynamics? , 2008, Comput. Math. Appl..

[12]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

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

[14]  J. S. Halow,et al.  Smoothed particle hydrodynamics: Applications to heat conduction , 2003 .

[15]  B. Ben Moussa,et al.  Convergence of SPH Method for Scalar Nonlinear Conservation Laws , 2000, SIAM J. Numer. Anal..

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

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

[18]  Mehrdad T. Manzari,et al.  On Discretization of Second-order Derivatives in Smoothed Particle Hydrodynamics , 2008 .

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

[20]  Mark A Fleming,et al.  Meshless methods: An overview and recent developments , 1996 .

[21]  J. P. Hughes,et al.  Accuracy of SPH viscous flow models , 2008 .

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

[23]  Bertrand Alessandrini,et al.  An improved SPH method: Towards higher order convergence , 2007, J. Comput. Phys..

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

[25]  L. Brookshaw,et al.  A Method of Calculating Radiative Heat Diffusion in Particle Simulations , 1985, Publications of the Astronomical Society of Australia.

[26]  Michael Griebel,et al.  Meshfree Methods for Partial Differential Equations IV , 2005 .

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