Les travaux de cette thèse portent sur des méthodes de volumes finis sur maillages quelconque pour la discrétisation de problèmes d'évolution non linéaires modélisant le transport de contaminants en milieu poreux et les écoulements diphasiques.Au Chapitre 1, nous étudions une famille de schémas numériques pour la discrétisation d'une équation parabolique dégénérée de convection-reaction-diffusion modélisant le transport de contaminants dans un milieu poreux qui peut être hétérogène et anisotrope. La discrétisation du terme de diffusion est basée sur une famille de méthodes qui regroupe les schémas de volumes finis hybrides, de différences finies mimétiques et de volumes finis mixtes. Le terme de convection est traité à l'aide d'une famille de méthodes qui s'appuient sur les inconnues hybrides associées aux interfaces du maillage. Cette famille contient à la fois les schémas centré et amont. Les schémas que nous étudions permettent une discrétisation localement conservative des termes d'ordre un et d'ordre deux sur des maillages arbitraires en dimensions d'espace deux et trois. Nous démontrons qu'il existe une solution unique du problème discret qui converge vers la solution du problème continu et nous présentons des résultats numériques en dimensions d'espace deux et trois, en nous appuyant sur des maillages adaptatifs.Au Chapitre 2, nous proposons un schéma de volumes finis hybrides pour la discrétisation d'un problème d'écoulement diphasique incompressible et immiscible en milieu poreux. On suppose que ce problème a la forme d'une équation parabolique dégénérée de convection-diffusion en saturation couplée à une équation uniformément elliptique en pression. On considère un schéma implicite en temps, où les flux diffusifs sont discrétisés par la méthode des volumes finis hybride, ce qui permet de pouvoir traiter le cas d'un tenseur de perméabilité anisotrope et hétérogène sur un maillage très général, et l'on s'appuie sur un schéma de Godunov pour la discrétisation des flux convectifs, qui peuvent être non monotones et discontinus par rapport aux variables spatiales. On démontre l'existence d'une solution discrète, dont une sous-suite converge vers une solution faible du problème continu. On présente finalement des cas test bidimensionnels.Le Chapitre 3 porte sur un problème d'écoulement diphasique, dans lequel la courbe de pression capillaire admet des discontinuité spatiales. Plus précisément on suppose que l'écoulement prend place dans deux régions du sol aux propriétés très différentes, et l'on suppose que la loi de pression capillaire est discontinue en espace à la frontière entre les deux régions, si bien que la saturation de l'huile et la pression globale sont discontinues à travers cette frontière avec des conditions de raccord non linéaires à l'interface. On discrétise le problème à l'aide d'un schéma, qui coïncide avec un schéma de volumes finis standard dans chacune des deux régions, et on démontre la convergence d'une solution approchée vers une solution faible du problème continu. Les test numériques présentés à la fin du chapitre montrent que le schéma permet de reproduire le phénomène de piégeage de la phase huile.

Ce travail porte sur l'étude d'une méthode d'éléments finis discontinus (ou méthode de type Galerkin Discontinu, DGTD) basée sur l'utilisation d'un maillage héxaédrique régulier, proposée pour la résolution des équations de Maxwell dans le domaine temporel, afin de l'adapter à la caractérisation de structures rayonnantes et de l'associer à des techniques de domaines fictifs.<br />On présente tout d'abord une méthode Galerkin Discontinu s'appuyant sur une formulation centrée pour approcher les flux numériques aux interfaces du maillage et sur un schéma en temps explicite de type saute-mouton. Ainsi, le schéma obtenu est non-diffusif, stable, peu dispersif, parfaitement adapté à l'utilisation de maillages localement raffinés de manière non-conforme. La méthode a été dotée de parois absorbantes performantes (modèle Unsplit-PML), permettant de prendre en compte facilement des objets à cheval entre le domaine de calcul et la couche absorbante. Nous avons ensuite utilisé la méthode pour effectuer des calculs d'impédances, de paramètres S et de T.O.S. sur des structures rayonnantes planaires. La comparaison entre la simulation et la mesure de ces structures montre le bon fonctionnement de la méthode.<br />Nous avons alors couplé une méthode de domaines fictifs avec la méthode DGTD afin de prendre en compte la présence d'objets métalliques à géométries complexes. La méthode des domaines fictifs utilise deux maillages de manière indépendante: un maillage cartésien, pour faire évoluer le champ électromagnétique dans l'espace libre, et un maillage surfacique qui permet de prendre en compte l'objet métallique. La convergence de la méthode (pour la méthode FDTD) est liée à une relation de compatibilité entre le maillage volumique et le maillage surfacique: le plus petit élément du maillage surfacique impose la taille des éléments du maillage volumique. Ainsi, pour des objets présentant de tout petits détails, cette condition n'est assurée que si le maillage volumique est de l'ordre du plus petit élément du maillage surfacique, ce qui peut devenir extrêmement contraignant sans le recours à des techniques de raffinement local, telle que celle rendue possible par la méthode Galerkin Discontinu et utilisée ici.

