Smoothedparticle hydrodynamics
Smoothedparticle hydrodynamics (SPH) is a computational method used for simulating the mechanics of continuum media, such as solid mechanics and fluid flows. It was developed by Gingold and Monaghan^{[1]} and Lucy^{[2]} in 1977, initially for astrophysical problems. It has been used in many fields of research, including astrophysics, ballistics, volcanology, and oceanography. It is a meshfree Lagrangian method (where the coordinates move with the fluid), and the resolution of the method can easily be adjusted with respect to variables such as density.
Contents
Method
History
Advantages
 By construction, SPH is a meshfree method, which makes it ideally suited to simulate problems dominated by complex boundary dynamics, like free surface flows, or large boundary displacement.
 On top of that, the lack of a mesh significantly simplifies the model implementation and its parallelization, even for manycore architectures.^{[3]}^{[4]}
 SPH can be easily extended to a wide variety of fields, and hybridized with some other models, as discussed in Advanced techniques and applications.
 As discussed in section on weakly compressible SPH, the method has great conservation features.
Limitations
Examples
Fluid dynamics
Smoothedparticle hydrodynamics is being increasingly used to model fluid motion as well. This is due to several benefits over traditional gridbased techniques. First, SPH guarantees conservation of mass without extra computation since the particles themselves represent mass. Second, SPH computes pressure from weighted contributions of neighboring particles rather than by solving linear systems of equations. Finally, unlike gridbased techniques, which must track fluid boundaries, SPH creates a free surface for twophase interacting fluids directly since the particles represent the denser fluid (usually water) and empty space represents the lighter fluid (usually air). For these reasons, it is possible to simulate fluid motion using SPH in real time. However, both gridbased and SPH techniques still require the generation of renderable free surface geometry using a polygonization technique such as metaballs and marching cubes, point splatting, or 'carpet' visualization. For gas dynamics it is more appropriate to use the kernel function itself to produce a rendering of gas column density (e.g., as done in the SPLASH visualisation package).
One drawback over gridbased techniques is the need for large numbers of particles to produce simulations of equivalent resolution. In the typical implementation of both uniform grids and SPH particle techniques, many voxels or particles will be used to fill water volumes that are never rendered. However, accuracy can be significantly higher with sophisticated gridbased techniques, especially those coupled with particle methods (such as particle level sets), since it is easier to enforce the incompressibility condition in these systems. SPH for fluid simulation is being used increasingly in realtime animation and games where accuracy is not as critical as interactivity.
Recent work in SPH for fluid simulation has increased performance, accuracy, and areas of application:
 B. Solenthaler, 2009, develops PredictiveCorrective SPH (PCISPH) to allow for better incompressibility constraints^{[5]}
 M. Ihmsen et al., 2010, introduce boundary handling and adaptive timestepping for PCISPH for accurate rigid body interactions^{[6]}
 K. Bodin et al., 2011, replace the standard equation of state pressure with a density constraint and apply a variational time integrator^{[7]}
 R. Hoetzlein, 2012, develops efficient GPUbased SPH for large scenes in Fluids v.3^{[8]}
 N. Akinci et al., 2012, introduce a versatile boundary handling and twoway SPHrigid coupling technique that is completely based on hydrodynamic forces; the approach is applicable to different types of SPH solvers ^{[9]}
 M. Macklin et al., 2013 simulates incompressible flows inside the Position Based Dynamics framework, for bigger timesteps ^{[10]}
 N. Akinci et al., 2013, introduce a versatile surface tension and twoway fluidsolid adhesion technique that allows simulating a variety of interesting physical effects that are observed in reality^{[11]}
 J. Kyle and E. Terrell, 2013, apply SPH to FullFilm Lubrication^{[12]}
 A. Mahdavi and N. Talebbeydokhti, 2015, propose a hybrid algorithm for implementation of solid boundary condition and simulate flow over a sharp crested weir^{[13]}
 S. Tavakkol et al., 2016, develop curvSPH, which makes the horizontal and vertical size of particles independent and generates uniform mass distribution along curved boundaries^{[14]}
