Meshfree methods explained

In the field of numerical analysis, meshfree methods are those that do not require connection between nodes of the simulation domain, i.e. a mesh, but are rather based on interaction of each node with all its neighbors. As a consequence, original extensive properties such as mass or kinetic energy are no longer assigned to mesh elements but rather to the single nodes. Meshfree methods enable the simulation of some otherwise difficult types of problems, at the cost of extra computing time and programming effort. The absence of a mesh allows Lagrangian simulations, in which the nodes can move according to the velocity field.

Motivation

Numerical methods such as the finite difference method, finite-volume method, and finite element method were originally defined on meshes of data points. In such a mesh, each point has a fixed number of predefined neighbors, and this connectivity between neighbors can be used to define mathematical operators like the derivative. These operators are then used to construct the equations to simulate—such as the Euler equations or the Navier–Stokes equations.

But in simulations where the material being simulated can move around (as in computational fluid dynamics) or where large deformations of the material can occur (as in simulations of plastic materials), the connectivity of the mesh can be difficult to maintain without introducing error into the simulation. If the mesh becomes tangled or degenerate during simulation, the operators defined on it may no longer give correct values. The mesh may be recreated during simulation (a process called remeshing), but this can also introduce error, since all the existing data points must be mapped onto a new and different set of data points. Meshfree methods are intended to remedy these problems. Meshfree methods are also useful for:

• Simulations where creating a useful mesh from the geometry of a complex 3D object may be especially difficult or require human assistance
• Simulations where nodes may be created or destroyed, such as in cracking simulations
• Simulations where the problem geometry may move out of alignment with a fixed mesh, such as in bending simulations
• Simulations containing nonlinear material behavior, discontinuities or singularities

Example

In a traditional finite difference simulation, the domain of a one-dimensional simulation would be some function

u(x, t) , represented as a mesh of data values u ^{n} _{i} at points x _{i} , wherei =0 , 1 , 2 ... n =0 , 1 , 2 ... x _{i+1} - x _{i=h \forall} i t _{n+1} - t _{n=k \forall} n We can define the derivatives that occur in the equation being simulated using some finite difference formulae on this domain, for example{\partial u\over \partial x} = {u ^{n\over} _{i -1} 2 h} and{\partial u\over \partial t} = {u ^{n +1} _{i} - u ^{n\over} _{i} k} Then we can use these definitions of u(x,t) and its spatial and temporal derivatives to write the equation being simulated in finite difference form, then simulate the equation with one of many finite difference methods.

In this simple example, the steps (here the spatial step

h and timestep k ) are constant along all the mesh, and the left and right mesh neighbors of the data value at x _{i} are the values at x _{i-1} and x _{i+1} , respectively. Generally in finite differences one can allow very simply for steps variable along the mesh, but all the original nodes should be preserved and they can move independently only by deforming the original elements. If even only two of all the nodes change their order, or even only one node is added to or removed from the simulation, that creates a defect in the original mesh and the simple finite difference approximation can no longer hold.

Smoothed-particle hydrodynamics (SPH), one of the oldest meshfree methods, solves this problem by treating data points as physical particles with mass and density that can move around over time, and carry some value

u _{i} with them. SPH then defines the value of u(x,t) between the particles byu(x,t _{n)} = \sum _{i} m _{i} (u ^{n} _{i})/(\rho _{i}) W(|x - x _{i|)}

where

m _{i} is the mass of particle i , \rho _{i} is the density of particle i , and W is a kernel function that operates on nearby data points and is chosen for smoothness and other useful qualities. By linearity, we can write the spatial derivative as{\partial u\over \partial x} = \sum _{i} m _{i} (u ^{n} _{i})/(\rho _{i}) {\partial W(|x - x _{i|)} \over \partial x}

Then we can use these definitions of

u(x,t) and its spatial derivatives to write the equation being simulated as an ordinary differential equation, and simulate the equation with one of many numerical methods. In physical terms, this means calculating the forces between the particles, then integrating these forces over time to determine their motion.

