# Hydrodynamique des particules lissées

L'hydrodynamique des particules lissées, en anglais Smoothed particle hydrodynamics (SPH), est une méthode de calcul utilisée pour simuler la mécanique des milieux continus, comme la mécanique des solides ou les écoulements de fluides. Elle a été développée par Gingold, Monaghan et Lucy en 1977, initialement pour des problèmes d'astrophysique. Elle a été utilisée dans de nombreux domaines de recherche, incluant l'astrophysique, la balistique, la volcanologie et océanologie. Il s'agit d'une méthode lagrangienne (où les coordonnées se déplacent avec le fluide) sans maillage, et la résolution peut facilement être ajustée en fonction de variables physiques telles que la densité.

## Méthodologie
### Avantages

* Par construction, la SPH est une méthode sans maillage, ce qui la rend adaptée aux simulation de problèmes caractérisés par une dynamique de frontière complexe, comme les écoulements de surface libre ou les grands déplacements de frontière.
* L'absence de maillage simplifie considérablement l'implémentation du modèle et sa parallélisation, même pour les architectures hautement multicœurs$^,$.
* La SPH peut être facilement étendue à une grande variété de domaines et hybridée avec d'autres modèles, comme indiqué dans Modélisations physiques.
* Comme discuté dans la section sur la SPH faiblement compressible, la méthode est particulièrement conservative.
* Les simulations SPH ont un coût par nombre de particules nettement inférieur à celui des simulations avec maillage par nombre de maille lorsque la métrique d'intérêt est liée à la densité du fluide (comme, la fonction de densité de probabilité des fluctuations de densité). Ceci s'explique car la SPH concentre l'effort de résolution là où se situe le besoin.

### Limitations

* La définition des conditions aux limites telles que les entrées et sorties et les parois est moins aisée dans la SPH qu'avec les méthodes s'appuyant sur un maillage. En fait, « le traitement des conditions aux limites est certainement l'un des points techniques les plus difficiles de la méthode SPH ». Ce défi est en partie dû au fait que dans la SPH, les particules proches de la frontière changent avec le temps. Néanmoins, les conditions aux limites des murs pour la SPH sont disponibles$^,$$^,$.
* Le coût de calcul des simulations SPH par nombre de particules est considérablement plus élevé que le coût des simulations basées sur la grille par nombre de cellules lorsque la métrique d'intérêt n'est pas directement liée à la densité (mais, par exemple, au spectre d'énergie cinétique). Par conséquent, en négligeant les problèmes d'accélération par calcul parallèle, la simulation des flux à densité constante (par exemple, l'aérodynamique externe) est plus efficace avec les méthodes à base de maillage.

## Exemples
### Dynamique des fluides

L'hydrodynamique des particules lissées est de plus en plus utilisée pour modéliser la dynamique des fluides. Cela est dû à plusieurs avantages de la méthode par rapport aux méthodes classiques à base de maillage. Premièrement, la SPH garantit la conservation de la masse sans calcul supplémentaire puisque la masse est incarnée par les particules. Deuxièmement, la SPH calcule la pression à partir des contributions pondérées des particules voisines plutôt qu'en résolvant des systèmes linéaires d'équations. Enfin, contrairement aux techniques à base de maillage lequel doit suivre les limites du fluide, la SPH crée spontanément une surface libre pour les fluides en interaction biphasique puisque les particules représentent le fluide dense (généralement l'eau) et l'espace vide représente le fluide léger (généralement l'air). Pour ces raisons, il est possible de suivre le mouvement du fluide en temps réel avec la SPH . Cependant, les méthodes à base de maillage et la SPH nécessitent toujours la génération d'une géométrie de surface libre pouvant être rendue à l'aide d'une technique de polygonisation telle que les métaballes ou les marching cubes, le splatting ou un plan d'interpolation. Pour la dynamique des gaz, il est plus approprié d'utiliser la fonction noyau elle-même pour produire un rendu de la densité de la colonne de gaz (par exemple, comme cela est fait dans le package de visualisation SPLASH).
Un inconvénient par rapport aux méthodes à base de maillage est le grand nombre de particules nécessaire pour atteindre une même résolution. Dans la mise en œuvre typique des méthodes à base de maillage et de particules SPH, de nombreux voxels ou particules sont utilisés pour remplir des volumes d'eau qui ne sont pas restitués. Mais la précision peut être considérablement plus élevée avec des techniques sophistiquées à base de maillage, en particulier celles recourant à des méthodes particulaires (telles que les ensembles de niveaux de particules), car il est plus facile d'appliquer la condition d'incompressibilité dans ces systèmes. La SPH est de plus en plus utilisée pour la simulation de fluides dans l'animation et les jeux en temps réel où l'interactivité prime sur la précision.
Des travaux récents en SPH appliquée à la simulation de fluides ont permis d'améliorer les performances, la précision et d'élargir les domaines d'application :
* B. Solenthaler, 2009, develops Predictive-Corrective SPH (PCISPH) to allow for better incompressibility constraints
* M. Ihmsen et al., 2010, introduce boundary handling and adaptive time-stepping for PCISPH for accurate rigid body interactions
* K. Bodin et al., 2011, replace the standard equation of state pressure with a density constraint and apply a variational time integrator
* R. Hoetzlein, 2012, develops efficient GPU-based SPH for large scenes in Fluids v.3
* N. Akinci et al., 2012, introduce a versatile boundary handling and two-way SPH-rigid coupling technique that is completely based on hydrodynamic forces; the approach is applicable to different types of SPH solvers
* M. Macklin et al., 2013 simulates incompressible flows inside the Position Based Dynamics framework, for bigger timesteps
* N. Akinci et al., 2013, introduce a versatile surface tension and two-way fluid-solid adhesion technique that allows simulating a variety of interesting physical effects that are observed in reality
* J. Kyle and E. Terrell, 2013, apply SPH to Full-Film Lubrication
* A. Mahdavi and N. Talebbeydokhti, 2015, propose a hybrid algorithm for implementation of solid boundary condition and simulate flow over a sharp crested weir
* 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
* W. Kostorz and A. Esmail-Yakas, 2020, propose a general, efficient and simple method for piecewise-planar SPH wall boundary conditions based on the semianalytical formulation of Kulasegaram et al.