Astrophysics
Smoothedparticle hydrodynamics's adaptive resolution, numerical conservation of physically conserved quantities, and ability to simulate phenomena covering many orders of magnitude make it ideal for computations in theoretical astrophysics.^{[15]}
Simulations of galaxy formation, star formation, stellar collisions,^{[16]} supernovae^{[17]} and meteor impacts are some of the wide variety of astrophysical and cosmological uses of this method.
SPH is used to model hydrodynamic flows, including possible effects of gravity. Incorporating other astrophysical processes which may be important, such as radiative transfer and magnetic fields is an active area of research in the astronomical community, and has had some limited success.^{[18]}^{[19]}
Solid mechanics
Libersky and Petschek^{[20]}^{[21]} extended SPH to Solid Mechanics. The main advantage of SPH in this application is the possibility of dealing with larger local distortion than gridbased methods. This feature has been exploited in many applications in Solid Mechanics: metal forming, impact, crack growth, fracture, fragmentation, etc.
Another important advantage of meshfree methods in general, and of SPH in particular, is that mesh dependence problems are naturally avoided given the meshfree nature of the method. In particular, mesh alignment is related to problems involving cracks and it is avoided in SPH due to the isotropic support of the kernel functions. However, classical SPH formulations suffer from tensile instabilities^{[22]} and lack of consistency.^{[23]} Over the past years, different corrections have been introduced to improve the accuracy of the SPH solution, leading to the RKPM by Liu et al..^{[24]} Randles and Libersky^{[25]} and Johnson and Beissel^{[26]} tried to solve the consistency problem in their study of impact phenomena.
Dyka et al.^{[27]}^{[28]} and Randles and Libersky^{[29]} introduced the stresspoint integration into SPH and Ted Belytschko et al.^{[30]} showed that the stresspoint technique removes the instability due to spurious singular modes, while tensile instabilities can be avoided by using a Lagrangian kernel. Many other recent studies can be found in the literature devoted to improve the convergence of the SPH method.
Recent improvements in understanding the convergence and stability of SPH have allowed for more widespread applications in Solid Mechanics. Other examples of applications and developments of the method include:
 Metal forming simulations.^{[31]}
 SPHbased method SPAM (Smoothed Particle Applied Mechanics) for impact fracture in solids by William G. Hoover.^{[32]}
 Modified SPH (SPH/MLSPH) for fracture and fragmentation.^{[33]}
 TaylorSPH (TSPH) for shock wave propagation in solids.^{[34]}
Numerical tools
Interpolations
The smoothedparticle hydrodynamics (SPH) method works by dividing the fluid into a set of discrete moving elements , referred to as particles. Their Lagrangian nature allows setting their position by integration of their velocity as:
These particles interact through a kernel function with characteristic radius known as the "smoothing length", typically represented in equations by . This means that the physical quantity of any particle can be obtained by summing the relevant properties of all the particles that lie within the range of the kernel, the latter being used as a weighting function . This can be understood in two steps. First an arbitrary field is written as a convolution with :
The error in making the above approximation is order . Secondly, the integral is approximated using a Riemann summation over the particles:
where the summation over includes all particles in the simulation. is the volume of particle , is the value of the quantity for particle and denotes position. For example, the density of particle can be expressed as:
where denotes the particle mass and the particle density, while is a short notation for . The error done in approximating the integral by a discrete sum depends on , on the particle size (i.e. , being the space dimension), and on the particle arrangement in space. The latter effect is still poorly known.^{[35]}
Kernel functions commonly used include the Gaussian function, the quintic spline and the Wendland kernel.^{[36]} The latter two kernels are compactly suppported (unlike the Gaussian, where there is a small contribution at any finite distance away), with support proportional to . This has the advantage of saving computational effort by not including the relatively minor contributions from distant particles.
Although the size of the smoothing length can be fixed in both space and time, this does not take advantage of the full power of SPH. By assigning each particle its own smoothing length and allowing it to vary with time, the resolution of a simulation can be made to automatically adapt itself depending on local conditions. For example, in a very dense region where many particles are close together, the smoothing length can be made relatively short, yielding high spatial resolution. Conversely, in lowdensity regions where individual particles are far apart and the resolution is low, the smoothing length can be increased, optimising the computation for the regions of interest.
Operators
For particles of constant mass, differentiating the interpolated density with respect to time yields
where is the gradient of with respect to . Comparing the above equation with the continuity equation in continuum mechanics shows that the righthand side is an approximation of ; hence one defines a discrete divergence operator as follows:
This operator gives an SPH approximation of at the particle for a given set of particles with given masses , positions and velocities .
Similarly, one can define a discrete gradient operator to approximate the pressure gradient at the position of particle :
where denote the set of particle pressures. There are several ways to define discrete operators in SPH; the above divergence and gradient formulae have the property to be skewadjoint, leading to nice conservation properties.^{[37]} On the other hand, while the divergence operator is zeroorder consistent, it can be seen that the approximate gradient is not so. Several techniques have been proposed to circumvent this issue, leading to renormalised operators (see e.g.^{[38]}).
Governing equations
The SPH operators can be used to discretize numbers of partial differential equations. For a compressible inviscid fluid, the Euler equations of mass conservation and momentum balance read:
All kinds of SPH divergence and gradient operators can practically be used for discretization purposes. Nevertheless, some perform better regarding physical and numerical effects. A frequently used form of the balance equations is based on the symmetric divergence operator and antisymmetric gradient:
Although there are several ways of discretizing the pressure gradient in the Euler equations, the above antisymmetric form is the most acknowledged one. It supports strict conservation of linear and angular momentum. This means that a force that is exerted on particle by particle equals the one that is exerted on particle by particle including the sign change of the effective direction, thanks to the antisymmetry property .
Variational principle
The above SPH governing equations can be derived from a Least action principle, starting from the Lagrangian of a particle system:
 ,
