Je suis ingénieur à l'INRIA rattaché à l'équipe projet CALVI. Je travaille actuellement dans l'équipe EDP de l'IRMA, le laboratoire de mathématique de l'université de Strasbourg.
|
|
CALVI est une équipe-projet de l'Institut National de Recherche en Informatique et en Automatique (INRIA) bi-localisé à Nancy et à Strasbourg et est une collaboration de trois organismes de recherche. Dans ce cadre, j'ai contribué fortement à la conception et l'implantation d'une plate-forme de simulation numérique. Elle consiste en un ensemble de logiciels qui permettent d'implanter, valider et comparer rapidement différentes méthodes numériques pour le système Vlasov-Maxwell.
Le principe de la plate-forme consiste à décomposer une méthode numérique en un ensemble de modules, qu'on peut regrouper en trois grandes classes : les cas-tests, les exports et les calculs. Les modules cas-tests contiennent la description de tests standards nécessaires à la validation des méthodes numériques. Les modules d'exports servent à exporter les résultats de la simulation sur disque dans le but de les visualiser. Enfin les modules de calculs contiennent réellement les méthodes numériques qui résolvent des équations particulières. En conséquence, une nouvelle méthode numérique peut être facilement introduite dans la plate-forme en créant les modules de calculs correspondants. Dès lors, les modules de cas-test et d'export peuvent directement être utilisés pour valider le nouveau schéma et le comparer avec les schémas existants.
EUFORIA (EU Fusion fOR Iter Applications) est un projet financé par l'Union européenne dans le cadre du septième programme cadre (7e PC) qui va fournir un cadre et une infrastructure logicielle pour la simulation du transport de coeur, du transport de bord et de la turbulence dans un Tokamak. Cette infrastrucutre fera le lien entre la communauté des grilles de calcul et du calcul haute performance avec la communauté de la modélisation pour la fusion.
Ma contribution dans ce projet est de fournir un ensemble d'outils de visualisation unifiés pour les codes intégrés à la plate-forme. Certains de ces outils permettent d'effectuer des post-traitements sur les données générées par les codes, d'autres sont directement intégrés au système de workflow (Kepler) pour permettre la visualisation et le suivi des résultats durant l'exécution du workflow. Les outils de visualisation développés dans ce cadre sont basés sur des logicels libres existants tel que Python, numpy et matplotlib d'une part et VTK et VisIt d'autre part.
Cette partie de ma thèse contribue à la mise au point et l'implantation efficace et parallèle d'une méthode numérique adaptative résolvant le système de Vlasov-Poisson en une dimension. En utilisant une décomposition en ondelettes, j'ai modifié la méthode semi-lagrangienne pour qu'elle utilise une grille adatative plutôt qu'une grille régulière. La complexité des algorithmes a ainsi été réduite et l'erreur commise est contrôlée par un seuil de compression fixé par l'utilisateur. Une structure de données dédiée a permis une implantation efficace et parallèle qui reflète la complexité des algorithmes et présente un très bon passage à l'échelle jusqu'à 64 processeurs en mémoire partagée.
![]() |
![]() |
![]() |
| Fonction de distribution (couleur = f(x, vx)) |
Discrétisation régulière 512 x 512 = 256K points |
Maillage de calcul adaptatif 43K points (17%) |
Cette partie de ma thèse consiste à élaborer un schéma de compression pour la visualisation de fonctions multidimensionnelles. En particulier, je me suis intéressé aux fonctions 4D à valeurs scalaires issues des codes résolvant le système de Vlasov-Poisson en quatre dimensions. La visualisation proposée utilise une méthode d'extraction de tranches 2D où l'utilisateur sélectionne et modifie interactivement les tranches visualisées. La taille de la donnée 4D ( > 2 Go) étant le problème principal, les données sont compressées à l'aide d'une base hiérarchique d'éléments finis. La localité des algorithmes de compression a permis leur intégration à différents codes de simulation parallèle pour qu'il puisse directement exporter la fonction 4D compressée. La réduction substantielle de la taille des données (facteur 10 à 100) ainsi que l'utilisation de structures de données dédiées permet au logiciel de visualisation plasmaViz de générer les coupes demandées par l'utilisateur en temps interactif ( > 10 images/s).

|
M. Haefele Simulation adaptative et visualisation haute performance de plasmas et de faisceaux de particules LSIIT, 2007 |
|
|
|
M. Haefele, F. Zara, G. Latu, J-M. Dichler A dedicated Compression Scheme for Large Multidimensional Functions Visualization To appear in Workshop on Super Visualization - Kos (Greece) 2008 |
|
M. Haefele, G. Latu, M. Gutnic A parallel Vlasov solver using a Wavelet based Adaptive Mesh Refinement. In ICPP'2005 7th workshop HPSEC Oslo, Norway, 2005 |
|
|
|
M. Gutnic, M. Haefele, I. Paun and E. Sonnendrücker Vlasov simulation on an adaptive phase space grid Computer Physic Communication, 164, 214-219, 2004 |
|
|
|
M. Gutnic, M. Haefele, I. Paun and E. Sonnendrücker Moments conservation in adaptive Vlasov solver International Computational Accelerator Physics Conference St Petersburg 2004. |
|
|
|
E. Sonnendrücker, M. Gutnic, M. Haefele and J.L. Lemaire Vlasov simulation of beams and halo. Particule Accelerator Conference Knoxville, U.S.A., 2005 |
|
|
|
S. Brandel, M. Haefele and D. Bechmann. A geological application in immersive virtual environments. Virtual Reality International Conference Laval, France, 2003. |
|