Ce travail est dédié à l´étude et la simulation numérique de la rupture dynamique des séismes en deux et trois dimensions d´espace par une méthode d´éléments finis discontinus. Après avoir transformé le système de l´élastodynamique en un système hyperbolique symétrique du premier ordre, nous proposons un schéma numérique basé sur des flux centrés et un schéma explicite en temps de type saute-mouton. A travers l´étude d´une énergie discrète du système, nous spécifions les conditions aux limites sur la faille afin de prendre en compte de manière faible la rupture en mode cisaillant que nous traitons. Nous montrons, qu´en l´absence de tractions tangentielles sur la faille, cette énergie est parfaitement conservée. Nous illustrons la capacité de notre méthode à travers divers cas tests sur des configurations complexes grâce à une implémentation parallèle.

Cette thèse aborde la problématique de la modélisation multiphysique en génie électrique, avec une application à l'étude des vibrations d'origine électromagnétique des cages de développantes. Cette étude comporte quatre parties : la construction de la densité de courant, le calcul des forces locales, le transfert de solutions entre maillages et la résolution des problèmes couplés. Un premier enjeu est de correctement représenter les courants, cette opération est effectuée en deux étapes : la construction de la densité de courant et l'annulation de la divergence. Si des structures complexes sont utilisées, l'imposition du courant ne peut pas toujours être réalisée à l'aide de méthodes analytiques. Une méthode basée sur une résolution électrocinétique ainsi qu'une méthode purement géométrique sont testées. Cette dernière donne des résultats plus proches de la densité de courant réelle. Parmi les nombreuses méthodes de calcul de forces, les méthodes des travaux virtuels et des forces de Laplace, considérées par la littérature comme les plus adaptées au calcul des forces locales, ont été étudiées. Nos travaux ont montré que bien que les forces de Laplace sont particulièrement précises, elles ne sont pas valables si la perméabilité n'est plus homogène. Ainsi, la méthode des travaux virtuels, applicable de manière universelle, est préférée. Afin de modéliser des problèmes multi-physiques complexes à l'aide de plusieurs codes de calculs dédiés, des méthodes de transferts entre maillages non conformes ont été développées. Les procédures d'interpolations, les méthodes localement conservatives et les projections orthogonales sont comparées. Les méthodes d'interpolations sont réputées rapides mais très diffusives tandis que les méthodes de projections sont considérées comme les plus précises. La méthode localement conservative peut être vue comme produisant des résultats comparables aux méthodes de projections, mais évite l'assemblage et la résolution de systèmes linéaires. La modélisation des problèmes multi-physiques est abordée à l'aide des méthodes de transferts de solutions. Pour une classe de problème donnée, l'assemblage d'un schéma de couplage n'est pas unique. Des tests sur des cas analytiques sont réalisés afin de déterminer, pour plusieurs types de couplages, les stratégies les plus appropriées.Ces travaux ont permis une application à la modélisation magnéto-mécanique des cages de développantes est présentée.