where is the particle specific internal energy. The EulerLagrange equation of variational mechanics reads, for each particle:
When applied to the above Lagrangian, it gives the following momentum equation:
 ,
where we used the thermodynamical property . Pluging the SPH density interpolation and differentiating explicitly leads to
which is the SPH momentum equation already mentioned, where we recognize the operator. This explains why linear momentum is conserved, and allows conservation of angular momentum and energy to be conserved as well.^{[39]}
Time integration
From the work done in the 80's and 90's on numerical integration of pointlike particles in large accelators, appropriate time integrators have been developed with accurate conservation properties on the long term; they are called symplectic integrators. The most popular in the SPH literature is the leapfrog scheme, which reads for each particle :
where is the time step, superscripts stand for time iterations while is the particle acceleration, given by the righthand side of the momentum equation.
Other symplectic integrators exist (see the reference textbook ^{[40]}). It is recommended to use a symplectic (even loworder) scheme instead of a high order nonsymplectic scheme, to avoid error accumulation after many iterations.
Integration of density has not been studied extensively (see below for more details).
Symplectic schemes are conservative but explicit, thus their numerical stability requires stability conditions, analogous to the CourantFriedrichsLewy condition (see below).
Boundary techniques
In case the SPH convolution shall be practised close to a boundary, i.e. closer than s · h, then the integral support is truncated. Indeed, when the convolution is affected by a boundary, the convolution shall be split in 2 integrals,
where B(r) is the compact support ball centered at r, with radius s · h, and Ω(r) denotes the part of the compact support inside the computational domain, Ω ∩ B(r). Hence, imposing boundary conditions in SPH is completely based on approximating the second integral on the right hand side. The same can be of course applied to the differential operators computation,
Several techniques has been introduced in the past to model boundaries in SPH.
Integral neglect
The most strightforward boundary model is neglecting the integral,
such that just the bulk interactions are taken into account,
This is a popular approach when freesurface is considered in monophase simulations.^{[41]}
The main benefit of this boundary condition is its obvious simplicity. However, several consistency issues shall be considered when this boundary technique is applied.^{[41]} That's in fact a heavy limitation on its potential applications.
Fluid Extension
Probably the most popular methodology, or at least the most traditional one, to impose boundary conditions in SPH, is Fluid Extension technique. Such technique is based on populating the compact support across the boundary with socalled ghost particles, conveniently imposing their field values.^{[42]}
Along this line, the integral neglect methodology can be considered as a particular case of fluid extensions, where the field, A, vanish outside the computational domain.
The main benefit of this methodology is the simplicity, provided that the boundary contribution is computed as part of the bulk interactions. Also, this methodology has been deeply analysed in the literature.^{[43]}^{[42]}^{[44]}
On the other hand, deploying ghost particles in the truncated domain is not a trivial task, such that modelling complex boundary shapes becomes cumbersome. The 2 most popular approaches to populate the empty domain with ghost particles are MirroredParticles ^{[45]} and FixedParticles.^{[42]}
Boundary Integral
The newest Boundary technique is the Boundary Integral methodology.^{[46]} In this methodology, the empty volume integral is replaced by a surface integral, and a renormalization:
with n_{j} the normal of the generic j^{th} boundary element. The surface term can be also solved considering a semianalytic expression.^{[46]}
Modelling Physics
Hydrodynamics
Weakly compressible approach
The weakly compressible SPH in fluid dynamics is based on the discretization of the NavierStokes equations. More precisely, the equations of compressible mass conservation and incompressible momentum balance. The governing equations are shown above. To close the system of differential equations, an appropriate equation of state is utilized.
Another way to determine the density is based on the SPH identity operator itself. Therefore, the density is defined as a state that can be determined from the present particle distribution utilizing the inherent description of the SPH which is introduced with the interpolation. To overcome undesired errors at the free surface through kernel truncation, the density formulation can again be integrated in time. ^{[46]}
Now that the system is described by a weakly compressible scheme, a coupling between pressure and mass density is needed to close the system of differential equations. Therefore, the socalled Cole equation ^{[47]} (sometimes mistakenly referred to as the "Tait equation") is introduced in SPH. It reads
where is the reference mass density and the artificial speed of sound. For water, is considered to give satisfactory values. The background pressure is added to avoid negative pressure values.
Real nearly incompressible fluids such as water are defined by very high speed of sounds of the order . Hence, pressure information travels fast compared to the actual bulk flow, which leads to very small Mach numbers . The momentum equation leads to the following relation:
where is the density change and the velocity vector. The relation shows that the relative change of density is negligible for very high Mach numbers. In theory, these fluids are assumed to be incompressible. If now a density variation of 1% is accepted, a Mach number of 0.1 has to be assumed. This gives a 10 times lower maximum time increment compared to other methods whose CFL condition relates to the maximum bulk velocity instead of the speed of sound. However, other methods often require multiple iterations per time step to assure the divergence free velocity field. Based on the latter discussion the SPH assumes the following values for the speed of sound:
where the maximum velocity vmax needs to be estimated, for e.g. by Torricelli’s law or an educated guess. Since only small density variations occur, a linear equation of state is sometimes sufficient:
Usually the weaklycompressible schemes are affected by a highfrequency spurious noise on the pressure and density fields. ^{[48]} This phenomenon is caused by the nonlinear interaction of acoustic waves and by fact that the scheme is explicit in time and centred in space^{[citation needed]}. Through the years, several techniques have been proposed to get rid of this problem. They can be classified in three different groups: the schemes that adopt density filters, the models that add a diffusive term in the continuity equation, the schemes that employ Riemann Solvers to model the particles' interaction.
These schemes apply a filter directly on the density field to remove the spurious numerical noise. The most used filters are the MLS (Mean Least Squares) and the Shepard filter ^{[48]} which can be applied at each time step or every n time steps (cit). The more frequent is the use of the filtering procedure, the more regular density and pressure fields are obtained. On the other hand, this leads to an increase of the computational costs. In long time simulations, the use of the filtering procedure may lead to an inconsistency between the global volume of fluid and the density field. Further, it does not ensure the enforcement of the dynamic freesurface boundary condition.
A different way to smooth out the density and pressure field is to add a diffusive term inside the continuity equation:
The first schemes that adopted such an approach were described in Ferrari et al. (2009) and in Molteni \& Colagrossi (2009?2010?) where the diffusive term was modelled as a Laplacian of the density field. A similar approach was also used in Fatehi \& Manzari (2010). In Antuono at al.(2010) a correction to the diffusive term of Molteni and Colagrossi (2009) was proposed to remove some inconsistencies close to the freesurface. In this case the adopted diffusive term was equivalent to a highorder differential operator on the density field (Antuono 2012). The scheme was called δSPH and preserved all the conservations properties of the classical SPH (link) (e.g., linear and angular momenta, total energy, see Antuono 2015) along with a smooth and regular representation of the density and pressure fields.
Incompressible approach
Viscosity modelling
In general, the description of hydrodynamic flows require a convenient treatment of diffusive processes to model the viscosity in the NavierStokes equations. It needs special consideration because it involves the laplacian differential operator. Since the direct computation does not provide satisfactory results, several approaches to model the diffusion have been proposed.
 Artificial viscosity
