Quadric SFDI for Laplacian Discretisation in Lagrangian Meshless Methods

Abstract

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 × 3 for 3D cases, compared with 6 × 6 or 10 × 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.

Keywords

Laplacian discretisation / Lagrangian meshless methods / QSFDI / Random/disordered particle distribution / Poisson’s equation / Patch tests

References

[1] Abbaszadeh M, Dehghan M. The interpolating element-free Galerkin method for solving Korteweg-de Vries-Rosenau-regularized long-wave equation with error analysis. Nonlinear Dyn, 2019, 96(2): 1345-1365

[2] Brookshaw L. A method of calculating radiative heat diffusion in particle simulations. Publ Astron Soc Aust, 1985, 6(2): 207-210

[3] Chen JK, Beraun JE, Carney TC. A corrective smoothed particle method for boundary value problems in heat conduction. Int J Numer Methods Eng, 1999, 46(2): 231-252

[4] Cummins SJ, Rudman M. An SPH projection method. J Comput Phys, 1999, 152(2): 584-607

[5] Dehghan M, Abbaszadeh M. Variational multiscale element-free Galerkin method combined with the moving Kriging interpolation for solving some partial differential equations with discontinuous solutions. Comput Appl Math, 2018, 37(3): 3869-3905

[6] Dehghan M, Abbaszadeh M. Error analysis and numerical simulation of magnetohydrodynamics (MHD) equation based on the interpolating element free Galerkin (IEFG) method. Appl Numer Math, 2019, 137: 252-273

[7] Fatehi R, Manzari MT. Error estimation in smoothed particle hydrodynamics and a new scheme for second derivatives. Comput Math Appl, 2011, 61(2): 482-498

[8] Gotoh H, Khayyer A. Current achievements and future perspectives for projection-based particle methods with applications in ocean engineering. J Ocean Eng Marine Energy, 2016, 2(3): 251-278

[9] Gotoh H, Khayyer A, Ikari H, Arikawa T, Shimosako K. On enhancement of incompressible SPH method for simulation of violent sloshing flows. Appl Ocean Res, 2014, 46: 104-115

[10] Hu XY, Adams NA. An incompressible multi-phase SPH method. J Comput Phys, 2007, 227(1): 264-278

[11] Ikari H, Khayyer A, Gotoh H. Corrected higher order Laplacian for enhancement of pressure calculation by projection-based particle methods with applications in ocean engineering. J Ocean Eng Marine Energy, 2015, 1(4): 361-376

[12] Khayyer A, Gotoh H. A higher order Laplacian model for enhancement and stabilization of pressure calculation by the MPS method. Appl Ocean Res, 2010, 32(1): 124-131

[13] Khayyer A, Gotoh H. A 3D higher order Laplacian model for enhancement and stabilization of pressure calculation in 3D MPS-based simulations. Appl Ocean Res, 2012, 37: 120-126

[14] Khayyer A, Gotoh H, Shao SD. Corrected incompressible SPH method for accurate water-surface tracking in breaking waves. Coast Eng, 2008, 55(3): 236-250

[15] Koshizuka S, Oka Y. Moving-particle semi-implicit method for fragmentation of incompressible fluid. Nucl Sci Eng, 1996, 123(3): 421-434

[16] Lee ES, Moulinec C, Xu R, Violeau D, Laurence D, Stansby PK. Comparisons of weakly compressible and truly incompressible algorithms for the SPH mesh free particle method. J Comput Phys, 2008, 227(18): 8417-8436

[17] Lind SJ, Stansby PK. High-order Eulerian incompressible smoothed particle hydrodynamics with transition to Lagrangian free surface motion. J Comput Phys, 2016, 326: 290-311

[18] Lind SJ, Xu R, Stansby PK, Rogers BD. Incompressible smoothed particle hydrodynamics for free-surface flows: a generalised diffusion-based algorithm for stability and validations for impulsive flows and propagating waves. J Comput Phys, 2012, 231(4): 1499-1523

[19] Lo EY, Shao S. Simulation of near-shore solitary wave mechanics by an incompressible SPH method. Appl Ocean Res, 2002, 24(5): 275-286

[20] Ma QW. MLPG method based on Rankine source solution for simulating nonlinear water waves. Comput Model Eng Sci, 2005, 9(2): 193-210

[21] Ma QW. Meshless local Petrov–Galerkin method for two-dimensional nonlinear water wave problems. J Comput Phys, 2005, 205(2): 611-625

[22] Ma QW. A new meshless interpolation scheme for MLPG_R method. Comput Model Eng Sci, 2008, 23(2): 75-90

[23] Ma QW, Zhou Y, Yan S. A review on approaches to solving Poisson’s equation in projection-based meshless methods for modelling strongly nonlinear water waves. J Ocean Eng Marine Energy, 2016, 2(3): 279-299

[24] Monaghan JJ. Simulating free surface flows with SPH. J Comput Phys, 1994, 110(2): 399-406

[25] Oger G, Doring M, Alessandrini B, Ferrant P. An improved SPH method: towards higher order convergence. J Comput Phys, 2007, 225(2): 1472-1492

[26] Quinlan NJ, Basa M, Lastiwka M. Truncation error in mesh-free particle methods. Int J Numer Methods Eng, 2006, 66(13): 2064-2085

[27] Schwaiger HF. An implicit corrected SPH formulation for thermal diffusion with linear free surface boundary conditions. Int J Numer Methods Eng, 2008, 75(6): 647-671

[28] Shao S, Ji C, Graham DI, Reeve DE, James PW, Chadwick AJ. Simulation of wave overtopping by an incompressible SPH model. Coast Eng, 2006, 53(9): 723-735

[29] Shao S, Lo EYM. Incompressible SPH method for simulating Newtonian and non-Newtonian flows with free surface. Adv Water Resour, 2003, 26(7): 787-800

[30] Tamai T, Koshizuka S. Least squares moving particle semi-implicit method. Comput Part Mech, 2014, 1(3): 277-305

[31] Tamai T, Murotani K, Koshizuka S. On the consistency and convergence of particle-based meshfree discretization schemes for the Laplace operator. Comput Fluids, 2017, 142: 79-85

[32] Xu R, Stansby P, Laurence D. Accuracy and stability in incompressible SPH (ISPH) based on the projection method and a new approach. J Comput Phys, 2009, 228(18): 6703-6025

[33] Zheng X, Ma QW, Duan WY. Incompressible SPH method based on Rankine source solution for violent water wave simulation. J Comput Phys, 2014, 276: 291-314

[34] Zheng X, Ma QW, Shao S. Study on SPH Viscosity Term Formulations. Appl Sci, 2018, 8(2): 249

[35] Zhou JT, Ma QW. MLPG method based on Rankine source solution for modelling 3D breaking waves. Comput Model Eng Sci, 2010, 56(2): 179

Funding

City, University of London