The advantage of SPH in this situation is that the formulae for

u(x,t) and its derivatives do not depend on any adjacency information about the particles; they can use the particles in any order, so it doesn't matter if the particles move around or even exchange places.

One disadvantage of SPH is that it requires extra programming to determine the nearest neighbors of a particle. Since the kernel function

W only returns nonzero results for nearby particles within twice the "smoothing length" (because we typically choose kernel functions with compact support), it would be a waste of effort to calculate the summations above over every particle in a large simulation. So typically SPH simulators require some extra code to speed up this nearest neighbor calculation.

History

One of the earliest meshfree methods is smoothed particle hydrodynamics, presented in 1977.^[1] Libersky et al.^[2] were the first to apply SPH in solid mechanics. The main drawbacks of SPH are inaccurate results near boundaries and tension instability that was first investigated by Swegle.^[3]

In the 1990s a new class of meshfree methods emerged based on the Galerkin method. This first method called the diffuse element method^[4] (DEM), pioneered by Nayroles et al., utilized the MLS approximation in the Galerkin solution of partial differential equations, with approximate derivatives of the MLS function. Thereafter Belytschko pioneered the Element Free Galerkin (EFG) method,^[5] which employed MLS with Lagrange multipliers to enforce boundary conditions, higher order numerical quadrature in the weak form, and full derivatives of the MLS approximation which gave better accuracy. Around the same time, the reproducing kernel particle method^[6] (RKPM) emerged, the approximation motivated in part to correct the kernel estimate in SPH: to give accuracy near boundaries, in non-uniform discretizations, and higher-order accuracy in general. Notably, in a parallel development, the Material point methods were developed around the same time^[7] which offer similar capabilities. Material point methods are widely used in the movie industry to simulate large deformation solid mechanics, such as snow in the movie Frozen.^[8] RKPM and other meshfree methods were extensively developed by Chen, Liu, and Li in the late 1990s for a variety of applications and various classes of problems.^[9] During the 1990s and thereafter several other varieties were developed including those listed below.

List of methods and acronyms

The following numerical methods are generally considered to fall within the general class of "meshfree" methods. Acronyms are provided in parentheses.

• Smoothed particle hydrodynamics (SPH) (1977)
• Diffuse element method (DEM) (1992)
• Dissipative particle dynamics (DPD) (1992)
• Element-free Galerkin method (EFG / EFGM) (1994)
• Reproducing kernel particle method (RKPM) (1995)
• Finite point method (FPM) (1996)
• Finite pointset method (FPM) (1998)
• hp-clouds
• Natural element method (NEM)
• Material point method (MPM)
• Meshless local Petrov Galerkin (MLPG) (1998)^[10]
• Generalized-strain mesh-free (GSMF) formulation (2016)^[11]
• Moving particle semi-implicit (MPS)
• Generalized finite difference method (GFDM)^{[12] [13]}
• Particle-in-cell (PIC)
• Moving particle finite element method (MPFEM)
• Finite cloud method (FCM)
• Boundary node method (BNM)
• Meshfree moving Kriging interpolation method (MK)
• Boundary cloud method (BCM)
• Method of fundamental solutions (MFS)
• Method of particular solution (MPS)
• Method of finite spheres (MFS)
• Discrete vortex method (DVM)
• Reproducing Kernel Particle Method (RKPM) (1995) ^[14]
• Generalized/Gradient Reproducing Kernel Particle Method (2011) ^[15]
• Finite mass method (FMM) (2000)^[16]
• Smoothed point interpolation method (S-PIM) (2005).
• Meshfree local radial point interpolation method (RPIM).
• Local radial basis function collocation Method (LRBFCM)^[17]
• Viscous vortex domains method (VVD)
• Cracking Particles Method (CPM) (2004)
• Discrete least squares meshless method (DLSM) (2006)
• Immersed Particle Method (IPM) (2006)
• Optimal Transportation Meshfree method (OTM) (2010)^[18]
• Repeated replacement method (RRM) (2012)^[19]
• Radial basis integral equation method^[20]
• Least-square collocation meshless method (2001)^[21]
• Exponential Basis Functions method (EBFs) (2010)^[22]