### Astrophysique

L'hydrodynamique des particules lissées est particulièrement adaptée aux calculs d'astrophysique théorique du fait d'une résolution adaptative, de la conservation numérique des grandeurs physiques constantes et de sa capacité à simuler des phénomènes couvrant différents ordres de grandeur.
Les simulations de formation de galaxies, de formation d'étoiles, de collisions stellaires, de supernovae et d'impacts d'astéroïdes font partie des divers champs d'application de la SPH en cosmologie.
La SPH est utilisée pour modéliser les écoulements hydrodynamiques en prenant en compte les effets de la gravité. L'intégration d'autres phénomènes physiques qui peuvent jouer un rôle important, tels que le transfert radiatif et les champs magnétiques, est un domaine de recherche actif chez les astrophysiciens mais au succès mitigé$^,$.

### Mécanique des solides déformables

Libersky et Petschek$^,$ ont étendu la SPH à la mécanique des solides déformables. Le principal avantage de la SPH dans cette application est la possibilité de traiter une distorsion locale plus importante qu'avec les méthodes à base de maillage. Cette fonctionnalité a été exploitée dans de nombreuses applications en mécanique : formage de métaux, impacts, propagation de fissures, rupture, fragmentation, etc.
Un autre avantage important des méthodes sans maillage, dont la SPH, est qu'elles s’affranchissent naturellement des problèmes de dépendance au maillage. En particulier, l'alignement de maillage est associé aux problèmes de fissures et est évité en SPH du fait du support isotrope des fonctions du noyau. Cependant, les formulations SPH classiques souffrent d'instabilités de traction et d'un manque de cohérence. Au cours des dernières années, différentes corrections ont été introduites pour améliorer la précision de la solution SPH, conduisant au RKPM par Liu et al.. Randles et Libersky ainsi que Johnson et Beissel ont tenté de résoudre le problème de cohérence dans leur étude des phénomènes d'impact.
Dyka et al.$^,$ et Randles et Libersky ont introduit l'intégration des points de contrainte dans la SPH et Ted Belytschko et al. ont montré que la technique du point de contrainte supprime l'instabilité due à des modes singuliers parasites, tandis que les instabilités de traction peuvent être évitées en utilisant un noyau lagrangien. De nombreuses autres études récentes se trouvent dans la littérature consacrée à l'amélioration de la convergence de la méthode SPH.