Introduced by Monaghan and Gingold ^{[49]} the artificial viscosity was used to deal with high Mach number fluid flows. It reads
Here, is controlling a volume viscosity while acts similar to the Neumann Richtmeyr artifical viscosity. The is defined by
The artificial viscosity also has shown to improve the overall stability of general flow simulations. Therefore, it is applied to inviscid problems in the following form
It is possible to not only stabilize inviscid simulations but also to model the physical viscosity by this approach. To do so
is substituted in the equation above , where is the number of spartial dimensions of the model. This approach introduces the bulk viscosity .
 Morris
For low Reynolds numbers the viscosity model by Morris ^{[50]} was proposed.
 LoShao
Additional physics
 Surface Tension
 Heat transfer
 Turbulence
Multiphase extensions
Astrophysics
Often in astrophysics, one wishes to model selfgravity in addition to pure hydrodynamics. The particlebased nature of SPH makes it ideal to combine with a particlebased gravity solver, for instance tree gravity code,^{[51]} particle mesh, or particleparticle particlemesh.
Solid mechanics
Others
Variants of the method
This section is empty. You can help by adding to it. (July 2018)

References
 ^ R.A. Gingold; J.J. Monaghan (1977). "Smoothed particle hydrodynamics: theory and application to nonspherical stars". Mon. Not. R. Astron. Soc. 181: 375–89.
 ^ L.B. Lucy (1977). "A numerical approach to the testing of the fission hypothesis". Astron. J. 82: 1013–1024.
 ^ Takahiro Harada; Seiichi Koshizuka; Yoichiro Kawaguchi (2007). Smoothed particle hydrodynamics on GPUs. Computer Graphics International. pp. 63–70.
 ^ Alejandro Crespo; Jose M. Dominguez; Anxo Barreiro; Moncho GomezGesteira; Benedict D. Rogers (2011). "GPUs, a new tool of acceleration in CFD: efficiency and reliability on smoothed particle hydrodynamics methods". PLOS One.
 ^ Solenthaler (2009). "PredictiveCorrective Incompressible SPH".
 ^ Imhsen (2010). "Boundary handling and adaptive timestepping for PCISPH". Workshop on Virtual Reality Interaction and Physical Simulation VRIPHYS.
 ^ Bodin (2011). "Constraint Fluids". IEEE Transactions on Visualization and Computer Graphics.