Related methods:

• Moving least squares (MLS) – provide general approximation method for arbitrary set of nodes
• Partition of unity methods (PoUM) – provide general approximation formulation used in some meshfree methods
• Continuous blending method (enrichment and coupling of finite elements and meshless methods) – see
• eXtended FEM, Generalized FEM (XFEM, GFEM) – variants of FEM (finite element method) combining some meshless aspects
• Smoothed finite element method (S-FEM) (2007)
• Gradient smoothing method (GSM) (2008)
• Local maximum-entropy (LME) – see
• Space-Time Meshfree Collocation Method (STMCM) – see,
• Meshfree Interface-Finite Element Method (MIFEM) (2015) - a hybrid finite element-meshfree method for numerical simulation of phase transformation and multiphase flow problems^[23]

Recent development

The primary areas of advancement in meshfree methods are to address issues with essential boundary enforcement, numerical quadrature, and contact and large deformations.^[24] The common weak form requires strong enforcement of the essential boundary conditions, yet meshfree methods in general lack the Kronecker delta property. This make essential boundary condition enforcement non-trivial, at least more difficult than the Finite element method, where they can be imposed directly. Techniques have been developed to overcome this difficulty and impose conditions strongly. Several methods have been developed to impose the essential boundary conditions weakly, including Lagrange multipliers, Nitche's method, and the penalty method.

As for quadrature, nodal integration is generally preferred which offers simplicity, efficiency, and keeps the meshfree method free of any mesh (as opposed to using Gauss quadrature, which necessitates a mesh to generate quadrature points and weights). Nodal integration however, suffers from numerical instability due to underestimation of strain energy associated with short-wavelength modes,^[25] and also yields inaccurate and non-convergent results due to under-integration of the weak form.^[26] One major advance in numerical integration has been the development of a stabilized conforming nodal integration (SCNI) which provides a nodal integration method which does not suffer from either of these problems.^[26] The method is based on strain-smoothing which satisfies the first order patch test. However, it was later realized that low-energy modes were still present in SCNI, and additional stabilization methods have been developed. This method has been applied to a variety of problems including thin and thick plates, poromechanics, convection-dominated problems, among others.^[24] More recently, a framework has been developed to pass arbitrary-order patch tests, based on a Petrov–Galerkin method.^[27]

One recent advance in meshfree methods aims at the development of computational tools for automation in modeling and simulations. This is enabled by the so-called weakened weak (W2) formulation based on the G space theory.^{[28] [29]} The W2 formulation offers possibilities to formulate various (uniformly) "soft" models that work well with triangular meshes. Because a triangular mesh can be generated automatically, it becomes much easier in re-meshing and hence enables automation in modeling and simulation. In addition, W2 models can be made soft enough (in uniform fashion) to produce upper bound solutions (for force-driving problems). Together with stiff models (such as the fully compatible FEM models), one can conveniently bound the solution from both sides. This allows easy error estimation for generally complicated problems, as long as a triangular mesh can be generated. Typical W2 models are the Smoothed Point Interpolation Methods (or S-PIM).^[30] The S-PIM can be node-based (known as NS-PIM or LC-PIM),^[31] edge-based (ES-PIM),^[32] and cell-based (CS-PIM).^[33] The NS-PIM was developed using the so-called SCNI technique.^[26] It was then discovered that NS-PIM is capable of producing upper bound solution and volumetric locking free.^[34] The ES-PIM is found superior in accuracy, and CS-PIM behaves in between the NS-PIM and ES-PIM. Moreover, W2 formulations allow the use of polynomial and radial basis functions in the creation of shape functions (it accommodates the discontinuous displacement functions, as long as it is in G1 space), which opens further rooms for future developments. The W2 formulation has also led to the development of combination of meshfree techniques with the well-developed FEM techniques, and one can now use triangular mesh with excellent accuracy and desired softness. A typical such a formulation is the so-called smoothed finite element method (or S-FEM).^[35] The S-FEM is the linear version of S-PIM, but with most of the properties of the S-PIM and much simpler.