La modélisation numérique des décharges plasma joue un rôle important dans la compréhension des mécanismes physiques ou chimiques ayant lieu dans les dispositifs assistés par plasma. Une grande partie de ces mécanismes est déjà prise en compte dans les codes actuels. En revanche, beaucoup d'entre eux ne permettent pas de travailler avec des géométries complexes. Cette limitation provient essentiellement de l'utilisation de maillages structurés, cartésiens. Ceux-ci ne sont pas bien adaptés aux géométries courbes. Les calculs en maillages structurés deviennent rapidement compliqués et spécifiques à une géométrie donnée. Notre travail concerne la modélisation de décharge pour un réacteur de traitement à la pression atmosphérique développé par Dow Corning. Sa configuration complexe ainsi que ses grandes dimensions nous ont incités à faire un nouveau code fonctionnant en maillages non structurés. Celui-ci doit être capable de s'adapter à la présence d'une pointe, d'arrondis et de multiples diélectriques mais aussi permettre le passage rapide à de nouvelles géométries. De plus ses grandes dimensions nécessitent l'utilisation de maillages raffinés uniquement aux endroits nécessaires (pointe, surfaces des diélectriques...). Le modèle mathématique utilisé est basé sur l'équation de Poisson couplée aux équations de transport de type dérive-diffusion. Plusieurs discrétisations numériques ont été testées dans des configurations physiques différentes. Nous présentons et validons les méthodes numériques choisies. Les résultats obtenus pour le réacteur Dow Corning sont alors exposés et commentés.

Ce manuscrit est dédié à la détection et la caractérisation de points d'intérêt dans les maillages. Nous montrons tout d'abord les limitations de la mesure de courbure sur des contours francs, mesure habituellement utilisée dans le domaine de l'analyse de maillages. Nous présentons ensuite une généralisation de l'opérateur SUSAN pour les maillages, nommé SUSAN-3D. La mesure de saillance proposée quantifie les variations locales de la surface et classe directement les points analysés en cinq catégories : saillant, crête, plat, vallée et creux. Les maillages considérés sont à variété uniforme avec ou sans bords et peuvent être réguliers ou irréguliers, denses ou non et bruités ou non. Nous étudions ensuite les performances de SUSAN-3D en les comparant à celles de deux opérateurs de courbure : l'opérateur de Meyer et l'opérateur de Stokely. Deux méthodes de comparaison des mesures de saillance et courbure sont proposées et utilisées sur deux types d'objets : des sphères et des cubes. Les sphères permettent l'étude de la précision sur des surfaces différentiables et les cubes sur deux types de contours non-différentiables : les arêtes et les coins. Nous montrons au travers de ces études les avantages de notre méthode qui sont une forte répétabilité de la mesure, une faible sensibilité au bruit et la capacité d'analyser les surfaces peu denses. Enfin, nous présentons une extension multi-échelle et une automatisation de la détermination des échelles d'analyse qui font de SUSAN-3D un opérateur générique et autonome d'analyse et de caractérisation pour les maillages

L'utilisation de modèles tridimensionnels dans les applications multimédia, prend de l'ampleur de jour en jour. Le développement des outils de modélisation, des scanners 3D, des cartes graphiques accélérées, du Web3D, etc. ouvre l'accès à des données tridimensionnelles de grande qualité. Les besoins, sans cesse croissants, concernant ce type de données, changent rapidement. S'il devient de plus en plus facile de créer de nouveaux modèles 3D, qu'en est-il du traitement et de l'analyse de ces modèles après leur création ? De nos jours, le concepteur d'objets 3D ne pose plus la question : " Comment créer un nouvel objet 3D ? ", mais plus vrai- semblablement " Comment retrouver un modèle 3D similaire à ceux en ma possession pour le réutiliser ? " et " Comment retrouver la structure d'un modèle 3D maillé sans connaissance a priori sur celui-ci ? " Cette habilitation a pour but d'apporter des éléments de réponse à ces deux questions. En réponse à la première question, nous avons développé un nouveau système bayésien pour retrouver des modèles 3D à partir d'une requête constituée d'une ou plusieurs vues 2D, ou d'un modèle 3D entier. Ce système a été testé dans un contexte applicatif industriel ainsi qu'avec un benchmark international. Chaque expérience a mis en évidence les excellents résultats de notre approche. La seconde question a été abordée sous l'angle de l'analyse topologique des maillages 3D grâce aux graphes de Reeb. Ce travail théorique a été appliqué à différents domaines comme la déformation automatique, l'indexation et la segmentation de maillages 3D. L'approche a toujours montré des résultats remarquables dans ces domaines. Finalement, la segmentation de maillages 3D, qui est une étape de pré-traitement fréquente avant d'autres analyses du maillage, a attiré notre attention. Nous avons proposé une métrique fiable et robuste pour la comparaison de segmentations et l'évaluation des performances des méthodes de segmentation de maillages 3D, ainsi qu'une approche de la segmentation par apprentissage qui surpasse les méthodes existantes. Pour terminer, de nouvelles pistes de recherche sur les maillages 3D sont ouvertes.

