Quadric SFDI for Laplacian Discretisation in Lagrangian Meshless Methods

In the Lagrangian meshless (particle) methods, such as the smoothed particle hydrodynamics (SPH), moving particle semi-implicit (MPS) method and meshless local Petrov-Galerkin method based on Rankine source solution (MLPG_R), the Laplacian discretisation is often required in order to solve the governing equations and/or estimate physical quantities (such as the viscous stresses). In some meshless applications, the Laplacians are also needed as stabilisation operators to enhance the pressure calculation. The particles in the Lagrangian methods move following the material velocity, yielding a disordered (random) particle distribution even though they may be distributed uniformly in the initial state. Different schemes have been developed for a direct estimation of second derivatives using finite difference, kernel integrations and weighted/moving least square method. Some of the schemes suffer from a poor convergent rate. Some have a better convergent rate but require inversions of high order matrices, yielding high computational costs. This paper presents a quadric semi-analytical finite-difference interpolation (QSFDI) scheme, which can achieve the same degree of the convergent rate as the best schemes available to date but requires the inversion of significant lower-order matrices, i.e. 3 x 3 for 3D cases, compared with 6 x 6 or 10 x 10 in the schemes with the best convergent rate. Systematic patch tests have been carried out for either estimating the Laplacian of given functions or solving Poisson's equations. The convergence, accuracy and robustness of the present schemes are compared with the existing schemes. It will show that the present scheme requires considerably less computational time to achieve the same accuracy as the best schemes available in literatures, particularly for estimating the Laplacian of given functions.