fbpx
Wikipedia

Meshfree methods

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.

20 points and their Voronoi cells

Motivation edit

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 edit

In a traditional finite difference simulation, the domain of a one-dimensional simulation would be some function  , represented as a mesh of data values   at points  , where

 
 
 
 

We can define the derivatives that occur in the equation being simulated using some finite difference formulae on this domain, for example

 

and

 

Then we can use these definitions of   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   and timestep  ) are constant along all the mesh, and the left and right mesh neighbors of the data value at   are the values at   and  , 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   with them. SPH then defines the value of   between the particles by

 

where   is the mass of particle  ,   is the density of particle  , and   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

 

Then we can use these definitions of   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   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   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 edit

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 edit

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

Related methods:

Recent development edit

The primary areas of advancement in meshfree methods are to address issues with essential boundary enforcement, numerical quadrature, and contact and large deformations.[25] 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,[26] and also yields inaccurate and non-convergent results due to under-integration of the weak form.[27] 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.[27] 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.[25] More recently, a framework has been developed to pass arbitrary-order patch tests, based on a Petrov–Galerkin method.[28]

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.[29][30] 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).[17] 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.[27] 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.[17][35]

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.[citation needed] 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 edit

References edit

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

Further reading edit

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

External links edit

  • The USACM blog on Meshfree Methods

