Précédent : Méthodes nodales Remonter : Actions de recherche Suivant : Actions industrielles

Méthodes numériques parallèles pour la propagation d'ondes

Participants : Olivier Anthore, Michel Kern

A la suite des travaux de G. Cohen et P. Joly (projet Ondes), nous nous sommes intéressés aux méthodes d'éléments finis d'ordre élevés pour la simulation des phénomènes d'ondes, dans le contexte de la décomposition de domaine. Dans la mesure où ces éléments finis ont été conçus pour conduire à des schémas explicites, il s'agit là d'une idée naturelle et assez simple à mettre en oeuvre. Nous avons mené trois actions sur ce sujet:

élastodynamique :

Pour son stage de DESS, Olivier Anthore a repris le code écrit par C. Tsogka qui traite la propagation d'ondes élastiques dans un milieu invariant par translation selon une direction et l'a réimplémenté en parallèle avec MPI.

acoustique :

M. Kern a poursuivi le travail débuté par S. Mbouayouéou sur les ondes acoustiques. Il s'agit d'un code 2D, utilisant MPI pour les communications. Le code a été essentiellement réécrit, et plusieurs cas tests montrent son bon comportement, et des performances correctes sur un nombre limité de processeurs. Le code a tourné sur le CRAY T3D du CEA/CENG (Grenoble), sur l'IBM SP/2 du CNUSC (Montpellier) et sur le SGI Power Challenge Array du Centre Charles Hermite (Nancy).

éléments finis en 3D :

Le parallélisme ne prend réllement son sens que pour les modèles 3D. La recherche de nouveaux éléments sur un tétraèdre a fait l'objet du stage de F. Hannoyer. Ces travaux sont décrits dans le rapport d'activité du projet Ondes.