^ Hoetzlein (2012). "Fluids v.3, A Large scale, Open Source Fluid Simulator. http://fluids3.com". External link in
title=
(help) 
^ Akinci (2012). "Versatile RigidFluid Coupling for Incompressible SPH http://www.nadir.tk/research". ACM TOG, SIGGRAPH proceedings. External link in
title=
(help) 
^ Macklin (2013). "Position Based Fluids http://blog.mmacklin.com/publications". ACM TOG, SIGGRAPH proceedings. External link in
title=
(help) 
^ Akinci (2013). "Versatile Surface Tension and Adhesion for SPH Fluids SPH http://www.nadir.tk/research". ACM TOG, SIGGRAPH proceedings. External link in
title=
(help)  ^ Journal of Tribology (2013). "Application of Smoothed Particle Hydrodynamics to FullFilm Lubrication".

^ Mahdavi and Talebbeydokhti (2015). "A hybrid solid boundary treatment algorithm for smoothed particle hydrodynamics https://www.researchgate.net/publication/282870566_A_hybrid_solid_boundary_treatment_algorithm_for_smoothed_particle_hydrodynamics". Scientia Iranica, Transaction A, Civil Engineering. 22 (4): 1457–1469. External link in
title=
(help)  ^ International Journal for Numerical Methods in Fluids (2016). "Curvilinear smoothed particle hydrodynamics".
 ^ https://arxiv.org/abs/0903.5075
 ^ Rosswog, Stephan (2015). "SPH Methods in the Modelling of Compact Objects". Living Rev Comput Astrophys. 1 (1). arXiv:1406.4224 . Bibcode:2015LRCA....1....1R. doi:10.1007/lrca20151.
 ^ https://arxiv.org/abs/astroph/0512532
 ^ http://www.astro.ex.ac.uk/people/mbate/Cluster/clusterRT.html
 ^ http://users.monash.edu.au/~dprice/pubs/spmhd/pricespmhd.pdf
 ^ Libersky, L.D.; Petschek, A.G. (1990). "Smooth Particle Hydrodynamics with Strength of Materials, Advances in the Free Lagrange Method". Lecture Notes in Physics. 395: 248–257. doi:10.1007/3540549609_58.
 ^ L.D. Libersky; A.G. Petschek; A.G. Carney; T.C. Hipp; J.R. Allahdadi; F.A. High (1993). "Strain Lagrangian hydrodynamics: a threedimensional SPH code for dynamic material response". J. Comput. Phys. 109: 67–75. Bibcode:1993JCoPh.109...67L. doi:10.1006/jcph.1993.1199.
 ^ J.W. Swegle; D.A. Hicks; S.W. Attaway (1995). "Smooth particle hydrodynamics stability analysis". J. Comput. Phys. 116: 123–134. Bibcode:1995JCoPh.116..123S. doi:10.1006/jcph.1995.1010.
 ^ T. Belytschko; Y. Krongauz; J. Dolbow; C. Gerlach (1998). "On the completeness of meshfree particle methods". Int. J. Numer. Methods Eng. 43: 785–819. Bibcode:1998IJNME..43..785B. doi:10.1002/(sici)10970207(19981115)43:5<785::aidnme420>3.0.co;29.
 ^ 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.
 ^ P.W. Randles; L.D. Libersky (1997). "Recent improvements in SPH modelling of hypervelocity impact". Int. J. Impact Eng. 20: 525–532. doi:10.1016/s0734743x(97)874416.
 ^ G.R. Johnson; S.R. Beissel (1996). "Normalized smoothing functions for SPH impact computations". Int. J. Numer. Methods Eng. 39: 2725–2741. Bibcode:1996IJNME..39.2725J. doi:10.1002/(sici)10970207(19960830)39:16<2725::aidnme973>3.0.co;29.
 ^ C.T. Dyka; R.P. Ingel (1995). "An approach for tension instability in Smoothed Particle Hydrodynamics". Comput. Struct. 57: 573–580. doi:10.1016/00457949(95)00059p.
 ^ C.T. Dyka; P.W. Randles; R.P. Ingel (1997). "Stress points for tension instability in SPH". Int. J. Numer. Methods Eng. 40: 2325–2341. Bibcode:1997IJNME..40.2325D. doi:10.1002/(sici)10970207(19970715)40:13<2325::aidnme161>3.0.co;28.
 ^ P.W. Randles; L.D. Libersky (2000). "Normalized SPH with stress points". Int. J. Numer. Methods Eng. 48: 1445–1462. Bibcode:2000IJNME..48.1445R. doi:10.1002/10970207(20000810)48:10<1445::aidnme831>3.0.co;29.
 ^ T. Belytschko; Y. Guo; W.K. Liu; S.P. Xiao (2000). "A unified stability analysis of meshless particle methods". Int. J. Numer. Methods Eng. 48: 1359–1400. Bibcode:2000IJNME..48.1359B. doi:10.1002/10970207(20000730)48:9<1359::aidnme829>3.0.co;2u.
 ^ J. Bonet; S. Kulasegaram (2000). "Correction and stabilization of smooth particle hydrodynamics methods with applications in metal forming simulations". Int. J. Numer. Methods Eng. 47: 1189–1214. Bibcode:2000IJNME..47.1189B. doi:10.1002/(sici)10970207(20000228)47:6<1189::aidnme830>3.0.co;2i.
 ^ W. G. Hoover; C. G. Hoover (2001). "SPAMbased recipes for continuum simulations". Computing in Science and Engineering. 3 (2): 78–85.
 ^ T. Rabczuk; J. Eibl; L. Stempniewski (2003). "Simulation of high velocity concrete fragmentation using SPH/MLSPH". Int. J. Numer. Methods Eng. 56: 1421–1444. Bibcode:2003IJNME..56.1421R. doi:10.1002/nme.617.
 ^ M.I. Herreros; M. Mabssout (2011). "A twosteps time discretization scheme using the SPH method for shock wave propagation". Comput. Methods Appl. Mech. Engrg. 200: 1833–1845. Bibcode:2011CMAME.200.1833H. doi:10.1016/j.cma.2011.02.006.
 ^ N.J. Quinlan; M. Basa; M. Lastiwka (2006). "Truncation error in meshfree particle methods". International Journal for Numerical Methods in Engineering (66).
 ^ H. Wendland (1995). "Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree". Advances in Computational Mathematics (4).
 ^ A. Mayrhofer; B.D. Rogers; D. Violeau; M. Ferrand (2013). "Investigation of wall bounded flows using SPH and the unified semianalytical wall boundary conditions". Computer Physics Communications (184).
 ^ J. Bonet; T.S. Lok (1999). "Variational and momentum preservation aspects of Smoothed Particle Hydrodynamics formulations". Computers Methods in Applied Mechanical Engineering (180).
 ^ J.J. Monaghan (2005). "Smoothed particle hydrodynamics". Reports on Progress in Physics (68).
 ^ E. Hairer; C. Lubich; G. Wanner (2006). Geometric Numerical Integration. Springer. ISBN 9783540306665.
 ^ ^{a} ^{b} Andrea Colagrossi; Matteo Antuono; David Le Touzè (2009). "Theoretical considerations on the freesurface role in the smoothedparticlehydrodynamics model". Physical Review E. 79.
 ^ ^{a} ^{b} ^{c} Bejamin Bouscasse; Andrea Colagrossi; Salvatore Marrone; Matteo Antuono (2013). "Nonlinear water wave interaction with floating bodies in SPH". Journal of Fluids and Structures. 42: 112–129.
 ^ Fabricio Macià; Matteo Antuono; Leo M González; Andrea Colagrossi (2011). "Theoretical analysis of the noslip boundary condition enforcement in SPH methods". Progress of theoretical physics. 125: 1091–1121.
 ^ Jose Luis CercosPita; Matteo Antuono; Andrea Colagrossi; Antonio Souto (2017). "SPH energy conservation for fluidsolid interactions". Computer Methods in Applied Mechanics and Engineering. 317: 771–791.
 ^ J. Campbell; R. Vignjevic; L. Libersky (2000). "A contact algorithm for smoothed particle hydrodynamics". Computer methods in applied mechanics and engineering. 184: 49–65.
 ^ ^{a} ^{b} ^{c} M. Ferrand, D.R. Laurence, B.D. Rogers, D. Violeau, C. Kassiotis (2013). "Unified semianalytical wall boundary conditions for inviscid, laminar or turbulent flows in the meshless SPH method". 71. Int. J. Numer. Meth. Fluids.
 ^ H. R. Cole (1948). Underwater Explosions. Princeton, New Jersey: Princeton University Press.
 ^ ^{a} ^{b} Colagrossi, Andrea; Landrini, Maurizio (2003). "Numerical simulation of interfacial flows by smoothed particle hydrodynamics". Journal of computational physics. 191 (2): 448–475. Bibcode:2003JCoPh.191..448C. doi:10.1016/S00219991(03)003243.
 ^ J.J. Monaghan; R.A. Gingold (1983). "Shock Simulation by the Particle Method". 52. Journal of Computational Physics: 347–389.
 ^ J.P. Morris; P.J. Fox; Y. Zhu (1997). "Modeling Low Reynolds Number Incompressible Flows Using SPH". 136. Journal of Computational Physics: 214–226.
 ^ Marios D. Dikaiakos; Joachim Stadel, PKDGRAV The Parallel kD Tree Gravity Code, retrieved February 1, 2017
