Les méthodes classiques de résolution des équations de l'élasticité linéaire sont les méthodes éléments finis. Ces méthodes produisent de très bons résultats et sont très largement analysées mathématiquement pour l'étude des déformations solides. Pour des problèmes de couplage solide/fluide, pour des situations réalistes en présence de discontinuités (modélisation des fronts de gel dans les sols humides), ou bien encore pour des domaines de calcul mieux adaptés aux maillages non conformes, il parait intéressant de disposer de solveurs Volumes Finis. Les méthodes Volumes Finis sont très largement utilisées en mécanique des fluides. Appliquées aux problèmes de convection, elles sont bien adaptées à la capture de solutions présentant des discontinuités et ne nécessitent pas de maillages conformes. De plus, elles présentent l'avantage de conserver au niveau discret les flux à travers les interfaces du maillage. C'est pourquoi sont développées et testées, dans cette thèse, plusieurs méthodes de volumes finis, qui permettent de traiter le problème de l'élasticité. On a, dans un premier temps, mis en œuvre la méthode LSGR (Least Squares Gradient Reconstruction), qui reconstruit des gradients par volumes à partir d'une formule de moindres carrés pondérés sur les volumes voisins. Elle est testée pour des maillages tétraédriques non structurés, et montre un ordre 1 de convergence. La méthode des Volumes Finis mixtes est ensuite présentée, basée sur la conservation d'un flux "pénalisé" à travers les interfaces. Cette pénalisation impose une contrainte sur le type de maillage utilisé, et des tests sont réalisés en 2d avec des maillages structurés et non structurés de quadrangles. On étend ensuite la méthode des Volumes Finis diamants à l'élasticité. Cette méthode détermine un gradient discret sur des sous volumes associés aux interfaces à partir de l'interpolation de la solution aux sommets du maillage. La convergence théorique est prouvée sous réserve de vérifier une condition de coercivité. Les résultats numériques, en 2d pour des maillages non structurés, conduisent à un ordre de convergence meilleur que celui prouvé. Enfin, la méthode DDFV (Discrete Duality Finite Volume), qui est une extension de la méthode Diamant, est présentée. Elle est basée sur une correspondance entre plusieurs maillages afin d'y construire des opérateurs discrets en "dualité discrète". On montre que la méthode est convergente d'ordre 1. Les illustrations numériques, réalisées en 2d et en 3d pour des maillages non structurés, montrent une convergence d'ordre 2, ce qui est fréquemment observé pour cette méthode.
| Identifer | oai:union.ndltd.org:CCSD/oai:tel.archives-ouvertes.fr:tel-00798769 |
| Date | 19 September 2012 |
| Creators | Martin, Benjamin |
| Publisher | École normale supérieure de Cachan - ENS Cachan |
| Source Sets | CCSD theses-EN-ligne, France |
| Language | French |
| Detected Language | French |
| Type | PhD thesis |
Page generated in 0.002 seconds