It is a general perception that meshfree methods are much more expensive than the FEM counterparts. The recent study has found however, some meshfree methods such as the S-PIM and S-FEM can be much faster than the FEM counterparts.

The S-PIM and S-FEM works well for solid mechanics problems. For CFD problems, the formulation can be simpler, via strong formulation. A Gradient Smoothing Methods (GSM) has also been developed recently for CFD problems, implementing the gradient smoothing idea in strong form.^{[36] [37]} The GSM is similar to [FVM], but uses gradient smoothing operations exclusively in nested fashions, and is a general numerical method for PDEs.

Nodal integration has been proposed as a technique to use finite elements to emulate a meshfree behaviour. However, the obstacle that must be overcome in using nodally integrated elements is that the quantities at nodal points are not continuous, and the nodes are shared among multiple elements.

See also

• Continuum mechanics
• Smoothed finite element method
• G space^[38]
• Weakened weak form^{[28] [29]}
• Boundary element method
• Immersed boundary method
• Stencil code
• Particle method

Further reading

• Garg . Sahil . Pant . Mohit . Meshfree Methods: A Comprehensive Review of Applications . International Journal of Computational Methods . 24 May 2018 . 15 . 4 . 1830001 . 10.1142/S0219876218300015 .
• Liu . M. B. . Liu . G. R. . Zong . Z. . An overview on smoothed particle hydrodynamics . International Journal of Computational Methods . 20 November 2011 . 05 . 1 . 135–188 . 10.1142/S021987620800142X .
• Book: Liu . G.R. . Liu . M.B. . 2003 . Smoothed Particle Hydrodynamics, a meshfree and Particle Method . World Scientific . 981-238-456-1 .
• Book: S.N. . Atluri. The Meshless Method (MLPG) for Domain & BIE Discretization . 2004 . Tech Science Press . 0-9657001-8-6 .
• Arroyo . M. . Ortiz . M. . Localmaximum-entropy approximation schemes: a seamless bridge between finite elements and meshfree methods . International Journal for Numerical Methods in Engineering . 26 March 2006 . 65 . 13 . 2167–2202 . 10.1002/nme.1534 . 2006IJNME..65.2167A . 10.1.1.68.2696 . 15974625 .
• Belytschko, T., Chen, J.S. (2007). Meshfree and Particle Methods, John Wiley and Sons Ltd.
• .
• Liu, G.R. 1st edn, 2002. Mesh Free Methods, CRC Press. .
• Li, S., Liu, W.K. (2004). Meshfree Particle Methods, Berlin: Springer Verlag.
• Huerta . Antonio . Fernández-Méndez . Sonia . Enrichment and coupling of the finite element and meshless methods . International Journal for Numerical Methods in Engineering . 20 August 2000 . 48 . 11 . 1615–1636 . 10.1002/1097-0207(20000820)48:11<1615::AID-NME883>3.0.CO;2-S . 2000IJNME..48.1615H . 2117/8264 . 122813651 . free .
• , also as electronic ed..
• Netuzhylov . Hennadiy . Zilian . Andreas . Space-time meshfree collocation method: Methodology and application to initial-boundary value problems . International Journal for Numerical Methods in Engineering . 15 October 2009 . 80 . 3 . 355–380 . 10.1002/nme.2638 . 2009IJNME..80..355N . 122969330 .
• Alhuri . Y. . Naji . A. . Ouazar . D. . Taik . A. . RBF Based Meshless Method for Large Scale Shallow Water Simulations: Experimental Validation . Mathematical Modelling of Natural Phenomena . 26 August 2010 . 5 . 7 . 4–10 . 10.1051/mmnp/20105701 . free .
• Sousa . Washington . de Oliveira . Rodrigo . Coulomb's Law Discretization Method: a New Methodology of Spatial Discretization for the Radial Point Interpolation Method . IEEE Antennas and Propagation Magazine . April 2015 . 57 . 2 . 277–293 . 10.1109/MAP.2015.2414571 . 2015IAPM...57..277S .
• Gross . B. J. . Trask . N. . Kuberry . P. . Atzberger . P. J. . Meshfree methods on manifolds for hydrodynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach . Journal of Computational Physics . 15 May 2020 . 409 . 109340 . 10.1016/j.jcp.2020.109340. 1905.10469 . 2020JCoPh.40909340G . 166228451 .
• Gross . B. J. . Kuberry . P. . Atzberger . P. J. . First-passage time statistics on surfaces of general shape: Surface PDE solvers using Generalized Moving Least Squares (GMLS) . Journal of Computational Physics . 15 March 2022 . 453 . 110932 . 10.1016/j.jcp.2021.110932 . 2102.02421 . 2022JCoPh.45310932G . 231802303 . en . 0021-9991.