Comprendre, analyser et modéliser l'environment 3D à partir d'images provenant de caméras et d'appareils photos est l'un des défis majeurs actuel de recherche en vision par ordinateur. Cette thèse s'interesse à plusieurs aspects géométriques et photometriques liés à la reconstruction de surface à partir de plusieurs caméras calibrées. La reconstruction 3D est vue comme un problème de rendu inverse, et vise à minimiser une fonctionnelle d'énergie afin d'optimiser un maillage triangulaire représentant la surface à reconstruire. L'énergie est définie via un modèle génératif faisant naturellement apparaître des attributs tels que la visibilité ou la photométrie. Ainsi, l'approche présentée peut indifférement s'adapter à divers cas d'application tels que la stéréovision multi-vues, la stéréo photométrique multi-vues ou encore le "shape from shading" multi-vues. Plusieurs approches sont proposées afin de résoudre les problèmes de correspondances de l'apparence pour des scènes non Lambertiennes, dont l'apparence varie en fonction du point de vue. La segmentation, la stéréo photométrique ou encore la réciprocité d'Helmholtz sont des éléments étudiés afin de contraindre la reconstruction. L'exploitation de ces contraintes dans le cadre de reconstruction multi-vues permet de reconstruire des modèles complets 3D avec une meilleure qualité.

