National audienceCet article représente un résumé des différents travaux effectués jusqu'à maintenant dans le cadre de ma thèse dans le département SIC du laboratoire XLIM et le département d2 (Fluide, Thermique, Combustion) du laboratoire P' de Poitiers sur les méthodes de reconstruction de volumes de particules 3d appliquées à la tomographie de vélocimétrie par image de particule (Tomo-PIV). Ces travaux s'inscrivent dans le cadre d'un projet européen intitulé AFDAR (Advanced Flow Diagnostics for Aeronautical Research project no. 265695). L'objectif de la thèse est de proposer une nouvelle technique de mesure tridimensionnelle de vitesse en mécanique des fluides par Tomo-PIV et de l'appliquer sur une étude expérimentale de la dissipation d'énergie turbulente derrière une grille

Alata, Olivier

Ben-Salah, Riadh

David, Laurent

Thomas, Lionel

Tremblais, Benoit

HAL-UJM

Nouvelle me´thode de reconstruction de volumes departicules 3D base´e sur les processus ponctuels marques: Application a` la TOMO-PIV en me´canique de fluideRiadh Ben-Salah, Olivier Alata, Benoit Tremblais, Lionel Thomas, LaurentDavidTo cite this version:Riadh Ben-Salah, Olivier Alata, Benoit Tremblais, Lionel Thomas, Laurent David. Nouvelleme´thode de reconstruction de volumes de particules 3D base´e sur les processus ponctuels mar-ques : Application a` la TOMO-PIV en me´canique de fluide. Workshop XLIM 2013, Sep 2013,Limoges, France. <hal-00913346>HAL Id: hal-00913346https://hal.archives-ouvertes.fr/hal-00913346Submitted on 3 Dec 2013HAL is a multi-disciplinary open accessarchive for the deposit and dissemination of sci-entific research documents, whether they are pub-lished or not. The documents may come fromteaching and research institutions in France orabroad, or from public or private research centers.L’archive ouverte pluridisciplinaire HAL, estdestine´e au de´poˆt et a` la diffusion de documentsscientifiques de niveau recherche, publie´s ou non,e´manant des e´tablissements d’enseignement et derecherche franc¸ais ou e´trangers, des laboratoirespublics ou prive´s.Workshop 26 Septembre 2013  NOUVELLE METHODE DE RECONSTRUCTION DE VOLUMES DE PARTICULES 3D BASEE SUR LES PROCESSUS PONCTUELS MARQUES : APPLICATION A LA TOMO-PIV EN MECANIQUE DE FLUIDE R. Ben Salah1, 2, O. Alata3, B. Tremblais2, L. Thomas1, L. David1  1  Institut Pprime – Département D2 – UPR CNRS 3346 2  Département XLIM-SIC – UMR CNRS 7252  3 Laboratoire Hubert Curien – UMR CNRS 5516  Résumé. Cet article représente un résumé des différents travaux effectués jusqu’à maintenant dans le cadre de ma thèse dans le département SIC du laboratoire XLIM et le département d2 (Fluide, Thermique, Combustion) du laboratoire P’ de Poitiers sur les méthodes de reconstruction de volumes de particules 3d appliquées à la tomographie de vélocimétrie par image de particule (Tomo-PIV). Ces travaux s’inscrivent dans le cadre d’un projet européen intitulé AFDAR (Advanced Flow Diagnostics for Aeronautical Research project no. 265695). L’objectif de la thèse est de proposer une nouvelle technique de mesure tridimensionnelle de vitesse en mécanique des fluides par Tomo-PIV et de l’appliquer sur une étude expérimentale de la dissipation d’énergie turbulente derrière une grille. I. INTRODUTION En 1914, le radiologue Karol Mayer a proposé le principe d’une technique nommée Tomographie. Les techniques tomographiques consistent à sonder la matière avec un type de rayonnement. Les rayonnements émergents vont jouer le rôle d’une sonde pour  effectuer des mesures et fournir une analyse de la matière explorée. Grâce à des calculs de reconstruction, elle fournit une cartographie du paramètre caractéristique du rayonnement étudié [1]. Depuis quelques années, la mécanique des fluides expérimentale s’est dotée d’une nouvelle méthode de mesure non-intrusive, qui permet d’acquérir des champs instantanés de vitesses. Cette méthode est la Vélocimétrie par Image de Particules (PIV). La PIV a pour origine la visualisation d’écoulements (flow visualisation) afin d’évaluer la vitesse ou d’avoir une idée sur les caractéristiques turbulentes d’un écoulement [2].  Récemment une nouvelle approche, appelée tomographie PIV (Tomo-PIV) fondée sur les systèmes tomographiques et la PIV, et étendue aux volumes, a été proposée. Elle permet la mesure instantanée des trois composantes de la vitesse dans un volume de mesure complet 3D. La PIV tomographique est basée sur le principe de la reconstruction tomographique en volume. Le principe de la Tomo-PIV consiste à ajouter des particules dans l’écoulement d'intérêt et à éclairer une portion du volume de l'écoulement par une source laser. La lumière diffusée par les particules est capturée par plusieurs caméras numériques de hautes résolutions suivant des angles d'observation différents. L’ensemble de projections de chaque caméra est fourni à un algorithme de reconstruction 3D. La plupart des solutions proposées reposent sur l'utilisation de méthodes algébriques de tomographie (ART, MART, SMART, MLOS...). II. POSITION DU PROBLÈME Les techniques de reconstruction tomographique proposées à ce jour ne prennent pas en compte la forme particulière des objets à reconstruire. En plus, la forte densité des particules dans les images, le bruit et les vibrations inhérentes aux installations expérimentales sont des éléments défavorables à une reconstruction de grande qualité. Vu la taille des données, le temps de traitement de ces méthodes et l’espace mémoire utilisé sont très élevés. Afin de résoudre ces problèmes, une solution se basant sur la parcimonie des volumes de particules à reconstruire pourrait être envisagée. Dans ce contexte, une méthode de reconstruction de volume de particules 3d basée sur les approches par processus “objet” tels que des processus ponctuels marqués semblent particulièrement bien adaptées à ces problèmes. III. PROCESSUS PONCTUEL MARQUE Un processus ponctuel marqué [3] [4] est un processus aléatoire dont les réalisations sont des configurations aléatoires d’objets. Chaque objet y est modélisé par un point (sa position), ainsi que par des marques (sa dimension, son orientation, sa forme, etc.) qui correspondent à ses propriétés. Ces processus fournissent donc une représentation naturellement parcimonieuse d’une configuration d’objets d’intérêt, présents dans une image. Ils permettent en effet de se détacher du modèle numérique induit par l’image, pour mieux se rapprocher du modèle physique. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Cam 1 Cam 2 Cam 3 Cam 4 t t + Δt t  + Δt Acquisition des projections Cam 1 Cam 4 Cam 3 Cam 2 Cam 1 Ecoulement Volume d’acquisition Particules traceurs Reconstruction tomographique Volume à T= t Volume à T= t + Δt Corrélation Champ de vitesse Figure 1: Principe de la Tomo-PIV Workshop 26 Septembre 2013 IV. APPROCHE PROPOSEE L’idée de cet algorithme consiste à reconstruire un volume de particules 3d en se basant sur l’énergie présente dans les projections acquises. Contrairement aux méthodes de reconstruction tomographique classiques, cette méthode part à partir d’un volume 3d inconnu, qui représente ici la population des particules 3d. Cette population sera obtenue à l'aide d'un algorithme de type Monte-Carlo par chaîne de Markov (MCMC) exploitant une dynamique de Metropolis-Hasting-Green: à chaque itération, il peut être envisagé un mouvement de naissance, de mort ou de déplacement d'une particule. Cette dynamique permet de simuler une densité de probabilité de type distribution de Gibbs dont l'énergie est composée de la somme d'une énergie d'attache aux données et d'une énergie dite "interne". Cette énergie permet de contraindre le comportement du processus d'une part en fonction de l'information que l'on cherche dans les données et d'autre part en fonction des interactions susceptibles d'apparaître entre les particules. V. TERME D’ATTACHE AUX DONNÉES Le terme d'attache aux données a été développé à partir de l'erreur quadratique moyenne obtenue entre la projection d'une population et les images acquises. Il doit permettre d'obtenir une population de particules cohérentes par rapport aux données acquises. Ce terme permettra au processus de converger vers la population la plus proche de la réalité. Il se base sur l’énergie provenant du voisinage des particules 2d sur chaque projection générée par le processus et la corrélation entre une particule 2d et ce qui lui correspond dans les projections réelles. La formule du terme d’attache aux données pour une particule est la suivante : �! � = �!,! � + �!,! � !!,!!~!!→!!!!! Avec, �!,! � = �!→! �!→! − 2�!,! !∈!!!→!!!!! Ce terme correspond au potentiel associé à une particule 3d (ξ). I! correspond à une projection (réelles), n correspond au nombre de projections, s correspond à un site (pixel) sur cette projection. ξ → s signifie que la projetée de la particule 3d ξ contient le site s. r!→! est la valeur projetée en s par la particule �. �!,! correspond à la valeur observée au site s de la projection �!. La deuxième partie du terme d’attache aux données �! �  correspond à l’information apportée par le voisinage d’une particule. Le terme �!,! �  s’écrit sous la forme suivante : �!,! � = 2�!→!�!!→! !∈!!!→!,!!→! VI. QUELQUES RÉSULTATS Les simulations ont été appliquées sur 4 projections simulées. Les images utilisées, de taille 20 x 20 pixels et le volume de particules de référence, de taille 20 x 20 x 15 voxels, ont été générés avec un programme développé au sein du laboratoire P’- d2. Les images et le volume contiennent 30 particules de forme gaussiennes isotropes et de diamètre 3 pixels. L’intensité centrale de toute les particules est la même. Le modèle de caméra utilisé est  le modèle Pinhole. Les simulations ont été réalisées avec et sans  terme d’attache aux données. La figure-2 montre l’évolution de l’erreur quadratique moyenne (EQM) par rapport au nombre de mouvements acceptés avec et sans attache aux données.     On peut remarquer que le terme d’attache aux données a permis de minimiser l’EQM entre la population simulée et la population générée par l’algorithme. VII. CONCLUSION ET PERSPECTIVE Les courbes obtenues nous ont montré l’efficacité du terme d’attache aux données avec les paramètres utilisés. Nous étudierons par la suite l'influence des différents paramètres du processus. Puis, nous validerons notre approche sur des données réelles. Enfin, une étude comparative avec les méthodes de reconstruction précédemment proposées en Tomo PIV sera réalisée. Ces travaux ont été financés par AFDAR et l’ANR VIVE3D. Leur support est vivement remercié. REFERENCES  [1] P. C. Seynaeve and J. I. Broos. “The history of tomography”, J Belge Radiol 78:284–288, Medline. 1995. [2] M. Raffel, C. Willert and J. Kompenhans, “Particle Image Velocimetry,a practical guide”.Ed. Springer Verlag, 1998. [3] C. J. Geyer and J. Moller. “Simulation and likelihood inference for spatial point processs”. Scandinavian Journal of Statistics series B, 21:359–373. 1994. [4] D. J. Daley, and D. Vere-Jones. “An Introduction to the Theory of Point Processes”. Springer Verlag, 2002. Figure 2: Evolution de l’EQM par rapport au nombre de mouvements acceptés  0 0.5 1 1.5 2 2.5 3 3.5x 1050.140.160.180.20.220.240.260.265M ouvement accepteEQM  Sans attache aux donnéesAvec attache aux données

Nouvelle méthode de reconstruction de volumes de particules 3D basée sur les processus ponctuels marques : Application à la TOMO-PIV en mécanique de fluide

HAL-UNILIM

https://hal.archives-ouvertes.fr/hal-00913346

Nouvelle méthode de reconstruction de volumes de particules 3D basée sur les processus ponctuels marques : Application à la TOMO-PIV en mécanique de fluide

Abstract

Similar works

Full text

Available Versions

HAL-UJM

HAL-UNILIM