Further reading
 Hoover, W. G. (2006). Smooth Particle Applied Mechanics: The State of the Art, World Scientific.
 Impact Modelling with SPH Stellingwerf, R. F., Wingate, C. A., Memorie della Societa Astronomia Italiana, Vol. 65, p. 1117 (1994).
 Amada, T., Imura, M., Yasumuro, Y., Manabe, Y. and Chihara, K. (2004) Particlebased fluid simulation on GPU, in proceedings of ACM Workshop on Generalpurpose Computing on Graphics Processors (August, 2004, Los Angeles, California).
 Desbrun, M. and Cani, MP. (1996). Smoothed Particles: a new paradigm for animating highly deformable bodies. In Proceedings of Eurographics Workshop on Computer Animation and Simulation (August 1996, Poitiers, France).
 Hegeman, K., Carr, N.A. and Miller, G.S.P. Particlebased fluid simulation on the GPU. In Proceedings of International Conference on Computational Science (Reading, UK, May 2006). Proceedings published as Lecture Notes in Computer Science v. 3994/2006 (SpringerVerlag).
 M. Kelager. (2006) Lagrangian Fluid Dynamics Using Smoothed Particle Hydrodynamics, M. Kelagar (MS Thesis, Univ. Copenhagen).
 Kolb, A. and Cuntz, N. (2005). Dynamic particle coupling for GPUbased fluid simulation. In Proceedings of the 18th Symposium on Simulation Techniques (2005) pp. 722–727.
 Liu, G.R. and Liu, M.B. Smoothed Particle Hydrodynamics: a meshfree particle method. Singapore: World Scientific (2003).
 Monaghan, J.J. (1992). Smoothed Particle Hydrodynamics. Annu. Rev. Astron. Astrophys. (1992). 30 : 54374.
 Muller, M., Charypar, D. and Gross, M. ] Particlebased Fluid Simulation for Interactive Applications, In Proceedings of Eurographics/SIGGRAPH Symposium on Computer Animation (2003), eds. D. Breen and M. Lin.
 Vesterlund, M. Simulation and Rendering of a Viscous Fluid Using Smoothed Particle Hydrodynamics, (MS Thesis, Umea University, Sweden).
 Violeau, D., Fluid Mechanics and the SPH method. Oxford University Press (2012).