Les motoristes aéronautiques souhaitent disposer de la représentation la plus fidèle possible du fonctionnement des propulseurs, dans une perspective d'amélioration continue de leurs performances. Les modèles numériques doivent donc intégrer au maximum les détails géométriques susceptibles d'influencer la physique de l'écoulement analysé. La prise en compte de tels effets technologiques s'avère difficile dans le contexte des solveurs structurés disponibles.Une stratégie hybride de prise en compte des effets technologiques fait coexister au sein d'un même domaine de calcul des zones structurées et non structurées. La flexibilité de génération d'un maillage non structuré permet une prise en compte aisée des détails géométriquement complexes tandis que la préservation de zones structurées dans une majeure partie du domaine de calcul permet de bénéficier de l'efficacité d'un solveur structuré. La présente thèse contribue au développement de cette stratégie hybride au sein du solveur elsA de l'ONERA en proposant des gains de précision et de robustesse par rapport à la version initialement développée pour établir la faisabilité et l'intérêt de l'approche. Après un état de l'art des techniques de discrétisation spatiale disponibles dans cette version initiale, différentes améliorations (techniques de moindres carrés, approche dite quasi-Green, méthode d'estimation des gradients aux faces) ont été analysées puis implémentées et validées sur des cas académiques. Le choix d'une stratégie hybride avec raccords coïncidents entre zones structurées et non-structurées conduit à des déformations de maillage dans la zone d'interface structuré / non-structuré qui ont exigé le développement supplémentaire de techniques d'amélioration de la robustesse (limiteurs physiques ou géométriques). Le solveur hybride rassemblant ces différentes fonctionnalités a permisde simuler avec succès des géométries d'aubes isolées dotées d'effets technologiques tels que congé de raccordement, trous de refroidissement, fentes de bord de fuite, cheminées internes d'alimentation. Enfin, une stratégie permettant l'utilisation de l'approche hybride en étage complet a été proposée et appliquée à la simulation hybride de l'interaction rotor/stator pour la configuration VKI-BRITE CT3, en stationnaire et en instationnaire, respectivement via une condition de plan de mélange et une condition de chorochronicité. / The aerospace engine manufacturers wish to rely on the most accurate description of their propulsion systems in order to continuously improve their performance levels. Therefore, numerical models must include as much as possible geometrical details likely to impact the physics of the flow under study. Taking into account such technological effects turns out to be a difficult task when working with available structured solvers. A hybrid strategy takes advantage of structured and unstructured zones within the same computational domain in order to efficiently describe technological effects. Geometrically complex local details are easily accounted for thanks to the flexibility of unstructured grid generation while keeping structured zones in the remainder of the flow domain allows to benefit from the tried and tested structured solver efficiency. The present work contributes to the development of such a hybrid strategy in ONERA elsA solver and enhances accuracy and robustness with respect to the solver initially developed to establish the feasibility and interest of hybridization. Following a review of the space discretization techniques available in the initial solver, several improvements (least square techniques, quasi-Green approach, computation of face gradients) have been analysed, then implemented and validated for academic test-cases. The choice of a hybrid strategy with coincident matching between structured and unstructured zones leads to highly deformed cells in the structured / unstructured interface region, requiring the development of supplementary robustness improvement techniques (physics- or geometry-based limiters). The hybrid solver gathering these various options allows to successfully compute isolated blade geometries including technological effects such as blade fillet, cooling holes, trailing edge cutbacks, internal coolant supply channel. Finally, a structured / unstructured strategy has been proposed and applied to the hybrid simulation of a rotor/stator interaction for the steady and unsteady

L’objectif de cette thèse est de répondre au besoin d’accélérer la restitution des résultats de calcul d’un code CFD pour la simulation d’écoulements hydrodynamiques quasi-incompressibles. Ce code présente l’originalité de résoudre explicitement les équations de Navier-Stokes sous l’hypothèse de faible compressibilité avec des schémas numériques d’ordre élevé. Les développements effectués visent à réduire les temps de calcul à précision équivalente.Une première partie est consacrée à l’implémentation d’une formulation purement incompressible avec une résolution implicite de la pression par un schéma de projection. La formulation incompressible autorise des pas de temps plus grand en s’affranchissant de la vitesse du son, mais au prix d’une algorithmique plus complexe et de la nécessité de résoudre un système linéaire. La comparaison des deux formulations,faiblement-compressible et incompressible, tend à montrer la pertinence du schéma de projection pour les écoulements laminaires instationnaires.Un deuxième axe de développement a consisté en la proposition d’une amélioration de la méthode de frontière immergée initialement présente dans le code.Si les résultats obtenus ne sont pas encore pleinement satisfaisants, ils montrent que la montée en ordre d’une méthode de frontière immergée peut être moins contraignante en formulation incompressible.Enfin la dernière partie traite de l’immersion rapide et robuste de géométries complexes telles qu’elles peuvent être rencontrées dans l’industrie. La localisation géométrique par arbre octal permet d’évaluer rapidement une fonction de distance signée indispensable pour la méthode de frontière immergée. / An original strategy to address hydrodynamic flow was recently proposed through a high-order weakly-compressible Cartesian grid approach. The method is based on a fully-explicit temporal scheme for solving the Navier-Stokes equations. The present thesis aims to reduce the computational time required to obtain the results without deteriorating the accuracy.A first part is dedicated to the implementation of a truly incompressible formulation with an implicit solution for the pressure field through a projection scheme. The incompressible solver allows larger time step size for time integration since the speed of sound tends to infinity. In return the algorithms are no longer straight forward and a linear system has to be solved through the Pressure Poisson Equation. Comparisons carried out between both formulations show that the projection scheme can be better adapted to efficiently simulate unsteady viscous flows. Then an improvement of the immersed boundary method has been proposed. Results are not fully satisfactory for now. However, it seems easier to develop a numerical scheme for the incompressible approach rather than the weakly-compressible one.Finally, the last part addresses the setup up of complex triangulations in immersed boundary simulations. A fast and robust procedure is developed for distance computation with an octree data structure.