meshfree, methods, field, numerical, analysis, meshfree, methods, those, that, require, connection, between, nodes, simulation, domain, mesh, rather, based, interaction, each, node, with, neighbors, consequence, original, extensive, properties, such, mass, kin. 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 20 points and their Voronoi cells Contents 1 Motivation 2 Example 3 History 4 List of methods and acronyms 5 Recent development 6 See also 7 References 8 Further reading 9 External linksMotivation editNumerical 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 singularitiesExample editIn a traditional finite difference simulation the domain of a one dimensional simulation would be some function u x t displaystyle u x t nbsp represented as a mesh of data values uin displaystyle u i n nbsp at points xi displaystyle x i nbsp where i 0 1 2 displaystyle i 0 1 2 nbsp n 0 1 2 displaystyle n 0 1 2 nbsp xi 1 xi h i displaystyle x i 1 x i h forall i nbsp tn 1 tn k n displaystyle t n 1 t n k forall n nbsp We can define the derivatives that occur in the equation being simulated using some finite difference formulae on this domain for example u x ui 1n ui 1n2h displaystyle partial u over partial x u i 1 n u i 1 n over 2h nbsp and u t uin 1 uink displaystyle partial u over partial t u i n 1 u i n over k nbsp Then we can use these definitions of u x t displaystyle u x t nbsp 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 displaystyle h nbsp and timestep k displaystyle k nbsp are constant along all the mesh and the left and right mesh neighbors of the data value at xi displaystyle x i nbsp are the values at xi 1 displaystyle x i 1 nbsp and xi 1 displaystyle x i 1 nbsp 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 ui displaystyle u i nbsp with them SPH then defines the value of u x t displaystyle u x t nbsp between the particles by u x tn imiuinriW x xi displaystyle u x t n sum i m i frac u i n rho i W x x i nbsp where mi displaystyle m i nbsp is the mass of particle i displaystyle i nbsp ri displaystyle rho i nbsp is the density of particle i displaystyle i nbsp and W displaystyle W nbsp 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 u x imiuinri W x xi x displaystyle partial u over partial x sum i m i frac u i n rho i partial W x x i over partial x nbsp Then we can use these definitions of u x t displaystyle u x t nbsp 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 displaystyle u x t nbsp 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 displaystyle W nbsp 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 editOne 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 editThe 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 17 Meshfree local radial point interpolation method RPIM 17 Local radial basis function collocation Method LRBFCM 18 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 19 Repeated replacement method RRM 2012 20 Radial basis integral equation method 21 Least square collocation meshless method 2001 22 Exponential Basis Functions method EBFs 2010 23 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 Huerta amp Fernandez Mendez 2000 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 Arroyo amp Ortiz 2006 Space Time Meshfree Collocation Method STMCM see Netuzhylov 2008 Netuzhylov amp Zilian 2009 Meshfree Interface Finite Element Method MIFEM 2015 a hybrid finite element meshfree method for numerical simulation of phase transformation and multiphase flow problems 24 Recent development editThe primary areas of advancement in meshfree methods are to address issues with essential boundary enforcement numerical quadrature and contact and large deformations 25 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 26 and also yields inaccurate and non convergent results due to under integration of the weak form 27 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 27 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 25 More recently a framework has been developed to pass arbitrary order patch tests based on a Petrov Galerkin method 28 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 29 30 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 17 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 27 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 17 35 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 citation needed 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 editContinuum mechanics Smoothed finite element method 35 G space 38 Weakened weak form 29 30 Boundary element method Immersed boundary method Stencil code Particle methodReferences edit Gingold R A Monaghan J J 1 December 1977 Smoothed particle hydrodynamics theory and application to non spherical stars Monthly Notices of the Royal Astronomical Society 181 3 375 389 Bibcode 1977MNRAS 181 375G doi 10 1093 mnras 181 3 375 Libersky Larry D Petschek Albert G Carney Theodore C Hipp Jim R Allahdadi Firooz A November 1993 High Strain Lagrangian Hydrodynamics Journal of Computational Physics 109 1 67 75 doi 10 1006 jcph 1993 1199 Swegle J W Hicks D L Attaway S W January 1995 Smoothed Particle Hydrodynamics Stability Analysis Journal of Computational Physics 116 1 123 134 Bibcode 1995JCoPh 116 123S doi 10 1006 jcph 1995 1010 Nayroles B Touzot G Villon P 1992 Generalizing the finite element method Diffuse approximation and diffuse elements Computational Mechanics 10 5 307 318 Bibcode 1992CompM 10 307N doi 10 1007 BF00364252 S2CID 121511161 Belytschko T Lu Y Y Gu L 30 January 1994 Element free Galerkin methods International Journal for Numerical Methods in Engineering 37 2 229 256 Bibcode 1994IJNME 37 229B doi 10 1002 nme 1620370205 Liu Wing Kam Jun Sukky Zhang Yi Fei 30 April 1995 Reproducing kernel particle methods International Journal for Numerical Methods in Fluids 20 8 9 1081 1106 Bibcode 1995IJNMF 20 1081L doi 10 1002 fld 1650200824 Sulsky D Chen Z Schreyer H L September 1994 A particle method for history dependent materials Computer Methods in Applied Mechanics and Engineering 118 1 2 179 196 doi 10 1016 0045 7825 94 90112 0 https www math ucla edu jteran papers SSCTS13 pdf bare URL PDF Liu W K Chen Y Jun S Chen J S Belytschko T Pan C Uras R A Chang C T March 1996 Overview and applications of the reproducing Kernel Particle methods Archives of Computational Methods in Engineering 3 1 3 80 doi 10 1007 BF02736130 S2CID 122241092 Atluri S N Zhu T 24 August 1998 A new Meshless Local Petrov Galerkin MLPG approach in computational mechanics Computational Mechanics 22 2 117 127 Bibcode 1998CompM 22 117A doi 10 1007 s004660050346 S2CID 3688083 Oliveira T Portela A December 2016 Weak form collocation A local meshless method in linear elasticity Engineering Analysis with Boundary Elements 73 144 160 doi 10 1016 j enganabound 2016 09 010 Chen Shang Ying Hsu Kuo Chin Fan Chia Ming 15 March 2021 Improvement of generalized finite difference method for stochastic subsurface flow modeling Journal of Computational Physics 429 110002 Bibcode 2021JCoPh 42910002C doi 10 1016 J JCP 2020 110002 S2CID 228828681 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 doi 10 1007 s00366 023 01897 6 ISSN 1435 5663 W K Liu S Jun Y F Zhang 1995 Reproducing kernel particle methods Int J Numer Methods Eng 20 8 9 1081 1106 Bibcode 1995IJNMF 20 1081L doi 10 1002 fld 1650200824 A Behzadan H M Shodja M Khezri 2011 A unified approach to the mathematical analysis of generalized RKPM gradient RKPM and GMLS Comput Methods Appl Mech Eng 200 5 8 540 576 Bibcode 2011CMAME 200 540B doi 10 1016 j cma 2010 07 017 Gauger Christoph Leinen Peter Yserentant Harry January 2000 The Finite Mass Method SIAM Journal on Numerical Analysis 37 6 1768 1799 doi 10 1137 S0036142999352564 a b c d Liu G R 2nd edn 2009 Mesh Free Methods CRC Press 978 1 4200 8209 9 Sarler B Vertnik R Meshfree Li B Habbal F Ortiz M 17 September 2010 Optimal transportation meshfree approximation schemes for fluid and plastic flows International Journal for Numerical Methods in Engineering 83 12 1541 1579 Bibcode 2010IJNME 83 1541L doi 10 1002 nme 2869 S2CID 18225521 Walker Wade A Langowski Jorg 6 July 2012 The Repeated Replacement Method A Pure Lagrangian Meshfree Method for Computational Fluid Dynamics PLOS ONE 7 7 e39999 Bibcode 2012PLoSO 739999W doi 10 1371 journal pone 0039999 PMC 3391243 PMID 22866175 Ooi E H Popov V May 2012 An efficient implementation of the radial basis integral equation method Engineering Analysis with Boundary Elements 36 5 716 726 doi 10 1016 j enganabound 2011 12 001 S2CID 122004658 Zhang Xiong Liu Xiao Hu Song Kang Zu Lu Ming Wan 30 July 2001 Least squares collocation meshless method International Journal for Numerical Methods in Engineering 51 9 1089 1100 Bibcode 2001IJNME 51 1089Z doi 10 1002 nme 200 S2CID 119952479 Boroomand B Soghrati S Movahedian B 2009 Exponential basis functions in solution of static and time harmonic elastic problems in a meshless style International Journal for Numerical Methods in Engineering 81 8 971 1018 doi 10 1002 nme 2718 S2CID 4943418 Ghoneim A March 2015 A meshfree interface finite element method for modelling isothermal solutal melting and solidification in binary systems Finite Elements in Analysis and Design 95 20 41 doi 10 1016 j finel 2014 10 002 a b Chen Jiun Shyan Hillman Michael Chi Sheng Wei April 2017 Meshfree Methods Progress Made after 20 Years Journal of Engineering Mechanics 143 4 04017001 doi 10 1061 ASCE EM 1943 7889 0001176 Belytschko Ted Guo Yong Kam Liu Wing Ping Xiao Shao 30 July 2000 A unified stability analysis of meshless particle methods International Journal for Numerical Methods in Engineering 48 9 1359 1400 Bibcode 2000IJNME 48 1359B doi 10 1002 1097 0207 20000730 48 9 lt 1359 AID NME829 gt 3 0 CO 2 U a b c Chen Jiun Shyan Wu Cheng Tang Yoon Sangpil You Yang 20 January 2001 A stabilized conforming nodal integration for Galerkin mesh free methods International Journal for Numerical Methods in Engineering 50 2 435 466 Bibcode 2001IJNME 50 435C doi 10 1002 1097 0207 20010120 50 2 lt 435 AID NME32 gt 3 0 CO 2 A Chen Jiun Shyan Hillman Michael Ruter Marcus 3 August 2013 An arbitrary order variationally consistent integration for Galerkin meshfree methods International Journal for Numerical Methods in Engineering 95 5 387 418 Bibcode 2013IJNME 95 387C doi 10 1002 nme 4512 S2CID 124640562 a b Liu G R 2009 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 81 9 1093 1126 doi 10 1002 nme 2719 S2CID 123009384 a b Liu G R 2009 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 81 9 1127 1156 doi 10 1002 nme 2720 S2CID 119378545 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 G R Liu G R Zhang Edge based Smoothed Point Interpolation Methods International Journal of Computational Methods 5 4 621 646 2008 Liu G R Zhang G Y 20 November 2011 A normed G space and weakened weak W2 formulation of a cell based smoothed point interpolation method International Journal of Computational Methods 06 1 147 179 doi 10 1142 S0219876209001796 Liu G R Zhang G Y 14 May 2008 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 74 7 1128 1161 Bibcode 2008IJNME 74 1128L doi 10 1002 nme 2204 S2CID 54088894 a b c Liu G R 2010 Smoothed Finite Element Methods CRC Press ISBN 978 1 4398 2027 8 page needed Liu G R Xu George X 10 December 2008 A gradient smoothing method GSM for fluid dynamics problems International Journal for Numerical Methods in Fluids 58 10 1101 1133 Bibcode 2008IJNMF 58 1101L doi 10 1002 fld 1788 S2CID 53983110 Zhang Jian Liu G R Lam K Y Li Hua Xu G November 2008 A gradient smoothing method GSM based on strong form governing equation for adaptive analysis of solid mechanics problems Finite Elements in Analysis and Design 44 15 889 909 doi 10 1016 j finel 2008 06 006 Liu G R 20 November 2011 On G space theory International Journal of Computational Methods 06 2 257 289 doi 10 1142 S0219876209001863 Further reading editGarg Sahil Pant Mohit 24 May 2018 Meshfree Methods A Comprehensive Review of Applications International Journal of Computational Methods 15 4 1830001 doi 10 1142 S0219876218300015 Liu M B Liu G R Zong Z 20 November 2011 An overview on smoothed particle hydrodynamics International Journal of Computational Methods 05 1 135 188 doi 10 1142 S021987620800142X Liu G R Liu M B 2003 Smoothed Particle Hydrodynamics a meshfree and Particle Method World Scientific ISBN 981 238 456 1 Atluri S N 2004 The Meshless Method MLPG for Domain amp BIE Discretization Tech Science Press ISBN 0 9657001 8 6 Arroyo M Ortiz M 26 March 2006 Localmaximum entropy approximation schemes a seamless bridge between finite elements and meshfree methods International Journal for Numerical Methods in Engineering 65 13 2167 2202 Bibcode 2006IJNME 65 2167A CiteSeerX 10 1 1 68 2696 doi 10 1002 nme 1534 S2CID 15974625 Belytschko T Chen J S 2007 Meshfree and Particle Methods John Wiley and Sons Ltd ISBN 0 470 84800 6 Belytschko T Huerta A Fernandez Mendez S Rabczuk T 2004 Meshless methods Encyclopedia of Computational Mechanics Vol 1 Chapter 10 John Wiley amp Sons ISBN 0 470 84699 2 Liu G R 1st edn 2002 Mesh Free Methods CRC Press ISBN 0 8493 1238 8 Li S Liu W K 2004 Meshfree Particle Methods Berlin Springer Verlag ISBN 3 540 22256 1 Huerta Antonio Fernandez Mendez Sonia 20 August 2000 Enrichment and coupling of the finite element and meshless methods International Journal for Numerical Methods in Engineering 48 11 1615 1636 Bibcode 2000IJNME 48 1615H doi 10 1002 1097 0207 20000820 48 11 lt 1615 AID NME883 gt 3 0 CO 2 S hdl 2117 8264 S2CID 122813651 Netuzhylov H 2008 A Space Time Meshfree Collocation Method for Coupled Problems on Irregularly Shaped Domains Dissertation TU Braunschweig CSE Computational Sciences in Engineering ISBN 978 3 00 026744 4 also as electronic ed Netuzhylov Hennadiy Zilian Andreas 15 October 2009 Space time meshfree collocation method Methodology and application to initial boundary value problems International Journal for Numerical Methods in Engineering 80 3 355 380 Bibcode 2009IJNME 80 355N doi 10 1002 nme 2638 S2CID 122969330 Alhuri Y Naji A Ouazar D Taik A 26 August 2010 RBF Based Meshless Method for Large Scale Shallow Water Simulations Experimental Validation Mathematical Modelling of Natural Phenomena 5 7 4 10 doi 10 1051 mmnp 20105701 Sousa Washington de Oliveira Rodrigo April 2015 Coulomb s Law Discretization Method a New Methodology of Spatial Discretization for the Radial Point Interpolation Method IEEE Antennas and Propagation Magazine 57 2 277 293 Bibcode 2015IAPM 57 277S doi 10 1109 MAP 2015 2414571 Gross B J Trask N Kuberry P Atzberger P J 15 May 2020 Meshfree methods on manifolds for hydrodynamic flows on curved surfaces A Generalized Moving Least Squares GMLS approach Journal of Computational Physics 409 109340 arXiv 1905 10469 Bibcode 2020JCoPh 40909340G doi 10 1016 j jcp 2020 109340 S2CID 166228451 Gross B J Kuberry P Atzberger P J 15 March 2022 First passage time statistics on surfaces of general shape Surface PDE solvers using Generalized Moving Least Squares GMLS Journal of Computational Physics 453 110932 arXiv 2102 02421 Bibcode 2022JCoPh 45310932G doi 10 1016 j jcp 2021 110932 ISSN 0021 9991 S2CID 231802303 External links editThe USACM blog on Meshfree Methods Retrieved from https en wikipedia org w index php title Meshfree methods amp oldid 1210952494, wikipedia, wiki, book, books, library,

article

, read, download, free, free download, mp3, video, mp4, 3gp, jpg, jpeg, gif, png, picture, music, song, movie, book, game, games.