Delaunay-admissiblité en dimensions 2 et 3
La méthode des éléments finis, largement utilisée en analyse numérique, requiert que le domaine considéré soit préalablement maillé, c'est-à-dire partitionné en un ensemble de polytopes généralement, mais pas nécessairement, simpliciaux. Parmi les méthodes permettant la génération de tels maill...
| Main Author: | |
|---|---|
| Language: | FRE |
| Published: |
Université Pierre et Marie Curie - Paris VI
2000
|
| Subjects: | |
| Online Access: | http://tel.archives-ouvertes.fr/tel-00607168 http://tel.archives-ouvertes.fr/docs/00/60/71/68/PDF/TU-0638.pdf |
| Summary: | La méthode des éléments finis, largement utilisée en analyse numérique, requiert que le domaine considéré soit préalablement maillé, c'est-à-dire partitionné en un ensemble de polytopes généralement, mais pas nécessairement, simpliciaux. Parmi les méthodes permettant la génération de tels maillages, la triangulation de Delaunay présente le double intérêt d'avoir un support théorique fondant la robustesse des algorithmes, ainsi que de produire des éléments de qualité, conditionnant fortement la précision des calculs ultérieurs. Elle présente cependant l'inconvénient de ne pas être à même de prendre en compte des considérations topologiques, lui interdisant de facto d'être utilisée en l'état pour produire des maillages. Un certain nombre de méthodes ont été proposées pour tenter de résoudre ce problème, mais aucune ne constitue une solution générale. Par ailleurs, les maillages qu'elles restituent ne possèdent plus la propriété de Delaunay. Ce travail étudie les conditions dans lesquelles une contrainte, en dimensions 2 et 3, apparaîtra dans toute triangulation de Delaunay du nuage de points auquel elle est associée. En particulier, des théorèmes de Delaunay-admissibilité a priori sont établis. A l'aide de ces résultats, des algorithmes de redéfinition de contraintes sont proposés, de telle sorte que les nouvelles discrétisations, recouvrements des anciennes, soient construites par toute triangulation de Delaunay. Ainsi, les contraintes étant satisfaites automatiquement, aucune opération de forçage a posteriori n'est requise, et les maillages produits sont de Delaunay. En raison du coût prohibitif de la convergence en dimension 3, deux algorithmes efficaces sont proposés, sans qu'aucune conjecture sur leur convergence ne soit formulée. A titre d'application de la méthode en dimension 3, l'interfaçage avec un mailleur de Delaunay contraint existant est étudié. En particulier, la pertinence de la méthode est illustrée grâce au déblocage d'une configuration que ce mailleur ne parvient pas à résoudre. D'autres applications possibles, ainsi que les développements en cours sont également évoqués. |
|---|