External links

• The USACM blog on Meshfree Methods

Notes and References

1. Gingold . R. A. . Monaghan . J. J. . Smoothed particle hydrodynamics: theory and application to non-spherical stars . Monthly Notices of the Royal Astronomical Society . 1 December 1977 . 181 . 3 . 375–389 . 10.1093/mnras/181.3.375 . 1977MNRAS.181..375G . free .
2. Libersky . Larry D. . Petschek . Albert G. . Carney . Theodore C. . Hipp . Jim R. . Allahdadi . Firooz A. . High Strain Lagrangian Hydrodynamics . Journal of Computational Physics . November 1993 . 109 . 1 . 67–75 . 10.1006/jcph.1993.1199 .
3. Swegle . J.W. . Hicks . D.L. . Attaway . S.W. . Smoothed Particle Hydrodynamics Stability Analysis . Journal of Computational Physics . January 1995 . 116 . 1 . 123–134 . 10.1006/jcph.1995.1010 . 1995JCoPh.116..123S .
4. Nayroles . B. . Touzot . G. . Villon . P. . Generalizing the finite element method: Diffuse approximation and diffuse elements . Computational Mechanics . 1992 . 10 . 5 . 307–318 . 10.1007/BF00364252 . 1992CompM..10..307N . 121511161 .
5. Belytschko . T. . Lu . Y. Y. . Gu . L. . Element-free Galerkin methods . International Journal for Numerical Methods in Engineering . 30 January 1994 . 37 . 2 . 229–256 . 10.1002/nme.1620370205 . 1994IJNME..37..229B .
6. Liu . Wing Kam . Jun . Sukky . Zhang . Yi Fei . Reproducing kernel particle methods . International Journal for Numerical Methods in Fluids . 30 April 1995 . 20 . 8–9 . 1081–1106 . 10.1002/fld.1650200824 . 1995IJNMF..20.1081L .
7. Sulsky . D. . Chen . Z. . Schreyer . H.L. . A particle method for history-dependent materials . Computer Methods in Applied Mechanics and Engineering . September 1994 . 118 . 1–2 . 179–196 . 10.1016/0045-7825(94)90112-0 .
8. https://www.math.ucla.edu/~jteran/papers/SSCTS13.pdf
9. Liu . W. K. . Chen . Y. . Jun . S. . Chen . J. S. . Belytschko . T. . Pan . C. . Uras . R. A. . Chang . C. T. . Overview and applications of the reproducing Kernel Particle methods . Archives of Computational Methods in Engineering . March 1996 . 3 . 1 . 3–80 . 10.1007/BF02736130 . 122241092 .
10. Atluri . S. N. . Zhu . T. . A new Meshless Local Petrov-Galerkin (MLPG) approach in computational mechanics . Computational Mechanics . 24 August 1998 . 22 . 2 . 117–127 . 10.1007/s004660050346 . 1998CompM..22..117A . 3688083 .
11. Oliveira . T. . Portela . A. . Weak-form collocation – A local meshless method in linear elasticity . Engineering Analysis with Boundary Elements . December 2016 . 73 . 144–160 . 10.1016/j.enganabound.2016.09.010 .
12. Chen . Shang-Ying . Hsu . Kuo-Chin . Fan . Chia-Ming . Improvement of generalized finite difference method for stochastic subsurface flow modeling . Journal of Computational Physics . 15 March 2021 . 429 . 110002 . 10.1016/J.JCP.2020.110002 . 2021JCoPh.42910002C . 228828681 .
13. Chen . Shang-Ying . Wei . Jian-Yu . Hsu . Kuo-Chin . 2023-10-01 . Data assimilation for real-time subsurface flow modeling with dynamically adaptive meshless node adjustments . Engineering with Computers . en . 10.1007/s00366-023-01897-6 . 1435-5663.
14. W.K. Liu . S. Jun . Y.F. Zhang. Reproducing kernel particle methods. Int. J. Numer. Methods Eng.. 1995. 20. 8–9. 1081–1106. 10.1002/fld.1650200824. 1995IJNMF..20.1081L .
15. A. Behzadan . H. M. Shodja . M. Khezri. A unified approach to the mathematical analysis of generalized RKPM, gradient RKPM, and GMLS. Comput. Methods. Appl. Mech. Eng.. 2011. 200. 5–8. 540–576. 10.1016/j.cma.2010.07.017. 2011CMAME.200..540B .
16. Gauger . Christoph . Leinen . Peter . Yserentant . Harry . The Finite Mass Method . SIAM Journal on Numerical Analysis . January 2000 . 37 . 6 . 1768–1799 . 10.1137/S0036142999352564 .
17. Sarler B, Vertnik R. Meshfree
18. Li . B. . Habbal . F. . Ortiz . M. . Optimal transportation meshfree approximation schemes for fluid and plastic flows . International Journal for Numerical Methods in Engineering . 17 September 2010 . 83 . 12 . 1541–1579 . 10.1002/nme.2869 . 2010IJNME..83.1541L . 18225521 .
19. Walker . Wade A. . Langowski . Jörg . The Repeated Replacement Method: A Pure Lagrangian Meshfree Method for Computational Fluid Dynamics . PLOS ONE . 6 July 2012 . 7 . 7 . e39999 . 10.1371/journal.pone.0039999 . 22866175 . 3391243 . 2012PLoSO...739999W . free .
20. Ooi . E.H. . Popov . V. . An efficient implementation of the radial basis integral equation method . Engineering Analysis with Boundary Elements . May 2012 . 36 . 5 . 716–726 . 10.1016/j.enganabound.2011.12.001 . 122004658 .
21. Zhang . Xiong . Liu . Xiao‐Hu . Song . Kang‐Zu . Lu . Ming‐Wan . Least‐squares collocation meshless method . International Journal for Numerical Methods in Engineering . 30 July 2001 . 51 . 9 . 1089–1100 . 10.1002/nme.200 . 2001IJNME..51.1089Z . 119952479 .
22. Boroomand . B. . Soghrati . S. . Movahedian . B. . Exponential basis functions in solution of static and time harmonic elastic problems in a meshless style . International Journal for Numerical Methods in Engineering . 2009 . 81 . 8 . 971–1018 . 10.1002/nme.2718 . 4943418 .
23. Ghoneim . A. . A meshfree interface-finite element method for modelling isothermal solutal melting and solidification in binary systems . Finite Elements in Analysis and Design . March 2015 . 95 . 20–41 . 10.1016/j.finel.2014.10.002 .
24. Chen . Jiun-Shyan . Jiun-Shyan Chen . Hillman . Michael . Chi . Sheng-Wei . Meshfree Methods: Progress Made after 20 Years . Journal of Engineering Mechanics . April 2017 . 143 . 4 . 04017001 . 10.1061/(ASCE)EM.1943-7889.0001176 .
25. Belytschko . Ted . Guo . Yong . Kam Liu . Wing . Ping Xiao . Shao . A unified stability analysis of meshless particle methods . International Journal for Numerical Methods in Engineering . 30 July 2000 . 48 . 9 . 1359–1400 . 10.1002/1097-0207(20000730)48:9<1359::AID-NME829>3.0.CO;2-U . 2000IJNME..48.1359B .
26. Chen . Jiun-Shyan . Jiun-Shyan Chen . Wu . Cheng-Tang . Yoon . Sangpil . You . Yang . A stabilized conforming nodal integration for Galerkin mesh-free methods . International Journal for Numerical Methods in Engineering . 20 January 2001 . 50 . 2 . 435–466 . 10.1002/1097-0207(20010120)50:2<435::AID-NME32>3.0.CO;2-A . 2001IJNME..50..435C .
27. Chen . Jiun-Shyan . Jiun-Shyan Chen . Hillman . Michael . Rüter . Marcus . An arbitrary order variationally consistent integration for Galerkin meshfree methods . International Journal for Numerical Methods in Engineering . 3 August 2013 . 95 . 5 . 387–418 . 10.1002/nme.4512 . 2013IJNME..95..387C . 124640562 .
28. Liu . G. R. . A G space theory and a weakened weak (W2) form for a unified formulation of compatible and incompatible methods: Part I theory . International Journal for Numerical Methods in Engineering . 2009 . 81 . 9 . 1093–1126 . 10.1002/nme.2719 . 123009384 .
29. Liu . G. R. . A G space theory and a weakened weak (W2) form for a unified formulation of compatible and incompatible methods: Part II applications to solid mechanics problems . International Journal for Numerical Methods in Engineering . 2009 . 81 . 9 . 1127–1156 . 10.1002/nme.2720 . 119378545 .
30. Liu, G.R. 2nd edn: 2009 Mesh Free Methods, CRC Press. 978-1-4200-8209-9
31. Liu GR, Zhang GY, Dai KY, Wang YY, Zhong ZH, Li GY and Han X, A linearly conforming point interpolation method (LC-PIM) for 2D solid mechanics problems, International Journal of Computational Methods, 2(4): 645–665, 2005.
32. G.R. Liu, G.R. Zhang. Edge-based Smoothed Point Interpolation Methods. International Journal of Computational Methods, 5(4): 621–646, 2008
33. Liu . G. R. . Zhang . G. Y. . A normed G space and weakened weak (W2) formulation of a cell-based smoothed point interpolation method . International Journal of Computational Methods . 20 November 2011 . 06 . 1 . 147–179 . 10.1142/S0219876209001796 .
34. Liu . G. R. . Zhang . G. Y. . Upper bound solution to elasticity problems: A unique property of the linearly conforming point interpolation method (LC-PIM) . International Journal for Numerical Methods in Engineering . 14 May 2008 . 74 . 7 . 1128–1161 . 10.1002/nme.2204 . 2008IJNME..74.1128L . 54088894 .
35. Liu, G.R., 2010 Smoothed Finite Element Methods, CRC Press, .
36. Liu . G. R. . Xu . George X. . A gradient smoothing method (GSM) for fluid dynamics problems . International Journal for Numerical Methods in Fluids . 10 December 2008 . 58 . 10 . 1101–1133 . 10.1002/fld.1788 . 2008IJNMF..58.1101L . 53983110 .
37. Zhang . Jian . Liu . G.R. . Lam . K.Y. . Li . Hua . Xu . G. . A gradient smoothing method (GSM) based on strong form governing equation for adaptive analysis of solid mechanics problems . Finite Elements in Analysis and Design . November 2008 . 44 . 15 . 889–909 . 10.1016/j.finel.2008.06.006 .
38. Liu . G. R. . On G space theory . International Journal of Computational Methods . 20 November 2011 . 06 . 2 . 257–289 . 10.1142/S0219876209001863 .