External links
 First large simulation of star formation using SPH
 SPHERIC (SPH European Research Interest Community)
 ITVO is the website of The Italian Theoretical Virtual Observatory created to query a database of numerical simulation archive.
 SPHC Image Gallery depicts a wide variety of test cases, experimental validations, and commercial applications of the SPH code SPHC.
 A derivation of the SPH model starting from NavierStokes equations
Software
 Algodoo is a 2D simulation framework for education using SPH
 AQUAgpusph is the free (GPLv3) SPH of the researchers, by the researchers, for the researchers
 dive.sph is a commercial webbased SPH engineering software for CFD purposes
 DualSPHysics is an open source SPH code based on SPHysics and using GPU computing
 FLUIDS v.1 is a simple, open source (Zlib), realtime 3D SPH implementation in C++ for liquids for CPU and GPU.
 Fluidix is a GPUbased particle simulation API available from OneZero Software
 GADGET [1] is a freely available (GPL) code for cosmological Nbody/SPH simulations
 GPUSPH SPH simulator with viscosity (GPLv3)
 Pasimodo is a program package for particlebased simulation methods, e.g. SPH
 Physics Abstraction Layer is an open source abstraction system that supports real time physics engines with SPH support
 Punto is a freely available visualisation tool for particle simulations
 pysph Open Source Framework for Smoothed Particle Hydrodynamics in Python (New BSD License)
 SimPARTIX is a commercial simulation package for SPH and DEM simulations from Fraunhofer IWM
 SPHflow
 SPHERA
 SPHysics is an open source SPH implementation in Fortran
 SPLASH is an open source (GPL) visualisation tool for SPH simulations
 SYMPLER: A freeware SYMbolic ParticLE simulatoR from the University of Freiburg.